map_fit_detail.h File Reference

Scalar helpers shared by the MAP/PH moment-matching headers. More...

#include <cmath>
#include "line/num/number.h"
#include "line/util/error.h"

Include dependency graph for map_fit_detail.h:

Go to the source code of this file.

Namespaces
namespace	line
namespace	line::mam

Detailed Description

Scalar helpers shared by the MAP/PH moment-matching headers.

These are not ports of any MATLAB function: they are the arithmetic primitives the closed forms need (square root, exponential, logarithm, integer power) resolved through argument-dependent lookup so that the same expression compiles for double and for Boost.Multiprecision types, plus a minimal complex arithmetic type.

The complex type exists because aph_fit's second fitting case (the Bobbio-Horvath-Telek chain K9..K22) takes square and cube roots of quantities that are negative for many feasible moment sets; the imaginary parts cancel in the final f. MATLAB evaluates that chain in complex arithmetic implicitly and the JAR does it explicitly with org.apache.commons.math3.complex.Complex, so the port must too. Roots use the principal branch, matching both references.

Definition in file map_fit_detail.h.