Commit 08a27a7c authored by Thomas Witkowski's avatar Thomas Witkowski
Browse files

Bugfix for integration of DOFVectors with dual mesh.

parent 2d8a6968
......@@ -758,6 +758,9 @@ namespace AMDiS {
TEST_EXIT_DBG(!quadFast || quadFast->getBasisFunctions() == feSpace->getBasisFcts())
("invalid basis functions");
if (smallElInfo->getMesh() == feSpace->getMesh())
return getVecAtQPs(smallElInfo, quad, quadFast, vecAtQPs);
const BasisFunction *basFcts = feSpace->getBasisFcts();
int nPoints =
quadFast ? quadFast->getQuadrature()->getNumPoints() : quad->getNumPoints();
......@@ -943,11 +946,11 @@ namespace AMDiS {
double det = dualElInfo.smallElInfo->getDet();
double tmp = 0.0;
vec1.getVecAtQPs(dualElInfo.smallElInfo, dualElInfo.largeElInfo,
quad, fastQuad, &(qp1[0]));
quad, NULL, &(qp1[0]));
vec2.getVecAtQPs(dualElInfo.smallElInfo, dualElInfo.largeElInfo,
quad, fastQuad, &(qp2[0]));
for (int iq = 0; iq < fastQuad->getNumPoints(); iq++)
tmp += fastQuad->getWeight(iq) * (*fct)(qp1[iq], qp2[iq]);
quad, NULL, &(qp2[0]));
for (int iq = 0; iq < quad->getNumPoints(); iq++)
tmp += quad->getWeight(iq) * (*fct)(qp1[iq], qp2[iq]);
value += tmp * det;
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment