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

add missing minus sign in inhomogeneousbulkproblem

parent 523c9fb7
No related branches found
No related tags found
No related merge requests found
......@@ -12,9 +12,9 @@ class InhomogeneousBulkProblem : public DGProblem<Coordinate, Scalar>
using Matrix = typename Base::Matrix;
//the exact solution at position pos
//1d: f(x) = x,
//2d: f(x,y) = x*y,
//3d: f(x,y,z) = x*y*z
//1d: p(x) = x,
//2d: p(x,y) = x*y,
//3d: p(x,y,z) = x*y*z
virtual Scalar exactSolution (const Coordinate& pos) const
{
return pos * Coordinate(1.0);
......@@ -29,7 +29,7 @@ class InhomogeneousBulkProblem : public DGProblem<Coordinate, Scalar>
//source term at position pos
virtual Scalar q (const Coordinate& pos) const
{
return pos * Coordinate(2.0);
return -pos * Coordinate(2.0);
}
//permeability tensor at position pos
......@@ -47,23 +47,23 @@ class InhomogeneousBulkProblem : public DGProblem<Coordinate, Scalar>
//returns the recommended quadrature order to compute an integral
//over x * q(x)
virtual int quadratureOrder_q () const
int quadratureOrder_q () const
{
return 10;
return 2;
}
//returns the recommended quadrature order to compute an integral
//over x * boundary(x)
virtual int quadratureOrderBoundary () const
//over x * boundary(x) and K(x) * boundary(x)
int quadratureOrderBoundary () const
{
return 10;
return 1 + dim;
}
//returns the recommended quadrature order to compute an integral
//over an entry K[i][j] of the permeability tensor
virtual int quadratureOrder_K () const
int quadratureOrder_K () const
{
return 10;
return 2;
}
};
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment