diff --git a/dune/mmdg/mmdg.hh b/dune/mmdg/mmdg.hh index cb595057f8d1daf35a474b4f973dfacd4660359b..f58bbf0a2a4c3ef759dc29bcd0b4ba034afe3737 100644 --- a/dune/mmdg/mmdg.hh +++ b/dune/mmdg/mmdg.hh @@ -589,7 +589,7 @@ private: const Scalar interfaceUpdate5 = weight * dEvaluation * ( (qpGlobal * iFrame[j]) * ( interfaceUpdate2 - 0.5 * interfaceUpdate3 - qpI * interfaceUpdate1 ) - - (KparDotTau_j * normal) * qpI * dEvaluation ); + - 0.5 * (KparDotTau_j * normal) * qpI * dEvaluation ); //quadrature for // int_intersct mu * jump(phi_iElem,i) * jump(phi_iElem,j) dr @@ -754,7 +754,7 @@ private: (*Base::A)[iElemIdxSLE + i + 1][iElemIdxSLE + i + 1] += weight * qpI * dEvaluation * (interfaceUpdate2 - 2.0 * interfaceUpdate3 - - qpI * qpI * interfaceUpdate1); + - qpI * interfaceUpdate1); for (int j = 0; j < i; j++) {