IpMa77SolverInterface.hpp

Go to the documentation of this file.

// 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.
//
// $Id: IpMa77SolverInterface.hpp 2085 2012-02-16 11:01:51Z stefan $
//
// Authors: Jonathan Hogg                           2009-07-29 
//          Carl Laird, Andreas Waechter     IBM    2004-03-17

#ifndef __IPMA77SOLVERINTERFACE_HPP__
#define __IPMA77SOLVERINTERFACE_HPP__

#include "IpSparseSymLinearSolverInterface.hpp"

extern "C" {
#include "hsl_ma77d.h"
}

namespace Ipopt
{

  class Ma77SolverInterface: public SparseSymLinearSolverInterface
  {
  private:
    int ndim_; // Number of dimensions
    double *val_; // Storage for variables
    int numneg_; // Number of negative pivots in last factorization
    void *keep_; // Stores pointer to factors (only understood by Fortran code!)
    bool pivtol_changed_; // indicates if pivtol has been changed

    /* Options */
    struct ma77_control control_;
    double umax_;

  public:

    Ma77SolverInterface() :
        val_(NULL), keep_(NULL), pivtol_changed_(false)
    {}
    ~Ma77SolverInterface();

    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_Full_Format_0_Offset; }

    bool ProvidesDegeneracyDetection() const { return false; }
    ESymSolverStatus DetermineDependentRows(const Index* ia,
                                            const Index* ja,
                                            std::list<Index>& c_deps)
    {
      return SYMSOLVER_FATAL_ERROR;
    }

    static void MetisOrder(const int dim, const Index *ptr, const Index *row, Index *perm);
  };

} // namespace Ipopt

#endif