//// 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 stencilgauss::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;}