erf.h

Go to the documentation of this file.

//
// some arbitrary maths bits
// stu meikle Wed Mar  8 10:02:15 2000
//

#ifndef ERF_H
#define ERF_H

#include <cmath>

namespace mimas {
  
  class mm_erf
    {
    private:
      double  *vals;
    public:
      mm_erf(void)
  {
    vals = NULL;
  }
  
      ~mm_erf(void)
  {
    if (vals)
      delete[] vals;
  }
  
      void  init(int l, double range, double acc)
  {
    vals = new double[l];
      
    double y;

    for (y=-range+acc/2.0;y<=range+acc;y+=acc)
      {
        vals[ (int)floor( y/acc ) + ( l / 2 ) ] = erf( y );
      }

  }
    
    };

};

#endif

---