LuSolver.h

#ifndef PSCF_LU_SOLVER_H
#define PSCF_LU_SOLVER_H
 
/*
* PSCF - Polymer Self-Consistent Field Theory
*
* Copyright 2016 - 2022, The Regents of the University of Minnesota
* Distributed under the terms of the GNU General Public License.
*/
 
#include <util/containers/Array.h>
#include <util/containers/Matrix.h>
#include <gsl/gsl_matrix.h>
#include <gsl/gsl_vector.h>
#include <gsl/gsl_permutation.h>
 
namespace Pscf 
{
 
   using namespace Util;
 
   class LuSolver
   {
   public:
 
      LuSolver();
 
      ~LuSolver();
 
      void allocate(int n);
 
      void computeLU(const Matrix<double>& A);
 
      void solve(Array<double>& b, Array<double>& x);
 
      void inverse (Matrix<double>& inv); 
 
   private:
 
      gsl_vector b_;
 
      gsl_vector x_;
 
      gsl_matrix* luPtr_;
 
      gsl_matrix* gMatInverse_;
 
      gsl_permutation* permPtr_;
 
      int signum_;
 
      int n_;
 
   };
 
}
#endif