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

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

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

revision 616 by elspeth, Wed Mar 22 02:46:56 2006 UTC revision 751 by bcumming, Mon Jun 26 01:46:34 2006 UTC
# Line 263  void Finley_Assemble_PDE(Finley_NodeFile Line 263  void Finley_Assemble_PDE(Finley_NodeFile
263                  Finley_checkPtr(index_row) || Finley_checkPtr(dVdv) || Finley_checkPtr(dSdV) || Finley_checkPtr(Vol) ))  {                  Finley_checkPtr(index_row) || Finley_checkPtr(dVdv) || Finley_checkPtr(dSdV) || Finley_checkPtr(Vol) ))  {
264    
265             /*  open loop over all colors: */             /*  open loop over all colors: */
266    #ifndef PASO_MPI
267             for (color=elements->minColor;color<=elements->maxColor;color++) {             for (color=elements->minColor;color<=elements->maxColor;color++) {
268                /*  open loop over all elements: */                /*  open loop over all elements: */
269                #pragma omp for private(e) schedule(static)                #pragma omp for private(e) schedule(static)
270                for(e=0;e<elements->numElements;e++){                for(e=0;e<elements->numElements;e++){
271                  if (elements->Color[e]==color) {                  if (elements->Color[e]==color) {
272    #else
273               for(e=0;e<elements->numElements;e++){
274    #endif  
275  //============================  //============================
276                    for (q=0;q<p.NN_row;q++) index_row[q]=p.label_row[elements->Nodes[INDEX2(p.row_node[q],e,p.NN)]];                    for (q=0;q<p.NN_row;q++) index_row[q]=p.label_row[elements->Nodes[INDEX2(p.row_node[q],e,p.NN)]];
277                    /* gather V-coordinates of nodes into V: */                    /* gather V-coordinates of nodes into V: */
# Line 320  void Finley_Assemble_PDE(Finley_NodeFile Line 324  void Finley_Assemble_PDE(Finley_NodeFile
324                                                                   getSampleData(Y,e),isExpanded(Y));                                                                   getSampleData(Y,e),isExpanded(Y));
325                         }                         }
326                         /* add  */                         /* add  */
327    #ifndef PASO_MPI
328                         Finley_Util_AddScatter(p.NN_row,index_row,p.numEqu,EM_F,getSampleData(F,0));                         Finley_Util_AddScatter(p.NN_row,index_row,p.numEqu,EM_F,getSampleData(F,0));
329    #else
330                           Finley_Util_AddScatter_upperBound(p.NN_row,index_row,p.numEqu,EM_F,getSampleData(F,0),p.degreeOfFreedomUpperBound);
331    #endif                                      
332                      }                      }
333                  }                  }
334                }                }
335    #ifndef PASO_MPI
336              }              }
337           }           }
338    #endif
339           /* clean up */           /* clean up */
340           THREAD_MEMFREE(EM_S);           THREAD_MEMFREE(EM_S);
341           THREAD_MEMFREE(EM_F);           THREAD_MEMFREE(EM_F);

Legend:
Removed from v.616  
changed lines
  Added in v.751

  ViewVC Help
Powered by ViewVC 1.1.26