Commit e58f7b86 authored by LEFEBVREJP email's avatar LEFEBVREJP email
Browse files

Merge branch 'clenshaw' into 'master'

Clenshaw

See merge request !75
parents 73150cfe 61034a1b
Pipeline #59111 passed with stages
in 14 minutes and 30 seconds
......@@ -40,6 +40,48 @@ SparseGrids::SparseGrids(const int input_dimension,
std::fill(mA.begin(), mA.end(), 0.0);
std::fill(mB.begin(), mB.end(), 1.0);
}
else if (type_1d_rule == "gauss-legendre-odd")
{
rule = TasGrid::rule_gausslegendreodd;
std::fill(mA.begin(), mA.end(), -1.0);
std::fill(mB.begin(), mB.end(), 1.0);
}
else if (type_1d_rule == "gauss-patterson")
{
rule = TasGrid::rule_gausspatterson;
std::fill(mA.begin(), mA.end(), -1.0);
std::fill(mB.begin(), mB.end(), 1.0);
}
else if (type_1d_rule == "gauss-chebyshev1")
{
rule = TasGrid::rule_gausschebyshev1;
std::fill(mA.begin(), mA.end(), -1.0);
std::fill(mB.begin(), mB.end(), 1.0);
}
else if (type_1d_rule == "gauss-chebyshev2")
{
rule = TasGrid::rule_gausschebyshev2;
std::fill(mA.begin(), mA.end(), -1.0);
std::fill(mB.begin(), mB.end(), 1.0);
}
else if (type_1d_rule == "gauss-chebyshev1-odd")
{
rule = TasGrid::rule_gausschebyshev1odd;
std::fill(mA.begin(), mA.end(), -1.0);
std::fill(mB.begin(), mB.end(), 1.0);
}
else if (type_1d_rule == "gauss-chebyshev2-odd")
{
rule = TasGrid::rule_gausschebyshev2odd;
std::fill(mA.begin(), mA.end(), -1.0);
std::fill(mB.begin(), mB.end(), 1.0);
}
else if (type_1d_rule == "clenshaw-curtis")
{
rule = TasGrid::rule_clenshawcurtis;
std::fill(mA.begin(), mA.end(), -1.0);
std::fill(mB.begin(), mB.end(), 1.0);
}
else
{
radix_insist(false, type_1d_rule+" is not supported");
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment