gauss_legendre.h

Go to the documentation of this file.

// SPDX-License-Identifier: BSD-3-Clause
// Copyright (c) 2026, Liam Pohlmann
 
#ifndef HUMMINGBIRD_QUADRATURE_GAUSS_LEGENDRE_H_
#define HUMMINGBIRD_QUADRATURE_GAUSS_LEGENDRE_H_
 
#include "quadrature/quadrature_base.h"
namespace hummingbird {
 
/**
 * @brief Defines a Gauss-Legendre quadrature set on the interval [-1,1]
 *
 */
class GaussLegendre : public QuadratureBase<double> {
 public:
  /**
   * @brief Construct a new GaussLegendre object
   *
   * @param n_points Number of quadrature points to be created
   */
  GaussLegendre(const unsigned int n_points);
 
 private:
  /**
   * @brief Compute the weight for abscissa value k for a n-th order
   * Gauss-Legendre polynomial
   *
   * @param k Abscissa number
   * @param n Order of Gauss-Legendre polynomial
   * @return Weight corresponding to the k-th abscissa values
   */
  double ComputeWeight(const size_t k, const size_t n) override;
};
}  // namespace hummingbird
 
#endif  // HUMMINGBIRD_QUADRATURE_GAUSS_LEGENDRE_H_