ReducedProblem.H

Go to the documentation of this file.

/*---------------------------------------------------------------------------*\
     ██╗████████╗██╗  ██╗ █████╗  ██████╗ █████╗       ███████╗██╗   ██╗
     ██║╚══██╔══╝██║  ██║██╔══██╗██╔════╝██╔══██╗      ██╔════╝██║   ██║
     ██║   ██║   ███████║███████║██║     ███████║█████╗█████╗  ██║   ██║
     ██║   ██║   ██╔══██║██╔══██║██║     ██╔══██║╚════╝██╔══╝  ╚██╗ ██╔╝
     ██║   ██║   ██║  ██║██║  ██║╚██████╗██║  ██║      ██║      ╚████╔╝
     ╚═╝   ╚═╝   ╚═╝  ╚═╝╚═╝  ╚═╝ ╚═════╝╚═╝  ╚═╝      ╚═╝       ╚═══╝
 
 * In real Time Highly Advanced Computational Applications for Finite Volumes
 * Copyright (C) 2017 by the ITHACA-FV authors
-------------------------------------------------------------------------------
License
    This file is part of ITHACA-FV
    ITHACA-FV is free software: you can redistribute it and/or modify
    it under the terms of the GNU Lesser General Public License as published by
    the Free Software Foundation, either version 3 of the License, or
    (at your option) any later version.
    ITHACA-FV 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 Lesser General Public License for more details.
    You should have received a copy of the GNU Lesser General Public License
    along with ITHACA-FV. If not, see <http://www.gnu.org/licenses/>.
Class
    reducedProblem
Description
    A general reduced problem class is implemented
SourceFiles
    reducedProblem.C
\*---------------------------------------------------------------------------*/
 
 
#ifndef reducedProblem_H
#define reducedProblem_H
 
#include "fvCFD.H"
#include "IOmanip.H"
#include "reductionProblem.H"
#include <Eigen/Eigen>
#include "newton_argument.H"
#include "Foam2Eigen.H"
 
 
/*---------------------------------------------------------------------------*\
                        Class reducedProblem Declaration
\*---------------------------------------------------------------------------*/
class reducedProblem
{
    public:
        // Constructors
        reducedProblem();
        virtual ~reducedProblem() {};
 
        explicit reducedProblem(reductionProblem& problem);
 
        // Variables
        scalar nu;
 
        reductionProblem* problem;
 
        // Functions
        virtual void solveOnline();
 
        static Eigen::MatrixXd solveLinearSys(List<Eigen::MatrixXd> LinSys,
                                              Eigen::MatrixXd x,
                                              Eigen::VectorXd& residual, const Eigen::MatrixXd& bc = Eigen::MatrixXd::Zero(0,
                                                      0), const std::string solverType = "fullPivLu");
        static Eigen::MatrixXd solveLinearSys(List<Eigen::MatrixXd> LinSys,
                                              Eigen::MatrixXd x,
                                              Eigen::VectorXd& residual, const std::string solverType);
};
 
 
class onlineInterp
{
    public:
        // Constructors
        onlineInterp();
        virtual ~onlineInterp() {};
 
        int Nmu_samples;
 
        //--------------------------------------------------------------------------
        Eigen::MatrixXd getInterpCoeffRBF(std::vector<SPLINTER::RBFSpline> rbfVec,
                                          Eigen::MatrixXd mu_interp);
 
        //--------------------------------------------------------------------------
        Eigen::MatrixXd getInterpCoeffSPL(std::vector<SPLINTER::BSpline> splVec,
                                          Eigen::MatrixXd mu_interp);
};
 
 
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
#endif