/[escript]/branches/domexper/dudley/src/Assemble_jacobeans.c
ViewVC logotype

Diff of /branches/domexper/dudley/src/Assemble_jacobeans.c

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 3220 by jfenwick, Thu Sep 16 05:20:34 2010 UTC revision 3221 by jfenwick, Wed Sep 29 01:00:21 2010 UTC
# Line 105  coordinates[INDEX2(P,nodes[INDEX2(2,e,nu Line 105  coordinates[INDEX2(P,nodes[INDEX2(2,e,nu
105              DTDXSET(0,0);              DTDXSET(0,0);
106              DTDXSET(1,0);              DTDXSET(1,0);
107              }              }
 //          volume[INDEX2(0,e,1)]=ABS(D)*(*quadweight);  
108          }          }
109          else    // numQuad==3          else    // numQuad==3
110          {          {
# Line 116  coordinates[INDEX2(P,nodes[INDEX2(2,e,nu Line 115  coordinates[INDEX2(P,nodes[INDEX2(2,e,nu
115                  DTDXSET(0,q);                  DTDXSET(0,q);
116                  DTDXSET(1,q);                  DTDXSET(1,q);
117              }              }
 //          volume[INDEX2(q,e,3)]=ABS(D)*(*quadweight);  
118              }              }
119          }          }
120          }          }
# Line 167  void Assemble_jacobeans_2D_M1D_E1D(doubl Line 165  void Assemble_jacobeans_2D_M1D_E1D(doubl
165          dTdX[INDEX4(0,1,0,e,numTest,DIM,numQuad)]=-1*dvdX01;          dTdX[INDEX4(0,1,0,e,numTest,DIM,numQuad)]=-1*dvdX01;
166          dTdX[INDEX4(1,0,0,e,numTest,DIM,numQuad)]=-1*dvdX00;          dTdX[INDEX4(1,0,0,e,numTest,DIM,numQuad)]=-1*dvdX00;
167          dTdX[INDEX4(1,1,0,e,numTest,DIM,numQuad)]=-1*dvdX01;          dTdX[INDEX4(1,1,0,e,numTest,DIM,numQuad)]=-1*dvdX01;
 //      volume[INDEX2(0,e,numQuad)]=sqrt(D)*(*quadweight);  
168          absD[e]=sqrt(D);          absD[e]=sqrt(D);
169          if (numQuad==2)          if (numQuad==2)
170          {          {
# Line 175  void Assemble_jacobeans_2D_M1D_E1D(doubl Line 172  void Assemble_jacobeans_2D_M1D_E1D(doubl
172              dTdX[INDEX4(0,1,1,e,numTest,DIM,numQuad)]=dvdX01;              dTdX[INDEX4(0,1,1,e,numTest,DIM,numQuad)]=dvdX01;
173              dTdX[INDEX4(1,0,1,e,numTest,DIM,numQuad)]=dvdX00;              dTdX[INDEX4(1,0,1,e,numTest,DIM,numQuad)]=dvdX00;
174              dTdX[INDEX4(1,1,1,e,numTest,DIM,numQuad)]=dvdX01;              dTdX[INDEX4(1,1,1,e,numTest,DIM,numQuad)]=dvdX01;
 //          volume[INDEX2(1,e,numQuad)]=sqrt(D)*(*quadweight);  
175          }          }
176          }          }
177      }      }
# Line 256  void Assemble_jacobeans_3D(double* coord Line 252  void Assemble_jacobeans_3D(double* coord
252              dTdX[INDEX4(s,1,q,e,numTest,DIM,numQuad)]=DTDV_3D[s][0]*dvdX01+DTDV_3D[s][1]*dvdX11+DTDV_3D[s][2]*dvdX21;              dTdX[INDEX4(s,1,q,e,numTest,DIM,numQuad)]=DTDV_3D[s][0]*dvdX01+DTDV_3D[s][1]*dvdX11+DTDV_3D[s][2]*dvdX21;
253              dTdX[INDEX4(s,2,q,e,numTest,DIM,numQuad)]=DTDV_3D[s][0]*dvdX02+DTDV_3D[s][1]*dvdX12+DTDV_3D[s][2]*dvdX22;              dTdX[INDEX4(s,2,q,e,numTest,DIM,numQuad)]=DTDV_3D[s][0]*dvdX02+DTDV_3D[s][1]*dvdX12+DTDV_3D[s][2]*dvdX22;
254              }              }
 //          volume[INDEX2(q,e,numQuad)]=ABS(D)*(*quadweight);  
255                  }                  }
256          }          }
257      }      }
# Line 331  void Assemble_jacobeans_3D_M2D_E2D(doubl Line 326  void Assemble_jacobeans_3D_M2D_E2D(doubl
326              dTdX[INDEX4(s,1,q,e,numTest,DIM,numQuad)]=DTDV[s][0]*dvdX01+DTDV[s][1]*dvdX11;              dTdX[INDEX4(s,1,q,e,numTest,DIM,numQuad)]=DTDV[s][0]*dvdX01+DTDV[s][1]*dvdX11;
327              dTdX[INDEX4(s,2,q,e,numTest,DIM,numQuad)]=DTDV[s][0]*dvdX02+DTDV[s][1]*dvdX12;              dTdX[INDEX4(s,2,q,e,numTest,DIM,numQuad)]=DTDV[s][0]*dvdX02+DTDV[s][1]*dvdX12;
328              }              }
 //          volume[INDEX2(q,e,numQuad)]=sqrt(D)*(*quadweight);  
329          }          }
330          }          }
331      }      }

Legend:
Removed from v.3220  
changed lines
  Added in v.3221

  ViewVC Help
Powered by ViewVC 1.1.26