/[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 2748 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_1D(Ass Line 70  void  Finley_Assemble_PDE_Single2_1D(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,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,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 91  void  Finley_Assemble_PDE_Single2_1D(Ass Line 85  void  Finley_Assemble_PDE_Single2_1D(Ass
85               for(e=0;e<elements->numElements;e++){               for(e=0;e<elements->numElements;e++){
86                  if (elements->Color[e]==color) {                  if (elements->Color[e]==color) {
87    
88                    A_p=getSampleDataRO(A,e,ABuff);                    A_p=getSampleDataRO(A,e);
89                    C_p=getSampleDataRO(C,e,CBuff);                    C_p=getSampleDataRO(C,e);
90                    B_p=getSampleDataRO(B,e,BBuff);                    B_p=getSampleDataRO(B,e);
91                    D_p=getSampleDataRO(D,e,DBuff);                    D_p=getSampleDataRO(D,e);
92                    X_p=getSampleDataRO(X,e,XBuff);                    X_p=getSampleDataRO(X,e);
93                    Y_p=getSampleDataRO(Y,e,YBuff);                    Y_p=getSampleDataRO(Y,e);
94    
95                    for (isub=0; isub<p.numSub; isub++) {                    for (isub=0; isub<p.numSub; isub++) {
96                        Vol=&(p.row_jac->volume[INDEX3(0,isub,e, p.numQuadSub,p.numSub)]);                        Vol=&(p.row_jac->volume[INDEX3(0,isub,e, p.numQuadSub,p.numSub)]);
# Line 268  void  Finley_Assemble_PDE_Single2_1D(Ass Line 262  void  Finley_Assemble_PDE_Single2_1D(Ass
262    
263        } /* end of pointer check */        } /* end of pointer check */
264     } /* end parallel region */     } /* end parallel region */
    freeSampleBuffer(ABuff);  
    freeSampleBuffer(BBuff);  
    freeSampleBuffer(CBuff);  
    freeSampleBuffer(DBuff);  
    freeSampleBuffer(XBuff);  
    freeSampleBuffer(YBuff);  
265  }  }
266  /*  /*
267   * $Log$   * $Log$

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

  ViewVC Help
Powered by ViewVC 1.1.26