120 lines
4.4 KiB
C++
120 lines
4.4 KiB
C++
/***************************************************************************
|
|
* Copyright (c) 2017 Lorenz Lechner *
|
|
* *
|
|
* This file is part of the FreeCAD CAx development system. *
|
|
* *
|
|
* This library is free software; you can redistribute it and/or *
|
|
* modify it under the terms of the GNU Library General Public *
|
|
* License as published by the Free Software Foundation; either *
|
|
* version 2 of the License, or (at your option) any later version. *
|
|
* *
|
|
* This library is distributed in the hope that it will be useful, *
|
|
* but WITHOUT ANY WARRANTY; without even the implied warranty of *
|
|
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the *
|
|
* GNU Library General Public License for more details. *
|
|
* *
|
|
* You should have received a copy of the GNU Library General Public *
|
|
* License along with this library; see the file COPYING.LIB. If not, *
|
|
* write to the Free Software Foundation, Inc., 59 Temple Place, *
|
|
* Suite 330, Boston, MA 02111-1307, USA *
|
|
* *
|
|
***************************************************************************/
|
|
|
|
|
|
|
|
#ifndef NURBS_H
|
|
#define NURBS_H
|
|
|
|
#include <Eigen/Geometry>
|
|
#include <Eigen/IterativeLinearSolvers>
|
|
#include <Eigen/SparseCore>
|
|
#include <tuple>
|
|
|
|
namespace nurbs{
|
|
|
|
using trip = Eigen::Triplet<double>;
|
|
using spMat = Eigen::SparseMatrix<double>;
|
|
|
|
struct NurbsBase2D
|
|
{
|
|
//
|
|
NurbsBase2D(){}
|
|
NurbsBase2D(Eigen::VectorXd u_knots, Eigen::VectorXd v_knots,
|
|
Eigen::VectorXd weights,
|
|
int degree_u=3, int degree_v=3);
|
|
int degree_u;
|
|
int degree_v;
|
|
Eigen::VectorXd u_knots;
|
|
Eigen::VectorXd v_knots;
|
|
Eigen::VectorXd weights;
|
|
|
|
std::vector<std::function<double(double)>> u_functions;
|
|
std::vector<std::function<double(double)>> v_functions;
|
|
|
|
std::vector<std::function<double(double)>> Du_functions;
|
|
std::vector<std::function<double(double)>> Dv_functions;
|
|
|
|
std::vector<std::function<double(double)>> DDu_functions;
|
|
std::vector<std::function<double(double)>> DDv_functions;
|
|
|
|
void computeFirstDerivatives();
|
|
void computeSecondDerivatives();
|
|
|
|
Eigen::VectorXd getInfluenceVector(Eigen::Vector2d u);
|
|
spMat getInfluenceMatrix(Eigen::Matrix<double, Eigen::Dynamic, 2> U);
|
|
|
|
Eigen::VectorXd getDuVector(Eigen::Vector2d u);
|
|
spMat getDuMatrix(Eigen::Matrix<double, Eigen::Dynamic, 2> U);
|
|
|
|
Eigen::VectorXd getDvVector(Eigen::Vector2d u);
|
|
spMat getDvMatrix(Eigen::Matrix<double, Eigen::Dynamic, 2> U);
|
|
|
|
Eigen::Matrix<double, Eigen::Dynamic, 2> getUVMesh(int num_u_points, int num_v_points);
|
|
|
|
std::tuple<NurbsBase2D, Eigen::MatrixXd> interpolateUBS(
|
|
Eigen::Matrix<double, Eigen::Dynamic, 3> poles,
|
|
int degree_u,
|
|
int degree_v,
|
|
int num_u_poles,
|
|
int num_v_poles,
|
|
int num_u_points,
|
|
int num_v_points);
|
|
};
|
|
|
|
struct NurbsBase1D
|
|
{
|
|
NurbsBase1D(){}
|
|
NurbsBase1D(Eigen::VectorXd u_knots, Eigen::VectorXd weights, int degree_u=3);
|
|
int degree_u;
|
|
Eigen::VectorXd u_knots;
|
|
Eigen::VectorXd weights;
|
|
std::vector<std::function<double(double)>> u_functions;
|
|
std::vector<std::function<double(double)>> Du_functions;
|
|
std::vector<std::function<double(double)>> DDu_functions;
|
|
|
|
void computeFirstDerivatives();
|
|
void computeSecondDerivatives();
|
|
|
|
Eigen::VectorXd getInfluenceVector(double u);
|
|
spMat getInfluenceMatrix(Eigen::VectorXd u);
|
|
|
|
Eigen::VectorXd getDuVector(double u);
|
|
spMat getDuMatrix(Eigen::VectorXd u);
|
|
|
|
static Eigen::VectorXd getKnotSequence(double u_min, double u_max, int deg, int num_poles);
|
|
static Eigen::VectorXd getWeightList(Eigen::VectorXd knots, int u_deg);
|
|
|
|
Eigen::VectorXd getUMesh(int num_u_points);
|
|
|
|
std::tuple<NurbsBase1D, Eigen::Matrix<double, Eigen::Dynamic, 3>> interpolateUBS(
|
|
Eigen::Matrix<double, Eigen::Dynamic, 3> poles,
|
|
int degree,
|
|
int num_u_poles,
|
|
int num_u_points);
|
|
|
|
};
|
|
|
|
}
|
|
|
|
#endif
|