MumpsSolver.hpp

Go to the documentation of this file.

//==============================================================================
//
//                       Seismo Virtual Laboratory
//             Module for Serial and Parallel Analysis of seismic 
//         wave propagation and soil-structure interaction simulation
//         Copyright (C) 2018-2021, The California Institute of Technology
//                         All Rights Reserved.
//
// Commercial use of this program without express permission of the California
// Institute of Technology, is strictly  prohibited. See  file "COPYRIGHT"  in
// main  directory  for  information on  usage  and  redistribution, and for a
// DISCLAIMER OF ALL WARRANTIES.
//
//==============================================================================
//
// Written by:
//   Danilo S. Kusanovic (dkusanov@caltech.edu)
//   Elnaz E. Seylabi    (elnaze@unr.edu)
//
// Supervised by:
//   Domniki M. Asimaki  (domniki@caltech.edu)
//
// References : 
//   [1] A Fully Asynchronous Multifrontal Solver Using Distributed Dynamic 
//       Scheduling, P. R. Amestoy and I. S. Duff and J. Koster and J.-Y. 
//       L'Excellent, SIAM Journal on Matrix Analysis and Applications, 
//       volume 23 (1) 15-41, 2003. 
//   [2] Hybrid scheduling for the parallel solution of linear systems, P. R. 
//       Amestoy and A. Guermouche and J.-Y. L'Excellent and S. Pralet, Journal 
//       of Parallel Computing, volume 32(2), 136-156, 2006.
//
// Description:
// This file contains the "MumpsSolver" solver declaration; this object 
// performs a (MU)ltifrontal (M)assively (P)arallel (S)parse direct Solver on 
// the matrix A, assuming that such matrix is symmetric positive definite (SPD), 
// general symmetric (SYM) or unsymmetric (USYM) matrix
//------------------------------------------------------------------------------

#ifndef _MUMPSSOLVER_HPP_
#define _MUMPSSOLVER_HPP_

#define JOB_INIT -1

#define JOB_END -2

#define USE_COMM_WORLD -987654

#define ICNTL(I) icntl[(I)-1]

#include <dmumps_c.h>
#include <Eigen/Dense>
#include <Eigen/SparseCore>

#include "LinearSystem.hpp"


class MumpsSolver : public LinearSystem {

    public:
        MumpsSolver(unsigned int option=1, bool flag = false);

        ~MumpsSolver();

        bool SolveSystem(Eigen::SparseMatrix<double> &A, Eigen::VectorXd &b);

        const Eigen::VectorXd& GetSolution();

    protected:
        unsigned int nz;

        unsigned int Option;

        bool Flag;

        bool Factored;

        bool Initialized;

        Eigen::VectorXd x;

        DMUMPS_STRUC_C id;
};

#endif