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

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

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

revision 3202 by jfenwick, Thu Sep 23 23:24:09 2010 UTC revision 3203 by jfenwick, Thu Sep 23 23:51:26 2010 UTC
# Line 65  void Dudley_Assemble_PDE_System2_3D(Asse Line 65  void Dudley_Assemble_PDE_System2_3D(Asse
65      double *F_p = (requireWrite(F), getSampleDataRW(F, 0)); /* use comma, to get around the mixed code and declarations thing */      double *F_p = (requireWrite(F), getSampleDataRW(F, 0)); /* use comma, to get around the mixed code and declarations thing */
66  //    double *S = p.row_jac->BasisFunctions->S;  //    double *S = p.row_jac->BasisFunctions->S;
67      const double* S=p.shapeFns;      const double* S=p.shapeFns;
68      dim_t len_EM_S = p.row_numShapesTotal * p.col_numShapesTotal * p.numEqu * p.numComp;      dim_t len_EM_S = p.row_numShapesTotal * p.row_numShapesTotal * p.numEqu * p.numComp;
69      dim_t len_EM_F = p.row_numShapesTotal * p.numEqu;      dim_t len_EM_F = p.row_numShapesTotal * p.numEqu;
70    
71  #pragma omp parallel private(color,EM_S, EM_F, Vol, DSDX, A_p, B_p, C_p, D_p, X_p, Y_p, A_q, B_q, C_q, D_q, X_q, Y_q,row_index,q, s,r,k,m,rtmp, rtmp0, rtmp1, rtmp2, rtmp00, rtmp01, rtmp02, rtmp10, rtmp11, rtmp12, rtmp20, rtmp21, rtmp22,add_EM_F, add_EM_S)  #pragma omp parallel private(color,EM_S, EM_F, Vol, DSDX, A_p, B_p, C_p, D_p, X_p, Y_p, A_q, B_q, C_q, D_q, X_q, Y_q,row_index,q, s,r,k,m,rtmp, rtmp0, rtmp1, rtmp2, rtmp00, rtmp01, rtmp02, rtmp10, rtmp11, rtmp12, rtmp20, rtmp21, rtmp22,add_EM_F, add_EM_S)
# Line 113  void Dudley_Assemble_PDE_System2_3D(Asse Line 113  void Dudley_Assemble_PDE_System2_3D(Asse
113                  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.numQuadTotal)]);
114                  for (s = 0; s < p.row_numShapes; s++)                  for (s = 0; s < p.row_numShapes; s++)
115                  {                  {
116                      for (r = 0; r < p.col_numShapes; r++)                      for (r = 0; r < p.row_numShapes; r++)
117                      {                      {
118                      for (k = 0; k < p.numEqu; k++)                      for (k = 0; k < p.numEqu; k++)
119                      {                      {
# Line 162  void Dudley_Assemble_PDE_System2_3D(Asse Line 162  void Dudley_Assemble_PDE_System2_3D(Asse
162                  {                  {
163                  for (s = 0; s < p.row_numShapes; s++)                  for (s = 0; s < p.row_numShapes; s++)
164                  {                  {
165                      for (r = 0; r < p.col_numShapes; r++)                      for (r = 0; r < p.row_numShapes; r++)
166                      {                      {
167                      rtmp00 = 0;                      rtmp00 = 0;
168                      rtmp01 = 0;                      rtmp01 = 0;
# Line 221  void Dudley_Assemble_PDE_System2_3D(Asse Line 221  void Dudley_Assemble_PDE_System2_3D(Asse
221                  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.numQuadTotal)]);
222                  for (s = 0; s < p.row_numShapes; s++)                  for (s = 0; s < p.row_numShapes; s++)
223                  {                  {
224                      for (r = 0; r < p.col_numShapes; r++)                      for (r = 0; r < p.row_numShapes; r++)
225                      {                      {
226                      for (k = 0; k < p.numEqu; k++)                      for (k = 0; k < p.numEqu; k++)
227                      {                      {
# Line 249  void Dudley_Assemble_PDE_System2_3D(Asse Line 249  void Dudley_Assemble_PDE_System2_3D(Asse
249                  {                  {
250                  for (s = 0; s < p.row_numShapes; s++)                  for (s = 0; s < p.row_numShapes; s++)
251                  {                  {
252                      for (r = 0; r < p.col_numShapes; r++)                      for (r = 0; r < p.row_numShapes; r++)
253                      {                      {
254                      rtmp0 = 0;                      rtmp0 = 0;
255                      rtmp1 = 0;                      rtmp1 = 0;
# Line 286  void Dudley_Assemble_PDE_System2_3D(Asse Line 286  void Dudley_Assemble_PDE_System2_3D(Asse
286                  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.numQuadTotal)]);
287                  for (s = 0; s < p.row_numShapes; s++)                  for (s = 0; s < p.row_numShapes; s++)
288                  {                  {
289                      for (r = 0; r < p.col_numShapes; r++)                      for (r = 0; r < p.row_numShapes; r++)
290                      {                      {
291                      for (k = 0; k < p.numEqu; k++)                      for (k = 0; k < p.numEqu; k++)
292                      {                      {
# Line 314  void Dudley_Assemble_PDE_System2_3D(Asse Line 314  void Dudley_Assemble_PDE_System2_3D(Asse
314                  {                  {
315                  for (s = 0; s < p.row_numShapes; s++)                  for (s = 0; s < p.row_numShapes; s++)
316                  {                  {
317                      for (r = 0; r < p.col_numShapes; r++)                      for (r = 0; r < p.row_numShapes; r++)
318                      {                      {
319                      rtmp0 = 0;                      rtmp0 = 0;
320                      rtmp1 = 0;                      rtmp1 = 0;
# Line 351  void Dudley_Assemble_PDE_System2_3D(Asse Line 351  void Dudley_Assemble_PDE_System2_3D(Asse
351                  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.numQuadTotal)]);
352                  for (s = 0; s < p.row_numShapes; s++)                  for (s = 0; s < p.row_numShapes; s++)
353                  {                  {
354                      for (r = 0; r < p.col_numShapes; r++)                      for (r = 0; r < p.row_numShapes; r++)
355                      {                      {
356                      for (k = 0; k < p.numEqu; k++)                      for (k = 0; k < p.numEqu; k++)
357                      {                      {
# Line 375  void Dudley_Assemble_PDE_System2_3D(Asse Line 375  void Dudley_Assemble_PDE_System2_3D(Asse
375                  {                  {
376                  for (s = 0; s < p.row_numShapes; s++)                  for (s = 0; s < p.row_numShapes; s++)
377                  {                  {
378                      for (r = 0; r < p.col_numShapes; r++)                      for (r = 0; r < p.row_numShapes; r++)
379                      {                      {
380                      rtmp = 0;                      rtmp = 0;
381                      for (q = 0; q < p.numQuadTotal; q++)                      for (q = 0; q < p.numQuadTotal; q++)
# Line 486  void Dudley_Assemble_PDE_System2_3D(Asse Line 486  void Dudley_Assemble_PDE_System2_3D(Asse
486                             p.row_DOF_UpperBound);                             p.row_DOF_UpperBound);
487              if (add_EM_S)              if (add_EM_S)
488                  Dudley_Assemble_addToSystemMatrix(Mat, p.row_numShapesTotal, row_index, p.numEqu,                  Dudley_Assemble_addToSystemMatrix(Mat, p.row_numShapesTotal, row_index, p.numEqu,
489                                    p.col_numShapesTotal, row_index, p.numComp, EM_S);                                    p.row_numShapesTotal, row_index, p.numComp, EM_S);
490              }       /* end color check */              }       /* end color check */
491          }       /* end element loop */          }       /* end element loop */
492          }           /* end color loop */          }           /* end color loop */

Legend:
Removed from v.3202  
changed lines
  Added in v.3203

  ViewVC Help
Powered by ViewVC 1.1.26