GBasis_3D.hpp

Go to the documentation of this file.

// _LICENSE_HEADER_
//
// Copyright (C) 2019 - 2025.
// Terms register on the GPL-3.0 license.
//
// This file can be redistributed and/or modified under the license terms.
//
// See top level LICENSE file for more details.
//
// This file can be used citing references in CITATION.cff file.
 
#ifndef __GBasis_3D_HPP
#define __GBasis_3D_HPP
 
#include "GBasis_Data.hpp"
#include "Monomials_3D.hpp"
 
namespace Polydim
{
namespace Utilities
{
class GBasis_3D final
{
  private:
    Polydim::Utilities::Monomials_3D monomials;
 
    std::vector<Eigen::Vector4i> VectorDecompositionIndices(const Polydim::Utilities::GBasis_Data &data,
                                                            const Eigen::VectorXi &expo) const;
 
  public:
    Polydim::Utilities::GBasis_Data Compute(const unsigned int polynomial_degree);
 
    std::vector<Eigen::MatrixXd> VanderGBigOPlus(const Polydim::Utilities::GBasis_Data &data, const Eigen::MatrixXd &vander) const;
 
    std::vector<std::vector<Eigen::MatrixXd>> VectorDecomposition(const Polydim::Utilities::GBasis_Data &data) const
    {
        std::vector<std::vector<Eigen::MatrixXd>> result(3);
        for (unsigned int i = 0; i < data.Dimension; i++)
        {
            result[i].resize(2);
            result[i][0] = data.VectorDecomposition[i][0];
            result[i][1].setZero(data.Nk, 2 * data.Nkm1 + data.DimFirstBasis);
            result[i][1] << data.VectorDecomposition[i][1], data.VectorDecomposition[i][2], data.VectorDecomposition[i][3];
        }
        return result;
    }
};
} // namespace Utilities
} // namespace Polydim
 
#endif