IpLapackSolverInterface.hpp

Go to the documentation of this file.

#ifndef __IPLAPACKSOLVERINTERFACE_HPP__

#define __IPLAPACKSSOLVERINTERFACE_HPP__

#include "IpSparseSymLinearSolverInterface.hpp"

namespace Ipopt
{

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

    virtual ~LapackSolverInterface();

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

    static void RegisterOptions(SmartPtr<RegisteredOptions> roptions);

  private:
    LapackSolverInterface(const LapackSolverInterface&);

    void operator=(const LapackSolverInterface&);

    Index n;

    Index nz;

    double* a;

    double* afact;

    int* irn_;
    int* jcn_;

    Index negevals_;

    bool isFactored;








    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);
    //MUMPS data structure
//    DMUMPS_STRUC_C mumps_data;

  };

} // namespace Ipopt
#endif