struct sFace{  struct EqrPlane	{		GLfloat a,b,c,d;	}eqr;    int vertIndex[3];     int neigh_face[3];	bool facelight;};struct s3DObject{   int numofFaces;   sFace *psFaces;};

void CShadow3DS::ImportShadowModel(GLuint Model_id, char *strFileName){   ImportModel(Model_id, strFileName); //步骤1,载入模型   s3DObject newsobject= {0};　//步骤2,这里是初始化shadow版本模型:遍历原模型各个对象,把该对象中有用的信息:面数,面索引等存入新结构s3DObject ,s3Face 中,并初始化各面的邻面索引为-1(无),非朝光,其余填充0   for(vector<t3DObject>::iterator p_r = GetModel().pObject.begin(); p_r!=GetModel().pObject.end(); p_r++)   {       newsobject.numofFaces= p_r->numOfFaces;	   newsobject.psFaces = new sFace [newsobject.numofFaces];	   memset(newsobject.psFaces, 0, sizeof(sFace) * newsobject.numofFaces);	   for(int j = 0 ; j< newsobject.numofFaces; j++)			   for(int i = 0 ; i<3; i++)		   {			   newsobject.psFaces[j].vertIndex[i] = p_r->pFaces[j].vertIndex[i];			   newsobject.psFaces[j].neigh_face[i] = NO_NEIGH_FACE;			   }//把遍历中每次得到的新结构压入vector<s3DObject >psObject中,它能满足shadow版本模型的各种数据存储需求	   psObject.push_back(newsobject);       }    SetConnectivity();//步骤3,找邻面//步骤4,对每个面,进行Face to Plane的转换   for(int iter_obj = 0; iter_obj < GetModel().numOfObjects; iter_obj++)	 for(int iter_face = 0 ; iter_face < GetModel().pObject[iter_obj].numOfFaces; iter_face++)				 FaceToPlane(iter_obj, iter_face); }

void CShadow3DS::RenderShadowModel(){   	ToModelCoord();//步骤1:转换光源到模型空间	glLoadIdentity();//把当前矩阵作为单位矩阵（接下来的变换以此为基准）   RenderModel();//步骤2:渲染原Model 	RenderShadow();//步骤345,判断并渲染阴影}//先看看步骤345:void CShadow3DS::RenderShadow(){....//转换光源到模型空间必要的反变换,后述 ShadowCast() {    FaceLightJudging();//步骤3.找出每帧中哪些面具有朝光性   .........//步骤4:接下来是蒙板值非0之判断,代码跟第一篇中三角形面片的ShadowCast()代码是一样的,功能也一样,只不过这里仅处理朝光的面,并且画ShadowVolume的方法(ShadowVolume() )不同   ShadowVolume() ;//具体的画法思路见"第二篇-第2个问题"  ..........//注意阴影锥各侧面（silhouette edge)分正面朝观察者与背面朝贯彻者观察者渲染两次.逆时针为正的保证由3DS模型自己提供(通过顶点索引)  ShadowVolume() ;   ........  }}

void CShadow3DS::ToModelCoord(){	CVector4D antiRotLight,antiTransLight;	CMatrix16 current_mvmatrix;//分别保存MV矩阵的旋转部分（包括缩放）和移动部分	antiRotLight = LightPos;    antiTransLight.set(0,0,0,1);//把当前矩阵作为单位矩阵（接下来的变换以此为基准），glLoadIdentity()其实也在压栈	glLoadIdentity(); //根据CLoad3DS中的模型变换，作其反变换  	glRotatef(-GetRot().y,0,1,0);  		glRotatef(-GetRot().z,0,0,1);	glRotatef(-GetRot().x,1,0,0);    glRotatef(-1, 1, 0, 0);	glScalef(1/GetSize(),1/GetSize(),1/GetSize());//存储旋转部分为：MV矩阵乘以光源实际坐标	glGetFloatv(GL_MODELVIEW_MATRIX,current_mvmatrix.mt);	antiRotLight = current_mvmatrix * antiRotLight;//继续	glTranslatef(-GetPos().x, -GetPos().y, -GetPos().z);//存储平移部分为：MV矩阵乘以（0，0，0，1）		glGetFloatv(GL_MODELVIEW_MATRIX,current_mvmatrix.mt);	 	antiTransLight = current_mvmatrix * antiTransLight;	 //结合这两部分到mvLightPos:在模型空间和视图空间之间玩无间道的家伙	     mvLightPos.x = antiRotLight.x + antiTransLight.x;    mvLightPos.y = antiRotLight.y + antiTransLight.y;    mvLightPos.z = antiRotLight.z + antiTransLight.z;	    mvLightPos.w = 1.0;	}//.......接下来在RenderShadow()中,不压栈地作正变换,使两者变换抵消,实际中因而只剩下RenderModel()中受glPushMatrix()保护的正变换,一切跟没做这些的时候一致,做这些处理只是为了得到上面那个mvLightPos//mvLightPos,应用于判断面朝光性判断和画ShadowVolume,它在模型空间中可以不断变化,但在视图空间(基于实际视觉)中永远与实际光源坐标重合,即相对于光源不动.(注意,模型在这个过程中与之相反.)