NurbsUtils.h

Go to the documentation of this file.

// This file is part of the MercuryDPM project (https://www.mercurydpm.org).
// Copyright (c), The MercuryDPM Developers Team. All rights reserved.
// License: BSD 3-Clause License; see the LICENSE file in the root directory.
 
#ifndef MERCURYDPM_NURBSUTILS_H
#define MERCURYDPM_NURBSUTILS_H
 
#include <vector>
#include <algorithm>
#include "Math/Vector.h"
 
namespace NurbsUtils
{
 
 
template <typename T>
class array2 {
public:
    array2(size_t nRows, size_t nCols, T fillValue = 0.0)
    : rows(nRows), cols(nCols) {
        data.resize(rows * cols, fillValue);
    }
 
    T operator()(size_t row, size_t col) const {
        return data[row*cols + col];
    }
 
    T& operator()(size_t row, size_t col) {
        return data[row*cols + col];
    }
 
private:
    size_t rows, cols;
    std::vector<T> data;
};
 
bool isKnotVectorMonotonic(const std::vector<double>& knots);
 
 
bool close(double a, double b, double eps = std::numeric_limits<double>::epsilon());
 
int findSpan(int degree, const std::vector<double>& knots, double u);
 
 
double bsplineOneBasis(int i, int deg, const std::vector<double>& U, double u);
 
void bsplineBasis(int deg, int span, const std::vector<double>& knots, double u, std::vector<double>& N);
 
void bsplineDerBasis(int deg, int span, const std::vector<double>& knots, double u,
                     int nDers, std::vector<std::vector<double>>& ders);
 
std::vector<Mdouble> createUniformKnotVector(unsigned int numberOfControlPoints, unsigned int degree, bool clampedAtStart, bool clampedAtEnd);
 
void normalizeKnotVector(std::vector<Mdouble>& knots);
 
void extendKnotVector(std::vector<Mdouble>& knots, unsigned int degree, unsigned int numStart, unsigned int numEnd, bool forceBothEndsUniform = false);
 
Vec3D evaluate(Mdouble u, Mdouble v, std::vector<Mdouble> knotsU, std::vector<Mdouble> knotsV,
               std::vector<std::vector<Vec3D>> controlPoints, std::vector<std::vector<Mdouble>> weights);
}
 
#endif //MERCURYDPM_NURBSUTILS_H