gauss.h

The Spectral Toolkit is a C++ spectral transform library written by Rodney James and Chuck Panaccion

//
// spectral toolkit 
// copyright (c) 2005 university corporation for atmospheric research
// licensed under the gnu general public license
//

#ifndef __gauss__
#define __gauss__

#include "spectral.h"

namespace spectral
{
  /// Base class for Gaussian quadratures.  Includes common methods
  /// and members for the Legendre and Lobatto quadratures.

  class gauss
  {  
  public:
    gauss(int n);
    ~gauss();
    /// Gaussian point function.  
    /// \param i index from \f$ 0, \ldots, n-1 \f$
    /// \return value of point at i
    inline real point(int i) { return(points[i]); }
    /// Gaussian weight function.  
    /// \param i index from \f$ 0, \ldots, n-1 \f$
    /// \return value of weight at i
    inline real weight(int i) { return(weights[i]); }
  protected:
    /// Internal machine precision used for iteration.  Set in constructor.
    real eps;
    /// Internal array to store Gaussian weights.
    real *weights;
    /// Internal array to store Gaussian points.
    real *points;
    /// Internal length (total).
    int n;
  };
}

#endif

// Local Variables:
// mode:C++
// End: