test/domains/compute_normal_test.cpp

#include "medusa/bits/domains/compute_normal_fwd.hpp"
 
#include "gtest/gtest.h"
 
#include <Eigen/Geometry>
 
namespace mm {
 
TEST(Domains, NormalTest3D) {
    Eigen::Matrix<double, 3, 2> J;
    J.col(1) = Vec3d(0.5, 0.5, 0);
    J.col(0) = Vec3d(1, 0.5, 0);
    Vec3d normal = -J.col(0).cross(J.col(1));
    normal.normalize();
    Vec3d normal2 = surface_fill_internal::compute_normal(J);
    EXPECT_LT((normal - normal2).norm(), 1e-15);
}
 
TEST(Domains, NormalTest1Dto3D) {
    Eigen::Matrix<double, 3, 1> J;
    J.col(0) = Vec3d(0.5, 0.5, 0);
    Vec3d normal = surface_fill_internal::compute_normal(J);
    EXPECT_LT(std::abs(normal.dot(J.col(0))), 1e-15);
}
 
TEST(Domains, NormalTest2D) {
    Eigen::Matrix<double, 2, 1> J;
    J << 0.1, 2.3;
    Vec2d normal(J(1), -J(0)); normal.normalize();
    Vec2d normal2 = surface_fill_internal::compute_normal(J);
    EXPECT_LT((normal - normal2).norm(), 1e-15);
}
 
}  // namespace mm