allemang.d/toddcox-visualize
1
0
Fork 0
mirror of https://github.com/allemangD/toddcox-visualize.git synced 2025-11-10 03:52:48 -05:00
Code Issues Packages Projects Releases Wiki Activity

WIP: Separate points, cells logic

Browse Source
This commit is contained in:
David Allemang
2022-02-24 16:56:23 -05:00
parent d99e824a06
commit 7a95155919
1 changed files with 43 additions and 34 deletions
Download Patch File Download Diff File Expand all files Collapse all files

77
src/geometrytest.cpp
View File

@@ -9,53 +9,62 @@
#include <fstream>
#include <iostream>
int main() {
std::vector<int> symbol = {5, 3, 3};
vec4 root{0.80, 0.02, 0.02, 0.02};
template<int N>
using ProjectiveNf = Eigen::Transform<float, N, Eigen::Projective>;
template<int N>
using VectorNf = Eigen::Vector<float, N>;
auto group = tc::schlafli(symbol);
auto gens = generators(group);
auto combos = combinations(gens, 2);
const auto &inds = merge<3>(hull<3>(group, combos, {
{0, 1},
}));
auto cosets = group.solve();
auto mirrors = mirror<4>(group);
template<int N>
Eigen::Matrix<float, N, Eigen::Dynamic> make_points(
const tc::Group &group,
const Eigen::Vector<float, N> &root
) {
auto mirrors = mirror<N>(group);
auto corners = plane_intersections(mirrors);
auto start = barycentric(corners, root);
using Projective = Eigen::Transform<float, 4, Eigen::Projective>;
Projective transform = Projective::Identity();
auto cosets = group.solve();
auto higher = cosets.path.walk<vec4, vec4>(start, mirrors, reflect<vec4>);
std::transform(
higher.begin(), higher.end(), higher.begin(),
[&](const vec4 &v) {
return (transform * v.homogeneous()).hnormalized();
}
);
auto verts = cosets.path.walk<VectorNf<N>, VectorNf<N>>(start, mirrors, reflect<vec4>);
// std::vector<vec4> lower(higher.size());
// std::transform(higher.begin(), higher.end(), lower.begin(), stereo<4>);
// const auto &verts = lower;
const auto &verts = higher;
using PT = Eigen::Matrix<float, 4, Eigen::Dynamic>;
using CT = Eigen::Matrix<unsigned, 3, Eigen::Dynamic>;
PT points(4, verts.size());
Eigen::Matrix<float, N, Eigen::Dynamic> points(root.size(), verts.size());
std::copy(verts.begin(), verts.end(), points.colwise().begin());
CT cells = inds;
return points;
}
template<int N>
Eigen::Matrix<unsigned, N, Eigen::Dynamic> make_cells(
const tc::Group &group,
const std::vector<std::vector<int>> &exclude = {}
) {
auto gens = generators(group);
auto combos = combinations(gens, N - 1);
Eigen::Matrix<unsigned, N, Eigen::Dynamic> cells = merge<N>(hull<N>(group, combos, exclude));
return cells;
}
int main() {
std::vector<int> symbol = {5, 3, 3};
auto group = tc::schlafli(symbol);
vec4 root{0.80, 0.02, 0.02, 0.02};
auto points = make_points(group, root);
auto cells = make_cells<3>(group, {{0, 1}});
ml::Mesh mesh(points, cells);
// auto transform = ProjectiveNf<vec4::RowsAtCompileTime>::Identity();
// transform.translation() << 0, 0.5, 0, 0;
// points = (transform * points.colwise().homogeneous()).colwise().hnormalized();
std::cout << points.rows() << " " << points.cols() << std::endl;
std::cout << cells.rows() << " " << cells.cols() << std::endl;
ml::Mesh mesh(points, cells);
ml::write(mesh, std::ofstream("dodeca.pak", std::ios::out | std::ios::binary));
return EXIT_SUCCESS;