基于bSpline插值的直纹面与旋转曲面.docVIP

基于bSpline插值的直纹面和旋转曲面（上）   此程序主要分为两部分，一部分是抽象旋转曲面和直纹面的类，另一部分是视图界面的类组织，二者通过在视图界面定义数据成员建立联系。   第一部分： 旋转曲面类：   class Revsurf:public bspline,public Mesh                         //旋转曲面 { public:        Revsurf(Point3* pt1,int num,int uu,int vv,int mm);             //构造函数        ~Revsurf(); public:        void writeMesh(char *fname);                             //将生成的点、法线、面片写入文件        void makefaces();                                       //生成面片及索引        void makenormals();                                     //生成法线及索引        int ctlpt_num;                                          //控制点的个数        int v;                                                  //经度方向        int u;                                                  //纬度方向        Point3 *ctlpt;                                           //保存控制点        void makepoints();                                       //生成点及索引 };   Revsurf::Revsurf(Point3 *pt1,int num,int uu,int vv,int mm)                                          //旋转曲面 :ctlpt(pt1),ctlpt_num(num),u(uu),v(vv),bspline(mm,num) {        numVerts=u*v;        numNorms=u*v;        numFaces=(u-1)*v;            pt=new Point3[numVerts];        norm=new Vector3[numNorms];        face=new Face[numFaces]; }   Revsurf::~Revsurf() { }   /******************************旋转曲面生成顶点*********************************/ void Revsurf::makepoints() {        int i=0,j=0,k=0;        float t=0,x=0,y=0,z=0;        Point3 *p;                                        //记录两条u线的点        p=new Point3[u];          //细节处理1e-5        for(t=knot[0]+1e-5;iu;t+=knot[ctlpt_num+m-1]/(float)(u-1)-1e-5)                           //采集两条u线的点        {               p[i].set(0,0,0);               for(k=0;kctlpt_num;k++)               {                      p[i].x+=ctlpt[k].x*getbsplineValue(k,m,t,knot);                      p[i].y+=ctlpt[k].y*getbsplineValue(k,m,t,knot);                      p[i].z+=ctlpt[k].z*getbsplineValue(k,m,t,knot);               }               i++;        }          glBegin(GL_LINE_STRIP