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

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

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

revision 3204 by jfenwick, Thu Sep 23 23:59:39 2010 UTC revision 3205 by jfenwick, Fri Sep 24 00:30:43 2010 UTC
# Line 94  void Dudley_Assemble_PDE_System2_2D(Asse Line 94  void Dudley_Assemble_PDE_System2_2D(Asse
94              X_p = getSampleDataRO(X, e);              X_p = getSampleDataRO(X, e);
95              Y_p = getSampleDataRO(Y, e);              Y_p = getSampleDataRO(Y, e);
96              double vol = p.row_jac->absD[e] * p.row_jac->quadweight;              double vol = p.row_jac->absD[e] * p.row_jac->quadweight;
97              DSDX = &(p.row_jac->DSDX[INDEX5(0, 0, 0, 0, e, p.numShapes, DIM, p.numQuadTotal, 1)]);              DSDX = &(p.row_jac->DSDX[INDEX5(0, 0, 0, 0, e, p.numShapes, DIM, p.numQuad, 1)]);
98              for (q = 0; q < len_EM_S; ++q)              for (q = 0; q < len_EM_S; ++q)
99                  EM_S[q] = 0;                  EM_S[q] = 0;
100              for (q = 0; q < len_EM_F; ++q)              for (q = 0; q < len_EM_F; ++q)
# Line 111  void Dudley_Assemble_PDE_System2_2D(Asse Line 111  void Dudley_Assemble_PDE_System2_2D(Asse
111                  add_EM_S = TRUE;                  add_EM_S = TRUE;
112                  if (extendedA)                  if (extendedA)
113                  {                  {
114                  A_q = &(A_p[INDEX6(0, 0, 0, 0, 0, 0, p.numEqu, DIM, p.numComp, DIM, p.numQuadTotal)]);                  A_q = &(A_p[INDEX6(0, 0, 0, 0, 0, 0, p.numEqu, DIM, p.numComp, DIM, p.numQuad)]);
115                  for (s = 0; s < p.numShapes; s++)                  for (s = 0; s < p.numShapes; s++)
116                  {                  {
117                      for (r = 0; r < p.numShapes; r++)                      for (r = 0; r < p.numShapes; r++)
# Line 121  void Dudley_Assemble_PDE_System2_2D(Asse Line 121  void Dudley_Assemble_PDE_System2_2D(Asse
121                          for (m = 0; m < p.numComp; m++)                          for (m = 0; m < p.numComp; m++)
122                          {                          {
123                          rtmp = 0;                          rtmp = 0;
124                          for (q = 0; q < p.numQuadTotal; q++)                          for (q = 0; q < p.numQuad; q++)
125                          {                          {
126                              rtmp +=                              rtmp +=
127                              vol * (DSDX[INDEX3(s, 0, q, p.numShapes, DIM)] *                              vol * (DSDX[INDEX3(s, 0, q, p.numShapes, DIM)] *
# Line 153  void Dudley_Assemble_PDE_System2_2D(Asse Line 153  void Dudley_Assemble_PDE_System2_2D(Asse
153                      rtmp01 = 0;                      rtmp01 = 0;
154                      rtmp10 = 0;                      rtmp10 = 0;
155                      rtmp11 = 0;                      rtmp11 = 0;
156                      for (q = 0; q < p.numQuadTotal; q++)                      for (q = 0; q < p.numQuad; q++)
157                      {                      {
158                          rtmp0 = vol * DSDX[INDEX3(s, 0, q, p.numShapes, DIM)];                          rtmp0 = vol * DSDX[INDEX3(s, 0, q, p.numShapes, DIM)];
159                          rtmp1 = vol * DSDX[INDEX3(s, 1, q, p.numShapes, DIM)];                          rtmp1 = vol * DSDX[INDEX3(s, 1, q, p.numShapes, DIM)];
# Line 186  void Dudley_Assemble_PDE_System2_2D(Asse Line 186  void Dudley_Assemble_PDE_System2_2D(Asse
186                  add_EM_S = TRUE;                  add_EM_S = TRUE;
187                  if (extendedB)                  if (extendedB)
188                  {                  {
189                  B_q = &(B_p[INDEX5(0, 0, 0, 0, 0, p.numEqu, DIM, p.numComp, p.numQuadTotal)]);                  B_q = &(B_p[INDEX5(0, 0, 0, 0, 0, p.numEqu, DIM, p.numComp, p.numQuad)]);
190                  for (s = 0; s < p.numShapes; s++)                  for (s = 0; s < p.numShapes; s++)
191                  {                  {
192                      for (r = 0; r < p.numShapes; r++)                      for (r = 0; r < p.numShapes; r++)
# Line 196  void Dudley_Assemble_PDE_System2_2D(Asse Line 196  void Dudley_Assemble_PDE_System2_2D(Asse
196                          for (m = 0; m < p.numComp; m++)                          for (m = 0; m < p.numComp; m++)
197                          {                          {
198                          rtmp = 0;                          rtmp = 0;
199                          for (q = 0; q < p.numQuadTotal; q++)                          for (q = 0; q < p.numQuad; q++)
200                          {                          {
201                              rtmp += vol * S[INDEX2(r, q, p.numShapes)] *                              rtmp += vol * S[INDEX2(r, q, p.numShapes)] *
202                              (DSDX[INDEX3(s, 0, q, p.numShapes, DIM)] *                              (DSDX[INDEX3(s, 0, q, p.numShapes, DIM)] *
# Line 218  void Dudley_Assemble_PDE_System2_2D(Asse Line 218  void Dudley_Assemble_PDE_System2_2D(Asse
218                      {                      {
219                      rtmp0 = 0;                      rtmp0 = 0;
220                      rtmp1 = 0;                      rtmp1 = 0;
221                      for (q = 0; q < p.numQuadTotal; q++)                      for (q = 0; q < p.numQuad; q++)
222                      {                      {
223                          rtmp = vol * S[INDEX2(r, q, p.numShapes)];                          rtmp = vol * S[INDEX2(r, q, p.numShapes)];
224                          rtmp0 += rtmp * DSDX[INDEX3(s, 0, q, p.numShapes, DIM)];                          rtmp0 += rtmp * DSDX[INDEX3(s, 0, q, p.numShapes, DIM)];
# Line 246  void Dudley_Assemble_PDE_System2_2D(Asse Line 246  void Dudley_Assemble_PDE_System2_2D(Asse
246                  add_EM_S = TRUE;                  add_EM_S = TRUE;
247                  if (extendedC)                  if (extendedC)
248                  {                  {
249                  C_q = &(C_p[INDEX5(0, 0, 0, 0, 0, p.numEqu, p.numComp, DIM, p.numQuadTotal)]);                  C_q = &(C_p[INDEX5(0, 0, 0, 0, 0, p.numEqu, p.numComp, DIM, p.numQuad)]);
250                  for (s = 0; s < p.numShapes; s++)                  for (s = 0; s < p.numShapes; s++)
251                  {                  {
252                      for (r = 0; r < p.numShapes; r++)                      for (r = 0; r < p.numShapes; r++)
# Line 256  void Dudley_Assemble_PDE_System2_2D(Asse Line 256  void Dudley_Assemble_PDE_System2_2D(Asse
256                          for (m = 0; m < p.numComp; m++)                          for (m = 0; m < p.numComp; m++)
257                          {                          {
258                          rtmp = 0;                          rtmp = 0;
259                          for (q = 0; q < p.numQuadTotal; q++)                          for (q = 0; q < p.numQuad; q++)
260                          {                          {
261                              rtmp += vol * S[INDEX2(s, q, p.numShapes)] *                              rtmp += vol * S[INDEX2(s, q, p.numShapes)] *
262                              (C_q[INDEX4(k, m, 0, q, p.numEqu, p.numComp, DIM)] *                              (C_q[INDEX4(k, m, 0, q, p.numEqu, p.numComp, DIM)] *
# Line 278  void Dudley_Assemble_PDE_System2_2D(Asse Line 278  void Dudley_Assemble_PDE_System2_2D(Asse
278                      {                      {
279                      rtmp0 = 0;                      rtmp0 = 0;
280                      rtmp1 = 0;                      rtmp1 = 0;
281                      for (q = 0; q < p.numQuadTotal; q++)                      for (q = 0; q < p.numQuad; q++)
282                      {                      {
283                          rtmp = vol * S[INDEX2(s, q, p.numShapes)];                          rtmp = vol * S[INDEX2(s, q, p.numShapes)];
284                          rtmp0 += rtmp * DSDX[INDEX3(r, 0, q, p.numShapes, DIM)];                          rtmp0 += rtmp * DSDX[INDEX3(r, 0, q, p.numShapes, DIM)];
# Line 306  void Dudley_Assemble_PDE_System2_2D(Asse Line 306  void Dudley_Assemble_PDE_System2_2D(Asse
306                  add_EM_S = TRUE;                  add_EM_S = TRUE;
307                  if (extendedD)                  if (extendedD)
308                  {                  {
309                  D_q = &(D_p[INDEX4(0, 0, 0, 0, p.numEqu, p.numComp, p.numQuadTotal)]);                  D_q = &(D_p[INDEX4(0, 0, 0, 0, p.numEqu, p.numComp, p.numQuad)]);
310                  for (s = 0; s < p.numShapes; s++)                  for (s = 0; s < p.numShapes; s++)
311                  {                  {
312                      for (r = 0; r < p.numShapes; r++)                      for (r = 0; r < p.numShapes; r++)
# Line 316  void Dudley_Assemble_PDE_System2_2D(Asse Line 316  void Dudley_Assemble_PDE_System2_2D(Asse
316                          for (m = 0; m < p.numComp; m++)                          for (m = 0; m < p.numComp; m++)
317                          {                          {
318                          rtmp = 0;                          rtmp = 0;
319                          for (q = 0; q < p.numQuadTotal; q++)                          for (q = 0; q < p.numQuad; q++)
320                          {                          {
321                              rtmp +=                              rtmp +=
322                              vol * S[INDEX2(s, q, p.numShapes)] *                              vol * S[INDEX2(s, q, p.numShapes)] *
# Line 336  void Dudley_Assemble_PDE_System2_2D(Asse Line 336  void Dudley_Assemble_PDE_System2_2D(Asse
336                      for (r = 0; r < p.numShapes; r++)                      for (r = 0; r < p.numShapes; r++)
337                      {                      {
338                      rtmp = 0;                      rtmp = 0;
339                      for (q = 0; q < p.numQuadTotal; q++)                      for (q = 0; q < p.numQuad; q++)
340                          rtmp +=                          rtmp +=
341                          vol * S[INDEX2(s, q, p.numShapes)] *                          vol * S[INDEX2(s, q, p.numShapes)] *
342                          S[INDEX2(r, q, p.numShapes)];                          S[INDEX2(r, q, p.numShapes)];
# Line 361  void Dudley_Assemble_PDE_System2_2D(Asse Line 361  void Dudley_Assemble_PDE_System2_2D(Asse
361                  add_EM_F = TRUE;                  add_EM_F = TRUE;
362                  if (extendedX)                  if (extendedX)
363                  {                  {
364                  X_q = &(X_p[INDEX4(0, 0, 0, 0, p.numEqu, DIM, p.numQuadTotal)]);                  X_q = &(X_p[INDEX4(0, 0, 0, 0, p.numEqu, DIM, p.numQuad)]);
365                  for (s = 0; s < p.numShapes; s++)                  for (s = 0; s < p.numShapes; s++)
366                  {                  {
367                      for (k = 0; k < p.numEqu; k++)                      for (k = 0; k < p.numEqu; k++)
368                      {                      {
369                      rtmp = 0;                      rtmp = 0;
370                      for (q = 0; q < p.numQuadTotal; q++)                      for (q = 0; q < p.numQuad; q++)
371                      {                      {
372                          rtmp +=                          rtmp +=
373                          vol * (DSDX[INDEX3(s, 0, q, p.numShapes, DIM)] *                          vol * (DSDX[INDEX3(s, 0, q, p.numShapes, DIM)] *
# Line 384  void Dudley_Assemble_PDE_System2_2D(Asse Line 384  void Dudley_Assemble_PDE_System2_2D(Asse
384                  {                  {
385                      rtmp0 = 0;                      rtmp0 = 0;
386                      rtmp1 = 0;                      rtmp1 = 0;
387                      for (q = 0; q < p.numQuadTotal; q++)                      for (q = 0; q < p.numQuad; q++)
388                      {                      {
389                      rtmp0 += vol * DSDX[INDEX3(s, 0, q, p.numShapes, DIM)];                      rtmp0 += vol * DSDX[INDEX3(s, 0, q, p.numShapes, DIM)];
390                      rtmp1 += vol * DSDX[INDEX3(s, 1, q, p.numShapes, DIM)];                      rtmp1 += vol * DSDX[INDEX3(s, 1, q, p.numShapes, DIM)];
# Line 404  void Dudley_Assemble_PDE_System2_2D(Asse Line 404  void Dudley_Assemble_PDE_System2_2D(Asse
404                  add_EM_F = TRUE;                  add_EM_F = TRUE;
405                  if (extendedY)                  if (extendedY)
406                  {                  {
407                  Y_q = &(Y_p[INDEX3(0, 0, 0, p.numEqu, p.numQuadTotal)]);                  Y_q = &(Y_p[INDEX3(0, 0, 0, p.numEqu, p.numQuad)]);
408                  for (s = 0; s < p.numShapes; s++)                  for (s = 0; s < p.numShapes; s++)
409                  {                  {
410                      for (k = 0; k < p.numEqu; k++)                      for (k = 0; k < p.numEqu; k++)
411                      {                      {
412                      rtmp = 0;                      rtmp = 0;
413                      for (q = 0; q < p.numQuadTotal; q++)                      for (q = 0; q < p.numQuad; q++)
414                          rtmp +=                          rtmp +=
415                          vol * S[INDEX2(s, q, p.numShapes)] * Y_q[INDEX2(k, q, p.numEqu)];                          vol * S[INDEX2(s, q, p.numShapes)] * Y_q[INDEX2(k, q, p.numEqu)];
416                      EM_F[INDEX2(k, s, p.numEqu)] += rtmp;                      EM_F[INDEX2(k, s, p.numEqu)] += rtmp;
# Line 421  void Dudley_Assemble_PDE_System2_2D(Asse Line 421  void Dudley_Assemble_PDE_System2_2D(Asse
421                  for (s = 0; s < p.numShapes; s++)                  for (s = 0; s < p.numShapes; s++)
422                  {                  {
423                      rtmp = 0;                      rtmp = 0;
424                      for (q = 0; q < p.numQuadTotal; q++)                      for (q = 0; q < p.numQuad; q++)
425                      rtmp += vol * S[INDEX2(s, q, p.numShapes)];                      rtmp += vol * S[INDEX2(s, q, p.numShapes)];
426                      for (k = 0; k < p.numEqu; k++)                      for (k = 0; k < p.numEqu; k++)
427                      EM_F[INDEX2(k, s, p.numEqu)] += rtmp * Y_p[k];                      EM_F[INDEX2(k, s, p.numEqu)] += rtmp * Y_p[k];

Legend:
Removed from v.3204  
changed lines
  Added in v.3205

  ViewVC Help
Powered by ViewVC 1.1.26