Medusa  1.1
Coordinate Free Mehless Method implementation
test/domains/NURBSShape_test.cpp
#include "medusa/bits/domains/NURBSShape.hpp"
#include "gtest/gtest.h"
namespace mm {
TEST(Domains, NURBSShapeCurve1PatchConstantFill) {
int p = 2;
Range<double> w({1, sqrt(3) / 2, 1, sqrt(3) / 2, 1, sqrt(3) / 2, 1, sqrt(3) / 2, 1});
Range<double> knots({0, 0, 0, PI / 2, PI / 2, PI, PI, 3 * PI / 2,
3 * PI / 2, 2 * PI, 2 * PI, 2 * PI});
Range<Vec3d> cp({Vec3d(1, 0, 0), Vec3d(1, 1, 0), Vec3d(0, 1, 0), Vec3d(-1, 1, 0),
Vec3d(-1, 0, 0), Vec3d(-1, -1, 0), Vec3d(0, -1, 0), Vec3d(1, -1, 0),
Vec3d(1, 0, 0)});
NURBSPatch<Vec3d, Vec1d> curve(cp, w, {knots}, {p});
NURBSShape<Vec3d, Vec1d> shape({curve});
shape.seed(123);
double h = 0.01;
DomainDiscretization<Vec3d> domain = shape.discretizeBoundaryWithStep(h);
// Calculate distances.
KDTree<Vec3d> tree(domain.positions());
Range<double> distances;
double avg_dist = 0;
for (const auto& pt : domain.positions()) {
distances.push_back(sqrt(tree.query(pt, 2).second[1]));
avg_dist += distances[distances.size() - 1];
}
avg_dist /= domain.positions().size();
auto err_avg = std::abs(avg_dist - h) / h;
EXPECT_LT(err_avg, 0.01);
auto err_min = std::abs(*std::min_element(distances.begin() + 1, distances.end() - 1) - h) / h;
EXPECT_LT(err_min, 0.01);
}
TEST(Domains, NURBSShapeCurve1PatchFunctionFill) {
int p = 2;
Range<double> w({1, sqrt(3) / 2, 1, sqrt(3) / 2, 1, sqrt(3) / 2, 1, sqrt(3) / 2, 1});
Range<double> knots({0, 0, 0, PI / 2, PI / 2, PI, PI, 3 * PI / 2,
3 * PI / 2, 2 * PI, 2 * PI, 2 * PI});
Range<Vec3d> cp({Vec3d(1, 0, 0), Vec3d(1, 1, 0), Vec3d(0, 1, 0), Vec3d(-1, 1, 0),
Vec3d(-1, 0, 0), Vec3d(-1, -1, 0), Vec3d(0, -1, 0), Vec3d(1, -1, 0),
Vec3d(1, 0, 0)});
NURBSPatch<Vec3d, Vec1d> curve(cp, w, {knots}, {p});
NURBSShape<Vec3d, Vec1d> shape({curve});
shape.seed(123);
double h = 0.01;
std::function<double(Vec3d)> func_h = [h](Vec3d x) {
return h * (pow((x(0) - 0.5), 2) + x(1) * x(1));
};
DomainDiscretization<Vec3d> domain = shape.discretizeBoundaryWithDensity(func_h, 0);
// Calculate distances.
KDTree<Vec3d> tree(domain.positions());
auto positions = domain.positions();
for (int i = 1; i < positions.size(); i++) {
auto pt = positions[i];
double dist = sqrt(tree.query(pt, 2).second[1]);
double h = func_h(pt);
double err = std::abs(dist - h) / h;
EXPECT_LT(err, 0.05);
}
}
TEST(Domains, NURBSShapeCurveMultiplePatchFill) {
int p = 3;
auto knots = Range<double>{0, 0, 0, 0, 1, 1, 1, 1};
Range<double> ws{1, 1, 1, 1};
Range<Vec2d> cps1{{1, 1}, {2, 3}, {4, 4}, {6, 1}};
Range<Vec2d> cps2{{1, 1}, {0, -3}, {8, -2}, {6, 1}};
NURBSPatch<Vec2d, Vec1d> p1(cps1, {ws}, {knots}, {p});
NURBSPatch<Vec2d, Vec1d> p2(cps2, {ws}, {knots}, {p});
Range<NURBSPatch<Vec2d, Vec1d>> patches{p1, p2};
NURBSShape<Vec2d, Vec1d> shape(patches);
shape.seed(1337);
double h = 0.01;
DomainDiscretization<Vec2d> domain = shape.discretizeBoundaryWithStep(h);
// Calculate distances.
KDTree<Vec2d> tree(domain.positions());
Range<double> distances;
double avg_dist = 0;
for (const auto& pt : domain.positions()) {
distances.push_back(sqrt(tree.query(pt, 2).second[1]));
avg_dist += distances[distances.size() - 1];
}
avg_dist /= domain.positions().size();
auto err_avg = std::abs(avg_dist - h) / h;
EXPECT_LT(err_avg, 0.01);
auto err_min = std::abs(*std::min_element(distances.begin() + 1, distances.end() - 1) - h) / h;
EXPECT_LT(err_min, 0.01);
EXPECT_GT(domain.size(), 1200);
}
TEST(Domains, NURBSShapeSurface1PatchConstantFill) {
auto p = std::array<int, 2>{1, 2};
auto knots = std::array<Range<double>, 2>{Range<double>({0, 0, 1, 1}),
Range<double>({0, 0, 0, 0.25, 0.25, 0.5, 0.5, 0.75, 0.75, 1, 1, 1})};
Range<Range<Vec<double, 4>>> wcp(2);
double x = 0.7071;
wcp[0].push_back(Vec<double, 4>(1, 0, 0, 1));
wcp[0].push_back(Vec<double, 4>(x, x, 0, x));
wcp[0].push_back(Vec<double, 4>(0, 1, 0, 1));
wcp[0].push_back(Vec<double, 4>(-x, x, 0, x));
wcp[0].push_back(Vec<double, 4>(-1, 0, 0, 1));
wcp[0].push_back(Vec<double, 4>(-x, -x, 0, x));
wcp[0].push_back(Vec<double, 4>(0, -1, 0, 1));
wcp[0].push_back(Vec<double, 4>(x, -x, 0, x));
wcp[0].push_back(Vec<double, 4>(1, 0, 0, 1));
wcp[1].push_back(Vec<double, 4>(1, 0, 1, 1));
wcp[1].push_back(Vec<double, 4>(x, x, x, x));
wcp[1].push_back(Vec<double, 4>(0, 1, 1, 1));
wcp[1].push_back(Vec<double, 4>(-x, x, x, x));
wcp[1].push_back(Vec<double, 4>(-1, 0, 1, 1));
wcp[1].push_back(Vec<double, 4>(-x, -x, x, x));
wcp[1].push_back(Vec<double, 4>(0, -1, 1, 1));
wcp[1].push_back(Vec<double, 4>(x, -x, x, x));
wcp[1].push_back(Vec<double, 4>(1, 0, 1, 1));
double h = 0.02;
NURBSPatch<Vec3d, Vec2d> surf(wcp, knots, p);
NURBSShape<Vec3d, Vec2d> nsh({surf});
nsh.seed(123);
DomainDiscretization<Vec3d> domain = nsh.discretizeBoundaryWithStep(h);
// Calculate distances.
KDTree<Vec3d> tree(domain.positions());
Range<double> distances;
double avg_dist = 0;
for (const auto& pt : domain.positions()) {
distances.push_back(sqrt(tree.query(pt, 2).second[1]));
avg_dist += distances[distances.size() - 1];
}
avg_dist /= domain.positions().size();
auto err_avg = std::abs(avg_dist - h) / h;
EXPECT_LT(err_avg, 0.05);
auto err_min = std::abs(*std::min_element(distances.begin() + 1, distances.end() - 1) - h) / h;
EXPECT_LT(err_min, 0.2);
}
TEST(Domains, NURBSShapeSurface1PatchFunctionFill) {
auto p = std::array<int, 2>{1, 2};
auto knots = std::array<Range<double>, 2>{Range<double>({0, 0, 1, 1}),
Range<double>({0, 0, 0, 0.25, 0.25, 0.5, 0.5, 0.75, 0.75, 1, 1, 1})};
Range<Range<Vec<double, 4>>> wcp(2);
double x = 0.7071;
wcp[0].push_back(Vec<double, 4>(1, 0, 0, 1));
wcp[0].push_back(Vec<double, 4>(x, x, 0, x));
wcp[0].push_back(Vec<double, 4>(0, 1, 0, 1));
wcp[0].push_back(Vec<double, 4>(-x, x, 0, x));
wcp[0].push_back(Vec<double, 4>(-1, 0, 0, 1));
wcp[0].push_back(Vec<double, 4>(-x, -x, 0, x));
wcp[0].push_back(Vec<double, 4>(0, -1, 0, 1));
wcp[0].push_back(Vec<double, 4>(x, -x, 0, x));
wcp[0].push_back(Vec<double, 4>(1, 0, 0, 1));
wcp[1].push_back(Vec<double, 4>(1, 0, 1, 1));
wcp[1].push_back(Vec<double, 4>(x, x, x, x));
wcp[1].push_back(Vec<double, 4>(0, 1, 1, 1));
wcp[1].push_back(Vec<double, 4>(-x, x, x, x));
wcp[1].push_back(Vec<double, 4>(-1, 0, 1, 1));
wcp[1].push_back(Vec<double, 4>(-x, -x, x, x));
wcp[1].push_back(Vec<double, 4>(0, -1, 1, 1));
wcp[1].push_back(Vec<double, 4>(x, -x, x, x));
wcp[1].push_back(Vec<double, 4>(1, 0, 1, 1));
double h = 0.02;
std::function<double(Vec3d)> func_h = [h](Vec3d x) {
return h * (pow((x(0) - 0.5), 2) + x(1) * x(1));
};
NURBSPatch<Vec3d, Vec2d> surf(wcp, knots, p);
NURBSShape<Vec3d, Vec2d> nsh({surf});
nsh.seed(123);
DomainDiscretization<Vec3d> domain = nsh.discretizeBoundaryWithDensity(func_h, 0);
// Calculate distances.
KDTree<Vec3d> tree(domain.positions());
auto positions = domain.positions();
for (int i = 1; i < positions.size(); i++) {
auto pt = positions[i];
double dist = sqrt(tree.query(pt, 2).second[1]);
double h = func_h(pt);
double err = std::abs(dist - h) / h;
EXPECT_LT(err, 0.2);
}
}
namespace tp_data {
const int numPatches = 32;
const int numVertices = 306;
int teapotPatches[numPatches][16] = {
{ 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16},
{ 4, 17, 18, 19, 8, 20, 21, 22, 12, 23, 24, 25, 16, 26, 27, 28},
{ 19, 29, 30, 31, 22, 32, 33, 34, 25, 35, 36, 37, 28, 38, 39, 40},
{ 31, 41, 42, 1, 34, 43, 44, 5, 37, 45, 46, 9, 40, 47, 48, 13},
{ 13, 14, 15, 16, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60},
{ 16, 26, 27, 28, 52, 61, 62, 63, 56, 64, 65, 66, 60, 67, 68, 69},
{ 28, 38, 39, 40, 63, 70, 71, 72, 66, 73, 74, 75, 69, 76, 77, 78},
{ 40, 47, 48, 13, 72, 79, 80, 49, 75, 81, 82, 53, 78, 83, 84, 57},
{ 57, 58, 59, 60, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96},
{ 60, 67, 68, 69, 88, 97, 98, 99, 92, 100, 101, 102, 96, 103, 104, 105},
{ 69, 76, 77, 78, 99, 106, 107, 108, 102, 109, 110, 111, 105, 112, 113, 114},
{ 78, 83, 84, 57, 108, 115, 116, 85, 111, 117, 118, 89, 114, 119, 120, 93},
{121, 122, 123, 124, 125, 126, 127, 128, 129, 130, 131, 132, 133, 134, 135, 136},
{124, 137, 138, 121, 128, 139, 140, 125, 132, 141, 142, 129, 136, 143, 144, 133},
{133, 134, 135, 136, 145, 146, 147, 148, 149, 150, 151, 152, 69, 153, 154, 155},
{136, 143, 144, 133, 148, 156, 157, 145, 152, 158, 159, 149, 155, 160, 161, 69},
{162, 163, 164, 165, 166, 167, 168, 169, 170, 171, 172, 173, 174, 175, 176, 177},
{165, 178, 179, 162, 169, 180, 181, 166, 173, 182, 183, 170, 177, 184, 185, 174},
{174, 175, 176, 177, 186, 187, 188, 189, 190, 191, 192, 193, 194, 195, 196, 197},
{177, 184, 185, 174, 189, 198, 199, 186, 193, 200, 201, 190, 197, 202, 203, 194},
{204, 204, 204, 204, 207, 208, 209, 210, 211, 211, 211, 211, 212, 213, 214, 215},
{204, 204, 204, 204, 210, 217, 218, 219, 211, 211, 211, 211, 215, 220, 221, 222},
{204, 204, 204, 204, 219, 224, 225, 226, 211, 211, 211, 211, 222, 227, 228, 229},
{204, 204, 204, 204, 226, 230, 231, 207, 211, 211, 211, 211, 229, 232, 233, 212},
{212, 213, 214, 215, 234, 235, 236, 237, 238, 239, 240, 241, 242, 243, 244, 245},
{215, 220, 221, 222, 237, 246, 247, 248, 241, 249, 250, 251, 245, 252, 253, 254},
{222, 227, 228, 229, 248, 255, 256, 257, 251, 258, 259, 260, 254, 261, 262, 263},
{229, 232, 233, 212, 257, 264, 265, 234, 260, 266, 267, 238, 263, 268, 269, 242},
{270, 270, 270, 270, 279, 280, 281, 282, 275, 276, 277, 278, 271, 272, 273, 274},
{270, 270, 270, 270, 282, 289, 290, 291, 278, 286, 287, 288, 274, 283, 284, 285},
{270, 270, 270, 270, 291, 298, 299, 300, 288, 295, 296, 297, 285, 292, 293, 294},
{270, 270, 270, 270, 300, 305, 306, 279, 297, 303, 304, 275, 294, 301, 302, 271}
};
double teapotVertices[numVertices][3] = {
{ 1.4000, 0.0000, 2.4000},
{ 1.4000, -0.7840, 2.4000},
{ 0.7840, -1.4000, 2.4000},
{ 0.0000, -1.4000, 2.4000},
{ 1.3375, 0.0000, 2.5312},
{ 1.3375, -0.7490, 2.5312},
{ 0.7490, -1.3375, 2.5312},
{ 0.0000, -1.3375, 2.5312},
{ 1.4375, 0.0000, 2.5312},
{ 1.4375, -0.8050, 2.5312},
{ 0.8050, -1.4375, 2.5312},
{ 0.0000, -1.4375, 2.5312},
{ 1.5000, 0.0000, 2.4000},
{ 1.5000, -0.8400, 2.4000},
{ 0.8400, -1.5000, 2.4000},
{ 0.0000, -1.5000, 2.4000},
{-0.7840, -1.4000, 2.4000},
{-1.4000, -0.7840, 2.4000},
{-1.4000, 0.0000, 2.4000},
{-0.7490, -1.3375, 2.5312},
{-1.3375, -0.7490, 2.5312},
{-1.3375, 0.0000, 2.5312},
{-0.8050, -1.4375, 2.5312},
{-1.4375, -0.8050, 2.5312},
{-1.4375, 0.0000, 2.5312},
{-0.8400, -1.5000, 2.4000},
{-1.5000, -0.8400, 2.4000},
{-1.5000, 0.0000, 2.4000},
{-1.4000, 0.7840, 2.4000},
{-0.7840, 1.4000, 2.4000},
{ 0.0000, 1.4000, 2.4000},
{-1.3375, 0.7490, 2.5312},
{-0.7490, 1.3375, 2.5312},
{ 0.0000, 1.3375, 2.5312},
{-1.4375, 0.8050, 2.5312},
{-0.8050, 1.4375, 2.5312},
{ 0.0000, 1.4375, 2.5312},
{-1.5000, 0.8400, 2.4000},
{-0.8400, 1.5000, 2.4000},
{ 0.0000, 1.5000, 2.4000},
{ 0.7840, 1.4000, 2.4000},
{ 1.4000, 0.7840, 2.4000},
{ 0.7490, 1.3375, 2.5312},
{ 1.3375, 0.7490, 2.5312},
{ 0.8050, 1.4375, 2.5312},
{ 1.4375, 0.8050, 2.5312},
{ 0.8400, 1.5000, 2.4000},
{ 1.5000, 0.8400, 2.4000},
{ 1.7500, 0.0000, 1.8750},
{ 1.7500, -0.9800, 1.8750},
{ 0.9800, -1.7500, 1.8750},
{ 0.0000, -1.7500, 1.8750},
{ 2.0000, 0.0000, 1.3500},
{ 2.0000, -1.1200, 1.3500},
{ 1.1200, -2.0000, 1.3500},
{ 0.0000, -2.0000, 1.3500},
{ 2.0000, 0.0000, 0.9000},
{ 2.0000, -1.1200, 0.9000},
{ 1.1200, -2.0000, 0.9000},
{ 0.0000, -2.0000, 0.9000},
{-0.9800, -1.7500, 1.8750},
{-1.7500, -0.9800, 1.8750},
{-1.7500, 0.0000, 1.8750},
{-1.1200, -2.0000, 1.3500},
{-2.0000, -1.1200, 1.3500},
{-2.0000, 0.0000, 1.3500},
{-1.1200, -2.0000, 0.9000},
{-2.0000, -1.1200, 0.9000},
{-2.0000, 0.0000, 0.9000},
{-1.7500, 0.9800, 1.8750},
{-0.9800, 1.7500, 1.8750},
{ 0.0000, 1.7500, 1.8750},
{-2.0000, 1.1200, 1.3500},
{-1.1200, 2.0000, 1.3500},
{ 0.0000, 2.0000, 1.3500},
{-2.0000, 1.1200, 0.9000},
{-1.1200, 2.0000, 0.9000},
{ 0.0000, 2.0000, 0.9000},
{ 0.9800, 1.7500, 1.8750},
{ 1.7500, 0.9800, 1.8750},
{ 1.1200, 2.0000, 1.3500},
{ 2.0000, 1.1200, 1.3500},
{ 1.1200, 2.0000, 0.9000},
{ 2.0000, 1.1200, 0.9000},
{ 2.0000, 0.0000, 0.4500},
{ 2.0000, -1.1200, 0.4500},
{ 1.1200, -2.0000, 0.4500},
{ 0.0000, -2.0000, 0.4500},
{ 1.5000, 0.0000, 0.2250},
{ 1.5000, -0.8400, 0.2250},
{ 0.8400, -1.5000, 0.2250},
{ 0.0000, -1.5000, 0.2250},
{ 1.5000, 0.0000, 0.1500},
{ 1.5000, -0.8400, 0.1500},
{ 0.8400, -1.5000, 0.1500},
{ 0.0000, -1.5000, 0.1500},
{-1.1200, -2.0000, 0.4500},
{-2.0000, -1.1200, 0.4500},
{-2.0000, 0.0000, 0.4500},
{-0.8400, -1.5000, 0.2250},
{-1.5000, -0.8400, 0.2250},
{-1.5000, 0.0000, 0.2250},
{-0.8400, -1.5000, 0.1500},
{-1.5000, -0.8400, 0.1500},
{-1.5000, 0.0000, 0.1500},
{-2.0000, 1.1200, 0.4500},
{-1.1200, 2.0000, 0.4500},
{ 0.0000, 2.0000, 0.4500},
{-1.5000, 0.8400, 0.2250},
{-0.8400, 1.5000, 0.2250},
{ 0.0000, 1.5000, 0.2250},
{-1.5000, 0.8400, 0.1500},
{-0.8400, 1.5000, 0.1500},
{ 0.0000, 1.5000, 0.1500},
{ 1.1200, 2.0000, 0.4500},
{ 2.0000, 1.1200, 0.4500},
{ 0.8400, 1.5000, 0.2250},
{ 1.5000, 0.8400, 0.2250},
{ 0.8400, 1.5000, 0.1500},
{ 1.5000, 0.8400, 0.1500},
{-1.6000, 0.0000, 2.0250},
{-1.6000, -0.3000, 2.0250},
{-1.5000, -0.3000, 2.2500},
{-1.5000, 0.0000, 2.2500},
{-2.3000, 0.0000, 2.0250},
{-2.3000, -0.3000, 2.0250},
{-2.5000, -0.3000, 2.2500},
{-2.5000, 0.0000, 2.2500},
{-2.7000, 0.0000, 2.0250},
{-2.7000, -0.3000, 2.0250},
{-3.0000, -0.3000, 2.2500},
{-3.0000, 0.0000, 2.2500},
{-2.7000, 0.0000, 1.8000},
{-2.7000, -0.3000, 1.8000},
{-3.0000, -0.3000, 1.8000},
{-3.0000, 0.0000, 1.8000},
{-1.5000, 0.3000, 2.2500},
{-1.6000, 0.3000, 2.0250},
{-2.5000, 0.3000, 2.2500},
{-2.3000, 0.3000, 2.0250},
{-3.0000, 0.3000, 2.2500},
{-2.7000, 0.3000, 2.0250},
{-3.0000, 0.3000, 1.8000},
{-2.7000, 0.3000, 1.8000},
{-2.7000, 0.0000, 1.5750},
{-2.7000, -0.3000, 1.5750},
{-3.0000, -0.3000, 1.3500},
{-3.0000, 0.0000, 1.3500},
{-2.5000, 0.0000, 1.1250},
{-2.5000, -0.3000, 1.1250},
{-2.6500, -0.3000, 0.9375},
{-2.6500, 0.0000, 0.9375},
{-2.0000, -0.3000, 0.9000},
{-1.9000, -0.3000, 0.6000},
{-1.9000, 0.0000, 0.6000},
{-3.0000, 0.3000, 1.3500},
{-2.7000, 0.3000, 1.5750},
{-2.6500, 0.3000, 0.9375},
{-2.5000, 0.3000, 1.1250},
{-1.9000, 0.3000, 0.6000},
{-2.0000, 0.3000, 0.9000},
{ 1.7000, 0.0000, 1.4250},
{ 1.7000, -0.6600, 1.4250},
{ 1.7000, -0.6600, 0.6000},
{ 1.7000, 0.0000, 0.6000},
{ 2.6000, 0.0000, 1.4250},
{ 2.6000, -0.6600, 1.4250},
{ 3.1000, -0.6600, 0.8250},
{ 3.1000, 0.0000, 0.8250},
{ 2.3000, 0.0000, 2.1000},
{ 2.3000, -0.2500, 2.1000},
{ 2.4000, -0.2500, 2.0250},
{ 2.4000, 0.0000, 2.0250},
{ 2.7000, 0.0000, 2.4000},
{ 2.7000, -0.2500, 2.4000},
{ 3.3000, -0.2500, 2.4000},
{ 3.3000, 0.0000, 2.4000},
{ 1.7000, 0.6600, 0.6000},
{ 1.7000, 0.6600, 1.4250},
{ 3.1000, 0.6600, 0.8250},
{ 2.6000, 0.6600, 1.4250},
{ 2.4000, 0.2500, 2.0250},
{ 2.3000, 0.2500, 2.1000},
{ 3.3000, 0.2500, 2.4000},
{ 2.7000, 0.2500, 2.4000},
{ 2.8000, 0.0000, 2.4750},
{ 2.8000, -0.2500, 2.4750},
{ 3.5250, -0.2500, 2.4938},
{ 3.5250, 0.0000, 2.4938},
{ 2.9000, 0.0000, 2.4750},
{ 2.9000, -0.1500, 2.4750},
{ 3.4500, -0.1500, 2.5125},
{ 3.4500, 0.0000, 2.5125},
{ 2.8000, 0.0000, 2.4000},
{ 2.8000, -0.1500, 2.4000},
{ 3.2000, -0.1500, 2.4000},
{ 3.2000, 0.0000, 2.4000},
{ 3.5250, 0.2500, 2.4938},
{ 2.8000, 0.2500, 2.4750},
{ 3.4500, 0.1500, 2.5125},
{ 2.9000, 0.1500, 2.4750},
{ 3.2000, 0.1500, 2.4000},
{ 2.8000, 0.1500, 2.4000},
{ 0.0000, 0.0000, 3.1500},
{ 0.0000, -0.0020, 3.1500},
{ 0.0020, 0.0000, 3.1500},
{ 0.8000, 0.0000, 3.1500},
{ 0.8000, -0.4500, 3.1500},
{ 0.4500, -0.8000, 3.1500},
{ 0.0000, -0.8000, 3.1500},
{ 0.0000, 0.0000, 2.8500},
{ 0.2000, 0.0000, 2.7000},
{ 0.2000, -0.1120, 2.7000},
{ 0.1120, -0.2000, 2.7000},
{ 0.0000, -0.2000, 2.7000},
{-0.0020, 0.0000, 3.1500},
{-0.4500, -0.8000, 3.1500},
{-0.8000, -0.4500, 3.1500},
{-0.8000, 0.0000, 3.1500},
{-0.1120, -0.2000, 2.7000},
{-0.2000, -0.1120, 2.7000},
{-0.2000, 0.0000, 2.7000},
{ 0.0000, 0.0020, 3.1500},
{-0.8000, 0.4500, 3.1500},
{-0.4500, 0.8000, 3.1500},
{ 0.0000, 0.8000, 3.1500},
{-0.2000, 0.1120, 2.7000},
{-0.1120, 0.2000, 2.7000},
{ 0.0000, 0.2000, 2.7000},
{ 0.4500, 0.8000, 3.1500},
{ 0.8000, 0.4500, 3.1500},
{ 0.1120, 0.2000, 2.7000},
{ 0.2000, 0.1120, 2.7000},
{ 0.4000, 0.0000, 2.5500},
{ 0.4000, -0.2240, 2.5500},
{ 0.2240, -0.4000, 2.5500},
{ 0.0000, -0.4000, 2.5500},
{ 1.3000, 0.0000, 2.5500},
{ 1.3000, -0.7280, 2.5500},
{ 0.7280, -1.3000, 2.5500},
{ 0.0000, -1.3000, 2.5500},
{ 1.3000, 0.0000, 2.4000},
{ 1.3000, -0.7280, 2.4000},
{ 0.7280, -1.3000, 2.4000},
{ 0.0000, -1.3000, 2.4000},
{-0.2240, -0.4000, 2.5500},
{-0.4000, -0.2240, 2.5500},
{-0.4000, 0.0000, 2.5500},
{-0.7280, -1.3000, 2.5500},
{-1.3000, -0.7280, 2.5500},
{-1.3000, 0.0000, 2.5500},
{-0.7280, -1.3000, 2.4000},
{-1.3000, -0.7280, 2.4000},
{-1.3000, 0.0000, 2.4000},
{-0.4000, 0.2240, 2.5500},
{-0.2240, 0.4000, 2.5500},
{ 0.0000, 0.4000, 2.5500},
{-1.3000, 0.7280, 2.5500},
{-0.7280, 1.3000, 2.5500},
{ 0.0000, 1.3000, 2.5500},
{-1.3000, 0.7280, 2.4000},
{-0.7280, 1.3000, 2.4000},
{ 0.0000, 1.3000, 2.4000},
{ 0.2240, 0.4000, 2.5500},
{ 0.4000, 0.2240, 2.5500},
{ 0.7280, 1.3000, 2.5500},
{ 1.3000, 0.7280, 2.5500},
{ 0.7280, 1.3000, 2.4000},
{ 1.3000, 0.7280, 2.4000},
{ 0.0000, 0.0000, 0.0000},
{ 1.5000, 0.0000, 0.1500},
{ 1.5000, 0.8400, 0.1500},
{ 0.8400, 1.5000, 0.1500},
{ 0.0000, 1.5000, 0.1500},
{ 1.5000, 0.0000, 0.0750},
{ 1.5000, 0.8400, 0.0750},
{ 0.8400, 1.5000, 0.0750},
{ 0.0000, 1.5000, 0.0750},
{ 1.4250, 0.0000, 0.0000},
{ 1.4250, 0.7980, 0.0000},
{ 0.7980, 1.4250, 0.0000},
{ 0.0000, 1.4250, 0.0000},
{-0.8400, 1.5000, 0.1500},
{-1.5000, 0.8400, 0.1500},
{-1.5000, 0.0000, 0.1500},
{-0.8400, 1.5000, 0.0750},
{-1.5000, 0.8400, 0.0750},
{-1.5000, 0.0000, 0.0750},
{-0.7980, 1.4250, 0.0000},
{-1.4250, 0.7980, 0.0000},
{-1.4250, 0.0000, 0.0000},
{-1.5000, -0.8400, 0.1500},
{-0.8400, -1.5000, 0.1500},
{ 0.0000, -1.5000, 0.1500},
{-1.5000, -0.8400, 0.0750},
{-0.8400, -1.5000, 0.0750},
{ 0.0000, -1.5000, 0.0750},
{-1.4250, -0.7980, 0.0000},
{-0.7980, -1.4250, 0.0000},
{ 0.0000, -1.4250, 0.0000},
{ 0.8400, -1.5000, 0.1500},
{ 1.5000, -0.8400, 0.1500},
{ 0.8400, -1.5000, 0.0750},
{ 1.5000, -0.8400, 0.0750},
{ 0.7980, -1.4250, 0.0000},
{ 1.4250, -0.7980, 0.0000}
};
} // namespace tp_data
TEST(Domains, NURBSShapeSurfaceMultiplePatchFill) {
Range<NURBSPatch<Vec3d, Vec2d>> patches;
std::array<int, 2> p{3, 3};
auto knots = std::array<Range<double>, 2>{Range<double>{0, 0, 0, 0, 1, 1, 1, 1},
Range<double>{0, 0, 0, 0, 1, 1, 1, 1}};
Range<Range<double>> w;
w.resize(4);
for (int i = 0; i < 4; i++) {
for (int j = 0; j < 4; j++) {
w[i].push_back(1);
}
}
for (int i = 0; i < tp_data::numPatches; i++) {
Range<Range<Vec3d>> cp;
cp.resize(4);
for (int j = 0; j < 4; j++) {
for (int k = 0; k < 4; k++) {
auto v = tp_data::teapotVertices[tp_data::teapotPatches[i][4 * j + k] - 1];
Vec3d vert(v[0], v[1], v[2]);
cp[j].push_back(vert);
}
}
patches.push_back(NURBSPatch<Vec3d, Vec2d>(cp, w, knots, p));
}
NURBSShape<Vec3d, Vec2d> shape(patches);
shape.seed(1337).boundaryProximity(1e-11);
double h = 0.05;
// Fill
DomainDiscretization<Vec3d> domain = shape.discretizeBoundaryWithStep(h);
// Calculate distances.
KDTree<Vec3d> tree(domain.positions());
Range<double> distances;
double avg_dist = 0;
for (const auto& pt : domain.positions()) {
distances.push_back(sqrt(tree.query(pt, 2).second[1]));
avg_dist += distances[distances.size() - 1];
}
avg_dist /= domain.positions().size();
auto err_avg = std::abs(avg_dist - h) / h;
EXPECT_LT(err_avg, 0.05);
auto err_min = std::abs(*std::min_element(distances.begin() + 1, distances.end() - 1) - h) / h;
EXPECT_LT(err_min, 0.7);
}
} // namespace mm
NURBSShape.hpp
mm
Root namespace for the whole library.
Definition: Gaussian.hpp:14
mm::Vec3d
Vec< double, 3 > Vec3d
Convenience typedef for 3d vector of doubles.
Definition: Vec_fwd.hpp:35
mm::PI
static const double PI
Mathematical constant pi in double precision.
Definition: Config.hpp:44