/[escript]/branches/doubleplusgood/dudley/src/Assemble_PDE_System2_2D.cpp
ViewVC logotype

Diff of /branches/doubleplusgood/dudley/src/Assemble_PDE_System2_2D.cpp

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

revision 3247 by caltinay, Wed Oct 6 05:53:06 2010 UTC revision 3251 by jfenwick, Thu Oct 7 04:02:30 2010 UTC
# Line 85  void Dudley_Assemble_PDE_System2_2D(Dudl Line 85  void Dudley_Assemble_PDE_System2_2D(Dudl
85          {          {
86              if (elements->Color[e] == color)              if (elements->Color[e] == color)
87              {              {
88                double vol = p.row_jac->absD[e] * p.row_jac->quadweight;
89    
90              A_p = getSampleDataRO(A, e);              A_p = getSampleDataRO(A, e);
91              B_p = getSampleDataRO(B, e);              B_p = getSampleDataRO(B, e);
# Line 92  void Dudley_Assemble_PDE_System2_2D(Dudl Line 93  void Dudley_Assemble_PDE_System2_2D(Dudl
93              D_p = getSampleDataRO(D, e);              D_p = getSampleDataRO(D, e);
94              X_p = getSampleDataRO(X, e);              X_p = getSampleDataRO(X, e);
95              Y_p = getSampleDataRO(Y, e);              Y_p = getSampleDataRO(Y, e);
             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.3247  
changed lines
  Added in v.3251

  ViewVC Help
Powered by ViewVC 1.1.26