PetscSolver.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] Kris Buschelman, Lisandro Dalcin, Alp Dener, Victor Eijkhout, William D.
//       Gropp, Dmitry Karpeyev, Dinesh Kaushik, Matthew G. Knepley, Dave A.
//       May, Lois Curfman McInnes, Richard Tran Mills, Todd Munson, Karl Rupp,
//       Patrick Sanan, Barry F. Smith, Stefano Zampini, Hong Zhang, and Hong
//       Zhang. PETSc Web page. https://www.mcs.anl.gov/petsc, 2019.
//       URL https://www.mcs.anl.gov/petsc
//
//   [2] Satish Balay, Shrirang Abhyankar, Mark F. Adams, Jed Brown, Peter Brune,
//       Kris Buschelman, Lisandro Dalcin, Alp Dener, Victor Eijkhout, William D.
//       Gropp, Dmitry Karpeyev, Dinesh Kaushik, Matthew G. Knepley, Dave A.
//       May, Lois Curfman McInnes, Richard Tran Mills, Todd Munson, Karl Rupp,
//       Patrick Sanan, Barry F. Smith, Stefano Zampini, Hong Zhang, and Hong
//       Zhang. PETSc users manual. Technical Report ANL-95/11 - Revision 3.12,
//       Argonne National Laboratory, 2019. URL https://www.mcs.anl.gov/petsc
//
// Description:
//------------------------------------------------------------------------------

#ifndef _PETSCSOLVER_HPP_
#define _PETSCSOLVER_HPP_

#include <Eigen/Dense>
#include <Eigen/SparseCore>

#include "LinearSystem.hpp"

class PetscSolver : public LinearSystem {

    public:
        PetscSolver(unsigned int dnz, unsigned int onz, double tol=1E-12, unsigned int kspnum=0);

        ~PetscSolver();

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

        const Eigen::VectorXd& GetSolution();

    protected:
        unsigned int d_nz;

        unsigned int o_nz;

        double Tolerance;

        Eigen::VectorXd x;
};

#endif