SensSchurDriver.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 __ASSCHURDRIVER_HPP__
#define __ASSCHURDRIVER_HPP__
 
#include "SensSchurData.hpp"
#include "SensPCalculator.hpp"
#include "IpVector.hpp"
#include "IpIteratesVector.hpp"
 
namespace Ipopt
{
 
class SIPOPTLIB_EXPORT SchurDriver: public AlgorithmStrategyObject
{
public:
 
   SchurDriver(
      SmartPtr<PCalculator> pcalc,
      SmartPtr<SchurData>   data_B
   )
      : pcalc_(pcalc),
        data_B_(data_B)
   { }
 
   virtual ~SchurDriver()
   { }
 
   virtual bool InitializeImpl(
      const OptionsList& /*options*/,
      const std::string& /*prefix*/
   )
   {
      return true;
   }
 
   virtual SmartPtr<const SchurData> data_A() const
   {
      return pcalc_->data_A();
   }
 
   virtual SmartPtr<const SchurData> data_B() const
   {
      return ConstPtr(data_B_);
   }
 
   virtual SmartPtr<SchurData> data_A_nonconst()
   {
      return pcalc_->data_A_nonconst();
   }
 
   virtual SmartPtr<SchurData> data_B_nonconst()
   {
      return data_B_;
   }
 
   virtual SmartPtr<const PCalculator> pcalc() const
   {
      return ConstPtr(pcalc_);
   }
 
   virtual SmartPtr<PCalculator> pcalc_nonconst()
   {
      return pcalc_;
   }
 
   /* Sets the Data for which this SchurMatrix will be built. */
 
   virtual bool SchurBuild() = 0;
 
   virtual bool SchurFactorize() = 0;
 
   virtual bool SchurSolve(
      SmartPtr<IteratesVector>       lhs,
      SmartPtr<const IteratesVector> rhs,
      SmartPtr<Vector>               delta_u,
      SmartPtr<IteratesVector>       sol = NULL
   ) = 0;
 
   /* Performs a backsolve on S and K; calls the latter with sol=K^(-1)*r_s=0
    virtual bool SchurSolve(SmartPtr<IteratesVector> lhs,
    SmartPtr<const IteratesVector> rhs,
    SmartPtr<Vector> delta_u) =0;
    */
private:
   SchurDriver()
   { }
 
   SmartPtr<PCalculator> pcalc_;
 
   SmartPtr<SchurData> data_B_;
};
 
}
 
#endif