gauss.cpp

来自「The Spectral Toolkit is a C++ spectral t」· C++ 代码 · 共 44 行

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

#include "gauss.h"

using namespace spectral;

/// Constructor for Gaussian quadrature.  Also sets machine epsilon.
/// \param M total length (not order) of stencil

gauss::gauss(int M)
{
  n=M;
  if(n<1)
    throw bad_parameter();
    
  points=new real[n];
  weights=new real[n];
    
  switch(sizeof(real))
    {
    case 4:
      eps=4.0e-07F;
      break;
    case 16:
      eps=2.0e-31L;
      break;
    default:
      eps=2.0e-15;
      break;
    }
}

/// Destructor for Gaussian quadrature.

gauss::~gauss()
{
  delete [] weights;
  delete [] points;
}