IpSumSymMatrix.hpp

Go to the documentation of this file.

// Copyright (C) 2004, 2008 International Business Machines and others.
// All Rights Reserved.
// This code is published under the Eclipse Public License.
//
// $Id: IpSumSymMatrix.hpp 2269 2013-05-05 11:32:40Z stefan $
//
// Authors:  Carl Laird, Andreas Waechter     IBM    2004-08-13

#ifndef __IPSUMSYMMATRIX_HPP__
#define __IPSUMSYMMATRIX_HPP__

#include "IpUtils.hpp"
#include "IpSymMatrix.hpp"

namespace Ipopt
{

  /* forward declarations */
  class SumSymMatrixSpace;

  class SumSymMatrix : public SymMatrix
  {
  public:


    SumSymMatrix(const SumSymMatrixSpace* owner_space);

    ~SumSymMatrix();

    void SetTerm(Index iterm, Number factor, const SymMatrix& matrix);

    void GetTerm(Index iterm, Number& factor, SmartPtr<const SymMatrix>& matrix) const;

    Index NTerms() const;

  protected:
    virtual void MultVectorImpl(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:
    SumSymMatrix();

    SumSymMatrix(const SumSymMatrix&);

    void operator=(const SumSymMatrix&);

    std::vector<Number> factors_;

    std::vector<SmartPtr<const SymMatrix> > matrices_;

    const SumSymMatrixSpace* owner_space_;
  };

  class SumSymMatrixSpace : public SymMatrixSpace
  {
  public:
    SumSymMatrixSpace(Index ndim, Index nterms)
        :
        SymMatrixSpace(ndim),
        nterms_(nterms)
    {}

    ~SumSymMatrixSpace()
    {}

    Index NTerms() const
    {
      return nterms_;
    }

    void SetTermSpace(Index term_idx, const SymMatrixSpace& space);

    SmartPtr<const SymMatrixSpace> GetTermSpace(Index term_idx) const;

    SumSymMatrix* MakeNewSumSymMatrix() const;

    virtual SymMatrix* MakeNewSymMatrix() const;

  private:
    Index nterms_;

    std::vector< SmartPtr<const SymMatrixSpace> > term_spaces_;
  };

} // namespace Ipopt
#endif