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

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

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

revision 3250 by caltinay, Wed Oct 6 05:53:06 2010 UTC revision 3251 by jfenwick, Thu Oct 7 04:02:30 2010 UTC
# Line 84  void Dudley_Assemble_PDE_Single2_3D(Dudl Line 84  void Dudley_Assemble_PDE_Single2_3D(Dudl
84          {          {
85              if (elements->Color[e] == color)              if (elements->Color[e] == color)
86              {              {
87                double vol = p.row_jac->absD[e] * p.row_jac->quadweight;
88    
89              A_p = getSampleDataRO(A, e);              A_p = getSampleDataRO(A, e);
90              B_p = getSampleDataRO(B, e);              B_p = getSampleDataRO(B, e);
# Line 92  void Dudley_Assemble_PDE_Single2_3D(Dudl Line 93  void Dudley_Assemble_PDE_Single2_3D(Dudl
93              X_p = getSampleDataRO(X, e);              X_p = getSampleDataRO(X, e);
94              Y_p = getSampleDataRO(Y, e);              Y_p = getSampleDataRO(Y, e);
95    
             double vol = p.row_jac->absD[e] * p.row_jac->quadweight;  
96              DSDX = &(p.row_jac->DSDX[INDEX5(0, 0, 0, 0, e, p.numShapes, DIM, p.numQuad, 1)]);              DSDX = &(p.row_jac->DSDX[INDEX5(0, 0, 0, 0, e, p.numShapes, DIM, p.numQuad, 1)]);
97              for (q = 0; q < len_EM_S; ++q)              for (q = 0; q < len_EM_S; ++q)
98                  EM_S[q] = 0;                  EM_S[q] = 0;

Legend:
Removed from v.3250  
changed lines
  Added in v.3251

  ViewVC Help
Powered by ViewVC 1.1.26