IpMa86SolverInterface.hpp

Go to the documentation of this file.

// Copyright (C) 2011, Science and Technology Facilities Council
// Copyright (C) 2009, Jonathan Hogg <jdh41.at.cantab.net>
// Copyright (C) 2004, 2007 International Business Machines and others.
// All Rights Reserved.
// This code is published under the Eclipse Public License.
//
// Authors: Jonathan Hogg                    STFC   2011-03-14
//          Jonathan Hogg                           2009-07-29
//          Carl Laird, Andreas Waechter     IBM    2004-03-17
 
#ifndef __IPMA86SOLVERINTERFACE_HPP__
#define __IPMA86SOLVERINTERFACE_HPP__
 
#include "IpSparseSymLinearSolverInterface.hpp"
extern "C"
{
#include "hsl_ma86d.h"
}
 
namespace Ipopt
{
 
class Ma86SolverInterface: public SparseSymLinearSolverInterface
{
private:
   enum order_opts
   {
      ORDER_AUTO,
      ORDER_AMD,
      ORDER_METIS
   };
 
   int ndim_;     
   double* val_;  
   int numneg_;   
   Index* order_; 
   void* keep_;   
   bool pivtol_changed_; 
 
   /* Options */
   struct ma86_control control_;
   double umax_;
   int ordering_;
 
public:
 
   Ma86SolverInterface()
      : val_(NULL),
        order_(NULL),
        keep_(NULL),
        pivtol_changed_(false)
   { }
 
   ~Ma86SolverInterface();
 
   static void RegisterOptions(
      SmartPtr<RegisteredOptions> roptions
   );
 
   bool InitializeImpl(
      const OptionsList& options,
      const std::string& prefix
   );
 
   ESymSolverStatus InitializeStructure(
      Index        dim,
      Index        nonzeros,
      const Index* ia,
      const Index* ja
   );
 
   double* GetValuesArrayPtr()
   {
      return val_;
   }
 
   ESymSolverStatus MultiSolve(
      bool         new_matrix,
      const Index* ia,
      const Index* ja,
      Index        nrhs,
      double*      rhs_vals,
      bool         check_NegEVals,
      Index        numberOfNegEVals
   );
 
   Index NumberOfNegEVals() const
   {
      return numneg_;
   }
 
   //* @name Options of Linear solver */
   bool IncreaseQuality();
 
   bool ProvidesInertia() const
   {
      return true;
   }
 
   EMatrixFormat MatrixFormat() const
   {
      return CSR_Format_1_Offset;
   }
 
   bool ProvidesDegeneracyDetection() const
   {
      return false;
   }
 
   ESymSolverStatus DetermineDependentRows(
      const Index*      /*ia*/,
      const Index*      /*ja*/,
      std::list<Index>& /*c_deps*/
   )
   {
      return SYMSOLVER_FATAL_ERROR;
   }
};
 
} // namespace Ipopt
 
#endif