gsIpOpt.h

Go to the documentation of this file.

 
#include <gsCore/gsLinearAlgebra.h>
#include <gsOptimizer/gsOptProblem.h>
#include <gsOptimizer/gsOptimizer.h>
 
#pragma once
 
namespace gismo
{
 
template <typename T> class gsIpOptTNLP;
template <typename T> class gsIpOptPrivate;
template <typename T>
class gsIpOpt : public gsOptimizer<T>
{
public:
    using Base = gsOptimizer<T>;
 
    friend class gsIpOptTNLP<T>;
 
public:
 
    gsIpOpt(gsOptProblem<T> * problem);
 
    virtual ~gsIpOpt();
 
    // // Currently not used
    // void defaultOptions()
    // {
    //     using Base::defaultOptions;
    // }
 
    void getOptions()
    {
        gsWarn<<"gsIpOpt.h has its options stored in filedata/options/ipopt.opt\n";
    }
 
public:
 
    const gsMatrix<T> & lambda() const { return m_lambda; }
 
public:
 
    void solve (const gsMatrix<T> & initialGuess);
 
    std::ostream &print(std::ostream &os) const
    {
        os << "Design variables:" << m_op->numDesignVars()
           << "\nNumber of constraints: " << m_op->numConstraints()
            //<< "\ndesign lower:" << m_op->desLowerBounds().transpose()
            //<< "\ndesign upper:" << m_op->desUpperBounds().transpose()
            //<< "\nconstr lower:" << m_op->conLowerBounds().transpose()
            //<< "\nconstr upper:" << m_op->conUpperBounds().transpose()
           << "\nNumber of active sensitivities: " << m_op->numConJacNonZero()
           << "\nSparsity of sensitivities: " <<std::fixed<<std::setprecision(2)
           << (100.0*m_op->numConJacNonZero())/(m_op->numDesignVars()*m_op->numConstraints()) << " %"
            //<< "\nm_op->conJacRows():" << m_op->conJacRows().transpose()
            //<< "\nm_op->conJacCols():" << m_op->conJacCols().transpose()
            //<< "\nm_op->curDesign():" << m_op->curDesign().transpose()
            ;
        return os;
    }
 
 
protected:
 
    gsMatrix<T> m_lambda;
 
protected:
 
// Members taken from Base
protected:
    using Base::m_op;
    using Base::m_numIterations;
    using Base::m_finalObjective;
    using Base::m_curDesign;
    using Base::m_options;
    using Base::m_verbose;
    using Base::m_maxIterations;
 
    using Base::defaultOptions;
    using Base::getOptions;
 
    // Statistics
 
    gsIpOptPrivate<T> * m_data;
 
private:
 
    gsIpOpt(const gsIpOpt & );
    gsIpOpt& operator=(const gsIpOpt & );
};
 
 
} // end namespace gismo
 
// note: statically compiled in header-only mode
#ifndef GISMO_BUILD_LIB
#include GISMO_HPP_HEADER(gsIpOpt.hpp)
#endif