/[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 3196 by jfenwick, Wed Sep 22 01:18:52 2010 UTC revision 3197 by jfenwick, Wed Sep 22 05:52:31 2010 UTC
# Line 23  This file is not to be included in .h fi Line 23  This file is not to be included in .h fi
23    
24  #include "paso/Common.h"    // I just want the types not all the includes that get dragged in - fix that  #include "paso/Common.h"    // I just want the types not all the includes that get dragged in - fix that
25    
26    // These are constructed from dsdv in ShapeFunction.c in finley
27  // The first two are just there for functions that want a pointer  // The first two are just there for functions that want a pointer
28  static const double DTDV_0D[1][1]={{0}};  static const double DTDV_0D[1][1]={{0}};
29  static const double DTDV_1D[2][2]={{0,0},{0,0}};  static const double DTDV_1D[2][2]={{-1.,1},{-1.,1.}};
30  static const double DTDV_2D[3][2]={{-1,-1}, {1,0}, {0,1}};  static const double DTDV_2D[3][2]={{-1,-1}, {1,0}, {0,1}};
31  static const double DTDV_3D[4][3]={{-1, -1, -1}, {1, 0, 0}, {0, 1, 0}, {0, 0, 1}};  static const double DTDV_3D[4][3]={{-1, -1, -1}, {1, 0, 0}, {0, 1, 0}, {0, 0, 1}};
32    
33    
34    // If these appear to be in a different order to finley it is because finley uses macros to hide Fortran array ordering
35    static const double DTDV_2D_alt[3*3][2]={{-1,1}, {0,-1.}, {0,1},
36    {-1,1}, {0,-1.}, {0,1},
37    {-1,1}, {0,-1.}, {0,1}
38    };
39    // the repetition is a hack
40    // Why didn't I just reorder DTDV_2D?   Well some code apparently depends on the order as written.
41    // should probably fix that
42    
43    
44  // [0] is reduced quadrature, [1] is full quadrature  // [0] is reduced quadrature, [1] is full quadrature
45  // in order the positions are POINT, LINE, TRI, TET  // in order the positions are POINT, LINE, TRI, TET
46  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}};

Legend:
Removed from v.3196  
changed lines
  Added in v.3197

  ViewVC Help
Powered by ViewVC 1.1.26