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

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

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

trunk/dudley/src/Assemble_PDE_Single2_1D.c revision 3911 by jfenwick, Thu Jun 14 01:01:03 2012 UTC branches/doubleplusgood/dudley/src/Assemble_PDE_Single2_1D.cpp revision 4261 by jfenwick, Wed Feb 27 06:09:33 2013 UTC
# Line 1  Line 1 
1    
2  /*******************************************************  /*****************************************************************************
3  *  *
4  * Copyright (c) 2003-2012 by University of Queensland  * Copyright (c) 2003-2013 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 = 1  */  /*      X = 1  */
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_1D(Dudley_Assemble_Parameters p, Dudley_ElementFile * elements,  void Dudley_Assemble_PDE_Single2_1D(Dudley_Assemble_Parameters p, Dudley_ElementFile * elements,
46                      Paso_SystemMatrix * Mat, escriptDataC * F,                      Paso_SystemMatrix * Mat, escriptDataC * F,
# Line 100  void Dudley_Assemble_PDE_Single2_1D(Dudl Line 102  void Dudley_Assemble_PDE_Single2_1D(Dudl
102                  EM_F[q] = 0;                  EM_F[q] = 0;
103              add_EM_F = FALSE;              add_EM_F = FALSE;
104              add_EM_S = FALSE;              add_EM_S = FALSE;
105                /**************************************************************/                /************************************************************************************/
106              /*   process A: */              /*   process A: */
107                /**************************************************************/                /************************************************************************************/
108              if (NULL != A_p)              if (NULL != A_p)
109              {              {
110                  add_EM_S = TRUE;                  add_EM_S = TRUE;
# Line 142  void Dudley_Assemble_PDE_Single2_1D(Dudl Line 144  void Dudley_Assemble_PDE_Single2_1D(Dudl
144                  }                  }
145                  }                  }
146              }              }
147                /**************************************************************/                /************************************************************************************/
148              /*   process B: */              /*   process B: */
149                /**************************************************************/                /************************************************************************************/
150              if (NULL != B_p)              if (NULL != B_p)
151              {              {
152                  add_EM_S = TRUE;                  add_EM_S = TRUE;
# Line 183  void Dudley_Assemble_PDE_Single2_1D(Dudl Line 185  void Dudley_Assemble_PDE_Single2_1D(Dudl
185                  }                  }
186                  }                  }
187              }              }
188                /**************************************************************/                /************************************************************************************/
189              /*   process C: */              /*   process C: */
190                /**************************************************************/                /************************************************************************************/
191              if (NULL != C_p)              if (NULL != C_p)
192              {              {
193                  add_EM_S = TRUE;                  add_EM_S = TRUE;
# Line 224  void Dudley_Assemble_PDE_Single2_1D(Dudl Line 226  void Dudley_Assemble_PDE_Single2_1D(Dudl
226                  }                  }
227                  }                  }
228              }              }
229                /************************************************************* */                /*********************************************************************************** */
230              /* process D */              /* process D */
231                /**************************************************************/                /************************************************************************************/
232              if (NULL != D_p)              if (NULL != D_p)
233              {              {
234                  add_EM_S = TRUE;                  add_EM_S = TRUE;
# Line 265  void Dudley_Assemble_PDE_Single2_1D(Dudl Line 267  void Dudley_Assemble_PDE_Single2_1D(Dudl
267                  }                  }
268                  }                  }
269              }              }
270                /**************************************************************/                /************************************************************************************/
271              /*   process X: */              /*   process X: */
272                /**************************************************************/                /************************************************************************************/
273              if (NULL != X_p)              if (NULL != X_p)
274              {              {
275                  add_EM_F = TRUE;                  add_EM_F = TRUE;
# Line 295  void Dudley_Assemble_PDE_Single2_1D(Dudl Line 297  void Dudley_Assemble_PDE_Single2_1D(Dudl
297                  }                  }
298                  }                  }
299              }              }
300               /**************************************************************/               /************************************************************************************/
301              /*   process Y: */              /*   process Y: */
302               /**************************************************************/               /************************************************************************************/
303              if (NULL != Y_p)              if (NULL != Y_p)
304              {              {
305                  add_EM_F = TRUE;                  add_EM_F = TRUE;
# Line 323  void Dudley_Assemble_PDE_Single2_1D(Dudl Line 325  void Dudley_Assemble_PDE_Single2_1D(Dudl
325                  }                  }
326                  }                  }
327              }              }
328                 /***********************************************************************************************/                 /*********************************************************************************************************************/
329              /* add the element matrices onto the matrix and right hand side                                */              /* add the element matrices onto the matrix and right hand side                                */
330                 /***********************************************************************************************/                 /*********************************************************************************************************************/
331              for (q = 0; q < p.row_numShapesTotal; q++)              for (q = 0; q < p.row_numShapesTotal; q++)
332                  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)]];
333    

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

  ViewVC Help
Powered by ViewVC 1.1.26