/[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 3202 by jfenwick, Thu Sep 16 02:57:17 2010 UTC revision 3203 by jfenwick, Thu Sep 23 23:51:26 2010 UTC
# Line 64  void Dudley_Assemble_PDE_Single2_1D(Asse Line 64  void Dudley_Assemble_PDE_Single2_1D(Asse
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;      dim_t len_EM_S = p.row_numShapesTotal * p.row_numShapesTotal;
68      dim_t len_EM_F = p.row_numShapesTotal;      dim_t len_EM_F = p.row_numShapesTotal;
69    
70  #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,rtmp,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,rtmp,add_EM_F, add_EM_S)
# Line 112  void Dudley_Assemble_PDE_Single2_1D(Asse Line 112  void Dudley_Assemble_PDE_Single2_1D(Asse
112                  A_q = &(A_p[INDEX4(0, 0, 0, 0, DIM, DIM, p.numQuadTotal)]);                  A_q = &(A_p[INDEX4(0, 0, 0, 0, DIM, DIM, p.numQuadTotal)]);
113                  for (s = 0; s < p.row_numShapes; s++)                  for (s = 0; s < p.row_numShapes; s++)
114                  {                  {
115                      for (r = 0; r < p.col_numShapes; r++)                      for (r = 0; r < p.row_numShapes; r++)
116                      {                      {
117                      rtmp = 0;                      rtmp = 0;
118                      for (q = 0; q < p.numQuadTotal; q++)                      for (q = 0; q < p.numQuadTotal; q++)
# Line 129  void Dudley_Assemble_PDE_Single2_1D(Asse Line 129  void Dudley_Assemble_PDE_Single2_1D(Asse
129                  {                  {
130                  for (s = 0; s < p.row_numShapes; s++)                  for (s = 0; s < p.row_numShapes; s++)
131                  {                  {
132                      for (r = 0; r < p.col_numShapes; r++)                      for (r = 0; r < p.row_numShapes; r++)
133                      {                      {
134                      rtmp = 0;                      rtmp = 0;
135                      for (q = 0; q < p.numQuadTotal; q++)                      for (q = 0; q < p.numQuadTotal; q++)
# Line 153  void Dudley_Assemble_PDE_Single2_1D(Asse Line 153  void Dudley_Assemble_PDE_Single2_1D(Asse
153                  B_q = &(B_p[INDEX3(0, 0, 0, DIM, p.numQuadTotal)]);                  B_q = &(B_p[INDEX3(0, 0, 0, DIM, p.numQuadTotal)]);
154                  for (s = 0; s < p.row_numShapes; s++)                  for (s = 0; s < p.row_numShapes; s++)
155                  {                  {
156                      for (r = 0; r < p.col_numShapes; r++)                      for (r = 0; r < p.row_numShapes; r++)
157                      {                      {
158                      rtmp = 0;                      rtmp = 0;
159                      for (q = 0; q < p.numQuadTotal; q++)                      for (q = 0; q < p.numQuadTotal; q++)
# Line 169  void Dudley_Assemble_PDE_Single2_1D(Asse Line 169  void Dudley_Assemble_PDE_Single2_1D(Asse
169                  {                  {
170                  for (s = 0; s < p.row_numShapes; s++)                  for (s = 0; s < p.row_numShapes; s++)
171                  {                  {
172                      for (r = 0; r < p.col_numShapes; r++)                      for (r = 0; r < p.row_numShapes; r++)
173                      {                      {
174                      rtmp = 0;                      rtmp = 0;
175                      for (q = 0; q < p.numQuadTotal; q++)                      for (q = 0; q < p.numQuadTotal; q++)
# Line 193  void Dudley_Assemble_PDE_Single2_1D(Asse Line 193  void Dudley_Assemble_PDE_Single2_1D(Asse
193                  C_q = &(C_p[INDEX3(0, 0, 0, DIM, p.numQuadTotal)]);                  C_q = &(C_p[INDEX3(0, 0, 0, DIM, p.numQuadTotal)]);
194                  for (s = 0; s < p.row_numShapes; s++)                  for (s = 0; s < p.row_numShapes; s++)
195                  {                  {
196                      for (r = 0; r < p.col_numShapes; r++)                      for (r = 0; r < p.row_numShapes; r++)
197                      {                      {
198                      rtmp = 0;                      rtmp = 0;
199                      for (q = 0; q < p.numQuadTotal; q++)                      for (q = 0; q < p.numQuadTotal; q++)
# Line 209  void Dudley_Assemble_PDE_Single2_1D(Asse Line 209  void Dudley_Assemble_PDE_Single2_1D(Asse
209                  {                  {
210                  for (s = 0; s < p.row_numShapes; s++)                  for (s = 0; s < p.row_numShapes; s++)
211                  {                  {
212                      for (r = 0; r < p.col_numShapes; r++)                      for (r = 0; r < p.row_numShapes; r++)
213                      {                      {
214                      rtmp = 0;                      rtmp = 0;
215                      for (q = 0; q < p.numQuadTotal; q++)                      for (q = 0; q < p.numQuadTotal; q++)
# Line 233  void Dudley_Assemble_PDE_Single2_1D(Asse Line 233  void Dudley_Assemble_PDE_Single2_1D(Asse
233                  D_q = &(D_p[INDEX2(0, 0, p.numQuadTotal)]);                  D_q = &(D_p[INDEX2(0, 0, p.numQuadTotal)]);
234                  for (s = 0; s < p.row_numShapes; s++)                  for (s = 0; s < p.row_numShapes; s++)
235                  {                  {
236                      for (r = 0; r < p.col_numShapes; r++)                      for (r = 0; r < p.row_numShapes; r++)
237                      {                      {
238                      rtmp = 0;                      rtmp = 0;
239                      for (q = 0; q < p.numQuadTotal; q++)                      for (q = 0; q < p.numQuadTotal; q++)
# Line 249  void Dudley_Assemble_PDE_Single2_1D(Asse Line 249  void Dudley_Assemble_PDE_Single2_1D(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                      rtmp = 0;                      rtmp = 0;
255                      for (q = 0; q < p.numQuadTotal; q++)                      for (q = 0; q < p.numQuadTotal; q++)
# Line 329  void Dudley_Assemble_PDE_Single2_1D(Asse Line 329  void Dudley_Assemble_PDE_Single2_1D(Asse
329                             p.row_DOF_UpperBound);                             p.row_DOF_UpperBound);
330              if (add_EM_S)              if (add_EM_S)
331                  Dudley_Assemble_addToSystemMatrix(Mat, p.row_numShapesTotal, row_index, p.numEqu,                  Dudley_Assemble_addToSystemMatrix(Mat, p.row_numShapesTotal, row_index, p.numEqu,
332                                    p.col_numShapesTotal, row_index, p.numComp, EM_S);                                    p.row_numShapesTotal, row_index, p.numComp, EM_S);
333    
334              }       /* end color check */              }       /* end color check */
335          }       /* end element loop */          }       /* end element loop */

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

  ViewVC Help
Powered by ViewVC 1.1.26