test/approximations/Multiquadric_test.cpp

#include "medusa/bits/approximations/Multiquadric.hpp"
 
#include "gtest/gtest.h"
 
namespace mm {
 
TEST(Approximations, MultiQuadraticEvaluate) {
    double s = 1/0.66;
    Multiquadric<double> mq(s);
    double r1 = 0;
    double r2 = 1.17;
    EXPECT_NEAR(1, mq(r1), 1e-15);
    EXPECT_NEAR(1.2286789653933200, mq(r2), 1e-15);
 
    EXPECT_NEAR(0.21780000000000000, mq(r1, 1), 1e-15);
    EXPECT_NEAR(0.17726355389365580, mq(r2, 1), 1e-15);
 
    EXPECT_NEAR(-0.0474368400000000, mq(r1, 2), 1e-15);
    EXPECT_NEAR(-0.025574107170419560, mq(r2, 2), 1e-15);
 
    EXPECT_NEAR(-0.01047023439346170, mq(r2, 6), 1e-15);
}
 
TEST(Approximations, DISABLED_MultiQuadraticUsageExample) {
    Multiquadric<double> m(1.324);
    double val = m(4.0);     // value at 4.0 (this is the squared radius)
    double der = m(0.0, 2);  // second derivative at 0.0
    std::cout << m << std::endl;
    std::cout << val << der;  // to not be unused
}
 
}  // namespace mm