jaseg/kicoil
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

Poke claude until CGAL's WASM floating-point rounding mode incompatibility is worked around

Browse source
This commit is contained in:
jaseg 2025-12-15 22:48:27 +01:00
parent d7dd0e363b
commit e1f795a04d
2 changed files with 20 additions and 8 deletions
Download patch file Download diff file Expand all files Collapse all files

2
cgal_skeleton_core/Makefile
View file

@ -21,7 +21,7 @@ WASI_SOURCES := \
exception_stubs.cpp \ exception_stubs.cpp \
INCLUDES := -I$(CACHEDIR)/$(CGAL)/include -I$(CACHEDIR)/$(BOOST) INCLUDES := -I$(CACHEDIR)/$(CGAL)/include -I$(CACHEDIR)/$(BOOST)
CXXFLAGS := -std=c++20 -g -Wall -Wextra -O2 CXXFLAGS := -std=c++20 -g -Wall -Wextra -O2 -DCGAL_ALWAYS_ROUND_TO_NEAREST
WASI_CXXFLAGS := -DCGAL_ALWAYS_ROUND_TO_NEAREST -DCGAL_CORE_USE_BOOST_BACKEND -DCGAL_DISABLE_GMP -DCGAL_USE_BOOST_MP -DFE_UPWARD=FE_TONEAREST -D_WASI_EMULATED_PROCESS_CLOCKS $(CXXFLAGS) WASI_CXXFLAGS := -DCGAL_ALWAYS_ROUND_TO_NEAREST -DCGAL_CORE_USE_BOOST_BACKEND -DCGAL_DISABLE_GMP -DCGAL_USE_BOOST_MP -DFE_UPWARD=FE_TONEAREST -D_WASI_EMULATED_PROCESS_CLOCKS $(CXXFLAGS)
LDFLAGS := LDFLAGS :=
HOST_LDFLAGS := $(LDFLAGS) -lgmp -lmpfr HOST_LDFLAGS := $(LDFLAGS) -lgmp -lmpfr

26
cgal_skeleton_core/skeleton_wrapper.cpp
View file

@ -1,12 +1,25 @@
#include <CGAL/Exact_predicates_inexact_constructions_kernel.h> #include <CGAL/Simple_cartesian.h>
#include <CGAL/Polygon_2.h> #include <CGAL/Polygon_2.h>
#include <CGAL/create_straight_skeleton_2.h> #include <CGAL/create_straight_skeleton_2.h>
#include <CGAL/number_utils.h>
#include <boost/multiprecision/cpp_int.hpp>
#include <boost/multiprecision/number.hpp>
#include <boost/multiprecision/rational_adaptor.hpp>
#include <iostream> #include <iostream>
#include <sstream> #include <sstream>
#include <string> #include <string>
#include <memory> #include <memory>
#include <iomanip>
typedef CGAL::Exact_predicates_inexact_constructions_kernel K; // Use exact rational arithmetic throughout to avoid filtered predicates
typedef boost::multiprecision::number<boost::multiprecision::rational_adaptor<boost::multiprecision::cpp_int_backend<>>> Exact_NT;
typedef CGAL::Simple_cartesian<Exact_NT> Exact_K;
// Input kernel with doubles for reading coordinates
typedef CGAL::Simple_cartesian<double> Input_K;
// Use exact kernel for computation
typedef Exact_K K;
typedef K::Point_2 Point; typedef K::Point_2 Point;
typedef CGAL::Polygon_2<K> Polygon_2; typedef CGAL::Polygon_2<K> Polygon_2;
typedef CGAL::Straight_skeleton_2<K> Ss; typedef CGAL::Straight_skeleton_2<K> Ss;
@ -14,7 +27,6 @@ typedef std::shared_ptr<Ss> SsPtr;
int main() int main()
{ {
// Read polygon from stdin
Polygon_2 poly; Polygon_2 poly;
std::string line; std::string line;
@ -26,7 +38,8 @@ int main()
std::istringstream iss(line); std::istringstream iss(line);
double x, y; double x, y;
if (iss >> x >> y) { if (iss >> x >> y) {
poly.push_back(Point(x, y)); // Convert double to exact rational
poly.push_back(Point(Exact_NT(x), Exact_NT(y)));
} else { } else {
std::cerr << "Error: Invalid input line: " << line << std::endl; std::cerr << "Error: Invalid input line: " << line << std::endl;
return EXIT_FAILURE; return EXIT_FAILURE;
@ -38,13 +51,11 @@ int main()
return EXIT_FAILURE; return EXIT_FAILURE;
} }
// Ensure counter-clockwise orientation
if (!poly.is_counterclockwise_oriented()) { if (!poly.is_counterclockwise_oriented()) {
poly.reverse_orientation(); poly.reverse_orientation();
} }
// Create interior straight skeleton SsPtr ss = CGAL::create_interior_straight_skeleton_2(poly.vertices_begin(), poly.vertices_end(), K());
SsPtr ss = CGAL::create_interior_straight_skeleton_2(poly.vertices_begin(), poly.vertices_end());
if (!ss) { if (!ss) {
std::cerr << "Error: Failed to create straight skeleton" << std::endl; std::cerr << "Error: Failed to create straight skeleton" << std::endl;
@ -64,6 +75,7 @@ int main()
auto v_target = he->opposite()->vertex(); auto v_target = he->opposite()->vertex();
// Get times (distance from boundary) // Get times (distance from boundary)
// Convert from exact rational to double for output
double t1 = CGAL::to_double(v_source->time()); double t1 = CGAL::to_double(v_source->time());
double t2 = CGAL::to_double(v_target->time()); double t2 = CGAL::to_double(v_target->time());