mirror of
https://github.com/OrcaSlicer/OrcaSlicer.git
synced 2026-10-09 08:41:14 +00:00
HanifKoh
8a6377f087
* Add Missing Includes Across src/libslic3r Every libslic3r source and header now directly includes the headers declaring what it uses, rather than relying on the precompiled header or transitive includes. Generated with clang-tidy misc-include-cleaner, with libslic3r headers spelled libslic3r/... so they resolve outside the library's private include paths. MultiMaterialSegmentation.hpp, Support/SupportParameters.hpp and Format/STEP.hpp are made self-contained by hand. * Make the libslic3r Headers Compile on Their Own Each now includes, or forward-declares, what it uses instead of relying on what its includers happened to include first. Left out: I18N.hpp, which errors on purpose when included from GUI code, and VoxelizeCSGMesh.hpp and SLA/bicubic.h, which nothing includes and which no longer compile at all. * Add the Includes Missing From the Hand-Fixed libslic3r Headers clang-tidy would not edit these headers while they failed to compile on their own, so the first pass skipped them. With the headers now self-contained, a second pass adds the rest. * Keep Windows Setup Ahead of the Added libslic3r Includes Print.cpp and Thread.cpp open with a _WIN32 block that has to come first; without the precompiled header, Print.cpp otherwise reaches windows.h through OCCT with NONLS defined and boost/regex fails. OpenVDBUtils.cpp and SLA/SupportTreeBuilder.cpp had includes inside #ifndef NOMINMAX, which libslic3r defines on Windows, so those were skipped there. .clang-tidy also ignores the MSVC STL and UCRT internals, Boost.Multiprecision's fwd.hpp and CPython's Windows include directory. * Re-Add libslic3r Includes After the Clipper2 2.0.1 Migration Rebasing onto main took main's version of the files the Clipper2 migration rewrote, so their added includes are restored here, along with includes for main's new code. Clipper2's individual headers are now ignored by clang-tidy: they only build the Z variant through clipper2_z.hpp, which defines USINGZ first, so including clipper.core.h and the like directly broke ClipperZUtils.cpp.
105 lines
2.8 KiB
C++
105 lines
2.8 KiB
C++
// PchipInterpolatorHelper.cpp
|
|
// OrcaSlicer
|
|
//
|
|
// Implementation file for the PchipInterpolatorHelper class
|
|
|
|
#include "PchipInterpolatorHelper.hpp"
|
|
#include <cstddef>
|
|
#include <iterator>
|
|
#include <stdexcept>
|
|
#include <cmath>
|
|
#include <algorithm>
|
|
#include <vector>
|
|
#include <utility>
|
|
|
|
/**
|
|
* @brief Constructs the PCHIP interpolator with given data points.
|
|
* @param x The x-coordinates of the data points.
|
|
* @param y The y-coordinates of the data points.
|
|
*/
|
|
PchipInterpolatorHelper::PchipInterpolatorHelper(const std::vector<double>& x, const std::vector<double>& y) {
|
|
setData(x, y);
|
|
}
|
|
|
|
/**
|
|
* @brief Sets the data points for the interpolator.
|
|
* @param x The x-coordinates of the data points.
|
|
* @param y The y-coordinates of the data points.
|
|
* @throw std::invalid_argument if x and y have different sizes or if they contain fewer than two points.
|
|
*/
|
|
void PchipInterpolatorHelper::setData(const std::vector<double>& x, const std::vector<double>& y) {
|
|
if (x.size() != y.size() || x.size() < 2) {
|
|
throw std::invalid_argument("Input vectors must have the same size and contain at least two points.");
|
|
}
|
|
x_ = x;
|
|
y_ = y;
|
|
sortData();
|
|
computePCHIP();
|
|
}
|
|
|
|
/**
|
|
* @brief Sorts the data points by x-coordinate.
|
|
*/
|
|
void PchipInterpolatorHelper::sortData() {
|
|
std::vector<std::pair<double, double>> data;
|
|
for (size_t i = 0; i < x_.size(); ++i) {
|
|
data.emplace_back(x_[i], y_[i]);
|
|
}
|
|
std::sort(data.begin(), data.end());
|
|
|
|
for (size_t i = 0; i < data.size(); ++i) {
|
|
x_[i] = data[i].first;
|
|
y_[i] = data[i].second;
|
|
}
|
|
}
|
|
|
|
/**
|
|
* @brief Computes the PCHIP coefficients.
|
|
*/
|
|
void PchipInterpolatorHelper::computePCHIP() {
|
|
size_t n = x_.size() - 1;
|
|
h_.resize(n);
|
|
delta_.resize(n);
|
|
d_.resize(n+1);
|
|
|
|
for (size_t i = 0; i < n; ++i) {
|
|
h_[i] = h(i);
|
|
delta_[i] = delta(i);
|
|
}
|
|
|
|
d_[0] = delta_[0];
|
|
d_[n] = delta_[n-1];
|
|
for (size_t i = 1; i < n; ++i) {
|
|
if (delta_[i-1] * delta_[i] > 0) {
|
|
double w1 = 2 * h_[i] + h_[i-1];
|
|
double w2 = h_[i] + 2 * h_[i-1];
|
|
d_[i] = (w1 + w2) / (w1 / delta_[i-1] + w2 / delta_[i]);
|
|
} else {
|
|
d_[i] = 0;
|
|
}
|
|
}
|
|
}
|
|
|
|
/**
|
|
* @brief Interpolates the value at a given point.
|
|
*/
|
|
double PchipInterpolatorHelper::interpolate(double xi) const {
|
|
if (xi <= x_.front()) return y_.front();
|
|
if (xi >= x_.back()) return y_.back();
|
|
|
|
auto it = std::lower_bound(x_.begin(), x_.end(), xi);
|
|
size_t i = std::distance(x_.begin(), it) - 1;
|
|
|
|
double h_i = h_[i];
|
|
double t = (xi - x_[i]) / h_i;
|
|
double t2 = t * t;
|
|
double t3 = t2 * t;
|
|
|
|
double h00 = 2 * t3 - 3 * t2 + 1;
|
|
double h10 = t3 - 2 * t2 + t;
|
|
double h01 = -2 * t3 + 3 * t2;
|
|
double h11 = t3 - t2;
|
|
|
|
return h00 * y_[i] + h10 * h_i * d_[i] + h01 * y_[i+1] + h11 * h_i * d_[i+1];
|
|
}
|