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

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

Parent Directory Parent Directory | Revision Log Revision Log


Revision 3196 - (hide annotations)
Wed Sep 22 01:18:52 2010 UTC (8 years, 6 months ago) by jfenwick
File MIME type: text/plain
File size: 1453 byte(s)
moving slowly
1 jfenwick 3184
2     /*******************************************************
3     *
4     * Copyright (c) 2003-2010 by University of Queensland
5     * Earth Systems Science Computational Center (ESSCC)
6     * http://www.uq.edu.au/esscc
7     *
8     * Primary Business: Queensland, Australia
9     * Licensed under the Open Software License version 3.0
10     * http://www.opensource.org/licenses/osl-3.0.php
11     *
12     *******************************************************/
13    
14    
15     /* Shape Function info
16     These tables are a much simplified version of content from finley's ShapeFunctions files
17    
18     This file is not to be included in .h files - only .c files should have any use for it
19     */
20    
21     #ifndef SHAPETABLE_DUDLEY
22     #define SHAPETABLE_DUDLEY
23    
24 jfenwick 3196 #include "paso/Common.h" // I just want the types not all the includes that get dragged in - fix that
25    
26    
27     // The first two are just there for functions that want a pointer
28     static const double DTDV_0D[1][1]={{0}};
29     static const double DTDV_1D[2][2]={{0,0},{0,0}};
30 jfenwick 3184 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}};
32    
33    
34 jfenwick 3196 // [0] is reduced quadrature, [1] is full quadrature
35     // in order the positions are POINT, LINE, TRI, TET
36     static const double QuadWeight[4][2]={{0, 0}, {1., 0.5}, {0.5, 1./6}, {1./6, 1./24}};
37 jfenwick 3184
38 jfenwick 3196 static const dim_t QuadNums[4][2] ={{0,0}, {1,2}, {1,3}, {1,4}};
39    
40     //shape functions at quadrature nodes
41     bool_t getQuadShape(dim_t sim, bool_t reduced, const double** shapearr);
42    
43    
44    
45    
46    
47    
48    
49 jfenwick 3184 #endif
50    
51    
52    

  ViewVC Help
Powered by ViewVC 1.1.26