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

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

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

trunk/finley/src/Assemble_PDE_Single2_2D.c revision 2881 by jfenwick, Thu Jan 28 02:03:15 2010 UTC branches/domexper/dudley/src/Assemble_PDE_Single2_2D.c revision 3086 by jfenwick, Thu Aug 5 05:07:58 2010 UTC
# Line 44  Line 44 
44    
45  /**************************************************************/  /**************************************************************/
46    
47  void  Finley_Assemble_PDE_Single2_2D(Assemble_Parameters p, Finley_ElementFile* elements,  void  Dudley_Assemble_PDE_Single2_2D(Assemble_Parameters p, Dudley_ElementFile* elements,
48                                       Paso_SystemMatrix* Mat, escriptDataC* F,                                       Paso_SystemMatrix* Mat, escriptDataC* F,
49                                       escriptDataC* A, escriptDataC* B, escriptDataC* C, escriptDataC* D, escriptDataC* X, escriptDataC* Y) {                                       escriptDataC* A, escriptDataC* B, escriptDataC* C, escriptDataC* D, escriptDataC* X, escriptDataC* Y) {
50    
# Line 75  void  Finley_Assemble_PDE_Single2_2D(Ass Line 75  void  Finley_Assemble_PDE_Single2_2D(Ass
75         EM_F=THREAD_MEMALLOC(len_EM_F,double);         EM_F=THREAD_MEMALLOC(len_EM_F,double);
76         row_index=THREAD_MEMALLOC(p.row_numShapesTotal,index_t);         row_index=THREAD_MEMALLOC(p.row_numShapesTotal,index_t);
77                                                                                                                                                                                                                                                                                                                                                                                                            
78         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) ) {
79    
80            for (color=elements->minColor;color<=elements->maxColor;color++) {            for (color=elements->minColor;color<=elements->maxColor;color++) {
81               /*  open loop over all elements: */               /*  open loop over all elements: */
# Line 280  void  Finley_Assemble_PDE_Single2_2D(Ass Line 280  void  Finley_Assemble_PDE_Single2_2D(Ass
280                        /***********************************************************************************************/                        /***********************************************************************************************/
281                                
282                        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)]];
283                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);
284                        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);
285    
286                    } /* end of isub */                    } /* end of isub */
287                  } /* end color check */                  } /* end color check */

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

  ViewVC Help
Powered by ViewVC 1.1.26