superLuMtSolver.h

Go to the documentation of this file.

/************************************************************************
**
** Copyright (C) 2014 by Carlos Augusto Teixera Mendes
** All rights reserved.
**
** This file is part of the "GeMA" software. It's use should respect
** the terms in the license agreement that can be found together 
** with this source code.
** It is provided AS IS, with NO WARRANTY OF ANY KIND, 
** INCLUDING THE WARRANTY OF DESIGN, MERCHANTABILITY AND FITNESS FOR 
** A PARTICULAR PURPOSE.
**
************************************************************************/
 
#ifndef _GEMA_SUPERLUMT_SOLVER_H_
#define _GEMA_SUPERLUMT_SOLVER_H_
 
#include <gmNumSolver.h>
#include <gmSparseMatrix.h>
#include <QObject>
 
#include <slu_mt_ddefs.h> // For SuperMatrix and needed additional definitions used by supermatrix.h
 
class GmRunProgress;
 
extern "C" {
  int_t sp_ienv(int_t ispec);
}
 
class SuperLuMtSolver : public QObject, public GmNumSolver
{
public:
 
  SuperLuMtSolver(GmSimulationData* simulation, QString id, QString description, const GmLogCategory& logger);
  virtual ~SuperLuMtSolver();
 
  // Comments on the base class
  virtual const char* pluginName() const { return "SuperLuMtSolver"; }
 
  // Comments on the base class
  virtual const char* pluginType() const { return "ns"; }
 
  virtual bool loadPrivateData(LuaTable& table);
 
  virtual GmSolverMatrix* newMatrix(int nlin, int ncol, const GmSolverMatrix* sharedLayout,
                                    int numThreads = -1) const;
 
  // Comments on the base class
  virtual bool setSolverTolerance(double tol) { Q_UNUSED(tol); return false; }
 
  virtual bool solveImpl(GmSolverMatrix* A, const GmVector& b, GmVector& x, bool xIsGuess = false) const;
 
  virtual void printParameters(const GmLogCategory& logger);
 
  virtual QString quickDescription() const;
 
  // Comments on the base class
  virtual int trackItemId() const { return _trackSolveItem; }
 
private:
  SuperMatrix* superLuMatrix       (GmSolverMatrix* A) const;
  void         releaseSuperLuMatrix(SuperMatrix* sA) const;
 
  GmSparseMatrixOptions  _matrixOptions;      
  superlumt_options_t    _sluOptions;         
  int                    _sp_ienvOptions[9];  
  bool                   _singleThreadedBlas; 
 
  static int _numSolvers;      
  int _trackSolveItem = -1;
};
 
#endif