gsOptimizer.h

Go to the documentation of this file.

 
#include <gsCore/gsLinearAlgebra.h>
#include <gsOptimizer/gsOptProblem.h>
#include <gsIO/gsOptionList.h>
 
#pragma once
 
namespace gismo
{
 
template <typename T>
class gsOptimizer
{
public:
 
    gsOptimizer()
    :
    m_op(nullptr)
    {};
 
    gsOptimizer(gsOptProblem<T> * problem)
    :
    m_op(problem)
    {
 
    }
 
    virtual ~gsOptimizer() { };
 
public:
 
    virtual void defaultOptions()
    {
        m_options.addInt("MaxIterations","Maximum iterations",100);
        m_options.addInt("Verbose","Verbosity level",0);
    }
 
    virtual void getOptions()
    {
        m_maxIterations = m_options.getInt("MaxIterations");
        m_verbose = m_options.getInt("Verbose");
    }
 
    virtual bool intermediateCallback() { return true;}
 
public:
 
    const gsMatrix<T> & currentDesign() const { return m_curDesign; }
    gsMatrix<T> & currentDesign() { return m_curDesign; }
 
    T currentObjValue() const
    {
        gsAsConstVector<T> tmp(m_curDesign.data(), m_op->numDesignVars());
        return m_op->evalObj(tmp);
    }
 
    T objective()    const { return m_finalObjective; }
 
    int iterations() const { return m_numIterations; }
 
    gsOptionList & options() { return m_options; }
 
public:
 
    virtual void solve (const gsMatrix<T> & initialGuess) = 0;
    virtual void solve ()
    {
        m_curDesign.resize(m_op->numDesignVars(),1);
        m_curDesign.setZero();
        this->solve(m_curDesign);
    }
 
    std::ostream &print(std::ostream &os) const
    {
        os << "gsOptimizer";
        return os;
    }
 
 
protected:
 
    gsOptProblem<T> * m_op;
 
    gsMatrix<T>     m_curDesign;
 
    gsOptionList m_options;
    index_t m_maxIterations, m_verbose;
 
protected:
 
    // Statistics
    index_t     m_numIterations;
    T           m_finalObjective;
 
private:
 
    gsOptimizer(const gsOptimizer & );
    gsOptimizer& operator=(const gsOptimizer & );
};
 
} // end namespace gismo
 
// note: must be statically compiled in header-only mode
//#ifndef GISMO_BUILD_LIB
//#include GISMO_HPP_HEADER(gsOptimizer.hpp)
//#endif