dense_svd_decomposition.h

Go to the documentation of this file.

//    |  /           |
//    ' /   __| _` | __|  _ \   __|
//    . \  |   (   | |   (   |\__ `
//   _|\_\_|  \__,_|\__|\___/ ____/
//                   Multi-Physics
//
//  License:         BSD License
//                   Kratos default license: kratos/license.txt
//
//  Main authors:    Ruben Zorrilla
//
 
#if !defined(KRATOS_DENSE_SVD_H_INCLUDED)
#define KRATOS_DENSE_SVD_H_INCLUDED
 
// External includes
 
// Project includes
#include "includes/define.h"
#include "includes/kratos_parameters.h"
 
namespace Kratos {
 
 
 
 
 
 
template<class TDenseSpaceType>
class DenseSingularValueDecomposition
{
public:
 
 
    KRATOS_CLASS_POINTER_DEFINITION(DenseSingularValueDecomposition);
 
    typedef typename TDenseSpaceType::DataType DataType;
    typedef typename TDenseSpaceType::VectorType VectorType;
    typedef typename TDenseSpaceType::MatrixType MatrixType;
 
 
    DenseSingularValueDecomposition() = default;
 
    virtual ~DenseSingularValueDecomposition() = default;
 
 
 
 
    static std::string Name()
    {
        return "dense_singular_value_decomposition";
    }
 
    virtual void Compute(
        MatrixType& rInputMatrix,
        Parameters Settings) = 0;
 
    virtual void Compute(
        MatrixType& rInputMatrix,
        VectorType& rVectorS,
        MatrixType& rMatrixU,
        MatrixType& rMatrixV,
        Parameters Settings) = 0;
 
    virtual void MatrixU(MatrixType& rMatrixU) = 0;
 
    virtual void MatrixV(MatrixType& rMatrixV) = 0;
 
    virtual void SingularValues(VectorType& rVectorS) = 0;
 
    virtual std::size_t NonZeroSingularValues() = 0;
 
    virtual void SetThreshold(const double RelTolerance) = 0;
 
    virtual std::size_t Rank() = 0;
 
    virtual void PrintInfo(std::ostream &rOStream) const = 0;
 
 
 
 
 
 
 
 
 
};
 
 
 
 
 
} // namespace Kratos
 
#endif // defined(KRATOS_DENSE_SVD_H_INCLUDED)