IpTSymLinearSolver.hpp

Go to the documentation of this file.

// Copyright (C) 2004, 2007 International Business Machines and others.
// All Rights Reserved.
// This code is published under the Eclipse Public License.
//
// $Id: IpTSymLinearSolver.hpp 2206 2013-04-13 15:29:15Z stefan $
//
// Authors:  Carl Laird, Andreas Waechter     IBM    2004-03-17

#ifndef __IPTSYMLINEARSOLVER_HPP__
#define __IPTSYMLINEARSOLVER_HPP__

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

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;

    bool ProvidesDegeneracyDetection() const;
    ESymSolverStatus DetermineDependentRows(Index n_rows, Index n_cols,
                                            Index n_jac_nz,
                                            Number* jac_c_vals,
                                            Index* jac_c_iRow,
                                            Index* jac_c_jCol,
                                            std::list<Index>& c_deps);

    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_;
    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