All AMAP(2) representations matching three moments and the autocorrelation decay rate (matlab/lib/m3a/m3a/amap2/amap2_fitall_gamma.m). More...

#include <vector>
#include "line/api/mam/amap2_assemble.h"
#include "line/api/mam/map_fit_detail.h"
#include "line/api/mam/map_moment.h"
#include "line/num/number.h"
#include "line/util/error.h"

Include dependency graph for amap2_fitall_gamma.h:

Namespaces
namespace	line
namespace	line::mam

Functions
template<class T>
std::vector< Map< T > >	line::mam::amap2_fitall_gamma (const T &M1, const T &M2, const T &M3, const T &GAMMA, const T &degentol, const T &r12tol)
	Every AMAP(2) matching (M1, M2, M3, GAMMA).
template<class T>
std::vector< Map< T > >	line::mam::amap2_fitall_gamma (const T &M1, const T &M2, const T &M3, const T &GAMMA)
	amap2_fitall_gamma with the MATLAB defaults degentol = 1e-8, r12tol = 1e-6.

Detailed Description

All AMAP(2) representations matching three moments and the autocorrelation decay rate (matlab/lib/m3a/m3a/amap2/amap2_fitall_gamma.m).

The two phase means h1, h2 come from the same quadratic as aph2_fitall. For GAMMA >= 0 the first canonical form leaves a further quadratic in the branching probability r2, with discriminant z = M1^2 G^2 + (2 M1 h1 + 2 M1 h2 - 4 h1 h2 - 2 M1^2) G

M1^2 - 2 M1 h1 - 2 M1 h2 + h1^2 + 2 h1 h2 + h2^2, so up to four AMAP(2)s can match the same characteristics; for GAMMA < 0 the second canonical form determines r1 and r2 uniquely per (h1, h2). Solutions whose probabilities fall outside [0, 1] beyond r12tol are discarded, the others are clamped into [0, 1].

Unlike aph2_fitall this returns an empty vector when the characteristics are infeasible: no approximate fitting is performed, exactly as in the reference.

Gated on transcendental arithmetic: the two moment discriminants and the SCV <= 1 lower bound of the third moment are square roots.

Definition in file amap2_fitall_gamma.h.

Namespaces

Functions

Detailed Description