escript
Revision_
|
00001 00002 /***************************************************************************** 00003 * 00004 * Copyright (c) 2003-2014 by University of Queensland 00005 * http://www.uq.edu.au 00006 * 00007 * Primary Business: Queensland, Australia 00008 * Licensed under the Open Software License version 3.0 00009 * http://www.opensource.org/licenses/osl-3.0.php 00010 * 00011 * Development until 2012 by Earth Systems Science Computational Center (ESSCC) 00012 * Development 2012-2013 by School of Earth Sciences 00013 * Development from 2014 by Centre for Geoscience Computing (GeoComp) 00014 * 00015 *****************************************************************************/ 00016 00017 /* Shape Function info 00018 These tables are a much simplified version of content from finley's ShapeFunctions files 00019 00020 This file is not to be included in .h files - only .c files should have any use for it 00021 */ 00022 00023 #ifndef SHAPETABLE_DUDLEY 00024 #define SHAPETABLE_DUDLEY 00025 00026 #include "esysUtils/types.h" 00027 00028 #include "ElementType.h" 00029 00030 /* These are constructed from dsdv in ShapeFunction.c in finley 00031 The first two are just there for functions that want a pointer 00032 */ 00033 static const double DTDV_0D[1][1] = { {0} }; 00034 static const double DTDV_1D[2][2] = { {-1., 1}, {-1., 1.} }; 00035 00036 /* The repetition here is a hack to prevent out of bounds access */ 00037 static const double DTDV_2D[3 * 3][2] = { {-1, 1}, {0, -1.}, {0, 1}, 00038 {-1, 1}, {0, -1.}, {0, 1}, 00039 {-1, 1}, {0, -1.}, {0, 1} 00040 }; 00041 static const double DTDV_3D[4][3] = { {-1, -1, -1}, {1, 0, 0}, {0, 1, 0}, {0, 0, 1} }; 00042 00043 /* Index by the following by Dudley_ElementTypeID 00044 * The number of local dimensions (as opposed to dimension of the embedding space) */ 00045 static const dim_t localDims[8] = { 0, 1, 2, 3, 0, 1, 2, 0 }; 00046 static const dim_t Dims[8] = { 0, 1, 2, 3, 1, 2, 3, 0 }; 00047 00048 /* the following lists are only used for face elements defined by numNodesOnFace>0 */ 00049 static const dim_t numNodesOnFaceMap[8] = { 1, 2, 3, 4, 1, 2, 4, -1 }; /* if the element is allowed as a face element, numNodesOnFace defines the number of nodes defining the face */ 00050 static const dim_t shiftNodesMap[8][4] = { {0}, {1, 0}, {1, 2, 0}, {-1}, {0, 1, 2}, {1, 0, 2}, {1, 2, 0, 3}, {0} }; /* defines a permutation of the nodes which rotates the nodes on the face */ 00051 static const dim_t reverseNodesMap[8][4] = { {-1}, {-1}, {0, 2, 1}, {-1}, {-1}, {-1}, {0, 2, 1, 3}, {0} }; /* reverses the order of the nodes on a face. the permutation has keep 0 fixed. */ 00052 /* shiftNodes={-1} or reverseNodes={-1} are ignored. */ 00053 00054 /* [0] is reduced quadrature, [1] is full quadrature */ 00055 /* in order the positions are POINT, LINE, TRI, TET */ 00056 static const double QuadWeight[4][2] = { {0, 0}, {1., 0.5}, {0.5, 1. / 6}, {1. / 6, 1. / 24} }; 00057 00058 /* number of quadrature points per element */ 00059 static const dim_t QuadNums[4][2] = { {0, 0}, {1, 2}, {1, 3}, {1, 4} }; 00060 00061 /*shape functions at quadrature nodes */ 00062 bool getQuadShape(dim_t sim, bool reduced, const double **shapearr); 00063 00064 const char *getElementName(Dudley_ElementTypeId id); 00065 00066 #endif