/[escript]/branches/domexper/dudley/src/Assemble.h
ViewVC logotype

Annotation of /branches/domexper/dudley/src/Assemble.h

Parent Directory Parent Directory | Revision Log Revision Log


Revision 3202 - (hide annotations)
Thu Sep 23 23:24:09 2010 UTC (8 years, 11 months ago) by jfenwick
File MIME type: text/plain
File size: 5307 byte(s)


1 jgs 82
2 ksteube 1312 /*******************************************************
3 ksteube 1811 *
4 jfenwick 2881 * Copyright (c) 2003-2010 by University of Queensland
5 ksteube 1811 * Earth Systems Science Computational Center (ESSCC)
6     * http://www.uq.edu.au/esscc
7     *
8     * Primary Business: Queensland, Australia
9     * Licensed under the Open Software License version 3.0
10     * http://www.opensource.org/licenses/osl-3.0.php
11     *
12     *******************************************************/
13 jgs 82
14 ksteube 1811
15 jgs 82 /**************************************************************/
16    
17 ksteube 1312 /* assemblage routines: header file */
18 jgs 82
19     /**************************************************************/
20 jgs 150
21 jfenwick 3082 #ifndef INC_DUDLEY_ASSEMBLE
22     #define INC_DUDLEY_ASSEMBLE
23 jgs 150
24     /**************************************************************/
25    
26 jgs 82 #include "ReferenceElements.h"
27 jfenwick 3082 #include "Dudley.h"
28 jgs 82 #include "ElementFile.h"
29     #include "NodeFile.h"
30 robwdcock 682 #include "escript/DataC.h"
31     #include "paso/SystemMatrix.h"
32 jgs 82
33     struct Assemble_Parameters {
34 jfenwick 3137 dim_t numQuadTotal; /* total number of quadrature nodes */
35 gross 2748 dim_t numDim; /* spatial dimension */
36     dim_t NN; /* leading dimension of element node table */
37     dim_t numElements; /* number of elements */
38 jgs 82
39 jgs 123 dim_t numEqu;
40 gross 798 index_t* row_DOF;
41     dim_t row_DOF_UpperBound;
42 jfenwick 3086 Dudley_ElementFile_Jacobeans* row_jac;
43 gross 2748 dim_t row_numShapesTotal;
44     dim_t row_numShapes;
45    
46 jgs 123 dim_t numComp;
47 gross 798 index_t * col_DOF;
48     dim_t col_DOF_UpperBound;
49 jfenwick 3086 Dudley_ElementFile_Jacobeans* col_jac;
50 gross 2748 dim_t col_numShapesTotal;
51     dim_t col_numShapes;
52 jfenwick 3202
53     const double* shapeFns;
54     bool_t reducedIntegrationOrder; // can we get rid of this one?
55    
56 jgs 82 };
57    
58     typedef struct Assemble_Parameters Assemble_Parameters;
59    
60    
61 jfenwick 3090 #define Dudley_Assemble_reducedIntegrationOrder(__in__) ( (getFunctionSpaceType(__in__) == DUDLEY_REDUCED_ELEMENTS) || (getFunctionSpaceType(__in__) == DUDLEY_REDUCED_FACE_ELEMENTS) )
62 gross 1062
63 jfenwick 3086 void Dudley_Assemble_PDE(Dudley_NodeFile*,Dudley_ElementFile*,Paso_SystemMatrix*,escriptDataC*,
64 jgs 82 escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*) ;
65 jfenwick 3086 void Dudley_Assemble_LumpedSystem(Dudley_NodeFile*,Dudley_ElementFile*, escriptDataC*, escriptDataC*);
66 gross 1204
67 jfenwick 3086 void Assemble_getAssembleParameters(Dudley_NodeFile*,Dudley_ElementFile*,Paso_SystemMatrix*,escriptDataC*,bool_t, Assemble_Parameters*);
68     void Dudley_Assemble_PDE_System2_3D(Assemble_Parameters, Dudley_ElementFile*,Paso_SystemMatrix*,escriptDataC*,
69 gross 798 escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*);
70 jfenwick 3086 void Dudley_Assemble_PDE_System2_2D(Assemble_Parameters, Dudley_ElementFile*,Paso_SystemMatrix*,escriptDataC*,
71 gross 798 escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*);
72 jfenwick 3086 void Dudley_Assemble_PDE_System2_1D(Assemble_Parameters, Dudley_ElementFile*,Paso_SystemMatrix*,escriptDataC*,
73 gross 798 escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*);
74 jfenwick 3156
75 jfenwick 3086 void Dudley_Assemble_PDE_Single2_3D(Assemble_Parameters, Dudley_ElementFile*,Paso_SystemMatrix*,escriptDataC*,
76 gross 798 escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*);
77 jfenwick 3086 void Dudley_Assemble_PDE_Single2_2D(Assemble_Parameters, Dudley_ElementFile*,Paso_SystemMatrix*,escriptDataC*,
78 gross 798 escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*);
79 jfenwick 3086 void Dudley_Assemble_PDE_Single2_1D(Assemble_Parameters, Dudley_ElementFile*,Paso_SystemMatrix*,escriptDataC*,
80 gross 798 escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*);
81    
82    
83 jfenwick 3086 void Dudley_Assemble_NodeCoordinates(Dudley_NodeFile*,escriptDataC*);
84     void Dudley_Assemble_setNormal(Dudley_NodeFile*, Dudley_ElementFile*, escriptDataC*);
85     void Dudley_Assemble_interpolate(Dudley_NodeFile*,Dudley_ElementFile*,escriptDataC*, escriptDataC*);
86     void Dudley_Assemble_gradient(Dudley_NodeFile*, Dudley_ElementFile*,escriptDataC*, escriptDataC*);
87     void Dudley_Assemble_integrate(Dudley_NodeFile*,Dudley_ElementFile*,escriptDataC*,double*) ;
88     void Dudley_Assemble_getSize(Dudley_NodeFile*,Dudley_ElementFile*, escriptDataC*);
89     void Dudley_Assemble_CopyNodalData(Dudley_NodeFile* nodes,escriptDataC* out,escriptDataC* in);
90     void Dudley_Assemble_CopyElementData(Dudley_ElementFile* elements,escriptDataC* out,escriptDataC* in);
91     void Dudley_Assemble_AverageElementData(Dudley_ElementFile* elements,escriptDataC* out,escriptDataC* in);
92     void Dudley_Assemble_addToSystemMatrix(Paso_SystemMatrix*,dim_t,index_t*, dim_t,dim_t,index_t*,dim_t, double*);
93 gross 776
94 jfenwick 3184 void Assemble_jacobeans_2D(double*, dim_t, dim_t, dim_t, index_t*, double*, double* abs_D, double* quadweight, index_t*);
95     void Assemble_jacobeans_2D_M1D_E1D(double*, dim_t, dim_t, dim_t, index_t*, double*, double* abs_D, double* quadweight,index_t*);
96     void Assemble_jacobeans_3D(double*, dim_t, dim_t, dim_t, index_t*, double*, double* abs_D, double* quadweight,index_t*);
97     void Assemble_jacobeans_3D_M2D_E2D(double*, dim_t, dim_t, dim_t, index_t*, double*, double* abs_D, double* quadweight,index_t*);
98 ksteube 1312
99 gross 2748
100 jfenwick 3086 void Dudley_Assemble_LumpedSystem(Dudley_NodeFile* nodes,Dudley_ElementFile* elements, escriptDataC* lumpedMat, escriptDataC* D);
101     #endif /* #ifndef INC_DUDLEY_ASSEMBLE */
102 jgs 82

Properties

Name Value
svn:eol-style native
svn:keywords Author Date Id Revision

  ViewVC Help
Powered by ViewVC 1.1.26