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