|
|
|
#include <array>
|
|
|
|
#include <bitset>
|
|
|
|
#include <iostream>
|
|
|
|
#include <booluarray.hpp>
|
|
|
|
|
|
|
|
|
|
|
|
#include <cstddef>
|
|
|
|
#include <iostream>
|
|
|
|
#include <iomanip>
|
|
|
|
#include <fstream>
|
|
|
|
#include <vector>
|
|
|
|
#include "bernstein.hpp"
|
|
|
|
#include "multiloop.hpp"
|
|
|
|
#include "quadrature_multipoly.hpp"
|
|
|
|
#include "binomial.hpp"
|
|
|
|
|
|
|
|
#include "real.hpp"
|
|
|
|
#include "uvector.hpp"
|
|
|
|
#include "vector"
|
|
|
|
#include "xarray.hpp"
|
|
|
|
#include <chrono>
|
|
|
|
#include <memory>
|
|
|
|
|
|
|
|
#include "organizer/primitive.hpp"
|
|
|
|
#include "organizer/organizer.hpp"
|
|
|
|
|
|
|
|
#include "organizer/blobtree.hpp"
|
|
|
|
|
|
|
|
|
|
|
|
using namespace algoim::Organizer;
|
|
|
|
using namespace algoim;
|
|
|
|
|
|
|
|
void casePolyhedron1()
|
|
|
|
{
|
|
|
|
// std::vector<std::shared_ptr<PrimitiveDesc>> ps;
|
|
|
|
|
|
|
|
// mesh
|
|
|
|
std::vector<uvector3> vertices = {uvector3(-0.8, -0.8, -0.8),
|
|
|
|
uvector3(-0.8, -0.8, 0.8),
|
|
|
|
uvector3(-0.8, 0.8, -0.8),
|
|
|
|
uvector3(-0.8, 0.8, 0.8),
|
|
|
|
uvector3(0.8, -0.8, -0.8),
|
|
|
|
uvector3(0.8, -0.8, 0.8),
|
|
|
|
uvector3(0.8, 0.8, -0.8),
|
|
|
|
uvector3(0.8, 0.8, 0.8)};
|
|
|
|
std::vector<int> indices = {3, 2, 0, 1, // left
|
|
|
|
4, 6, 7, 5, // right
|
|
|
|
6, 2, 3, 7, // top
|
|
|
|
1, 0, 4, 5, // bottom
|
|
|
|
7, 3, 1, 5, // front
|
|
|
|
2, 6, 4, 0}; // back
|
|
|
|
std::vector<int> scan = {4, 8, 12, 16, 20, 24};
|
|
|
|
|
|
|
|
// ps.emplace_back(std::make_shared<MeshDesc>(vertices, indices, scan));
|
|
|
|
// ps.emplace_back(std::make_shared<SphereDesc>(0.8, 0., 1.));
|
|
|
|
// ps.emplace_back(std::make_shared<PowerTensor>(makeSphere(0.2)));
|
|
|
|
// ps.emplace_back(std::make_shared<MeshDesc>(MeshDesc(vertices, indices, scan)));
|
|
|
|
basicTask(std::make_shared<MeshDesc>(MeshDesc(vertices, indices, scan)));
|
|
|
|
}
|
|
|
|
|
|
|
|
void casePolyhedron2()
|
|
|
|
{
|
|
|
|
// std::vector<std::shared_ptr<PrimitiveDesc>> ps;
|
|
|
|
|
|
|
|
// mesh
|
|
|
|
std::vector<uvector3> vertices = {uvector3(-0.8, -0.8, -0.8),
|
|
|
|
uvector3(-0.8, -0.8, 0.8),
|
|
|
|
uvector3(-0.8, 0.8, -0.8),
|
|
|
|
uvector3(-0.8, 0.8, 0.8),
|
|
|
|
uvector3(0.8, -0.8, -0.8),
|
|
|
|
uvector3(0.8, -0.8, 0.8),
|
|
|
|
uvector3(0.8, 0.8, -0.8),
|
|
|
|
uvector3(0.8, 0.8, 0.8)};
|
|
|
|
std::vector<int> indices = {
|
|
|
|
3,
|
|
|
|
2,
|
|
|
|
0,
|
|
|
|
1, // left
|
|
|
|
6,
|
|
|
|
2,
|
|
|
|
3,
|
|
|
|
7 // top
|
|
|
|
};
|
|
|
|
std::vector<int> scan = {4, 8};
|
|
|
|
|
|
|
|
// ps.emplace_back(std::make_shared<MeshDesc>(vertices, indices, scan));
|
|
|
|
// ps.emplace_back(std::make_shared<SphereDesc>(0.8, 0., 1.));
|
|
|
|
// ps.emplace_back(std::make_shared<PowerTensor>(makeSphere(0.2)));
|
|
|
|
// ps.emplace_back(std::make_shared<MeshDesc>(MeshDesc(vertices, indices, scan)));
|
|
|
|
basicTask(std::make_shared<MeshDesc>(MeshDesc(vertices, indices, scan)));
|
|
|
|
}
|
|
|
|
|
|
|
|
void casePolyhedron3()
|
|
|
|
{
|
|
|
|
std::vector<std::shared_ptr<PrimitiveDesc>> primitiveDescriptions;
|
|
|
|
// std::vector<uvector3> vertices = {uvector3(-1.6, 0, 0),
|
|
|
|
// uvector3(-1.6, 0, 1.6),
|
|
|
|
// uvector3(-1.6, 1.6, 0),
|
|
|
|
// uvector3(-1.6, 1.6, 1.6),
|
|
|
|
// uvector3(1.6, 0, 0),
|
|
|
|
// uvector3(1.6, 0, 1.6),
|
|
|
|
// uvector3(1.6, 1.6, 0),
|
|
|
|
// uvector3(1.6, 1.6, 1.6)};
|
|
|
|
std::vector<uvector3> vertices = {uvector3(-0.8, -0.8, -0.8),
|
|
|
|
uvector3(-0.8, -0.8, 0.8),
|
|
|
|
uvector3(-0.8, 0.8, -0.8),
|
|
|
|
uvector3(-0.8, 0.8, 0.8),
|
|
|
|
uvector3(0.8, -0.8, -0.8),
|
|
|
|
uvector3(0.8, -0.8, 0.8),
|
|
|
|
uvector3(0.8, 0.8, -0.8),
|
|
|
|
uvector3(0.8, 0.8, 0.8)};
|
|
|
|
std::vector<int> indices = {3,
|
|
|
|
2,
|
|
|
|
0,
|
|
|
|
1, // left
|
|
|
|
6,
|
|
|
|
2,
|
|
|
|
3,
|
|
|
|
7, // top
|
|
|
|
// 7,
|
|
|
|
// 3,
|
|
|
|
// 1,
|
|
|
|
// 5, // front
|
|
|
|
2,
|
|
|
|
6,
|
|
|
|
4,
|
|
|
|
0}; // back
|
|
|
|
std::vector<int> scan = {4, 8, 12};
|
|
|
|
basicTask(std::make_shared<MeshDesc>(MeshDesc(vertices, indices, scan)));
|
|
|
|
}
|
|
|
|
|
|
|
|
void casePolyhedronSphere()
|
|
|
|
{
|
|
|
|
// PI * r^3 / 6
|
|
|
|
auto phi0 = std::make_shared<SphereDesc>(SphereDesc(0.8, uvector3(-0.8, 0.8, -0.8), 1.));
|
|
|
|
std::vector<uvector3> vertices = {uvector3(-0.8, -0.8, -0.8),
|
|
|
|
uvector3(-0.8, -0.8, 0.8),
|
|
|
|
uvector3(-0.8, 0.8, -0.8),
|
|
|
|
uvector3(-0.8, 0.8, 0.8),
|
|
|
|
uvector3(0.8, -0.8, -0.8),
|
|
|
|
uvector3(0.8, -0.8, 0.8),
|
|
|
|
uvector3(0.8, 0.8, -0.8),
|
|
|
|
uvector3(0.8, 0.8, 0.8)};
|
|
|
|
std::vector<int> indices = {3,
|
|
|
|
2,
|
|
|
|
0,
|
|
|
|
1, // left
|
|
|
|
6,
|
|
|
|
2,
|
|
|
|
3,
|
|
|
|
7, // top
|
|
|
|
2,
|
|
|
|
6,
|
|
|
|
4,
|
|
|
|
0}; // back
|
|
|
|
std::vector<int> scan = {4, 8, 12};
|
|
|
|
std::vector<std::shared_ptr<PrimitiveDesc>> primitiveDescriptions(2);
|
|
|
|
primitiveDescriptions[0] = phi0;
|
|
|
|
primitiveDescriptions[1] = std::make_shared<MeshDesc>(MeshDesc(vertices, indices, scan));
|
|
|
|
basicTask(primitiveDescriptions);
|
|
|
|
}
|
|
|
|
|
|
|
|
void testDeCasteljau()
|
|
|
|
{
|
|
|
|
auto phiDesc = std::make_shared<SphereDesc>(SphereDesc(0.8, uvector3(0), 1.));
|
|
|
|
tensor3 tensor(nullptr, 3), tensor01(nullptr, 3), tensorBernstein(nullptr, 3), transformedTensorBernstein(nullptr, 3);
|
|
|
|
algoim_spark_alloc(real, tensor, tensor01, tensorBernstein, transformedTensorBernstein);
|
|
|
|
makeSphere(*phiDesc, tensor);
|
|
|
|
uvector3 xmax = 1, xmin = -1, range = xmax - xmin;
|
|
|
|
Organizer::detail::powerTransformation(range, xmin, tensor, tensor01);
|
|
|
|
Organizer::detail::power2BernsteinTensor(tensor01, tensorBernstein);
|
|
|
|
|
|
|
|
// bernstein::deCasteljau(tensorBernstein, uvector3(0), uvector3(0.2690598923241497 + 0.0000001),
|
|
|
|
// transformedTensorBernstein);
|
|
|
|
bernstein::deCasteljau(tensorBernstein,
|
|
|
|
uvector3(0.2690598923241497 + 0.0000001),
|
|
|
|
uvector3(1 - 0.2690598923241497 - 0.0000001),
|
|
|
|
transformedTensorBernstein);
|
|
|
|
uvector3 testX = uvector3(0.5);
|
|
|
|
|
|
|
|
int sign = bernstein::uniformSign(transformedTensorBernstein);
|
|
|
|
real evalX = evalPower(tensor, testX);
|
|
|
|
|
|
|
|
std::cout << "evalX original rep = " << evalX << std::endl;
|
|
|
|
evalX = bernstein::evalBernsteinPoly(tensorBernstein, testX);
|
|
|
|
std::cout << "evalX bernstein within 0-1 = " << evalX << std::endl;
|
|
|
|
evalX = bernstein::evalBernsteinPoly(transformedTensorBernstein, testX);
|
|
|
|
std::cout << "evalX bernstein after Decasteljau = " << evalX << std::endl;
|
|
|
|
std::cout << "sign = " << sign << std::endl;
|
|
|
|
}
|
|
|
|
|
|
|
|
void testBlob()
|
|
|
|
{
|
|
|
|
orgzer::Blob blob = {5, 3, 4, 5, 6};
|
|
|
|
|
|
|
|
// std::cout << blob.isPrimitive << std::endl;
|
|
|
|
// std::cout << blob.nodeOp << std::endl;
|
|
|
|
// std::cout << blob.ignoreMod << std::endl;
|
|
|
|
// std::cout << blob.isLeft << std::endl;
|
|
|
|
// std::cout << blob.ancestor << std::endl;
|
|
|
|
std::cout << sizeof(blob) << std::endl;
|
|
|
|
std::cout << sizeof(int) << std::endl;
|
|
|
|
}
|
|
|
|
|
|
|
|
void testPrimitive()
|
|
|
|
{
|
|
|
|
// casePolyhedron1();
|
|
|
|
// casePolyhedronSphere();
|
|
|
|
// testDeCasteljau();
|
|
|
|
testBlob();
|
|
|
|
}
|