poly_levelset_test.cpp

/*
 * Unit tests for the Regression module PolyLevelset submodule
 * author : Sachin (sthekke@mpi-cbg.de)
 * date : 18.01.2023
 *
 */
 
 
#define BOOST_TEST_DYN_LINK
#include <boost/test/unit_test.hpp>
#include "poly_levelset.hpp"
#include "Vector/vector_dist.hpp"
#include "DMatrix/EMatrix.hpp"
 
BOOST_AUTO_TEST_SUITE( Regression_test )
 
 
 
BOOST_AUTO_TEST_CASE ( PolyLevelset_Sphere )
{
    Box<3,double> domain({-2.0,-2.0,-2.0},{2.0,2.0,2.0});
 
    // Here we define the boundary conditions of our problem
    size_t bc[3]={PERIODIC,PERIODIC,PERIODIC};
 
    // extended boundary around the domain, and the processor domain
    Ghost<3,double> g(0.01);
 
    using vectorType = vector_dist<3,double, aggregate<double, double> >;
 
    vectorType vd(1024,domain,bc,g);
 
    constexpr int mean_curvature = 0;
    constexpr int gauss_curvature = 1;
 
    // Initialize points on sphere
    auto it = vd.getDomainIterator();
    while (it.isNext())
    {
        auto key = it.get();
        double theta = ((double)rand() / RAND_MAX) * M_PI;
        double phi = ((double)rand() / RAND_MAX) * 2.0 * M_PI;
 
        vd.getPos(key)[0] = cos(theta) * sin(phi);
        vd.getPos(key)[1] = cos(theta) * cos(phi);
        vd.getPos(key)[2] = sin(theta);
 
        vd.template getProp<mean_curvature>(key) = 0.0;
        vd.template getProp<gauss_curvature>(key) = 0.0;
        
        ++it;
    }    
    vd.map();
    
    // auto model = PolyLevelset<3>(vd, 1e-4);
/*
    double max_err = -1.0;
    auto it2 = vd.getDomainIterator();
    while (it2.isNext())
    {
        auto key = it2.get();
        Point<3, double> pos = {vd.getPos(key)[0], vd.getPos(key)[1], vd.getPos(key)[2]};
        // vd.template getProp<mean_curvature>(key) = model.estimate_mean_curvature_at(pos);
        // vd.template getProp<gauss_curvature>(key) = model->estimate_gauss_curvature_at(vd.getPos(key));
 
        double val = vd.getProp<mean_curvature>(key);
        double actual = 1.0;
        double err = std::abs(actual - val);
        if (err > max_err) max_err = err;
 
        ++it2;
    }
 
 
    double tolerance = 1e-4;
    bool check;
    if (std::abs(max_err) < tolerance)
        check = true;
    else
        check = false;
    std::cout<<"Max err (poly level) = "<<max_err<<"\n";
    
    BOOST_TEST( check );
    */
    // if(model)
    //     delete model;
 
}
 
 
BOOST_AUTO_TEST_SUITE_END()