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

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

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

revision 2769 by gross, Tue Nov 17 07:32:59 2009 UTC revision 2770 by jfenwick, Wed Nov 25 01:24:51 2009 UTC
# Line 70  void  Finley_Assemble_PDE_Single2_3D(Ass Line 70  void  Finley_Assemble_PDE_Single2_3D(Ass
70      dim_t len_EM_S=p.row_numShapesTotal*p.col_numShapesTotal;      dim_t len_EM_S=p.row_numShapesTotal*p.col_numShapesTotal;
71      dim_t len_EM_F=p.row_numShapesTotal;      dim_t len_EM_F=p.row_numShapesTotal;
72    
     void* ABuff=allocSampleBuffer(A);  
     void* BBuff=allocSampleBuffer(B);  
     void* CBuff=allocSampleBuffer(C);  
     void* DBuff=allocSampleBuffer(D);  
     void* XBuff=allocSampleBuffer(X);  
     void* YBuff=allocSampleBuffer(Y);  
73      #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, rtmp00, rtmp01, rtmp02, rtmp10, rtmp11, rtmp12, rtmp20, rtmp21, rtmp22, rtmp0, rtmp1, rtmp2,add_EM_F, add_EM_S, isub)      #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, rtmp00, rtmp01, rtmp02, rtmp10, rtmp11, rtmp12, rtmp20, rtmp21, rtmp22, rtmp0, rtmp1, rtmp2,add_EM_F, add_EM_S, isub)
74      {      {
75         EM_S=THREAD_MEMALLOC(len_EM_S,double);         EM_S=THREAD_MEMALLOC(len_EM_S,double);
# Line 90  void  Finley_Assemble_PDE_Single2_3D(Ass Line 84  void  Finley_Assemble_PDE_Single2_3D(Ass
84               for(e=0;e<elements->numElements;e++){               for(e=0;e<elements->numElements;e++){
85                  if (elements->Color[e]==color) {                  if (elements->Color[e]==color) {
86    
87                     A_p=getSampleDataRO(A,e,ABuff);                     A_p=getSampleDataRO(A,e);
88                     B_p=getSampleDataRO(B,e,BBuff);                     B_p=getSampleDataRO(B,e);
89                     C_p=getSampleDataRO(C,e,CBuff);                     C_p=getSampleDataRO(C,e);
90                     D_p=getSampleDataRO(D,e,DBuff);                     D_p=getSampleDataRO(D,e);
91                     X_p=getSampleDataRO(X,e,XBuff);                     X_p=getSampleDataRO(X,e);
92                     Y_p=getSampleDataRO(Y,e,YBuff);                     Y_p=getSampleDataRO(Y,e);
93    
94                     for (isub=0; isub<p.numSub; isub++) {                     for (isub=0; isub<p.numSub; isub++) {
95    
# Line 338  void  Finley_Assemble_PDE_Single2_3D(Ass Line 332  void  Finley_Assemble_PDE_Single2_3D(Ass
332    
333        } /* end of pointer check */        } /* end of pointer check */
334     } /* end parallel region */     } /* end parallel region */
    freeSampleBuffer(ABuff);  
    freeSampleBuffer(BBuff);  
    freeSampleBuffer(CBuff);  
    freeSampleBuffer(DBuff);  
    freeSampleBuffer(XBuff);  
    freeSampleBuffer(YBuff);  
335  }  }

Legend:
Removed from v.2769  
changed lines
  Added in v.2770

  ViewVC Help
Powered by ViewVC 1.1.26