SensPCalculator.hpp

Go to the documentation of this file.

// Copyright 2009, 2011 Hans Pirnay
// All Rights Reserved.
// This code is published under the Eclipse Public License.
//
// Date   : 2009-05-06
 
#ifndef __ASPCALCULATOR_HPP__
#define __ASPCALCULATOR_HPP__
 
#include "IpAlgStrategy.hpp"
#include "SensSimpleBacksolver.hpp"
#include "SensSchurData.hpp"
 
namespace Ipopt
{
 
class SIPOPTLIB_EXPORT PCalculator: public AlgorithmStrategyObject
{
public:
   PCalculator(
      SmartPtr<SensBacksolver> backsolver,
      SmartPtr<SchurData>      A_data
   )
      : backsolver_(backsolver),
        data_A_init(ConstPtr(A_data->MakeNewSchurDataCopy())),
        data_A_(A_data)
   { }
 
   virtual ~PCalculator()
   { }
 
   /* Overloaded from AlgorithmStrategyObject */
   virtual bool InitializeImpl(
      const OptionsList& /*options*/,
      const std::string& /*prefix*/
   )
   {
      return true;
   }
 
   virtual bool ComputeP() = 0;
 
   virtual bool GetSchurMatrix(
      const SmartPtr<const SchurData>& B,
      SmartPtr<Matrix>&                S
   ) = 0;
 
   virtual void PrintImpl(
      const Journalist&  jnlst,
      EJournalLevel      level,
      EJournalCategory   category,
      const std::string& name,
      Index              indent,
      const std::string& prefix
   ) const = 0;
 
   void Print(
      const Journalist&  jnlst,
      EJournalLevel      level,
      EJournalCategory   category,
      const std::string& name,
      Index              indent = 0,
      const std::string& prefix = ""
   ) const
   {
      if( jnlst.ProduceOutput(level, category) )
      {
         PrintImpl(jnlst, level, category, name, indent, prefix);
      }
   }
 
   void Print(
      SmartPtr<const Journalist> jnlst,
      EJournalLevel              level,
      EJournalCategory           category,
      const std::string&         name,
      Index                      indent,
      const std::string&         prefix
   ) const
   {
      if( IsValid(jnlst) && jnlst->ProduceOutput(level, category) )
      {
         PrintImpl(*jnlst, level, category, name, indent, prefix);
      }
   }
 
   SmartPtr<SensBacksolver> Solver() const
   {
      return backsolver_;
   }
 
   SmartPtr<const SchurData> data_A() const
   {
      return ConstPtr(data_A_);
   }
 
   SmartPtr<SchurData> data_A_nonconst() const
   {
      return data_A_;
   }
 
   void reset_data_A()
   {
      data_A_ = data_A_init->MakeNewSchurDataCopy();
   }
 
private:
 
   SmartPtr<SensBacksolver> backsolver_;
 
   SmartPtr<const SchurData> data_A_init;
   SmartPtr<SchurData> data_A_;
};
 
}
 
#endif