IpTransposeMatrix.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.
//
// Authors:  Andreas Waechter              IBM    2008-08-25
 
#ifndef __IPTRANSPOSEMATRIX_HPP__
#define __IPTRANSPOSEMATRIX_HPP__
 
#include "IpMatrix.hpp"
 
namespace Ipopt
{
 
/* forward declarations */
class TransposeMatrixSpace;
 
class TransposeMatrix: public Matrix
{
public:
 
   TransposeMatrix(
      const TransposeMatrixSpace* owner_space
   );
 
   ~TransposeMatrix()
   { }
 
   SmartPtr<const Matrix> OrigMatrix() const
   {
      return ConstPtr(orig_matrix_);
   }
 
protected:
   virtual void MultVectorImpl(
      Number        alpha,
      const Vector& x,
      Number        beta,
      Vector&       y
   ) const
   {
      DBG_ASSERT(IsValid(orig_matrix_));
      orig_matrix_->TransMultVector(alpha, x, beta, y);
   }
 
   virtual void TransMultVectorImpl(
      Number        alpha,
      const Vector& x,
      Number        beta,
      Vector&       y
   ) const
   {
      DBG_ASSERT(IsValid(orig_matrix_));
      orig_matrix_->MultVector(alpha, x, beta, y);
   }
 
   virtual bool HasValidNumbersImpl() const
   {
      DBG_ASSERT(IsValid(orig_matrix_));
      return orig_matrix_->HasValidNumbers();
   }
 
   virtual void ComputeRowAMaxImpl(
      Vector& rows_norms,
      bool    init
   ) const
   {
      DBG_ASSERT(IsValid(orig_matrix_));
      orig_matrix_->ComputeColAMax(rows_norms, init);
   }
 
   virtual void ComputeColAMaxImpl(
      Vector& rows_norms,
      bool    init
   ) const
   {
      DBG_ASSERT(IsValid(orig_matrix_));
      orig_matrix_->ComputeRowAMax(rows_norms, init);
   }
 
   virtual void PrintImpl(
      const Journalist&  jnlst,
      EJournalLevel      level,
      EJournalCategory   category,
      const std::string& name,
      Index              indent,
      const std::string& prefix
   ) const;
 
private:
 
   TransposeMatrix();
 
   TransposeMatrix(
      const TransposeMatrix&
   );
 
   void operator=(
      const TransposeMatrix&
   );
 
   SmartPtr<Matrix> orig_matrix_;
};
 
class TransposeMatrixSpace: public MatrixSpace
{
public:
 
   TransposeMatrixSpace(
      const MatrixSpace* orig_matrix_space
   )
      : MatrixSpace(orig_matrix_space->NCols(), orig_matrix_space->NRows()),
        orig_matrix_space_(orig_matrix_space)
   { }
 
   virtual ~TransposeMatrixSpace()
   { }
 
   virtual Matrix* MakeNew() const
   {
      return MakeNewTransposeMatrix();
   }
 
   TransposeMatrix* MakeNewTransposeMatrix() const
   {
      return new TransposeMatrix(this);
   }
 
   Matrix* MakeNewOrigMatrix() const
   {
      return orig_matrix_space_->MakeNew();
   }
 
private:
 
   TransposeMatrixSpace();
 
   TransposeMatrixSpace(
      const TransposeMatrixSpace&
   );
 
   void operator=(
      const TransposeMatrixSpace&
   );
 
   SmartPtr<const MatrixSpace> orig_matrix_space_;
};
 
} // namespace Ipopt
#endif