/[escript]/trunk/finley/src/Assemble_LumpedSystem.c
ViewVC logotype

Diff of /trunk/finley/src/Assemble_LumpedSystem.c

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

revision 1645 by phornby, Tue Jul 15 05:26:38 2008 UTC revision 1650 by phornby, Tue Jul 15 08:58:22 2008 UTC
# Line 112  void Finley_Assemble_LumpedSystem(Finley Line 112  void Finley_Assemble_LumpedSystem(Finley
112      expandedD=isExpanded(D);      expandedD=isExpanded(D);
113      S=p.row_jac->ReferenceElement->S;      S=p.row_jac->ReferenceElement->S;
114    
115      #pragma omp parallel private(color, EM_lumpedMat, row_index, Vol, D_p, s, q, k, rtmp, m_t, diagS)      #pragma omp parallel private(color, EM_lumpedMat, row_index, Vol, D_p, s, q, k, rtmp)
116      {      {
117         EM_lumpedMat=THREAD_MEMALLOC(len_EM_lumpedMat,double);         EM_lumpedMat=THREAD_MEMALLOC(len_EM_lumpedMat,double);
118         row_index=THREAD_MEMALLOC(p.row_NN,index_t);         row_index=THREAD_MEMALLOC(p.row_NN,index_t);
# Line 132  void Finley_Assemble_LumpedSystem(Finley Line 132  void Finley_Assemble_LumpedSystem(Finley
132                             *           Number of PDEs: 1                             *           Number of PDEs: 1
133                             *  D_p varies over element: True                             *  D_p varies over element: True
134                             */                             */
135                              #pragma omp parallel private(m_t)
136                            m_t=0; /* mass of the element: m_t */                            m_t=0; /* mass of the element: m_t */
137                            for (q=0;q<p.numQuad;q++) {                            for (q=0;q<p.numQuad;q++) {
138                                m_t+=Vol[q]*D_p[q];                                m_t+=Vol[q]*D_p[q];
139                            }                                                      }                          
140                              #pragma omp parallel private(diagS)
141                            diagS=0; /* diagonal sum: S */                            diagS=0; /* diagonal sum: S */
142                            for (s=0;s<p.row_NS;s++) {                            for (s=0;s<p.row_NS;s++) {
143                                rtmp=0;                                rtmp=0;

Legend:
Removed from v.1645  
changed lines
  Added in v.1650

  ViewVC Help
Powered by ViewVC 1.1.26