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

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

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

revision 3143 by jfenwick, Thu Sep 2 05:03:17 2010 UTC revision 3144 by jfenwick, Fri Sep 3 00:49:02 2010 UTC
# Line 153  void Dudley_Assemble_LumpedSystem(Dudley Line 153  void Dudley_Assemble_LumpedSystem(Dudley
153                                 #endif                                 #endif
154                                 for (q=0;q<p.row_numShapesTotal;q++)                                 for (q=0;q<p.row_numShapesTotal;q++)
155                     {                     {
156                   row_index[q]=p.row_DOF[elements->Nodes[INDEX2(p.row_node[INDEX2(q,0, p.row_numShapesTotal)],e,p.NN)]];                   row_index[q]=p.row_DOF[elements->Nodes[INDEX2(q,e,p.NN)]];
157                     }                     }
158                                 Dudley_Util_AddScatter(p.row_numShapesTotal,row_index,p.numEqu,EM_lumpedMat,lumpedMat_p, p.row_DOF_UpperBound);                                 Dudley_Util_AddScatter(p.row_numShapesTotal,row_index,p.numEqu,EM_lumpedMat,lumpedMat_p, p.row_DOF_UpperBound);
159                         } /* end color check */                                 } /* end color check */        
# Line 191  void Dudley_Assemble_LumpedSystem(Dudley Line 191  void Dudley_Assemble_LumpedSystem(Dudley
191                                         EM_lumpedMat[INDEX2(0,s,p.numEqu)]=rtmp*D_p[0];                                         EM_lumpedMat[INDEX2(0,s,p.numEqu)]=rtmp*D_p[0];
192                                     }                                     }
193                                 #endif                                 #endif
194                                 for (q=0;q<p.row_numShapesTotal;q++) row_index[q]=p.row_DOF[elements->Nodes[INDEX2(p.row_node[INDEX2(q, 0, p.row_numShapesTotal)],e,p.NN)]];                                 for (q=0;q<p.row_numShapesTotal;q++) row_index[q]=p.row_DOF[elements->Nodes[INDEX2(q,e,p.NN)]];
195                                 Dudley_Util_AddScatter(p.row_numShapesTotal,row_index,p.numEqu,EM_lumpedMat,lumpedMat_p, p.row_DOF_UpperBound);                                 Dudley_Util_AddScatter(p.row_numShapesTotal,row_index,p.numEqu,EM_lumpedMat,lumpedMat_p, p.row_DOF_UpperBound);
196                         } /* end color check */                         } /* end color check */
197                      } /* end element loop */                      } /* end element loop */
# Line 233  void Dudley_Assemble_LumpedSystem(Dudley Line 233  void Dudley_Assemble_LumpedSystem(Dudley
233                                         }                                         }
234                    }                    }
235                    #endif                    #endif
236                    for (q=0;q<p.row_numShapesTotal;q++) row_index[q]=p.row_DOF[elements->Nodes[INDEX2(p.row_node[INDEX2(q,0,p.row_numShapesTotal)],e,p.NN)]];                    for (q=0;q<p.row_numShapesTotal;q++) row_index[q]=p.row_DOF[elements->Nodes[INDEX2(q,e,p.NN)]];
237                    Dudley_Util_AddScatter(p.row_numShapesTotal,row_index,p.numEqu,EM_lumpedMat,lumpedMat_p, p.row_DOF_UpperBound);                    Dudley_Util_AddScatter(p.row_numShapesTotal,row_index,p.numEqu,EM_lumpedMat,lumpedMat_p, p.row_DOF_UpperBound);
238                         } /* end color check */                         } /* end color check */
239                      } /* end element loop */                      } /* end element loop */
# Line 271  void Dudley_Assemble_LumpedSystem(Dudley Line 271  void Dudley_Assemble_LumpedSystem(Dudley
271                                       }                                       }
272                   }                   }
273                    #endif                    #endif
274                    for (q=0;q<p.row_numShapesTotal;q++) row_index[q]=p.row_DOF[elements->Nodes[INDEX2(p.row_node[INDEX2(q,0,p.row_numShapesTotal)],e,p.NN)]];                    for (q=0;q<p.row_numShapesTotal;q++) row_index[q]=p.row_DOF[elements->Nodes[INDEX2(q,e,p.NN)]];
275                    Dudley_Util_AddScatter(p.row_numShapesTotal,row_index,p.numEqu,EM_lumpedMat,lumpedMat_p, p.row_DOF_UpperBound);                    Dudley_Util_AddScatter(p.row_numShapesTotal,row_index,p.numEqu,EM_lumpedMat,lumpedMat_p, p.row_DOF_UpperBound);
276                         } /* end color check */                         } /* end color check */
277                      } /* end element loop */                      } /* end element loop */

Legend:
Removed from v.3143  
changed lines
  Added in v.3144

  ViewVC Help
Powered by ViewVC 1.1.26