rbf_shape_functions_utility.h

Go to the documentation of this file.

//    |  /           |
//    ' /   __| _` | __|  _ \   __|
//    . \  |   (   | |   (   |\__ `
//   _|\_\_|  \__,_|\__|\___/ ____/
//                   Multi-Physics
//
//  License:         BSD License
//                   Kratos default license: kratos/license.txt
//
//  Main authors:    Sebastian Ares de Parga
//
 
#if !defined(KRATOS_rbf_shape_functions_utility_H_INCLUDED)
#define  KRATOS_rbf_shape_functions_utility_H_INCLUDED
 
 
// System includes
 
// External includes
 
// Project includes
#include "containers/array_1d.h"
#include "includes/define.h"
#include "spaces/ublas_space.h"
#include "utilities/dense_householder_qr_decomposition.h"
 
namespace Kratos
{
 
class KRATOS_API(KRATOS_CORE) RBFShapeFunctionsUtility
{
 
public:
 
 
    using DenseSpace = UblasSpace<double, Matrix, Vector>;
 
    using DenseQRPointerType = typename DenseQRDecomposition<DenseSpace>::Pointer;
 
    using KratosCoreQRType = DenseHouseholderQRDecomposition<DenseSpace>;
 
 
 
 
    static double EvaluateRBF(
        const double x,
        const double h);
 
    static void CalculateShapeFunctions(
        const Matrix& rPoints,
        const array_1d<double,3>& rX,
        const double h,
        Vector& rN,
        DenseQRPointerType pDenseQR = nullptr);
 
    static void CalculateShapeFunctions(
        const Matrix& rPoints,
        const array_1d<double,3>& rX,
        Vector& rN,
        DenseQRPointerType pDenseQR = nullptr);
 
    static double CalculateShapeFunctionsAndInterpolation(
        const Matrix& rPoints,
        const array_1d<double,3>& rX,
        const double h,
        Vector& rN,
        Vector& rY);
 
private:
 
    RBFShapeFunctionsUtility(){};
 
    static double CalculateInverseMultiquadricShapeParameter(const Matrix& rPoints);
 
};
 
}  // namespace Kratos.
 
#endif // KRATOS_rbf_shape_functions_utility_H_INCLUDED  defined