ESYS13  Revision_
Assemble.h
Go to the documentation of this file.
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