Skip to content
Snippets Groups Projects
Commit 376ca5e1 authored by Hörl, Maximilian's avatar Hörl, Maximilian
Browse files

fix wrong quadrature order in Poisson problem

parent 1269124d
No related branches found
No related tags found
No related merge requests found
...@@ -56,7 +56,7 @@ class DGProblem ...@@ -56,7 +56,7 @@ class DGProblem
} }
//returns the recommended quadrature order to compute an integral //returns the recommended quadrature order to compute an integral
//over x * boundary(x) //over x * boundary(x) and K(x) * boundary(x)
virtual int quadratureOrderBoundary () const virtual int quadratureOrderBoundary () const
{ {
return 1; return 1;
......
...@@ -8,6 +8,8 @@ template<class Coordinate, class Scalar = double> ...@@ -8,6 +8,8 @@ template<class Coordinate, class Scalar = double>
class Poisson : public DGProblem<Coordinate, Scalar> class Poisson : public DGProblem<Coordinate, Scalar>
{ {
public: public:
static constexpr int dim = Coordinate::dimension;
//the exact solution at position pos //the exact solution at position pos
//1d: f(x) = x, //1d: f(x) = x,
//2d: f(x,y) = x*y, //2d: f(x,y) = x*y,
...@@ -29,6 +31,13 @@ class Poisson : public DGProblem<Coordinate, Scalar> ...@@ -29,6 +31,13 @@ class Poisson : public DGProblem<Coordinate, Scalar>
{ {
return true; return true;
}; };
//returns the recommended quadrature order to compute an integral
//over x * boundary(x) and K(x) * boundary(x)
int quadratureOrderBoundary () const
{
return dim;
}
}; };
#endif #endif
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment