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

Diff of /branches/domexper/dudley/src/ShapeTable.h

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 3250 by caltinay, Wed Oct 6 05:53:06 2010 UTC revision 3251 by jfenwick, Thu Oct 7 04:02:30 2010 UTC
# Line 20  This file is not to be included in .h fi Line 20  This file is not to be included in .h fi
20  #ifndef SHAPETABLE_DUDLEY  #ifndef SHAPETABLE_DUDLEY
21  #define SHAPETABLE_DUDLEY  #define SHAPETABLE_DUDLEY
22    
23  #include "esysUtils/types.h"    // I just want the types not all the includes that get dragged in - fix that  #include "esysUtils/types.h"    
24    
25  #include "ElementType.h"  #include "ElementType.h"
26    
27  // These are constructed from dsdv in ShapeFunction.c in finley  /* These are constructed from dsdv in ShapeFunction.c in finley
28  // The first two are just there for functions that want a pointer     The first two are just there for functions that want a pointer
29    */
30  static const double DTDV_0D[1][1] = { {0} };  static const double DTDV_0D[1][1] = { {0} };
31  static const double DTDV_1D[2][2] = { {-1., 1}, {-1., 1.} };  static const double DTDV_1D[2][2] = { {-1., 1}, {-1., 1.} };
32    
33  // The repetition here is a hack to prevent out of bounds access  /* The repetition here is a hack to prevent out of bounds access */
34  static const double DTDV_2D[3 * 3][2] = { {-1, 1}, {0, -1.}, {0, 1},  static const double DTDV_2D[3 * 3][2] = { {-1, 1}, {0, -1.}, {0, 1},
35  {-1, 1}, {0, -1.}, {0, 1},  {-1, 1}, {0, -1.}, {0, 1},
36  {-1, 1}, {0, -1.}, {0, 1}  {-1, 1}, {0, -1.}, {0, 1}
# Line 45  static const dim_t Dims[8] = { 0, 1, 2, Line 46  static const dim_t Dims[8] = { 0, 1, 2,
46  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 */  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 */
47  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 */  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 */
48  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. */  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. */
   
49                            /* shiftNodes={-1} or reverseNodes={-1} are ignored. */                            /* shiftNodes={-1} or reverseNodes={-1} are ignored. */
50    
51  // [0] is reduced quadrature, [1] is full quadrature  /* [0] is reduced quadrature, [1] is full quadrature */
52  // in order the positions are POINT, LINE, TRI, TET  /* in order the positions are POINT, LINE, TRI, TET */
53  static const double QuadWeight[4][2] = { {0, 0}, {1., 0.5}, {0.5, 1. / 6}, {1. / 6, 1. / 24} };  static const double QuadWeight[4][2] = { {0, 0}, {1., 0.5}, {0.5, 1. / 6}, {1. / 6, 1. / 24} };
54    
55  // number of quadrature points per element  /* number of quadrature points per element */
56  static const dim_t QuadNums[4][2] = { {0, 0}, {1, 2}, {1, 3}, {1, 4} };  static const dim_t QuadNums[4][2] = { {0, 0}, {1, 2}, {1, 3}, {1, 4} };
57    
58  //shape functions at quadrature nodes  /*shape functions at quadrature nodes */
59  bool_t getQuadShape(dim_t sim, bool_t reduced, const double **shapearr);  bool_t getQuadShape(dim_t sim, bool_t reduced, const double **shapearr);
60    
61  const char *getElementName(Dudley_ElementTypeId id);  const char *getElementName(Dudley_ElementTypeId id);

Legend:
Removed from v.3250  
changed lines
  Added in v.3251

  ViewVC Help
Powered by ViewVC 1.1.26