Salome HOME
Debug on GENERAL_24 and GENERAL_48
authorageay <ageay>
Mon, 5 Aug 2013 06:44:52 +0000 (06:44 +0000)
committerageay <ageay>
Mon, 5 Aug 2013 06:44:52 +0000 (06:44 +0000)
src/INTERP_KERNEL/SplitterTetra.cxx

index b5e4999af0cedbfe9162f212a56d979eeb64472a..ddc4dbf27b0be18b25f496959a67f2d3579d7af3 100644 (file)
@@ -66,10 +66,11 @@ namespace INTERP_KERNEL
               tmp2[0]=0.; tmp2[1]=0.; tmp2[2]=0.;
               for(int j=0;j<4;j++,conn+=4)
                 {
-                  tmp2[0]+=coords[3*nodalConnBg[GENERAL_24_SUB_NODES_WO[4*i+j]]+0];
-                  tmp2[1]+=coords[3*nodalConnBg[GENERAL_24_SUB_NODES_WO[4*i+j]]+1];
-                  tmp2[2]+=coords[3*nodalConnBg[GENERAL_24_SUB_NODES_WO[4*i+j]]+3];
-                  conn[0]=nodalConnBg[GENERAL_24_SUB_NODES_WO[4*i+j]];
+                  int tmp3(nodalConnBg[GENERAL_24_SUB_NODES_WO[4*i+j]]);
+                  tmp2[0]+=coords[3*tmp3+0];
+                  tmp2[1]+=coords[3*tmp3+1];
+                  tmp2[2]+=coords[3*tmp3+2];
+                  conn[0]=tmp3;
                   conn[1]=nodalConnBg[GENERAL_24_SUB_NODES_WO[4*i+(j+1)%4]];
                   conn[2]=-(i+1); conn[3]=-7;
                 }
@@ -92,9 +93,9 @@ namespace INTERP_KERNEL
             }
           for(int i=0;i<7;i++,tmp2+=3)
             {
-              tmp2[0]=(tmp[3*nodalConnBg[(GENERAL_48_SUB_NODES[2*i+24]-8)]+0]+tmp[3*nodalConnBg[(GENERAL_48_SUB_NODES[2*i+25]-8)]+0])/2.;
-              tmp2[1]=(tmp[3*nodalConnBg[(GENERAL_48_SUB_NODES[2*i+24]-8)]+1]+tmp[3*nodalConnBg[(GENERAL_48_SUB_NODES[2*i+25]-8)]+1])/2.;
-              tmp2[2]=(tmp[3*nodalConnBg[(GENERAL_48_SUB_NODES[2*i+24]-8)]+2]+tmp[3*nodalConnBg[(GENERAL_48_SUB_NODES[2*i+25]-8)]+2])/2.;
+              tmp2[0]=(tmp[3*(GENERAL_48_SUB_NODES[2*i+24]-8)+0]+tmp[3*(GENERAL_48_SUB_NODES[2*i+25]-8)+0])/2.;
+              tmp2[1]=(tmp[3*(GENERAL_48_SUB_NODES[2*i+24]-8)+1]+tmp[3*(GENERAL_48_SUB_NODES[2*i+25]-8)+1])/2.;
+              tmp2[2]=(tmp[3*(GENERAL_48_SUB_NODES[2*i+24]-8)+2]+tmp[3*(GENERAL_48_SUB_NODES[2*i+25]-8)+2])/2.;
             }
           int *conn(&tetrasNodalConn[0]);
           std::vector<double> dummy;