IpExpandedMultiVectorMatrix.hpp

Go to the documentation of this file.

// Copyright (C) 2009 International Business Machines and others.
// All Rights Reserved.
// This code is published under the Eclipse Public License.
//
// Authors:  Andreas Waechter             IBM    2009-11-05
 
#ifndef __IPEXPANDEDMULTIVECTORMATRIX_HPP__
#define __IPEXPANDEDMULTIVECTORMATRIX_HPP__
 
#include "IpUtils.hpp"
#include "IpMatrix.hpp"
#include "IpExpansionMatrix.hpp"
 
namespace Ipopt
{
 
class ExpandedMultiVectorMatrixSpace;
 
class ExpandedMultiVectorMatrix: public Matrix
{
public:
 
   ExpandedMultiVectorMatrix(
      const ExpandedMultiVectorMatrixSpace* owner_space
   );
 
   virtual ~ExpandedMultiVectorMatrix()
   { }
 
   SmartPtr<ExpandedMultiVectorMatrix> MakeNewExpandedMultiVectorMatrix() const;
 
   void SetVector(
      Index                  i,
      SmartPtr<const Vector> vec
   );
 
   inline SmartPtr<const Vector> GetVector(
      Index i
   ) const
   {
      DBG_ASSERT(i < NRows());
      return vecs_[i];
   }
 
   SmartPtr<const VectorSpace> RowVectorSpace() const;
 
   SmartPtr<const ExpandedMultiVectorMatrixSpace> ExpandedMultiVectorMatrixOwnerSpace() const;
 
   SmartPtr<const ExpansionMatrix> GetExpansionMatrix() const;
 
protected:
   virtual void MultVectorImpl(
      Number        alpha,
      const Vector& x,
      Number        beta,
      Vector&       y
   ) const;
 
   virtual void TransMultVectorImpl(
      Number        alpha,
      const Vector& x,
      Number        beta,
      Vector&       y
   ) const;
 
   virtual bool HasValidNumbersImpl() const;
 
   virtual void ComputeRowAMaxImpl(
      Vector& rows_norms,
      bool    init
   ) const;
 
   virtual void ComputeColAMaxImpl(
      Vector& cols_norms,
      bool    init
   ) const;
 
   virtual void PrintImpl(
      const Journalist&  jnlst,
      EJournalLevel      level,
      EJournalCategory   category,
      const std::string& name,
      Index              indent,
      const std::string& prefix
   ) const;
 
private:
 
   ExpandedMultiVectorMatrix();
 
   ExpandedMultiVectorMatrix(
      const ExpandedMultiVectorMatrix&
   );
 
   void operator=(
      const ExpandedMultiVectorMatrix&
   );
 
   const ExpandedMultiVectorMatrixSpace* owner_space_;
 
   std::vector<SmartPtr<const Vector> > vecs_;
 
};
 
class ExpandedMultiVectorMatrixSpace: public MatrixSpace
{
public:
 
   ExpandedMultiVectorMatrixSpace(
      Index                           nrows,
      const VectorSpace&              vec_space,
      SmartPtr<const ExpansionMatrix> exp_matrix
   );
 
   virtual ~ExpandedMultiVectorMatrixSpace()
   { }
 
   ExpandedMultiVectorMatrix* MakeNewExpandedMultiVectorMatrix() const
   {
      return new ExpandedMultiVectorMatrix(this);
   }
 
   virtual Matrix* MakeNew() const
   {
      return MakeNewExpandedMultiVectorMatrix();
   }
 
   SmartPtr<const VectorSpace> RowVectorSpace() const
   {
      return vec_space_;
   }
 
   SmartPtr<const ExpansionMatrix> GetExpansionMatrix() const
   {
      return exp_matrix_;
   }
 
private:
   SmartPtr<const VectorSpace> vec_space_;
 
   SmartPtr<const ExpansionMatrix> exp_matrix_;
};
 
inline SmartPtr<ExpandedMultiVectorMatrix> ExpandedMultiVectorMatrix::MakeNewExpandedMultiVectorMatrix() const
{
   return owner_space_->MakeNewExpandedMultiVectorMatrix();
}
 
inline SmartPtr<const VectorSpace> ExpandedMultiVectorMatrix::RowVectorSpace() const
{
   return owner_space_->RowVectorSpace();
}
 
inline SmartPtr<const ExpansionMatrix> ExpandedMultiVectorMatrix::GetExpansionMatrix() const
{
   return owner_space_->GetExpansionMatrix();
}
 
inline SmartPtr<const ExpandedMultiVectorMatrixSpace> ExpandedMultiVectorMatrix::ExpandedMultiVectorMatrixOwnerSpace() const
{
   return owner_space_;
}
 
} // namespace Ipopt
#endif