moment_raw_from_upfactorial.h File Reference

Raw moments from rising-factorial moments. More...

#include <vector>
#include "line/num/number.h"
#include "line/util/error.h"
#include "line/util/matrix.h"
#include "line/util/population.h"
#include "line/api/moment/moment_apply.h"
#include "line/api/moment/moment_stirling2.h"

Include dependency graph for moment_raw_from_upfactorial.h:

Go to the source code of this file.

Namespaces
namespace	line
namespace	line::moment

Functions
template<class T>
std::vector< T >	line::moment::moment_raw_from_upfactorial (const std::vector< T > &fp)
	m_i = sum_j (-1)^(i-j) S(i,j) fp_j.

Detailed Description

Raw moments from rising-factorial moments.

Templated port of matlab/src/api/moment/moment_raw_from_upfactorial.m. Every operation is integer or rational, so the exact instantiation returns the transform with no rounding at all. That matters more here than almost anywhere else in the API: the alternating binomial sums of the moment conversions cancel catastrophically in double arithmetic once the order grows.

Definition in file moment_raw_from_upfactorial.h.