IpWsmpSolverInterface.hpp

Go to the documentation of this file.

// Copyright (C) 2005, 2010 International Business Machines and others.
// All Rights Reserved.
// This code is published under the Eclipse Public License.
//
// $Id: IpWsmpSolverInterface.hpp 1861 2010-12-21 21:34:47Z andreasw $
//
// Authors:  Carl Laird, Andreas Waechter     IBM    2005-03-17


#ifndef __IPWSMPSOLVERINTERFACE_HPP__
#define __IPWSMPSOLVERINTERFACE_HPP__

#include "IpSparseSymLinearSolverInterface.hpp"

//#define PARDISO_MATCHING_PREPROCESS

namespace Ipopt
{

  class WsmpSolverInterface: public SparseSymLinearSolverInterface
  {
  public:
    WsmpSolverInterface();

    virtual ~WsmpSolverInterface();

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


    virtual ESymSolverStatus InitializeStructure(Index dim, Index nonzeros,
        const Index *ia,
        const Index *ja);

    virtual double* GetValuesArrayPtr();

    virtual ESymSolverStatus MultiSolve(bool new_matrix,
                                        const Index* ia,
                                        const Index* ja,
                                        Index nrhs,
                                        double* rhs_vals,
                                        bool check_NegEVals,
                                        Index numberOfNegEVals);

    virtual Index NumberOfNegEVals() const;

    //* @name Options of Linear solver */

    virtual bool IncreaseQuality();

    virtual bool ProvidesInertia() const
    {
      return true;
    }
    EMatrixFormat MatrixFormat() const
    {
      return CSR_Format_1_Offset;
    }

    static void RegisterOptions(SmartPtr<RegisteredOptions> roptions);

    virtual bool ProvidesDegeneracyDetection() const;

    virtual ESymSolverStatus DetermineDependentRows(const Index* ia,
        const Index* ja,
        std::list<Index>& c_deps);

  private:
    WsmpSolverInterface(const WsmpSolverInterface&);

    void operator=(const WsmpSolverInterface&);

    Index dim_;

    Index nonzeros_;

    double* a_;

#ifdef PARDISO_MATCHING_PREPROCESS

    ipfint* ia2;
    ipfint* ja2;
    double* a2_;
    ipfint* perm2;
    double* scale2;

#endif


    Index wsmp_num_threads_;
    Number wsmp_pivtol_;
    Number wsmp_pivtolmax_;
    Index wsmp_scaling_;
    Number wsmp_singularity_threshold_;
    Index wsmp_write_matrix_iteration_;
    bool skip_inertia_check_;
    bool wsmp_no_pivoting_;

    Index matrix_file_number_;

    Index negevals_;

    bool initialized_;
    bool printed_num_threads_;
    bool pivtol_changed_;
    bool have_symbolic_factorization_;
    Index factorizations_since_recomputed_ordering_;

    ipfint* IPARM_;
    double* DPARM_;
    ipfint* PERM_;
    ipfint* INVP_;
    ipfint* MRP_;

    ESymSolverStatus SymbolicFactorization(const Index* ia, const Index* ja);

    ESymSolverStatus InternalSymFact(const Index* ia, const Index* ja,
                                     Index numberOfNegEVals);

    ESymSolverStatus Factorization(const Index* ia,
                                   const Index* ja,
                                   bool check_NegEVals,
                                   Index numberOfNegEVals);

    ESymSolverStatus Solve(const Index* ia,
                           const Index* ja,
                           Index nrhs,
                           double *rhs_vals);
  };

} // namespace Ipopt
#endif