Commit 8518668b authored by Praetorius, Simon's avatar Praetorius, Simon
Browse files

small error in NavierStokes_TH_MultiPhase

parent 9bbfbc0b
......@@ -80,7 +80,7 @@ void NavierStokes_TH_MultiPhase::fillOperators()
/// < (1/tau)*rho*u_i^old , psi >
Operator *opTimeOld = new Operator(getFeSpace(i), getFeSpace(i));
if (density1==density2)
opTimeOld->addTerm(new VecAtQP_ZOT(prob->getSolution()->getDOFVector(i), new AMDiS::Factor(density1)));
opTimeOld->addTerm(new VecAtQP_ZOT(prob->getSolution()->getDOFVector(i), new AMDiS::Factor<double>(density1)));
else
opTimeOld->addTerm(new Vec2AtQP_ZOT(densityPhase, prob->getSolution()->getDOFVector(i), NULL));
prob->addVectorOperator(*opTimeOld, i, getInvTau(), getInvTau());
......
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