/[escript]/trunk/finley/src/Assemble_LumpedSystem.c
ViewVC logotype

Diff of /trunk/finley/src/Assemble_LumpedSystem.c

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

revision 2989 by gross, Thu Mar 18 01:45:55 2010 UTC revision 2990 by gross, Fri Mar 19 05:57:22 2010 UTC
# Line 102  void Finley_Assemble_LumpedSystem(Finley Line 102  void Finley_Assemble_LumpedSystem(Finley
102      }      }
103    }    }
104    if (Finley_noError()) {    if (Finley_noError()) {
      printf("SADSADSADAS 2\n");  
105      requireWrite(lumpedMat);      requireWrite(lumpedMat);
106      lumpedMat_p=getSampleDataRW(lumpedMat,0);      lumpedMat_p=getSampleDataRW(lumpedMat,0);
107      len_EM_lumpedMat=p.row_numShapesTotal*p.numEqu;      len_EM_lumpedMat=p.row_numShapesTotal*p.numEqu;
# Line 161  void Finley_Assemble_LumpedSystem(Finley Line 160  void Finley_Assemble_LumpedSystem(Finley
160                      } /* end element loop */                      } /* end element loop */
161                    } /* end color loop */                    } /* end color loop */
162               } else  {  /* with constant D */                 } else  {  /* with constant D */  
163  printf("SADSADSADAS\n");  
164                   for (color=elements->minColor;color<=elements->maxColor;color++) {                   for (color=elements->minColor;color<=elements->maxColor;color++) {
165                      /*  open loop over all elements: */                      /*  open loop over all elements: */
166                      #pragma omp for private(e) schedule(static)                      #pragma omp for private(e) schedule(static)
# Line 189  printf("SADSADSADAS\n"); Line 188  printf("SADSADSADAS\n");
188                                         for (q=0;q<p.numQuadSub;q++) rtmp+=Vol[q]*S[INDEX2(s,q,p.row_numShapes)];                                         for (q=0;q<p.numQuadSub;q++) rtmp+=Vol[q]*S[INDEX2(s,q,p.row_numShapes)];
189                                         EM_lumpedMat[INDEX2(0,s,p.numEqu)]=rtmp*D_p[0];                                         EM_lumpedMat[INDEX2(0,s,p.numEqu)]=rtmp*D_p[0];
190                                     }                                     }
 for (s=0;s<p.row_numShapes;s++) printf(" %d %d : %e\n",isub,s,EM_lumpedMat[INDEX2(0,s,p.numEqu)]);  
191                                 #endif                                 #endif
192                                 for (q=0;q<p.row_numShapesTotal;q++) row_index[q]=p.row_DOF[elements->Nodes[INDEX2(p.row_node[INDEX2(q,isub,p.row_numShapesTotal)],e,p.NN)]];                                 for (q=0;q<p.row_numShapesTotal;q++) row_index[q]=p.row_DOF[elements->Nodes[INDEX2(p.row_node[INDEX2(q,isub,p.row_numShapesTotal)],e,p.NN)]];
193                                 Finley_Util_AddScatter(p.row_numShapesTotal,row_index,p.numEqu,EM_lumpedMat,lumpedMat_p, p.row_DOF_UpperBound);                                 Finley_Util_AddScatter(p.row_numShapesTotal,row_index,p.numEqu,EM_lumpedMat,lumpedMat_p, p.row_DOF_UpperBound);

Legend:
Removed from v.2989  
changed lines
  Added in v.2990

  ViewVC Help
Powered by ViewVC 1.1.26