【文章內容簡介】 m_p1y))。 if(dis) goto loop。 }loop: return (dis) ? : 。}/* */double evaluate_point( Vec3 pos,Octree *o ){ return demo9( pos,o )。}int evaluate1_point( Octree *o ){ int tCnt=o()+1。 return(tCnt)。}/* *//* returns 1 if the octree should be split, 0 otherwise *//* this implementation checks whether all 8 corner values are the same (if any corner 1..7 is different to corner 0 then the function returns 1)*/int octree_needs_to_be_split( Octree* o ){ /*int i。 double v = ovalue[0]。 for(i=1。 i 8。 i++) if( ovalue[i] != v) return 1。 /* if we got here, then all corners have the same value */ //return 0。 if(odensity40) return 1。 else return 0。}void isoface(Octree* o){ char bf[2]。 //CPtrList vexlist,norlist。 float v1,v2,v3。 ifstream file()。 if(!o()) o()。 if(!o()) o()。 if(!()) { while(!()) { filebfv1v2v3。 if(bf[0]==39。v39。amp。amp。bf[1]==0) o(new pvex_nor(v1,v2,v3))。 if(bf[0]==39。v39。amp。amp。bf[1]==39。n39。) o(new pvex_nor(v1,v2,v3))。 } //ovex=vexlist。 } ()。}/* Linearly interpolate the position where an isosurface cuts an edge between two vertices, each with their own scalar value*/Vec3 VertexInterp(double isolevel,Vec3 p1,Vec3 p2,double valp1,double valp2){ double mu。 Vec3 p= makeVec3(0,0,0)。 if (fabs(isolevelvalp1) ) return(p1)。 if (fabs(isolevelvalp2) ) return(p2)。 if (fabs(valp1valp2) ) return(p1)。 mu = (isolevel valp1) / (valp2 valp1)。 p[0] = p1[0] + mu * (p2[0] p1[0])。 p[1] = p1[1] + mu * (p2[1] p1[1])。 p[2] = p1[2] + mu * (p2[2] p1[2])。 return(p)。}/* Given a grid cell and an isolevel, calculate the triangular facets required to represent the isosurface through the cell. Return the number of triangular facets, the array triangles will be loaded up with the vertices at most 5 triangular facets. 0 will be returned if the grid cell is either totally above of totally below the isolevel.*/int Polygonise(GRIDCELL grid,double isolevel,TRIANGLE *triangles){ int i,ntriang。 int cubeindex。 Vec3 vertlist[12]。 int edgeTable[256]={ 0x0 , 0x109, 0x203, 0x30a, 0x406, 0x50f, 0x605, 0x70c, 0x80c, 0x905, 0xa0f, 0xb06, 0xc0a, 0xd03, 0xe09, 0xf00, 0x190, 0x99 , 0x393, 0x29a, 0x596, 0x49f, 0x795, 0x69c, 0x99c, 0x895, 0xb9f, 0xa96, 0xd9a, 0xc93, 0xf99, 0xe90, 0x230, 0x339, 0x33 , 0x13a, 0x636, 0x73f, 0x435, 0x53c, 0xa3c, 0xb35, 0x83f, 0x936, 0xe3a, 0xf33, 0xc39, 0xd30, 0x3a0, 0x2a9, 0x1a3, 0xaa , 0x7a6, 0x6af, 0x5a5, 0x4ac, 0xbac, 0xaa5, 0x9af, 0x8a6, 0xfaa, 0xea3, 0xda9, 0xca0, 0x460, 0x569, 0x663, 0x76a, 0x66 , 0x16f, 0x265, 0x36c, 0xc6c, 0xd65, 0xe6f, 0xf66, 0x86a, 0x963, 0xa69, 0xb60, 0x5f0, 0x4f9, 0x7f3, 0x6fa, 0x1f6, 0xff , 0x3f5, 0x2fc, 0xdfc, 0xcf5, 0xfff, 0xef6, 0x9fa, 0x8f3, 0xbf9, 0xaf0, 0x650, 0x759, 0x453, 0x55a, 0x256, 0x35f, 0x55 , 0x15c, 0xe5c, 0xf55, 0xc5f, 0xd56, 0xa5a, 0xb53, 0x859, 0x950, 0x7c0, 0x6c9, 0x5c3, 0x4ca, 0x3c6, 0x2cf, 0x1c5, 0xcc , 0xfcc, 0xec5, 0xdcf, 0xcc6, 0xbca, 0xac3, 0x9c9, 0x8c0, 0x8c0, 0x9c9, 0xac3, 0xbca, 0xcc6, 0xdcf, 0xec5, 0xfcc, 0xcc , 0x1c5, 0x2cf, 0x3c6, 0x4ca, 0x5c3, 0x6c9, 0x7c0, 0x950, 0x859, 0xb53, 0xa5a, 0xd56, 0xc5f, 0xf55, 0xe5c, 0x15c, 0x55 , 0x35f, 0x256, 0x55a, 0x453, 0x759, 0x650, 0xaf0, 0xbf9, 0x8f3, 0x9fa, 0xef6, 0xfff, 0xcf5, 0xdfc, 0x2fc, 0x3f5, 0xff , 0x1f6, 0x6fa, 0x7f3, 0x4f9, 0x5f0, 0xb60, 0xa69, 0x963, 0x86a, 0xf66, 0xe6f, 0xd65, 0xc6c, 0x36c, 0x265, 0x16f, 0x66 , 0x76a, 0x663, 0x569, 0x460, 0xca0, 0xda9, 0xea3, 0xfaa, 0x8a6, 0x9af, 0xaa5, 0xbac, 0x4ac, 0x5a5, 0x6af, 0x7a6, 0xaa , 0x1a3, 0x2a9, 0x3a0, 0xd30, 0xc39, 0xf33, 0xe3a, 0x936, 0x83f, 0xb35, 0xa3c, 0x53c, 0x435, 0x73f, 0x636, 0x13a, 0x33 , 0x339, 0x230, 0xe90, 0xf99, 0xc93, 0xd9a, 0xa96, 0xb9f, 0x895, 0x99c, 0x69c, 0x795, 0x49f, 0x596, 0x29a, 0x393, 0x99 , 0x190, 0xf00, 0xe09, 0xd03, 0xc0a, 0xb06, 0xa0f, 0x905, 0x80c, 0x70c, 0x605, 0x50f, 0x406, 0x30a, 0x203, 0x109, 0x0 }。int triTable[256][16] ={{1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1},{0, 8, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1},{0, 1, 9, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1},{1, 8, 3, 9, 8, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1},{1, 2, 10, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1},{0, 8, 3, 1, 2, 10, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1},{9, 2, 10, 0, 2, 9, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1},{2, 8, 3, 2, 10, 8, 10, 9, 8, 1, 1, 1, 1, 1, 1, 1},{3, 11, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1},{0, 11, 2, 8, 11, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1},{1, 9, 0, 2, 3, 11, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1},{1, 11, 2, 1, 9, 11, 9, 8, 11, 1, 1, 1, 1, 1, 1, 1},{3, 10, 1, 11, 10, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1},{0, 10, 1, 0, 8, 10, 8, 11, 10, 1, 1, 1, 1, 1, 1, 1},{3, 9, 0, 3, 11, 9, 11, 10, 9, 1, 1, 1, 1, 1, 1, 1},{9, 8, 10, 10, 8, 11, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1},{4, 7, 8, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1},{4, 3, 0, 7, 3, 4, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1},{0, 1, 9, 8, 4, 7, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1},{4, 1, 9, 4, 7, 1, 7, 3, 1, 1, 1, 1, 1, 1, 1, 1},{1, 2, 10, 8, 4, 7, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1},{3, 4, 7, 3, 0, 4, 1, 2, 10, 1, 1, 1, 1, 1, 1, 1},{9, 2, 10, 9, 0, 2, 8, 4, 7, 1, 1, 1, 1, 1, 1, 1},{2, 10, 9, 2, 9, 7, 2, 7, 3, 7, 9, 4, 1, 1, 1, 1},{8, 4, 7, 3, 11, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1},{11, 4, 7, 11, 2, 4, 2, 0, 4, 1, 1, 1, 1, 1, 1, 1},{9, 0, 1, 8, 4, 7, 2, 3, 11, 1, 1, 1,