/[escript]/branches/doubleplusgood/dudley/src/Assemble_PDE_System2_3D.cpp
ViewVC logotype

Diff of /branches/doubleplusgood/dudley/src/Assemble_PDE_System2_3D.cpp

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

revision 3980 by jfenwick, Thu Jun 14 01:01:03 2012 UTC revision 3981 by jfenwick, Fri Sep 21 02:47:54 2012 UTC
# Line 1  Line 1 
1    
2  /*******************************************************  /*****************************************************************************
3  *  *
4  * Copyright (c) 2003-2012 by University of Queensland  * Copyright (c) 2003-2012 by University of Queensland
5  * Earth Systems Science Computational Center (ESSCC)  * http://www.uq.edu.au
 * http://www.uq.edu.au/esscc  
6  *  *
7  * Primary Business: Queensland, Australia  * Primary Business: Queensland, Australia
8  * Licensed under the Open Software License version 3.0  * Licensed under the Open Software License version 3.0
9  * http://www.opensource.org/licenses/osl-3.0.php  * http://www.opensource.org/licenses/osl-3.0.php
10  *  *
11  *******************************************************/  * Development until 2012 by Earth Systems Science Computational Center (ESSCC)
12    * Development since 2012 by School of Earth Sciences
13    *
14    *****************************************************************************/
15    
16  /**************************************************************/  /************************************************************************************/
17    
18  /*    assembles the system of numEq PDEs into the stiffness matrix S right hand side F  */  /*    assembles the system of numEq PDEs into the stiffness matrix S right hand side F  */
19  /*    the shape functions for test and solution must be identical */  /*    the shape functions for test and solution must be identical */
# Line 30  Line 32 
32  /*      X = p.numEqu x 3  */  /*      X = p.numEqu x 3  */
33  /*      Y = p.numEqu   */  /*      Y = p.numEqu   */
34    
35  /**************************************************************/  /************************************************************************************/
36    
37  #include "Assemble.h"  #include "Assemble.h"
38  #include "Util.h"  #include "Util.h"
# Line 38  Line 40 
40  #include <omp.h>  #include <omp.h>
41  #endif  #endif
42    
43  /**************************************************************/  /************************************************************************************/
44    
45  void Dudley_Assemble_PDE_System2_3D(Dudley_Assemble_Parameters p, Dudley_ElementFile * elements,  void Dudley_Assemble_PDE_System2_3D(Dudley_Assemble_Parameters p, Dudley_ElementFile * elements,
46                      Paso_SystemMatrix * Mat, escriptDataC * F,                      Paso_SystemMatrix * Mat, escriptDataC * F,
# Line 102  void Dudley_Assemble_PDE_System2_3D(Dudl Line 104  void Dudley_Assemble_PDE_System2_3D(Dudl
104              add_EM_F = FALSE;              add_EM_F = FALSE;
105              add_EM_S = FALSE;              add_EM_S = FALSE;
106    
107                /**************************************************************/                /************************************************************************************/
108              /*   process A: */              /*   process A: */
109                /**************************************************************/                /************************************************************************************/
110              if (NULL != A_p)              if (NULL != A_p)
111              {              {
112                  add_EM_S = TRUE;                  add_EM_S = TRUE;
# Line 210  void Dudley_Assemble_PDE_System2_3D(Dudl Line 212  void Dudley_Assemble_PDE_System2_3D(Dudl
212                  }                  }
213                  }                  }
214              }              }
215                /**************************************************************/                /************************************************************************************/
216              /*   process B: */              /*   process B: */
217                /**************************************************************/                /************************************************************************************/
218              if (NULL != B_p)              if (NULL != B_p)
219              {              {
220                  add_EM_S = TRUE;                  add_EM_S = TRUE;
# Line 275  void Dudley_Assemble_PDE_System2_3D(Dudl Line 277  void Dudley_Assemble_PDE_System2_3D(Dudl
277                  }                  }
278                  }                  }
279              }              }
280                /**************************************************************/                /************************************************************************************/
281              /*   process C: */              /*   process C: */
282                /**************************************************************/                /************************************************************************************/
283              if (NULL != C_p)              if (NULL != C_p)
284              {              {
285                  add_EM_S = TRUE;                  add_EM_S = TRUE;
# Line 340  void Dudley_Assemble_PDE_System2_3D(Dudl Line 342  void Dudley_Assemble_PDE_System2_3D(Dudl
342                  }                  }
343                  }                  }
344              }              }
345                /************************************************************* */                /*********************************************************************************** */
346              /* process D */              /* process D */
347                /**************************************************************/                /************************************************************************************/
348              if (NULL != D_p)              if (NULL != D_p)
349              {              {
350                  add_EM_S = TRUE;                  add_EM_S = TRUE;
# Line 392  void Dudley_Assemble_PDE_System2_3D(Dudl Line 394  void Dudley_Assemble_PDE_System2_3D(Dudl
394                  }                  }
395                  }                  }
396              }              }
397                /**************************************************************/                /************************************************************************************/
398              /*   process X: */              /*   process X: */
399                /**************************************************************/                /************************************************************************************/
400              if (NULL != X_p)              if (NULL != X_p)
401              {              {
402                  add_EM_F = TRUE;                  add_EM_F = TRUE;
# Line 441  void Dudley_Assemble_PDE_System2_3D(Dudl Line 443  void Dudley_Assemble_PDE_System2_3D(Dudl
443                  }                  }
444                  }                  }
445              }              }
446               /**************************************************************/               /************************************************************************************/
447              /*   process Y: */              /*   process Y: */
448               /**************************************************************/               /************************************************************************************/
449              if (NULL != Y_p)              if (NULL != Y_p)
450              {              {
451                  add_EM_F = TRUE;                  add_EM_F = TRUE;
# Line 474  void Dudley_Assemble_PDE_System2_3D(Dudl Line 476  void Dudley_Assemble_PDE_System2_3D(Dudl
476                  }                  }
477              }              }
478    
479                 /***********************************************************************************************/                 /*********************************************************************************************************************/
480              /* add the element matrices onto the matrix and right hand side                                */              /* add the element matrices onto the matrix and right hand side                                */
481                 /***********************************************************************************************/                 /*********************************************************************************************************************/
482              for (q = 0; q < p.numShapes; q++)              for (q = 0; q < p.numShapes; q++)
483                  row_index[q] = p.row_DOF[elements->Nodes[INDEX2(q, e, p.NN)]];                  row_index[q] = p.row_DOF[elements->Nodes[INDEX2(q, e, p.NN)]];
484    

Legend:
Removed from v.3980  
changed lines
  Added in v.3981

  ViewVC Help
Powered by ViewVC 1.1.26