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

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

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

revision 3911 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 = 3  */  /*      X = 3  */
33  /*      Y = scalar   */  /*      Y = scalar   */
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_Single2_3D(Dudley_Assemble_Parameters p, Dudley_ElementFile * elements,  void Dudley_Assemble_PDE_Single2_3D(Dudley_Assemble_Parameters p, Dudley_ElementFile * elements,
46                      Paso_SystemMatrix * Mat, escriptDataC * F,                      Paso_SystemMatrix * Mat, escriptDataC * F,
# Line 101  void Dudley_Assemble_PDE_Single2_3D(Dudl Line 103  void Dudley_Assemble_PDE_Single2_3D(Dudl
103              add_EM_F = FALSE;              add_EM_F = FALSE;
104              add_EM_S = FALSE;              add_EM_S = FALSE;
105    
106                /**************************************************************/                /************************************************************************************/
107              /*   process A: */              /*   process A: */
108                /**************************************************************/                /************************************************************************************/
109              if (NULL != A_p)              if (NULL != A_p)
110              {              {
111                  add_EM_S = TRUE;                  add_EM_S = TRUE;
# Line 193  void Dudley_Assemble_PDE_Single2_3D(Dudl Line 195  void Dudley_Assemble_PDE_Single2_3D(Dudl
195                  }                  }
196                  }                  }
197              }              }
198                /**************************************************************/                /************************************************************************************/
199              /*   process B: */              /*   process B: */
200                /**************************************************************/                /************************************************************************************/
201              if (NULL != B_p)              if (NULL != B_p)
202              {              {
203                  add_EM_S = TRUE;                  add_EM_S = TRUE;
# Line 242  void Dudley_Assemble_PDE_Single2_3D(Dudl Line 244  void Dudley_Assemble_PDE_Single2_3D(Dudl
244                  }                  }
245                  }                  }
246              }              }
247                /**************************************************************/                /************************************************************************************/
248              /*   process C: */              /*   process C: */
249                /**************************************************************/                /************************************************************************************/
250              if (NULL != C_p)              if (NULL != C_p)
251              {              {
252                  add_EM_S = TRUE;                  add_EM_S = TRUE;
# Line 291  void Dudley_Assemble_PDE_Single2_3D(Dudl Line 293  void Dudley_Assemble_PDE_Single2_3D(Dudl
293                  }                  }
294                  }                  }
295              }              }
296                /************************************************************* */                /*********************************************************************************** */
297              /* process D */              /* process D */
298                /**************************************************************/                /************************************************************************************/
299              if (NULL != D_p)              if (NULL != D_p)
300              {              {
301                  add_EM_S = TRUE;                  add_EM_S = TRUE;
# Line 327  void Dudley_Assemble_PDE_Single2_3D(Dudl Line 329  void Dudley_Assemble_PDE_Single2_3D(Dudl
329                  }                  }
330                  }                  }
331              }              }
332                /**************************************************************/                /************************************************************************************/
333              /*   process X: */              /*   process X: */
334                /**************************************************************/                /************************************************************************************/
335              if (NULL != X_p)              if (NULL != X_p)
336              {              {
337                  add_EM_F = TRUE;                  add_EM_F = TRUE;
# Line 368  void Dudley_Assemble_PDE_Single2_3D(Dudl Line 370  void Dudley_Assemble_PDE_Single2_3D(Dudl
370                  }                  }
371                  }                  }
372              }              }
373               /**************************************************************/               /************************************************************************************/
374              /*   process Y: */              /*   process Y: */
375               /**************************************************************/               /************************************************************************************/
376              if (NULL != Y_p)              if (NULL != Y_p)
377              {              {
378                  add_EM_F = TRUE;                  add_EM_F = TRUE;
# Line 396  void Dudley_Assemble_PDE_Single2_3D(Dudl Line 398  void Dudley_Assemble_PDE_Single2_3D(Dudl
398                  }                  }
399                  }                  }
400              }              }
401                 /***********************************************************************************************/                 /*********************************************************************************************************************/
402              /* add the element matrices onto the matrix and right hand side                                */              /* add the element matrices onto the matrix and right hand side                                */
403                 /***********************************************************************************************/                 /*********************************************************************************************************************/
404    
405              for (q = 0; q < p.numShapes; q++)              for (q = 0; q < p.numShapes; q++)
406                  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)]];

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

  ViewVC Help
Powered by ViewVC 1.1.26