12345678910111213141516171819202122232425262728 |
- #ifndef COEFFICIENTS_HPP
- #define COEFFICIENTS_HPP
- #include <mpfr.h>
- class Coefficients{
- private:
- static constexpr size_t prec=4096;
- static constexpr size_t imax=1000;
- static constexpr size_t ncoeffs=((imax+2)*(imax+1))/2;
- double coeffs[ncoeffs];
- size_t pos(size_t i,size_t j);
- public:
- Coefficients();
- double get(size_t i,size_t j);
- };
- inline size_t
- Coefficients::pos(size_t i,size_t j){
- return i+j*(j+1)/2;
- }
- inline double
- Coefficients::get(size_t i,size_t j){
- return (i<j)?coeffs[pos(i,j)]:coeffs[pos(j,i)];
- }
- #endif
|