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

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

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

revision 3221 by jfenwick, Wed Sep 29 01:00:21 2010 UTC revision 3224 by jfenwick, Wed Sep 29 05:19:37 2010 UTC
# Line 124  void Dudley_Assemble_PDE_Single2_1D(Asse Line 124  void Dudley_Assemble_PDE_Single2_1D(Asse
124                      EM_S[INDEX4(0, 0, s, r, p.numEqu, p.numComp, p.row_numShapesTotal)] += rtmp;                      EM_S[INDEX4(0, 0, s, r, p.numEqu, p.numComp, p.row_numShapesTotal)] += rtmp;
125                      }                      }
126                  }                  }
127                  } else                  }
128                    else
129                  {                  {
130                  for (s = 0; s < p.row_numShapes; s++)                  for (s = 0; s < p.row_numShapes; s++)
131                  {                  {
# Line 164  void Dudley_Assemble_PDE_Single2_1D(Asse Line 165  void Dudley_Assemble_PDE_Single2_1D(Asse
165                      EM_S[INDEX4(0, 0, s, r, p.numEqu, p.numComp, p.row_numShapesTotal)] += rtmp;                      EM_S[INDEX4(0, 0, s, r, p.numEqu, p.numComp, p.row_numShapesTotal)] += rtmp;
166                      }                      }
167                  }                  }
168                  } else                  }
169                    else
170                  {                  {
171                  for (s = 0; s < p.row_numShapes; s++)                  for (s = 0; s < p.row_numShapes; s++)
172                  {                  {
# Line 204  void Dudley_Assemble_PDE_Single2_1D(Asse Line 206  void Dudley_Assemble_PDE_Single2_1D(Asse
206                      EM_S[INDEX4(0, 0, s, r, p.numEqu, p.numComp, p.row_numShapesTotal)] += rtmp;                      EM_S[INDEX4(0, 0, s, r, p.numEqu, p.numComp, p.row_numShapesTotal)] += rtmp;
207                      }                      }
208                  }                  }
209                  } else                  }
210                    else
211                  {                  {
212                  for (s = 0; s < p.row_numShapes; s++)                  for (s = 0; s < p.row_numShapes; s++)
213                  {                  {
# Line 244  void Dudley_Assemble_PDE_Single2_1D(Asse Line 247  void Dudley_Assemble_PDE_Single2_1D(Asse
247                      EM_S[INDEX4(0, 0, s, r, p.numEqu, p.numComp, p.row_numShapesTotal)] += rtmp;                      EM_S[INDEX4(0, 0, s, r, p.numEqu, p.numComp, p.row_numShapesTotal)] += rtmp;
248                      }                      }
249                  }                  }
250                  } else                  }
251                    else
252                  {                  {
253                  for (s = 0; s < p.row_numShapes; s++)                  for (s = 0; s < p.row_numShapes; s++)
254                  {                  {
# Line 279  void Dudley_Assemble_PDE_Single2_1D(Asse Line 283  void Dudley_Assemble_PDE_Single2_1D(Asse
283                          X_q[INDEX2(0, q, DIM)];                          X_q[INDEX2(0, q, DIM)];
284                      EM_F[INDEX2(0, s, p.numEqu)] += rtmp;                      EM_F[INDEX2(0, s, p.numEqu)] += rtmp;
285                  }                  }
286                  } else                  }
287                    else
288                  {                  {
289                  for (s = 0; s < p.row_numShapes; s++)                  for (s = 0; s < p.row_numShapes; s++)
290                  {                  {
# Line 306  void Dudley_Assemble_PDE_Single2_1D(Asse Line 311  void Dudley_Assemble_PDE_Single2_1D(Asse
311                      rtmp += vol * S[INDEX2(s, q, p.row_numShapes)] * Y_q[q];                      rtmp += vol * S[INDEX2(s, q, p.row_numShapes)] * Y_q[q];
312                      EM_F[INDEX2(0, s, p.numEqu)] += rtmp;                      EM_F[INDEX2(0, s, p.numEqu)] += rtmp;
313                  }                  }
314                  } else                  }
315                    else
316                  {                  {
317                  for (s = 0; s < p.row_numShapes; s++)                  for (s = 0; s < p.row_numShapes; s++)
318                  {                  {

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

  ViewVC Help
Powered by ViewVC 1.1.26