diff --git a/dumux/decoupled/2p2c/fvpressure2p2c.hh b/dumux/decoupled/2p2c/fvpressure2p2c.hh index d5d05a72d12e0ab7eb28cebdc28d163d24d59909..06adf310c2faed323e39dd928d2e0a8c64f8e930 100644 --- a/dumux/decoupled/2p2c/fvpressure2p2c.hh +++ b/dumux/decoupled/2p2c/fvpressure2p2c.hh @@ -953,12 +953,12 @@ void FVPressure2P2C<TypeTag>::updateMaterialLawsInElement(const Element& element fluidState.saturation(wPhaseIdx)); // TODO: get right criterion, do output for evaluation //converge criterion - if (abs(oldPc-pc)<10) + if (fabs(oldPc-pc)<10) maxiter = 1; iterout = iter; } if(iterout !=0) - Dune::dinfo << iterout << "times iteration of pc was applied at Idx " << globalIdx << ", pc delta still " << abs(oldPc-pc) << std::endl; + Dune::dinfo << iterout << "times iteration of pc was applied at Idx " << globalIdx << ", pc delta still " << fabs(oldPc-pc) << std::endl; } // initialize phase properties not stored in fluidstate cellData.setViscosity(wPhaseIdx, FluidSystem::viscosity(fluidState, wPhaseIdx)); diff --git a/dumux/material/components/iapws/common.hh b/dumux/material/components/iapws/common.hh index 200dbc8850a8804d803486b8f74d0703f1a10798..db52883e08c0f10a7a3c9c58516fe21778bd8481 100644 --- a/dumux/material/components/iapws/common.hh +++ b/dumux/material/components/iapws/common.hh @@ -213,7 +213,7 @@ public: * rhobar + thcond_b2 * exp(thcond_B1 * ((rhobar + thcond_B2)*(rhobar + thcond_B2))); - Scalar DTbar = abs(Tbar - 1) + thcond_c4; + Scalar DTbar = fabs(Tbar - 1) + thcond_c4; Scalar DTbarpow = pow(DTbar, 3./5); Scalar Q = 2. + thcond_c5 / DTbarpow;