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

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

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

revision 2748 by gross, Tue Nov 17 07:32:59 2009 UTC revision 2881 by jfenwick, Thu Jan 28 02:03:15 2010 UTC
# Line 1  Line 1 
1    
2  /*******************************************************  /*******************************************************
3  *  *
4  * Copyright (c) 2003-2009 by University of Queensland  * Copyright (c) 2003-2010 by University of Queensland
5  * Earth Systems Science Computational Center (ESSCC)  * Earth Systems Science Computational Center (ESSCC)
6  * http://www.uq.edu.au/esscc  * http://www.uq.edu.au/esscc
7  *  *
# Line 69  void  Finley_Assemble_PDE_Single2_2D(Ass Line 69  void  Finley_Assemble_PDE_Single2_2D(Ass
69      dim_t len_EM_S=p.row_numShapesTotal*p.col_numShapesTotal;      dim_t len_EM_S=p.row_numShapesTotal*p.col_numShapesTotal;
70      dim_t len_EM_F=p.row_numShapesTotal;      dim_t len_EM_F=p.row_numShapesTotal;
71    
     void* ABuff=allocSampleBuffer(A);  
     void* BBuff=allocSampleBuffer(B);  
     void* CBuff=allocSampleBuffer(C);  
     void* DBuff=allocSampleBuffer(D);  
     void* XBuff=allocSampleBuffer(X);  
     void* YBuff=allocSampleBuffer(Y);  
72      #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,rtmp00, rtmp01, rtmp10, rtmp11, rtmp, rtmp0, rtmp1,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,rtmp00, rtmp01, rtmp10, rtmp11, rtmp, rtmp0, rtmp1,add_EM_F, add_EM_S, isub)
73      {      {
74         EM_S=THREAD_MEMALLOC(len_EM_S,double);         EM_S=THREAD_MEMALLOC(len_EM_S,double);
# Line 89  void  Finley_Assemble_PDE_Single2_2D(Ass Line 83  void  Finley_Assemble_PDE_Single2_2D(Ass
83               for(e=0;e<elements->numElements;e++){               for(e=0;e<elements->numElements;e++){
84                  if (elements->Color[e]==color) {                  if (elements->Color[e]==color) {
85    
86                     A_p=getSampleDataRO(A,e,ABuff);                     A_p=getSampleDataRO(A,e);
87                     B_p=getSampleDataRO(B,e,BBuff);                     B_p=getSampleDataRO(B,e);
88                     C_p=getSampleDataRO(C,e,CBuff);                     C_p=getSampleDataRO(C,e);
89                     D_p=getSampleDataRO(D,e,DBuff);                     D_p=getSampleDataRO(D,e);
90                     X_p=getSampleDataRO(X,e,XBuff);                     X_p=getSampleDataRO(X,e);
91                     Y_p=getSampleDataRO(Y,e,YBuff);                     Y_p=getSampleDataRO(Y,e);
92    
93                     for (isub=0; isub<p.numSub; isub++) {                     for (isub=0; isub<p.numSub; isub++) {
94    
# Line 300  void  Finley_Assemble_PDE_Single2_2D(Ass Line 294  void  Finley_Assemble_PDE_Single2_2D(Ass
294    
295        } /* end of pointer check */        } /* end of pointer check */
296     } /* end parallel region */     } /* end parallel region */
    freeSampleBuffer(ABuff);  
    freeSampleBuffer(BBuff);  
    freeSampleBuffer(CBuff);  
    freeSampleBuffer(DBuff);  
    freeSampleBuffer(XBuff);  
    freeSampleBuffer(YBuff);  
297  }  }

Legend:
Removed from v.2748  
changed lines
  Added in v.2881

  ViewVC Help
Powered by ViewVC 1.1.26