tetramod.c

su 的源代码库

                              tetra[it5].ifacet[1]=if11+nxhz-1;
                              tetra[it5].ifacet[2]=if4;        
                              tetra[it5].ifacet[3]=if8;

                              tetra[it6].ifacet[0]=if8;  /*lower right tetra*/
                              tetra[it6].ifacet[1]=if10+1;
                              tetra[it6].ifacet[2]=if12+nxhz-1;        
                              tetra[it6].ifacet[3]=if2+ntrap;

			      /*centered normal for hz (c1): left hz facet*/
		 	      tm_normal( 
                                    xhz[ihz][iyhz][ixhz],
                                    xhz[ihz][iyhz][ixhz+1],
				    xhz[ihz][iyhz+1][ixhz],
			            yhz[ihz][iyhz][ixhz],
                                    yhz[ihz][iyhz][ixhz+1],
       			            yhz[ihz][iyhz+1][ixhz],
			            zhz[ihz][iyhz][ixhz],
                                    zhz[ihz][iyhz][ixhz+1],
       			      	    zhz[ihz][iyhz+1][ixhz],
				    &facet[if1].cn[0],
				    &facet[if1].cn[1],
				    &facet[if1].cn[2]);

			      /*centered normal for hz (c1): right hz facet*/
		 	      tm_normal( 
                                    xhz[ihz][iyhz+1][ixhz+1],
                                    xhz[ihz][iyhz][ixhz+1],
				    xhz[ihz][iyhz+1][ixhz],
			            yhz[ihz][iyhz+1][ixhz+1],
                                    yhz[ihz][iyhz][ixhz+1],
       			            yhz[ihz][iyhz+1][ixhz],
			            zhz[ihz][iyhz+1][ixhz+1],
                                    zhz[ihz][iyhz][ixhz+1],
       			      	    zhz[ihz][iyhz+1][ixhz],
				    &facet[if2].cn[0],
				    &facet[if2].cn[1],
				    &facet[if2].cn[2]);

			      /*centered normal for if3: upper diagonal facet*/
		 	      tm_normal( 
                                    xhz[ihz][iyhz+1][ixhz],
                                    xhz[ihz][iyhz][ixhz+1],
				    xhz[ihz0][iyhz][ixhz+1],
			            yhz[ihz][iyhz+1][ixhz],
                                    yhz[ihz][iyhz][ixhz+1],
       			            yhz[ihz0][iyhz][ixhz+1],
			            zhz[ihz][iyhz+1][ixhz],
                                    zhz[ihz][iyhz][ixhz+1],
       			      	    zhz[ihz+1][iyhz][ixhz+1],
				    &facet[if3].cn[0],
				    &facet[if3].cn[1],
				    &facet[if3].cn[2]);

			      /*centered normal for if4: lower diagonal facet*/
		 	      tm_normal( 
                                    xhz[ihz][iyhz+1][ixhz],
                                    xhz[ihz0][iyhz+1][ixhz],
				    xhz[ihz0][iyhz][ixhz+1],
			            yhz[ihz][iyhz+1][ixhz],
                                    yhz[ihz0][iyhz+1][ixhz],
       			            yhz[ihz0][iyhz][ixhz+1],
			            zhz[ihz][iyhz+1][ixhz],
                                    zhz[ihz+1][iyhz+1][ixhz],
       			      	    zhz[ihz+1][iyhz][ixhz+1],
				    &facet[if4].cn[0],
				    &facet[if4].cn[1],
				    &facet[if4].cn[2]);

			      /*centered normal for if5: left slant1 facet*/
		 	      tm_normal( 
                                    xhz[ihz][iyhz][ixhz],
                                    xhz[ihz][iyhz+1][ixhz],
				    xhz[ihz0][iyhz][ixhz+1],
			            yhz[ihz][iyhz][ixhz],
                                    yhz[ihz][iyhz+1][ixhz],
       			            yhz[ihz0][iyhz][ixhz+1],
			            zhz[ihz][iyhz][ixhz],
                                    zhz[ihz][iyhz+1][ixhz],
       			      	    zhz[ihz+1][iyhz][ixhz+1],
				    &facet[if5].cn[0],
				    &facet[if5].cn[1],
				    &facet[if5].cn[2]);

			      /*centered normal for if6: right slant1 facet*/
		 	      tm_normal( 
                                    xhz[ihz][iyhz+1][ixhz+1],
                                    xhz[ihz][iyhz+1][ixhz],
				    xhz[ihz0][iyhz][ixhz+1],
			            yhz[ihz][iyhz+1][ixhz+1],
                                    yhz[ihz][iyhz+1][ixhz],
       			            yhz[ihz0][iyhz][ixhz+1],
			            zhz[ihz][iyhz+1][ixhz+1],
                                    zhz[ihz][iyhz+1][ixhz],
       			      	    zhz[ihz+1][iyhz][ixhz+1],
				    &facet[if6].cn[0],
				    &facet[if6].cn[1],
				    &facet[if6].cn[2]);

			      /*centered normal for if7: left slant 2 facet*/
		 	      tm_normal( 
                                    xhz[ihz][iyhz+1][ixhz],
                                    xhz[ihz0][iyhz][ixhz],
				    xhz[ihz0][iyhz][ixhz+1],
			            yhz[ihz][iyhz+1][ixhz],
                                    yhz[ihz0][iyhz][ixhz],
       			            yhz[ihz0][iyhz][ixhz+1],
			            zhz[ihz][iyhz+1][ixhz],
                                    zhz[ihz+1][iyhz][ixhz],
       			      	    zhz[ihz+1][iyhz][ixhz+1],
				    &facet[if7].cn[0],
				    &facet[if7].cn[1],
				    &facet[if7].cn[2]);

			      /*centered normal for if8: right slant 2 facet*/
		 	      tm_normal( 
                                    xhz[ihz][iyhz+1][ixhz+1],
                                    xhz[ihz0][iyhz+1][ixhz],
				    xhz[ihz0][iyhz][ixhz+1],
			            yhz[ihz][iyhz+1][ixhz+1],
                                    yhz[ihz0][iyhz+1][ixhz],
       			            yhz[ihz0][iyhz][ixhz+1],
			            zhz[ihz][iyhz+1][ixhz+1],
                                    zhz[ihz+1][iyhz+1][ixhz],
       			      	    zhz[ihz+1][iyhz][ixhz+1],
				    &facet[if8].cn[0],
				    &facet[if8].cn[1],
				    &facet[if8].cn[2]);

                              /*control point of left facet on hz*/
			      facet[if1].ip[0]=i000;
			      facet[if1].ip[1]=i100;
			      facet[if1].ip[2]=i010;
                        
			      /*control point of right facet on hz*/
			      facet[if2].ip[0]=i010;
			      facet[if2].ip[1]=i110;
			      facet[if2].ip[2]=i100;

                              /*upper diagonal facet*/
			      facet[if3].ip[0]=i010;
			      facet[if3].ip[1]=i100;
			      facet[if3].ip[2]=i101;
                        
			      /*lower diagonal facet*/
			      facet[if4].ip[0]=i010;
			      facet[if4].ip[1]=i011;
			      facet[if4].ip[2]=i101;

			      /*left slant1 facet*/
			      facet[if5].ip[0]=i000;
			      facet[if5].ip[1]=i010;
			      facet[if5].ip[2]=i101;

			      /*right slant1 facet*/
			      facet[if6].ip[0]=i010;
			      facet[if6].ip[1]=i110;
			      facet[if6].ip[2]=i101;

			      /*left slant2 facet*/
			      facet[if7].ip[0]=i001;
			      facet[if7].ip[1]=i010;
			      facet[if7].ip[2]=i101;

			      /*right slant2 facet*/
			      facet[if8].ip[0]=i011;
			      facet[if8].ip[1]=i110;
			      facet[if8].ip[2]=i101;

                        } else {

			      /*left top tetra*/
                              tetra[it1].ip[0]=i000;
                              tetra[it1].ip[1]=i010;
                              tetra[it1].ip[2]=i110;
                              tetra[it1].ip[3]=i111;

		      	      /*left middle tetra*/
                              tetra[it2].ip[0]=i000;    
                              tetra[it2].ip[1]=i010;
                              tetra[it2].ip[2]=i011;
                              tetra[it2].ip[3]=i111;

			      /*left bottom tetra*/
                              tetra[it3].ip[0]=i001;
                              tetra[it3].ip[1]=i011;
                              tetra[it3].ip[2]=i111;
                              tetra[it3].ip[3]=i000;

			      /*right top tetra*/
                              tetra[it4].ip[0]=i000;
                              tetra[it4].ip[1]=i100;
                              tetra[it4].ip[2]=i110;
                              tetra[it4].ip[3]=i111;

			      /*right middle tetra*/
                              tetra[it5].ip[0]=i000;
                              tetra[it5].ip[1]=i001;
                              tetra[it5].ip[2]=i111;
                              tetra[it5].ip[3]=i100;

			      /*right bottom tetra*/
                              tetra[it6].ip[0]=i100;
                              tetra[it6].ip[1]=i001;
                              tetra[it6].ip[2]=i111;
                              tetra[it6].ip[3]=i101;

			      /*upper left tetra*/
                              tetra[it1].ifacet[0]=if1;
                              tetra[it1].ifacet[1]=if3;
                              tetra[it1].ifacet[2]=if11+nxhz-1;        
                              tetra[it1].ifacet[3]=if5;

			      /*middle left tetra*/
                              tetra[it2].ifacet[0]=if5;
                              tetra[it2].ifacet[1]=if9;
                              tetra[it2].ifacet[2]=if12+nxhz-1;        
                              tetra[it2].ifacet[3]=if7;

			      /*lower left tetra*/
                              tetra[it3].ifacet[0]=if7;
                              tetra[it3].ifacet[1]=if10;
                              tetra[it3].ifacet[2]=if4;        
                              tetra[it3].ifacet[3]=if1+ntrap;

			      /*upper right tetra*/
                              tetra[it4].ifacet[0]=if2;
                              tetra[it4].ifacet[1]=if3;
                              tetra[it4].ifacet[2]=if9+1;        
                              tetra[it4].ifacet[3]=if6;

			      /*middle right tetra*/
                              tetra[it5].ifacet[0]=if6;
                              tetra[it5].ifacet[1]=if4;
                              tetra[it5].ifacet[2]=if11;        
                              tetra[it5].ifacet[3]=if8;

			      /*lower right tetra*/
                              tetra[it6].ifacet[0]=if8;
                              tetra[it6].ifacet[1]=if10+1;
                              tetra[it6].ifacet[2]=if12;        
                              tetra[it6].ifacet[3]=if2+ntrap;

			      /*centered normal for hz (c1): left hz facet*/
		 	      tm_normal( 
                                    xhz[ihz][iyhz][ixhz],
                                    xhz[ihz][iyhz+1][ixhz+1],
				    xhz[ihz][iyhz+1][ixhz],
			            yhz[ihz][iyhz][ixhz],
                                    yhz[ihz][iyhz+1][ixhz+1],
       			            yhz[ihz][iyhz+1][ixhz],
			            zhz[ihz][iyhz][ixhz],
                                    zhz[ihz][iyhz+1][ixhz+1],
       			      	    zhz[ihz][iyhz+1][ixhz],
				    &facet[if1].cn[0],
				    &facet[if1].cn[1],
				    &facet[if1].cn[2]);

			      /*centered normal for hz (c1): right hz facet*/
		 	      tm_normal( 
                                    xhz[ihz][iyhz+1][ixhz+1],
                                    xhz[ihz][iyhz][ixhz+1],
				    xhz[ihz][iyhz][ixhz],
			            yhz[ihz][iyhz+1][ixhz+1],
                                    yhz[ihz][iyhz][ixhz+1],
       			            yhz[ihz][iyhz][ixhz],
			            zhz[ihz][iyhz+1][ixhz+1],
                                    zhz[ihz][iyhz][ixhz+1],
       			      	    zhz[ihz][iyhz][ixhz],
				    &facet[if2].cn[0],
				    &facet[if2].cn[1],
				    &facet[if2].cn[2]);

			      /*centered normal for if3: upper diagonal facet*/
		 	      tm_normal( 
                                    xhz[ihz][iyhz][ixhz],
                                    xhz[ihz][iyhz+1][ixhz+1],
				    xhz[ihz0][iyhz+1][ixhz+1],
			            yhz[ihz][iyhz][ixhz],
                                    yhz[ihz][iyhz+1][ixhz+1],
       			            yhz[ihz0][iyhz+1][ixhz+1],
			            zhz[ihz][iyhz][ixhz],
                                    zhz[ihz][iyhz+1][ixhz+1],
       			      	    zhz[ihz+1][iyhz+1][ixhz+1],
				    &facet[if3].cn[0],
				    &facet[if3].cn[1],
				    &facet[if3].cn[2]);

			      /*centered normal for if4: lower diagonal facet*/
		 	      tm_normal( 
                                    xhz[ihz][iyhz][ixhz],
                                    xhz[ihz0][iyhz][ixhz],
				    xhz[ihz0][iyhz+1][ixhz+1],
			            yhz[ihz][iyhz][ixhz],
                                    yhz[ihz0][iyhz][ixhz],
       			            yhz[ihz0][iyhz+1][ixhz+1],
			            zhz[ihz][iyhz][ixhz],
                                    zhz[ihz+1][iyhz][ixhz],
       			      	    zhz[ihz+1][iyhz+1][ixhz+1],
				    &facet[if4].cn[0],
				    &facet[if4].cn[1],
				    &facet[if4].cn[2]);

			      /*centered normal for if5: left slant1 facet*/
		 	      tm_normal( 
                                    xhz[ihz][iyhz][ixhz],
                                    xhz[ihz][iyhz+1][ixhz],
				    xhz[ihz0][iyhz+1][ixhz+1],
			            yhz[ihz][iyhz][ixhz],
                                    yhz[ihz][iyhz+1][ixhz],
       			            yhz[ihz0][iyhz+1][ixhz+1],
			            zhz[ihz][iyhz][ixhz],
                                    zhz[ihz][iyhz+1][ixhz],
       			      	    zhz[ihz+1][iyhz+1][ixhz+1],
				    &facet[if5].cn[0],
				    &facet[if5].cn[1],
				    &facet[if5].cn[2]);

			      /*centered normal for if6: right slant1 facet*/
		 	      tm_normal( 
                                    xhz[ihz][iyhz][ixhz],
                                    xhz[ihz][iyhz][ixhz+1],
				    xhz[ihz0][iyhz+1][ixhz+1],
			            yhz[ihz][iyhz][ixhz],
                                    yhz[ihz][iyhz][ixhz+1],
       			            yhz[ihz0][iyhz+1][ixhz+1],
			            zhz[ihz][iyhz][ixhz],
                                    zhz[ihz][iyhz][ixhz+1],
       			      	    zhz[ihz+1][iyhz+1][ixhz+1],
				    &facet[if6].cn[0],
				    &facet[if6].cn[1],
				    &facet[if6].cn[2]);

			      /*centered normal for if7: left slant 2 facet*/
		 	      tm_normal(