@@ -110,7 +110,6 @@ namespace
110
110
{
111
111
const int offset = comp * nPieces;
112
112
const auto [nuConst, nuVar] = cadet::evaluateCubicPiecewisePolynomialSplit<ph_t , typename cadet::DoubleActivePromoter<param_t , cp_state_t >::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);
113
- LOG (Debug) << " cpQNuPowers " << comp << " : at " << static_cast <double >(pH) << " :" << static_cast <double >(nuConst) << " + " << static_cast <double >(nuVar) << " = " << static_cast <double >(nuConst) + static_cast <double >(nuVar);
114
113
const cp_state_t nu_i_0_over_nu0 = static_cast <param_t >(nuConst) / static_cast <param_t >(p.nu [0 ]);
115
114
const cp_state_t nu_i_pH_over_nu0 = nuVar / static_cast <param_t >(p.nu [0 ]);
116
115
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
127
126
else
128
127
{
129
128
const int offset = comp * nPieces;
130
- // return cadet::evaluateCubicPiecewisePolynomialSplit<ph_t, result_t>(pH, p.nuBreaks.data() + comp * (nPieces + 1), p.nu.data() + offset, p.nuLin.data() + offset, p.nuQuad.data() + offset, p.nuCube.data() + offset, nPieces);
131
- const auto [nuConst, nuVar] = cadet::evaluateCubicPiecewisePolynomialSplit<ph_t , result_t >(pH, p.nuBreaks .data () + comp * (nPieces + 1 ), p.nu .data () + offset, p.nuLin .data () + offset, p.nuQuad .data () + offset, p.nuCube .data () + offset, nPieces);
132
- LOG (Debug) << " evalNu " << comp << " : at " << static_cast <double >(pH) << " : " << static_cast <double >(nuConst) << " + " << static_cast <double >(nuVar) << " = " << static_cast <double >(nuConst) + static_cast <double >(nuVar);
133
- return std::make_tuple (nuConst, nuVar);
129
+ return cadet::evaluateCubicPiecewisePolynomialSplit<ph_t , result_t >(pH, p.nuBreaks .data () + comp * (nPieces + 1 ), p.nu .data () + offset, p.nuLin .data () + offset, p.nuQuad .data () + offset, p.nuCube .data () + offset, nPieces);
134
130
}
135
131
}
136
132
0 commit comments