LSST Applications g0265f82a02+0e5473021a,g02d81e74bb+0dd8ce4237,g1470d8bcf6+3ea6592b6f,g2079a07aa2+86d27d4dc4,g2305ad1205+5ca4c0b359,g295015adf3+d10818ec9d,g2a9a014e59+6f9be1b9cd,g2bbee38e9b+0e5473021a,g337abbeb29+0e5473021a,g3ddfee87b4+703ba97ebf,g487adcacf7+4fa16da234,g50ff169b8f+96c6868917,g52b1c1532d+585e252eca,g591dd9f2cf+ffa42b374e,g5a732f18d5+53520f316c,g64a986408d+0dd8ce4237,g858d7b2824+0dd8ce4237,g8a8a8dda67+585e252eca,g99cad8db69+d39438377f,g9ddcbc5298+9a081db1e4,ga1e77700b3+15fc3df1f7,ga8c6da7877+f1d96605c8,gb0e22166c9+60f28cb32d,gb6a65358fc+0e5473021a,gba4ed39666+c2a2e4ac27,gbb8dafda3b+e5339d463f,gc120e1dc64+da31e9920e,gc28159a63d+0e5473021a,gcf0d15dbbd+703ba97ebf,gdaeeff99f8+f9a426f77a,ge6526c86ff+889fc9d533,ge79ae78c31+0e5473021a,gee10cc3b42+585e252eca,gf18bd8381d+7268b93478,gff1a9f87cc+0dd8ce4237,w.2024.16
LSST Data Management Base Package
Loading...
Searching...
No Matches
Public Member Functions | List of all members
lsst::geom::polynomials::BinomialMatrix Class Reference

A class that computes binomial coefficients up to a certain power. More...

#include <BinomialMatrix.h>

Public Member Functions

 BinomialMatrix (int nMax)
 Construct an object that can compute binomial coefficients with \(n\) up to and including the given value.
 
double operator() (int n, int k) const
 Return the binomial coefficient.
 

Detailed Description

A class that computes binomial coefficients up to a certain power.

The binomial coefficient is defined as:

\[ \left(\begin{array}{ c } n \\ k \end{array}\right) = \frac{n!}{k!(n-k)!} \]

with both \(n\) and \(k\) nonnegative integers and \(k \le n\)

This class uses recurrence relations to avoid computing factorials directly, making it both more efficient and numerically stable.

Definition at line 45 of file BinomialMatrix.h.

Constructor & Destructor Documentation

◆ BinomialMatrix()

lsst::geom::polynomials::BinomialMatrix::BinomialMatrix ( int nMax)
explicit

Construct an object that can compute binomial coefficients with \(n\) up to and including the given value.

Definition at line 27 of file BinomialMatrix.cc.

27 : _matrix(Eigen::MatrixXd::Zero(nMax + 1, nMax + 1)) {
28 _matrix.col(0).setConstant(1);
29 _matrix.diagonal().setConstant(1);
30 for (int i = 2; i <= nMax; ++i) {
31 _matrix(i, 0) = 1.0;
32 _matrix(i, i) = 1.0;
33 for (int j = 1; j < i; ++j) {
34 _matrix(i, j) = _matrix(i - 1, j - 1) *
35 (static_cast<double>(i) / static_cast<double>(j));
36 }
37 }
38}

Member Function Documentation

◆ operator()()

double lsst::geom::polynomials::BinomialMatrix::operator() ( int n,
int k ) const
inline

Return the binomial coefficient.

No error checking is performed; the behavior of this method is undefined if the given values do not satisfy

n <= nMax && k <= n && n >=0 && k >= 0

Definition at line 63 of file BinomialMatrix.h.

63 {
64 return _matrix(n, k);
65 }

The documentation for this class was generated from the following files: