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