/[escript]/branches/ripleygmg_from_3668/ripley/src/Rectangle.cpp
ViewVC logotype

Diff of /branches/ripleygmg_from_3668/ripley/src/Rectangle.cpp

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

revision 3760 by caltinay, Mon Jan 9 05:21:18 2012 UTC revision 3761 by caltinay, Mon Jan 9 06:50:33 2012 UTC
# Line 2173  void Rectangle::assemblePDESingle(Paso_S Line 2173  void Rectangle::assemblePDESingle(Paso_S
2173    
2174                  } // end k0 loop                  } // end k0 loop
2175              } // end k1 loop              } // end k1 loop
2176          } // end of coloring          } // end of colouring
2177      } // end of parallel region      } // end of parallel region
2178  }  }
2179    
# Line 2194  void Rectangle::assemblePDESystem(Paso_S Line 2194  void Rectangle::assemblePDESystem(Paso_S
2194          numComp=mat->logical_col_block_size;          numComp=mat->logical_col_block_size;
2195      }      }
2196    
2197      /* GENERATOR SNIP_PDE_SYSTEM_PRE TOP */      /*** GENERATOR SNIP_PDE_SYSTEM_PRE TOP */
2198      const double w0 = -0.1555021169820365539*h1/h0;      const double w0 = -0.1555021169820365539*h1/h0;
2199      const double w1 = 0.041666666666666666667;      const double w1 = 0.041666666666666666667;
2200      const double w10 = -0.041666666666666666667*h0/h1;      const double w10 = -0.041666666666666666667*h0/h1;
# Line 2285  void Rectangle::assemblePDESystem(Paso_S Line 2285  void Rectangle::assemblePDESystem(Paso_S
2285                      vector<double> EM_S(4*4*numEq*numComp, 0);                      vector<double> EM_S(4*4*numEq*numComp, 0);
2286                      vector<double> EM_F(4*numEq, 0);                      vector<double> EM_F(4*numEq, 0);
2287                      const index_t e = k0 + m_NE0*k1;                      const index_t e = k0 + m_NE0*k1;
2288                      /* GENERATOR SNIP_PDE_SYSTEM TOP */                      /*** GENERATOR SNIP_PDE_SYSTEM TOP */
2289                      ///////////////                      ///////////////
2290                      // process A //                      // process A //
2291                      ///////////////                      ///////////////
# Line 2947  void Rectangle::assemblePDESystem(Paso_S Line 2947  void Rectangle::assemblePDESystem(Paso_S
2947                          double *F_p=rhs.getSampleDataRW(0);                          double *F_p=rhs.getSampleDataRW(0);
2948                          for (index_t i=0; i<rowIndex.size(); i++) {                          for (index_t i=0; i<rowIndex.size(); i++) {
2949                              if (rowIndex[i]<getNumDOF()) {                              if (rowIndex[i]<getNumDOF()) {
2950                                  F_p[rowIndex[i]]+=EM_F[i];                                  for (index_t eq=0; eq<numEq; eq++) {
2951                                  //cout << "F[" << rowIndex[i] << "]=" << F_p[rowIndex[i]] << endl;                                      F_p[INDEX2(eq,rowIndex[i],numEq)]+=EM_F[INDEX2(eq,i,numEq)];
2952                                        //cout << "F[" << INDEX2(eq,rowIndex[i],numEq) << "]=" << F_p[INDEX2(eq,rowIndex[i],numEq)] << endl;
2953                                    }
2954                              }                              }
2955                          }                          }
2956                          //cout << "---"<<endl;                          //cout << "---"<<endl;
# Line 2960  void Rectangle::assemblePDESystem(Paso_S Line 2962  void Rectangle::assemblePDESystem(Paso_S
2962    
2963                  } // end k0 loop                  } // end k0 loop
2964              } // end k1 loop              } // end k1 loop
2965          } // end of coloring          } // end of colouring
2966      } // end of parallel region      } // end of parallel region
2967  }  }
2968    

Legend:
Removed from v.3760  
changed lines
  Added in v.3761

  ViewVC Help
Powered by ViewVC 1.1.26