00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026 #ifndef U2D_GAUSSPOINTS_H
00027 #define U2D_GAUSSPOINTS_H
00028
00029 #include "u2d_mesh.h"
00030 #include "boost/numeric/mtl/mtl.hpp"
00031
00032 using namespace mtl;
00033
00039 class u2d_DomainGaussPoints {
00040 public:
00041 u2d_DomainGaussPoints(int _domain_size);
00042 u2d_DomainGaussPoints(const u2d_DomainGaussPoints& orig);
00043 virtual ~u2d_DomainGaussPoints();
00044
00045
00046 dense2D<double> dw_x;
00047 dense2D<double> dw_y;
00048 dense_vector<double> jacobian_weight;
00049 private:
00050 int domain_size;
00051
00052 };
00053
00059 class u2d_BoundaryGaussPoints {
00060 public:
00061 u2d_BoundaryGaussPoints(int _boundary_size);
00062 u2d_BoundaryGaussPoints(const u2d_BoundaryGaussPoints& orig);
00063 virtual ~u2d_BoundaryGaussPoints();
00064
00065
00066
00067 dense_vector<double> norm_x;
00068 dense_vector<double> norm_y;
00069 dense2D<double> dw_x;
00070 dense2D<double> dwp_x;
00071 dense_vector<double> jacobian_weight;
00072
00073 private:
00074 int boundary_size;
00075
00076 };
00077
00078
00079 #endif
00080
00081