ESYS13
Revision_
|
00001 00002 /******************************************************* 00003 * 00004 * Copyright (c) 2003-2012 by University of Queensland 00005 * Earth Systems Science Computational Center (ESSCC) 00006 * http://www.uq.edu.au/esscc 00007 * 00008 * Primary Business: Queensland, Australia 00009 * Licensed under the Open Software License version 3.0 00010 * http://www.opensource.org/licenses/osl-3.0.php 00011 * 00012 *******************************************************/ 00013 00014 00015 /**************************************************************/ 00016 00017 /* assemblage routines: header file */ 00018 00019 /**************************************************************/ 00020 00021 #ifndef INC_FINLEY_ASSEMBLE 00022 #define INC_FINLEY_ASSEMBLE 00023 00024 /**************************************************************/ 00025 00026 #include "ReferenceElements.h" 00027 #include "Finley.h" 00028 #include "ElementFile.h" 00029 #include "NodeFile.h" 00030 #include "escript/DataC.h" 00031 #include "paso/SystemMatrix.h" 00032 00033 struct Finley_Assemble_Parameters { 00034 dim_t numQuadTotal; /* total number of quadrature nodes = numQuadSub * numQuadSub */ 00035 dim_t numQuadSub; /* number of quadrature nodes per subelements */ 00036 dim_t numSides; /* number of sides */ 00037 dim_t numSub; /* number of subelements */ 00038 dim_t numDim; /* spatial dimension */ 00039 dim_t NN; /* leading dimension of element node table */ 00040 dim_t numElements; /* number of elements */ 00041 00042 dim_t numEqu; 00043 index_t* row_DOF; 00044 dim_t row_DOF_UpperBound; 00045 Finley_ElementFile_Jacobeans* row_jac; 00046 index_t* row_node; 00047 dim_t row_numShapesTotal; 00048 dim_t row_numShapes; 00049 00050 dim_t numComp; 00051 index_t * col_DOF; 00052 dim_t col_DOF_UpperBound; 00053 Finley_ElementFile_Jacobeans* col_jac; 00054 index_t* col_node; 00055 dim_t col_numShapesTotal; 00056 dim_t col_numShapes; 00057 }; 00058 00059 typedef struct Finley_Assemble_Parameters Finley_Assemble_Parameters; 00060 00061 00062 #define Finley_Assemble_reducedIntegrationOrder(__in__) ( (getFunctionSpaceType(__in__) == FINLEY_REDUCED_ELEMENTS) || (getFunctionSpaceType(__in__) == FINLEY_REDUCED_FACE_ELEMENTS) || (getFunctionSpaceType(__in__) == FINLEY_REDUCED_CONTACT_ELEMENTS_1) || (getFunctionSpaceType(__in__) == FINLEY_REDUCED_CONTACT_ELEMENTS_2) ) 00063 00064 void Finley_Assemble_PDE(Finley_NodeFile*,Finley_ElementFile*,Paso_SystemMatrix*,escriptDataC*, 00065 escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*) ; 00066 00067 void Finley_Assemble_getAssembleParameters(Finley_NodeFile*,Finley_ElementFile*,Paso_SystemMatrix*,escriptDataC*,bool_t, Finley_Assemble_Parameters*); 00068 void Finley_Assemble_PDE_System2_3D(Finley_Assemble_Parameters, Finley_ElementFile*,Paso_SystemMatrix*,escriptDataC*, 00069 escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*); 00070 void Finley_Assemble_PDE_System2_2D(Finley_Assemble_Parameters, Finley_ElementFile*,Paso_SystemMatrix*,escriptDataC*, 00071 escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*); 00072 void Finley_Assemble_PDE_System2_1D(Finley_Assemble_Parameters, Finley_ElementFile*,Paso_SystemMatrix*,escriptDataC*, 00073 escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*); 00074 void Finley_Assemble_PDE_System2_C(Finley_Assemble_Parameters , Finley_ElementFile*, Paso_SystemMatrix*, escriptDataC*, escriptDataC*, escriptDataC*); 00075 void Finley_Assemble_PDE_Single2_3D(Finley_Assemble_Parameters, Finley_ElementFile*,Paso_SystemMatrix*,escriptDataC*, 00076 escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*); 00077 void Finley_Assemble_PDE_Single2_2D(Finley_Assemble_Parameters, Finley_ElementFile*,Paso_SystemMatrix*,escriptDataC*, 00078 escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*); 00079 void Finley_Assemble_PDE_Single2_1D(Finley_Assemble_Parameters, Finley_ElementFile*,Paso_SystemMatrix*,escriptDataC*, 00080 escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*, escriptDataC*); 00081 void Finley_Assemble_PDE_Single2_C(Finley_Assemble_Parameters p, Finley_ElementFile*, Paso_SystemMatrix*, escriptDataC*, escriptDataC*, escriptDataC*); 00082 void Finley_Assemble_PDE_Points(Finley_Assemble_Parameters p, Finley_ElementFile *, Paso_SystemMatrix *, escriptDataC *, escriptDataC *, escriptDataC *); 00083 00084 00085 void Finley_Assemble_NodeCoordinates(Finley_NodeFile*,escriptDataC*); 00086 void Finley_Assemble_setNormal(Finley_NodeFile*, Finley_ElementFile*, escriptDataC*); 00087 void Finley_Assemble_interpolate(Finley_NodeFile*,Finley_ElementFile*,escriptDataC*, escriptDataC*); 00088 void Finley_Assemble_gradient(Finley_NodeFile*, Finley_ElementFile*,escriptDataC*, escriptDataC*); 00089 void Finley_Assemble_integrate(Finley_NodeFile*,Finley_ElementFile*,escriptDataC*,double*) ; 00090 void Finley_Assemble_getSize(Finley_NodeFile*,Finley_ElementFile*, escriptDataC*); 00091 void Finley_Assemble_CopyNodalData(Finley_NodeFile* nodes,escriptDataC* out,escriptDataC* in); 00092 void Finley_Assemble_CopyElementData(Finley_ElementFile* elements,escriptDataC* out,escriptDataC* in); 00093 void Finley_Assemble_AverageElementData(Finley_ElementFile* elements,escriptDataC* out,escriptDataC* in); 00094 void Finley_Assemble_addToSystemMatrix(Paso_SystemMatrix* in, const dim_t NN_Equa, const index_t* Nodes_Equa, const dim_t num_Equa, 00095 const dim_t NN_Sol, const index_t* Nodes_Sol, const dim_t num_Sol, const double* array); 00096 00097 void Finley_Assemble_jacobeans_1D(double*, dim_t, double*, dim_t, dim_t, dim_t, index_t*, double*, dim_t, double*, double*, double*, index_t*); 00098 void Finley_Assemble_jacobeans_2D(double*, dim_t, double*, dim_t, dim_t, dim_t, index_t*, double*, dim_t, double*, double*, double*, index_t*); 00099 void Finley_Assemble_jacobeans_2D_M1D_E2D(double*, dim_t, double*, dim_t, dim_t, dim_t, index_t*, double*, dim_t, double*, double*, double*, index_t*); 00100 void Finley_Assemble_jacobeans_2D_M1D_E2D_C(double*, dim_t, double*, dim_t, dim_t, dim_t, index_t*, double*, dim_t, double*, double*, double*, index_t*); 00101 void Finley_Assemble_jacobeans_2D_M1D_E1D(double*, dim_t, double*, dim_t, dim_t, dim_t, index_t*, double*, dim_t, double*, double*, double*, index_t*); 00102 void Finley_Assemble_jacobeans_2D_M1D_E1D_C(double*, dim_t, double*, dim_t, dim_t, dim_t, index_t*, double*, dim_t, double*, double*, double*, index_t*); 00103 void Finley_Assemble_jacobeans_3D(double*, dim_t, double*, dim_t, dim_t, dim_t, index_t*, double*, dim_t, double*, double*, double*, index_t*); 00104 void Finley_Assemble_jacobeans_3D_M2D_E3D(double*, dim_t, double*, dim_t, dim_t, dim_t, index_t*, double*, dim_t, double*, double*, double*, index_t*); 00105 void Finley_Assemble_jacobeans_3D_M2D_E3D_C(double*, dim_t, double*, dim_t, dim_t, dim_t, index_t*, double*, dim_t, double*, double*, double*, index_t*); 00106 void Finley_Assemble_jacobeans_3D_M2D_E2D(double*, dim_t, double*, dim_t, dim_t, dim_t, index_t*, double*, dim_t, double*, double*, double*, index_t*); 00107 void Finley_Assemble_jacobeans_3D_M2D_E2D_C(double*, dim_t, double*, dim_t, dim_t, dim_t, index_t*, double*, dim_t, double*, double*, double*, index_t*); 00108 00109 00110 void Finley_Assemble_LumpedSystem(Finley_NodeFile* nodes,Finley_ElementFile* elements, escriptDataC* lumpedMat, escriptDataC* D, const bool_t useHRZ); 00111 #endif /* #ifndef INC_FINLEY_ASSEMBLE */ 00112