3
#include <libciomr/libciomr.h>
9
StatCombData::StatCombData(int imax) :
13
throw std::runtime_error("ERROR: StatCombData::StatCombData -- argument out of range");
15
Fact_ = new PSI_FLOAT[imax+1];
17
for(int i=1; i<=imax; i++)
18
Fact_[i] = i*Fact_[i-1];
20
Fact2_ = new PSI_FLOAT[imax+1];
23
for(int i=2; i<=imax; i++)
24
Fact2_[i] = i*Fact_[i-2];
26
BinomC_ = block_matrix(imax+1,imax+1);
27
for(int n=0; n<=imax; n++)
28
for(int m=0; m<=n; m++)
29
BinomC_[n][m] = Fact_[n]/(Fact_[m]*Fact_[n-m]);
32
StatCombData::~StatCombData()