IpInexactCq.hpp

Go to the documentation of this file.

// Copyright (C) 2008 International Business Machines and others.
// All Rights Reserved.
// This code is published under the Eclipse Public License.
//
// $Id: IpInexactCq.hpp 1861 2010-12-21 21:34:47Z andreasw $
//
// Authors:  Andreas Waechter           IBM    2008-08-31
//               derived from IpIpoptCalculatedQuantities.hpp

#ifndef __IPINEXACTCQ_HPP__
#define __IPINEXACTCQ_HPP__

#include "IpIpoptCalculatedQuantities.hpp"
#include "IpInexactData.hpp"

namespace Ipopt
{

  class InexactCq : public IpoptAdditionalCq
  {
  public:

    InexactCq(IpoptNLP* ip_nlp,
              IpoptData* ip_data,
              IpoptCalculatedQuantities* ip_cq);

    virtual ~InexactCq();

    bool Initialize(const Journalist& jnlst,
                    const OptionsList& options,
                    const std::string& prefix);

    static void RegisterOptions(const SmartPtr<RegisteredOptions>& roptions);

    SmartPtr<const Vector> curr_jac_cdT_times_curr_cdminuss();

    SmartPtr<const Vector> curr_scaling_slacks();

    SmartPtr<const Vector> curr_slack_scaled_d_minus_s();

    Number curr_scaled_Ac_norm();

    Number curr_scaled_A_norm2();

    Number slack_scaled_norm(const Vector& x, const Vector &s);

    SmartPtr<const Vector> curr_W_times_vec_x(const Vector& vec_x);

    SmartPtr<const Vector> curr_W_times_vec_s(const Vector& vec_s);

    SmartPtr<const Vector> curr_Wu_x();

    SmartPtr<const Vector> curr_Wu_s();

    Number curr_uWu();

    SmartPtr<const Vector> curr_jac_times_normal_c();

    SmartPtr<const Vector> curr_jac_times_normal_d();

  private:
    InexactCq();

    InexactCq(const InexactCq&);

    void operator=(const InexactCq&);

    IpoptNLP* ip_nlp_;
    IpoptData* ip_data_;
    IpoptCalculatedQuantities* ip_cq_;

    InexactData& InexData()
    {
      InexactData& inexact_data =
        static_cast<InexactData&>(ip_data_->AdditionalData());
      DBG_ASSERT(dynamic_cast<InexactData*>(&ip_data_->AdditionalData()));
      return inexact_data;
    }

    CachedResults<SmartPtr<const Vector> > curr_jac_cdT_times_curr_cdminuss_cache_;
    CachedResults<SmartPtr<const Vector> > curr_scaling_slacks_cache_;
    CachedResults<SmartPtr<const Vector> > curr_slack_scaled_d_minus_s_cache_;
    CachedResults<Number> curr_scaled_Ac_norm_cache_;
    CachedResults<Number> slack_scaled_norm_cache_;
    CachedResults<SmartPtr<const Vector> > curr_W_times_vec_x_cache_;
    CachedResults<SmartPtr<const Vector> > curr_W_times_vec_s_cache_;
    CachedResults<SmartPtr<const Vector> > curr_Wu_x_cache_;
    CachedResults<SmartPtr<const Vector> > curr_Wu_s_cache_;
    CachedResults<Number> curr_uWu_cache_;
    CachedResults<SmartPtr<const Vector> > curr_jac_times_normal_c_cache_;
    CachedResults<SmartPtr<const Vector> > curr_jac_times_normal_d_cache_;

    Number slack_scale_max_;
  };

} // namespace Ipopt

#endif