/[escript]/branches/doubleplusgood/dudley/src/Assemble_PDE_System2_2D.cpp
ViewVC logotype

Diff of /branches/doubleplusgood/dudley/src/Assemble_PDE_System2_2D.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_System2_2D(Ass Line 70  void  Finley_Assemble_PDE_System2_2D(Ass
70      dim_t len_EM_S=p.row_numShapesTotal*p.col_numShapesTotal*p.numEqu*p.numComp;      dim_t len_EM_S=p.row_numShapesTotal*p.col_numShapesTotal*p.numEqu*p.numComp;
71      dim_t len_EM_F=p.row_numShapesTotal*p.numEqu;      dim_t len_EM_F=p.row_numShapesTotal*p.numEqu;
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,k,m,rtmp, rtmp0, rtmp1, rtmp00, rtmp10, rtmp01, rtmp11,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,k,m,rtmp, rtmp0, rtmp1, rtmp00, rtmp10, rtmp01, rtmp11,add_EM_F, add_EM_S, isub)
74      {      {
75    
# Line 91  void  Finley_Assemble_PDE_System2_2D(Ass Line 85  void  Finley_Assemble_PDE_System2_2D(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             B_p=getSampleDataRO(B,e,BBuff);             B_p=getSampleDataRO(B,e);
90                     C_p=getSampleDataRO(C,e,CBuff);                       C_p=getSampleDataRO(C,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    
# Line 110  void  Finley_Assemble_PDE_System2_2D(Ass Line 104  void  Finley_Assemble_PDE_System2_2D(Ass
104                        /**************************************************************/                        /**************************************************************/
105                        /*   process A: */                        /*   process A: */
106                        /**************************************************************/                        /**************************************************************/
107                        A_p=getSampleDataRO(A,e,ABuff);                        A_p=getSampleDataRO(A,e);
108                        if (NULL!=A_p) {                        if (NULL!=A_p) {
109                           add_EM_S=TRUE;                           add_EM_S=TRUE;
110                           if (extendedA) {                           if (extendedA) {
# Line 163  void  Finley_Assemble_PDE_System2_2D(Ass Line 157  void  Finley_Assemble_PDE_System2_2D(Ass
157                        /**************************************************************/                        /**************************************************************/
158                        /*   process B: */                        /*   process B: */
159                        /**************************************************************/                        /**************************************************************/
160                        B_p=getSampleDataRO(B,e,BBuff);                        B_p=getSampleDataRO(B,e);
161                        if (NULL!=B_p) {                        if (NULL!=B_p) {
162                           add_EM_S=TRUE;                           add_EM_S=TRUE;
163                           if (extendedB) {                           if (extendedB) {
# Line 206  void  Finley_Assemble_PDE_System2_2D(Ass Line 200  void  Finley_Assemble_PDE_System2_2D(Ass
200                        /**************************************************************/                        /**************************************************************/
201                        /*   process C: */                        /*   process C: */
202                        /**************************************************************/                        /**************************************************************/
203                        C_p=getSampleDataRO(C,e,CBuff);                        C_p=getSampleDataRO(C,e);
204                        if (NULL!=C_p) {                        if (NULL!=C_p) {
205                          add_EM_S=TRUE;                          add_EM_S=TRUE;
206                          if (extendedC) {                          if (extendedC) {
# Line 249  void  Finley_Assemble_PDE_System2_2D(Ass Line 243  void  Finley_Assemble_PDE_System2_2D(Ass
243                        /************************************************************* */                        /************************************************************* */
244                        /* process D */                        /* process D */
245                        /**************************************************************/                        /**************************************************************/
246                        D_p=getSampleDataRO(D,e,DBuff);                        D_p=getSampleDataRO(D,e);
247                        if (NULL!=D_p) {                        if (NULL!=D_p) {
248                          add_EM_S=TRUE;                          add_EM_S=TRUE;
249                          if (extendedD) {                          if (extendedD) {
# Line 284  void  Finley_Assemble_PDE_System2_2D(Ass Line 278  void  Finley_Assemble_PDE_System2_2D(Ass
278                        /**************************************************************/                        /**************************************************************/
279                        /*   process X: */                        /*   process X: */
280                        /**************************************************************/                        /**************************************************************/
281                        X_p=getSampleDataRO(X,e,XBuff);                        X_p=getSampleDataRO(X,e);
282                        if (NULL!=X_p) {                        if (NULL!=X_p) {
283                          add_EM_F=TRUE;                          add_EM_F=TRUE;
284                          if (extendedX) {                          if (extendedX) {
# Line 314  void  Finley_Assemble_PDE_System2_2D(Ass Line 308  void  Finley_Assemble_PDE_System2_2D(Ass
308                       /**************************************************************/                       /**************************************************************/
309                       /*   process Y: */                       /*   process Y: */
310                       /**************************************************************/                       /**************************************************************/
311                        Y_p=getSampleDataRO(Y,e,YBuff);                        Y_p=getSampleDataRO(Y,e);
312                        if (NULL!=Y_p) {                        if (NULL!=Y_p) {
313                          add_EM_F=TRUE;                          add_EM_F=TRUE;
314                          if (extendedY) {                          if (extendedY) {
# Line 353  void  Finley_Assemble_PDE_System2_2D(Ass Line 347  void  Finley_Assemble_PDE_System2_2D(Ass
347    
348        } /* end of pointer check */        } /* end of pointer check */
349     } /* end parallel region */     } /* end parallel region */
    freeSampleBuffer(ABuff);  
    freeSampleBuffer(BBuff);  
    freeSampleBuffer(CBuff);  
    freeSampleBuffer(DBuff);  
    freeSampleBuffer(XBuff);  
    freeSampleBuffer(YBuff);  
350  }  }
351  /*  /*
352   * $Log$   * $Log$

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

  ViewVC Help
Powered by ViewVC 1.1.26