IpTSymLinearSolver.hpp

Go to the documentation of this file.

// Copyright (C) 2004, 2006 International Business Machines and others.
// All Rights Reserved.
// This code is published under the Common Public License.
//
// $Id: IpTSymLinearSolver.hpp 759 2006-07-07 03:07:08Z andreasw $
//
// Authors:  Carl Laird, Andreas Waechter     IBM    2004-03-17

#ifndef __IPTSYMLINEARSOLVER_HPP__
#define __IPTSYMLINEATSOLVER_HPP__

#include "IpSymLinearSolver.hpp"
#include "IpSparseSymLinearSolverInterface.hpp"
#include "IpTSymScalingMethod.hpp"
#include "IpSymMatrix.hpp"
#include "IpTripletToCSRConverter.hpp"
#include "IpTripletToDenseConverter.hpp"
#include <vector>

namespace Ipopt
{

  class TSymLinearSolver: public SymLinearSolver
  {
  public:
    TSymLinearSolver(SmartPtr<SparseSymLinearSolverInterface> solver_interface,
                     SmartPtr<TSymScalingMethod> scaling_method);

    virtual ~TSymLinearSolver();

    bool InitializeImpl(const OptionsList& options,
                        const std::string& prefix);

    virtual ESymSolverStatus MultiSolve(const SymMatrix &A,
                                        std::vector<SmartPtr<const Vector> >& rhsV,
                                        std::vector<SmartPtr<Vector> >& solV,
                                        bool check_NegEVals,
                                        Index numberOfNegEVals);

    virtual Index NumberOfNegEVals() const;

    //* @name Options of Linear solver */

    virtual bool IncreaseQuality();

    virtual bool ProvidesInertia() const;

    static void RegisterOptions(SmartPtr<RegisteredOptions> roptions);

  private:
    TSymLinearSolver();

    TSymLinearSolver(const TSymLinearSolver&);

    void operator=(const TSymLinearSolver&);

    TaggedObject::Tag atag_;

    Index dim_;

    Index nonzeros_triplet_;
    Index nonzeros_compressed_;

    bool have_structure_;
    bool linear_scaling_on_demand_;
    bool initialized_;

    SmartPtr<SparseSymLinearSolverInterface> solver_interface_;
    SmartPtr<TSymScalingMethod> scaling_method_;
    double* scaling_factors_;
    bool use_scaling_;
    bool just_switched_on_scaling_;

    Index* airn_;
    Index* ajcn_;
    SmartPtr<TripletToCSRConverter> triplet_to_csr_converter_;
    SmartPtr<TripletToDenseConverter> triplet_to_dense_converter_;
    SparseSymLinearSolverInterface::EMatrixFormat matrix_format_;

    bool warm_start_same_structure_;

    ESymSolverStatus InitializeStructure(const SymMatrix& symT_A);

    void GiveMatrixToSolver(bool new_matrix, const SymMatrix& sym_A);
  };

} // namespace Ipopt
#endif