/[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 3187 by jfenwick, Thu Sep 16 02:57:17 2010 UTC revision 3203 by jfenwick, Thu Sep 23 23:51:26 2010 UTC
# Line 63  void Dudley_Assemble_PDE_System2_2D(Asse Line 63  void Dudley_Assemble_PDE_System2_2D(Asse
63      bool_t extendedX = isExpanded(X);      bool_t extendedX = isExpanded(X);
64      bool_t extendedY = isExpanded(Y);      bool_t extendedY = isExpanded(Y);
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      dim_t len_EM_S = p.row_numShapesTotal * p.col_numShapesTotal * p.numEqu * p.numComp;      const double* S = p.shapeFns;
68        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, rtmp00, rtmp10, rtmp01, rtmp11,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, rtmp00, rtmp10, rtmp01, rtmp11,add_EM_F, add_EM_S)
# Line 113  void Dudley_Assemble_PDE_System2_2D(Asse Line 114  void Dudley_Assemble_PDE_System2_2D(Asse
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.numQuadTotal)]);
115                  for (s = 0; s < p.row_numShapes; s++)                  for (s = 0; s < p.row_numShapes; s++)
116                  {                  {
117                      for (r = 0; r < p.col_numShapes; r++)                      for (r = 0; r < p.row_numShapes; r++)
118                      {                      {
119                      for (k = 0; k < p.numEqu; k++)                      for (k = 0; k < p.numEqu; k++)
120                      {                      {
# Line 146  void Dudley_Assemble_PDE_System2_2D(Asse Line 147  void Dudley_Assemble_PDE_System2_2D(Asse
147                  {                  {
148                  for (s = 0; s < p.row_numShapes; s++)                  for (s = 0; s < p.row_numShapes; s++)
149                  {                  {
150                      for (r = 0; r < p.col_numShapes; r++)                      for (r = 0; r < p.row_numShapes; r++)
151                      {                      {
152                      rtmp00 = 0;                      rtmp00 = 0;
153                      rtmp01 = 0;                      rtmp01 = 0;
# Line 188  void Dudley_Assemble_PDE_System2_2D(Asse Line 189  void Dudley_Assemble_PDE_System2_2D(Asse
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.numQuadTotal)]);
190                  for (s = 0; s < p.row_numShapes; s++)                  for (s = 0; s < p.row_numShapes; s++)
191                  {                  {
192                      for (r = 0; r < p.col_numShapes; r++)                      for (r = 0; r < p.row_numShapes; r++)
193                      {                      {
194                      for (k = 0; k < p.numEqu; k++)                      for (k = 0; k < p.numEqu; k++)
195                      {                      {
# Line 213  void Dudley_Assemble_PDE_System2_2D(Asse Line 214  void Dudley_Assemble_PDE_System2_2D(Asse
214                  {                  {
215                  for (s = 0; s < p.row_numShapes; s++)                  for (s = 0; s < p.row_numShapes; s++)
216                  {                  {
217                      for (r = 0; r < p.col_numShapes; r++)                      for (r = 0; r < p.row_numShapes; r++)
218                      {                      {
219                      rtmp0 = 0;                      rtmp0 = 0;
220                      rtmp1 = 0;                      rtmp1 = 0;
# Line 248  void Dudley_Assemble_PDE_System2_2D(Asse Line 249  void Dudley_Assemble_PDE_System2_2D(Asse
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.numQuadTotal)]);
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                      for (k = 0; k < p.numEqu; k++)                      for (k = 0; k < p.numEqu; k++)
255                      {                      {
# Line 273  void Dudley_Assemble_PDE_System2_2D(Asse Line 274  void Dudley_Assemble_PDE_System2_2D(Asse
274                  {                  {
275                  for (s = 0; s < p.row_numShapes; s++)                  for (s = 0; s < p.row_numShapes; s++)
276                  {                  {
277                      for (r = 0; r < p.col_numShapes; r++)                      for (r = 0; r < p.row_numShapes; r++)
278                      {                      {
279                      rtmp0 = 0;                      rtmp0 = 0;
280                      rtmp1 = 0;                      rtmp1 = 0;
# Line 308  void Dudley_Assemble_PDE_System2_2D(Asse Line 309  void Dudley_Assemble_PDE_System2_2D(Asse
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.numQuadTotal)]);
310                  for (s = 0; s < p.row_numShapes; s++)                  for (s = 0; s < p.row_numShapes; s++)
311                  {                  {
312                      for (r = 0; r < p.col_numShapes; r++)                      for (r = 0; r < p.row_numShapes; r++)
313                      {                      {
314                      for (k = 0; k < p.numEqu; k++)                      for (k = 0; k < p.numEqu; k++)
315                      {                      {
# Line 332  void Dudley_Assemble_PDE_System2_2D(Asse Line 333  void Dudley_Assemble_PDE_System2_2D(Asse
333                  {                  {
334                  for (s = 0; s < p.row_numShapes; s++)                  for (s = 0; s < p.row_numShapes; s++)
335                  {                  {
336                      for (r = 0; r < p.col_numShapes; r++)                      for (r = 0; r < p.row_numShapes; r++)
337                      {                      {
338                      rtmp = 0;                      rtmp = 0;
339                      for (q = 0; q < p.numQuadTotal; q++)                      for (q = 0; q < p.numQuadTotal; q++)
# Line 438  void Dudley_Assemble_PDE_System2_2D(Asse Line 439  void Dudley_Assemble_PDE_System2_2D(Asse
439                             p.row_DOF_UpperBound);                             p.row_DOF_UpperBound);
440              if (add_EM_S)              if (add_EM_S)
441                  Dudley_Assemble_addToSystemMatrix(Mat, p.row_numShapesTotal, row_index, p.numEqu,                  Dudley_Assemble_addToSystemMatrix(Mat, p.row_numShapesTotal, row_index, p.numEqu,
442                                    p.col_numShapesTotal, row_index, p.numComp, EM_S);                                    p.row_numShapesTotal, row_index, p.numComp, EM_S);
443    
444              }       /* end color check */              }       /* end color check */
445          }       /* end element loop */          }       /* end element loop */

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

  ViewVC Help
Powered by ViewVC 1.1.26