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.
//
// Authors:  Carl Laird, Andreas Waechter     IBM    2005-03-17
 
#ifndef __IPWSMPSOLVERINTERFACE_HPP__
#define __IPWSMPSOLVERINTERFACE_HPP__
 
#include "IpSparseSymLinearSolverInterface.hpp"
#include "IpPardisoSolverInterface.hpp"  // for IPOPT_DECL_SMAT_REORDERING_PARDISO_WSMP
#include "IpTypes.h"
 
namespace Ipopt
{
 
class WsmpSolverInterface: public SparseSymLinearSolverInterface
{
public:
 
   WsmpSolverInterface(
#ifdef PARDISO_MATCHING_PREPROCESS
      SmartPtr<LibraryLoader> pardisoloader_  
#endif
   );
 
   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
   );
 
   static void GetVersion(
      int& V,
      int& R,
      int& M
   );
 
   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
   Index* ia2;
   Index* ja2;
   double* a2_;
   Index* 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_;
 
 
   Index* IPARM_;
   double* DPARM_;
   Index* PERM_;
   Index* INVP_;
   Index* MRP_;
 
#ifdef PARDISO_MATCHING_PREPROCESS
   SmartPtr<LibraryLoader> pardisoloader;
   IPOPT_DECL_SMAT_REORDERING_PARDISO_WSMP(*smat_reordering_pardiso_wsmp);
#endif
 
   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