dev_fdaPDE/Descent_Direction_Factory.h at master · simonepanzeri/dev_fdaPDE · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
//
// Created by simonepanzeri on 02/12/2021.
//

#ifndef DEV_FDAPDE_DESCENT_DIRECTION_FACTORY_H
#define DEV_FDAPDE_DESCENT_DIRECTION_FACTORY_H

#include "FdaPDE.h"

//! @brief A Factory class: a class for the choice of the direction for the optimization algorithm.
template<UInt ORDER, UInt mydim, UInt ndim>
class DescentDirection_factory
{
public:
    //! A method that builds a pointer to the right object for the direction choice, taking as parameters a string and others objects needed for constructor.
    static std::unique_ptr<DirectionBase<ORDER,  mydim,  ndim>>
    createDirectionSolver(const DataProblem<ORDER, mydim, ndim>& dp,
                          const FunctionalProblem<ORDER, mydim, ndim>& fp, const std::string& d)
    {
        if (d=="Gradient")
            return make_unique<DirectionGradient<ORDER,mydim,ndim>>(fp);
        else if (d=="BFGS")
            return make_unique<DirectionBFGS<ORDER,mydim,ndim>>(fp, dp.getNumNodes());
        else{
            std::cout << "Unknown direction option - using gradient direction" << std::endl;
            //Rprintf("Unknown direction option - using gradient direction");

            return make_unique<DirectionGradient<ORDER,mydim,ndim>>(fp);
        }
    }

};

//! ####################################################################################################################
//! ######################################## SPACE-TIME PROBLEM ########################################################
//! ####################################################################################################################

//! @brief A Factory class: a class for the choice of the direction for the optimization algorithm (spatio-temporal setting).
template<UInt ORDER, UInt mydim, UInt ndim>
class DescentDirection_factory_time{
public:
    //! A method that builds a pointer to the right object for the direction choice, taking as parameters a string and others objects needed for constructor.
    static std::unique_ptr<DirectionBase<ORDER, mydim, ndim, FunctionalProblem_time<ORDER, mydim, ndim>>>
    createDirectionSolver(const DataProblem_time<ORDER, mydim, ndim>& dp,
                          const FunctionalProblem_time<ORDER, mydim, ndim>& fp, const std::string& d)
    {
        if (d=="Gradient") {
            std::cout << "Gradient direction" << std::endl;
            return make_unique<DirectionGradient<ORDER, mydim, ndim, FunctionalProblem_time<ORDER, mydim, ndim>>>(fp);
        }
        else if (d=="BFGS") {
            std::cout << "BFGS direction" << std::endl;
            return make_unique<DirectionBFGS<ORDER, mydim, ndim, FunctionalProblem_time<ORDER, mydim, ndim>>>(fp, dp.getNumNodes()*dp.getSplineNumber());
        }
        else{
            std::cout << "Unknown direction option - using gradient direction" << std::endl;
            //Rprintf("Unknown direction option - using gradient direction");
            return make_unique<DirectionGradient<ORDER, mydim, ndim, FunctionalProblem_time<ORDER, mydim, ndim>>>(fp);
        }
    }

};


#endif //DEV_FDAPDE_DESCENT_DIRECTION_FACTORY_H