23 const string str_vlevel = params.
get_string(
"verbose_level");
32 err += params.
fetch_int(
"number_of_poles", Np);
33 err += params.
fetch_int(
"exponent_numerator", n_exp);
34 err += params.
fetch_int(
"exponent_denominator", d_exp);
49 const double x_min,
const double x_max)
88 std::vector<double>& pole)
90 if (res.size() !=
m_Np) {
95 if (pole.size() !=
m_Np) {
101 for (
int i = 0; i <
m_Np; ++i) {
119 char filename[FILENAME_MAX];
121 snprintf(filename, FILENAME_MAX,
122 "parameter_rational.%1d_%1d%s_%02d_%010.8f_%07.3f",
130 int Np, n_exp, d_exp;
134 std::fstream parameterfile;
135 parameterfile.open(filename, std::ios::in);
136 if (!parameterfile.is_open()) {
137 vout.
crucial(
m_vl,
"Error at %s: failed to open parameter file. %s(%d)\n",
142 parameterfile >> Np >> n_exp >> d_exp;
143 parameterfile >> x_min >> x_max;
146 for (
int i = 0; i < Np; i++) {
149 parameterfile.close();
171 assert(fabs(
m_x_min - x_min) < 1.e-12);
174 assert(fabs(
m_x_min - x_min) < 1.e-8);
175 assert(fabs(
m_x_max - x_max) < 1.e-8);
184 for (
int k = 0; k <
m_Np; ++k) {
void get_parameters(double &norm, std::vector< double > &res, std::vector< double > &pole)
static const std::string class_name
void general(const char *format,...)
void set_parameters(const Parameters ¶ms)
int fetch_double(const string &key, double &value) const
std::vector< double > m_pole
int fetch_int(const string &key, int &value) const
Bridge::VerboseLevel m_vl
void crucial(const char *format,...)
int non_zero(const double v)
std::vector< double > m_res
string get_string(const string &key) const
static VerboseLevel set_verbose_level(const std::string &str)