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

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

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

revision 3085 by jfenwick, Tue Aug 3 04:28:03 2010 UTC revision 3086 by jfenwick, Thu Aug 5 05:07:58 2010 UTC
# Line 45  Line 45 
45    
46  /**************************************************************/  /**************************************************************/
47    
48  void  Finley_Assemble_PDE_System2_2D(Assemble_Parameters p, Finley_ElementFile* elements,  void  Dudley_Assemble_PDE_System2_2D(Assemble_Parameters p, Dudley_ElementFile* elements,
49                                       Paso_SystemMatrix* Mat, escriptDataC* F,                                       Paso_SystemMatrix* Mat, escriptDataC* F,
50                                       escriptDataC* A, escriptDataC* B, escriptDataC* C, escriptDataC* D, escriptDataC* X, escriptDataC* Y) {                                       escriptDataC* A, escriptDataC* B, escriptDataC* C, escriptDataC* D, escriptDataC* X, escriptDataC* Y) {
51    
# Line 77  void  Finley_Assemble_PDE_System2_2D(Ass Line 77  void  Finley_Assemble_PDE_System2_2D(Ass
77         EM_F=THREAD_MEMALLOC(len_EM_F,double);         EM_F=THREAD_MEMALLOC(len_EM_F,double);
78         row_index=THREAD_MEMALLOC(p.row_numShapesTotal,index_t);         row_index=THREAD_MEMALLOC(p.row_numShapesTotal,index_t);
79                                                                                                                                                                                                                                                                                                                                                                                                            
80         if (!Finley_checkPtr(EM_S) && !Finley_checkPtr(EM_F) && !Finley_checkPtr(row_index) ) {         if (!Dudley_checkPtr(EM_S) && !Dudley_checkPtr(EM_F) && !Dudley_checkPtr(row_index) ) {
81    
82            for (color=elements->minColor;color<=elements->maxColor;color++) {            for (color=elements->minColor;color<=elements->maxColor;color++) {
83               /*  open loop over all elements: */               /*  open loop over all elements: */
# Line 333  void  Finley_Assemble_PDE_System2_2D(Ass Line 333  void  Finley_Assemble_PDE_System2_2D(Ass
333                         /***********************************************************************************************/                         /***********************************************************************************************/
334                         for (q=0;q<p.row_numShapesTotal;q++) row_index[q]=p.row_DOF[elements->Nodes[INDEX2(p.row_node[INDEX2(q,isub,p.row_numShapesTotal)],e,p.NN)]];                         for (q=0;q<p.row_numShapesTotal;q++) row_index[q]=p.row_DOF[elements->Nodes[INDEX2(p.row_node[INDEX2(q,isub,p.row_numShapesTotal)],e,p.NN)]];
335    
336                         if (add_EM_F) Finley_Util_AddScatter(p.row_numShapesTotal,row_index,p.numEqu,EM_F,F_p, p.row_DOF_UpperBound);                         if (add_EM_F) Dudley_Util_AddScatter(p.row_numShapesTotal,row_index,p.numEqu,EM_F,F_p, p.row_DOF_UpperBound);
337                         if (add_EM_S) Finley_Assemble_addToSystemMatrix(Mat,p.row_numShapesTotal,row_index,p.numEqu,p.col_numShapesTotal,row_index,p.numComp,EM_S);                         if (add_EM_S) Dudley_Assemble_addToSystemMatrix(Mat,p.row_numShapesTotal,row_index,p.numEqu,p.col_numShapesTotal,row_index,p.numComp,EM_S);
338    
339                     } /* end of isub */                     } /* end of isub */
340                  } /* end color check */                  } /* end color check */

Legend:
Removed from v.3085  
changed lines
  Added in v.3086

  ViewVC Help
Powered by ViewVC 1.1.26