|
00001 // Copyright (C) 2003 Johan Hoffman and Anders Logg. 00002 // Licensed under the GNU GPL Version 2. 00003 00004 #ifndef __TETRAHEDRON_MIDPOINT_QUADRATURE_H 00005 #define __TETRAHEDRON_MIDPOINT_QUADRATURE_H 00006 00007 #include <dolfin/Quadrature.h> 00008 00009 namespace dolfin { 00010 00013 class TetrahedronMidpointQuadrature : public Quadrature { 00014 public: 00015 00016 TetrahedronMidpointQuadrature() : Quadrature(6) { 00017 00018 // Volume of tetrahedron 00019 m = 1 / 6.0; 00020 00021 // Quadrature points 00022 points[0] = Point(0.5, 0.0, 0.0); 00023 points[1] = Point(0.5, 0.5, 0.0); 00024 points[2] = Point(0.0, 0.5, 0.0); 00025 points[3] = Point(0.0, 0.0, 0.5); 00026 points[4] = Point(0.5, 0.0, 0.5); 00027 points[5] = Point(0.0, 0.5, 0.5); 00028 00029 // Quadrature weights 00030 weights[0] = m / 6.0; 00031 weights[1] = m / 6.0; 00032 weights[2] = m / 6.0; 00033 weights[3] = m / 6.0; 00034 weights[4] = m / 6.0; 00035 weights[5] = m / 6.0; 00036 00037 } 00038 00039 void show() const {}; 00040 00041 }; 00042 00043 } 00044 00045 #endif
Documentation automatically generated with Doxygen on 10 Sep 2004.