diff --git a/src/libcadet/model/binding/GeneralizedIonExchangeBinding.cpp b/src/libcadet/model/binding/GeneralizedIonExchangeBinding.cpp index 34ccd7995..3daec43da 100644 --- a/src/libcadet/model/binding/GeneralizedIonExchangeBinding.cpp +++ b/src/libcadet/model/binding/GeneralizedIonExchangeBinding.cpp @@ -110,7 +110,6 @@ namespace { const int offset = comp * nPieces; const auto [nuConst, nuVar] = cadet::evaluateCubicPiecewisePolynomialSplit::type>(pH, p.nuBreaks.data() + comp * (nPieces + 1), p.nu.data() + offset, p.nuLin.data() + offset, p.nuQuad.data() + offset, p.nuCube.data() + offset, nPieces); - LOG(Debug) << "cpQNuPowers " << comp << ": at " << static_cast(pH) << ":" << static_cast(nuConst) << " + " << static_cast(nuVar) << " = " << static_cast(nuConst) + static_cast(nuVar); const cp_state_t nu_i_0_over_nu0 = static_cast(nuConst) / static_cast(p.nu[0]); const cp_state_t nu_i_pH_over_nu0 = nuVar / static_cast(p.nu[0]); return {pow(cpBase, nu_i_0_over_nu0) * pow(cpVar, nu_i_pH_over_nu0), pow(qBase, nu_i_0_over_nu0) * pow(qVar, nu_i_pH_over_nu0)}; @@ -127,10 +126,7 @@ namespace else { const int offset = comp * nPieces; -// return cadet::evaluateCubicPiecewisePolynomialSplit(pH, p.nuBreaks.data() + comp * (nPieces + 1), p.nu.data() + offset, p.nuLin.data() + offset, p.nuQuad.data() + offset, p.nuCube.data() + offset, nPieces); - const auto [nuConst, nuVar] = cadet::evaluateCubicPiecewisePolynomialSplit(pH, p.nuBreaks.data() + comp * (nPieces + 1), p.nu.data() + offset, p.nuLin.data() + offset, p.nuQuad.data() + offset, p.nuCube.data() + offset, nPieces); - LOG(Debug) << "evalNu " << comp << ": at " << static_cast(pH) << ": " << static_cast(nuConst) << " + " << static_cast(nuVar) << " = " << static_cast(nuConst) + static_cast(nuVar); - return std::make_tuple(nuConst, nuVar); + return cadet::evaluateCubicPiecewisePolynomialSplit(pH, p.nuBreaks.data() + comp * (nPieces + 1), p.nu.data() + offset, p.nuLin.data() + offset, p.nuQuad.data() + offset, p.nuCube.data() + offset, nPieces); } }