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

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

Parent Directory Parent Directory | Revision Log Revision Log


Revision 3220 - (show annotations)
Wed Sep 29 00:33:16 2010 UTC (8 years, 6 months ago) by jfenwick
File MIME type: text/plain
File size: 5227 byte(s)
Removing references to Quadrature.? and ShapeFns

1
2 /*******************************************************
3 *
4 * Copyright (c) 2003-2010 by University of Queensland
5 * 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
14
15 /**************************************************************/
16
17 /* assemblage routines: header file */
18
19 /**************************************************************/
20
21 #ifndef INC_DUDLEY_ASSEMBLE
22 #define INC_DUDLEY_ASSEMBLE
23
24 /**************************************************************/
25
26 //#include "ReferenceElements.h"
27 #include "Dudley.h"
28 #include "ElementFile.h"
29 #include "NodeFile.h"
30 #include "escript/DataC.h"
31 #include "paso/SystemMatrix.h"
32
33 struct Assemble_Parameters {
34 dim_t numQuad; /* number of quadrature nodes */
35 dim_t numDim; /* spatial dimension */
36 dim_t NN; /* leading dimension of element node table */
37 dim_t numElements; /* number of elements */
38
39 dim_t numEqu;
40 index_t* row_DOF;
41 dim_t row_DOF_UpperBound;
42 Dudley_ElementFile_Jacobeans* row_jac;
43 // dim_t row_numShapesTotal;
44 // dim_t row_numShapes;
45 dim_t numShapes;
46
47 dim_t numComp;
48 index_t * col_DOF;
49 dim_t col_DOF_UpperBound;
50
51 const double* shapeFns;
52 bool_t reducedIntegrationOrder; // can we get rid of this one?
53
54 };
55
56 typedef struct Assemble_Parameters Assemble_Parameters;
57
58
59 #define Dudley_Assemble_reducedIntegrationOrder(__in__) ( (getFunctionSpaceType(__in__) == DUDLEY_REDUCED_ELEMENTS) || (getFunctionSpaceType(__in__) == DUDLEY_REDUCED_FACE_ELEMENTS) )
60
61 void Dudley_Assemble_PDE(Dudley_NodeFile*,Dudley_ElementFile*,Paso_SystemMatrix*,escriptDataC*,
62 escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*) ;
63 void Dudley_Assemble_LumpedSystem(Dudley_NodeFile*,Dudley_ElementFile*, escriptDataC*, escriptDataC*);
64
65 void Assemble_getAssembleParameters(Dudley_NodeFile*,Dudley_ElementFile*,Paso_SystemMatrix*,escriptDataC*,bool_t, Assemble_Parameters*);
66 void Dudley_Assemble_PDE_System2_3D(Assemble_Parameters, Dudley_ElementFile*,Paso_SystemMatrix*,escriptDataC*,
67 escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*);
68 void Dudley_Assemble_PDE_System2_2D(Assemble_Parameters, Dudley_ElementFile*,Paso_SystemMatrix*,escriptDataC*,
69 escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*);
70 void Dudley_Assemble_PDE_System2_1D(Assemble_Parameters, Dudley_ElementFile*,Paso_SystemMatrix*,escriptDataC*,
71 escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*);
72
73 void Dudley_Assemble_PDE_Single2_3D(Assemble_Parameters, Dudley_ElementFile*,Paso_SystemMatrix*,escriptDataC*,
74 escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*);
75 void Dudley_Assemble_PDE_Single2_2D(Assemble_Parameters, Dudley_ElementFile*,Paso_SystemMatrix*,escriptDataC*,
76 escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*);
77 void Dudley_Assemble_PDE_Single2_1D(Assemble_Parameters, Dudley_ElementFile*,Paso_SystemMatrix*,escriptDataC*,
78 escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*);
79
80
81 void Dudley_Assemble_NodeCoordinates(Dudley_NodeFile*,escriptDataC*);
82 void Dudley_Assemble_setNormal(Dudley_NodeFile*, Dudley_ElementFile*, escriptDataC*);
83 void Dudley_Assemble_interpolate(Dudley_NodeFile*,Dudley_ElementFile*,escriptDataC*, escriptDataC*);
84 void Dudley_Assemble_gradient(Dudley_NodeFile*, Dudley_ElementFile*,escriptDataC*, escriptDataC*);
85 void Dudley_Assemble_integrate(Dudley_NodeFile*,Dudley_ElementFile*,escriptDataC*,double*) ;
86 void Dudley_Assemble_getSize(Dudley_NodeFile*,Dudley_ElementFile*, escriptDataC*);
87 void Dudley_Assemble_CopyNodalData(Dudley_NodeFile* nodes,escriptDataC* out,escriptDataC* in);
88 void Dudley_Assemble_CopyElementData(Dudley_ElementFile* elements,escriptDataC* out,escriptDataC* in);
89 void Dudley_Assemble_AverageElementData(Dudley_ElementFile* elements,escriptDataC* out,escriptDataC* in);
90 void Dudley_Assemble_addToSystemMatrix(Paso_SystemMatrix*,dim_t,index_t*, dim_t,dim_t,index_t*,dim_t, double*);
91
92 void Assemble_jacobeans_2D(double*, dim_t, dim_t, dim_t, index_t*, double*, double* abs_D, double* quadweight, index_t*);
93 void Assemble_jacobeans_2D_M1D_E1D(double*, dim_t, dim_t, dim_t, index_t*, double*, double* abs_D, double* quadweight,index_t*);
94 void Assemble_jacobeans_3D(double*, dim_t, dim_t, dim_t, index_t*, double*, double* abs_D, double* quadweight,index_t*);
95 void Assemble_jacobeans_3D_M2D_E2D(double*, dim_t, dim_t, dim_t, index_t*, double*, double* abs_D, double* quadweight,index_t*);
96
97
98 void Dudley_Assemble_LumpedSystem(Dudley_NodeFile* nodes,Dudley_ElementFile* elements, escriptDataC* lumpedMat, escriptDataC* D);
99 #endif /* #ifndef INC_DUDLEY_ASSEMBLE */
100

Properties

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

  ViewVC Help
Powered by ViewVC 1.1.26