math_Sign_Zolotarev.h

Go to the documentation of this file.

 
#ifndef MATH_SIGN_ZOLOTAREV_INCLUDED
#define MATH_SIGN_ZOLOTAREV_INCLUDED
 
#include <cassert>
 
#include "bridge_defs.h"
#include "Parameters/commonParameters.h"
 
#include "IO/bridgeIO.h"
using Bridge::vout;
 
 
class Math_Sign_Zolotarev
{
 public:
  static const std::string class_name;
 
 protected:
  Bridge::VerboseLevel m_vl;
 
 private:
  int m_Np;
  double m_bmax;
  std::vector<double> m_cl;
  std::vector<double> m_bl;
 
 public:
  Math_Sign_Zolotarev(const int Np, const double bmax)
    : m_vl(CommonParameters::Vlevel()),
    m_Np(Np), m_bmax(bmax)
  {
    set_sign_parameters();
  }
 
  void get_sign_parameters(std::vector<double>& cl, std::vector<double>& bl);
 
  double sign(const double x);
 
 private:
  void set_sign_parameters();
  void poly_Zolotarev(const double bmax, double& UK);
  void Jacobi_elliptic(const double uu, const double emmc,
                       double& sn, double& cn, double& dn);
};
#endif