talbring/SU2/linear__solvers__structure_8hpp_source.html

 #pragma once

 #include "./mpi_structure.hpp"

 #include <climits>
 #include <limits>
 #include <cmath>
 #include <vector>
 #include <iostream>
 #include <cstdlib>
 #include <iomanip>
 #include <string>

 #include "option_structure.hpp"
 #include "vector_structure.hpp"
 #include "matrix_structure.hpp"
 #include "config_structure.hpp"
 #include "geometry_structure.hpp"

 using namespace std;

 class CSysSolve {

 private:

   su2double Sign(const su2double & x, const su2double & y) const;

   void ApplyGivens(const su2double & s, const su2double & c, su2double & h1, su2double & h2);

   void GenerateGivens(su2double & dx, su2double & dy, su2double & s, su2double & c);

   void SolveReduced(const int & n, const vector<vector<su2double> > & Hsbg,
                     const vector<su2double> & rhs, vector<su2double> & x);

   void ModGramSchmidt(int i, vector<vector<su2double> > & Hsbg, vector<CSysVector> & w);

   void WriteHeader(const string & solver, const su2double & restol, const su2double & resinit);

   void WriteHistory(const int & iter, const su2double & res, const su2double & resinit);

 public:

   unsigned long CG_LinSolver(const CSysVector & b, CSysVector & x, CMatrixVectorProduct & mat_vec,
                                   CPreconditioner & precond, su2double tol,
                                   unsigned long m, bool monitoring);

   unsigned long FGMRES_LinSolver(const CSysVector & b, CSysVector & x, CMatrixVectorProduct & mat_vec,
                       CPreconditioner & precond, su2double tol,
                       unsigned long m, su2double *residual, bool monitoring);

   unsigned long BCGSTAB_LinSolver(const CSysVector & b, CSysVector & x, CMatrixVectorProduct & mat_vec,
                         CPreconditioner & precond, su2double tol,
                         unsigned long m, su2double *residual, bool monitoring);

   unsigned long Solve(CSysMatrix & Jacobian, CSysVector & LinSysRes, CSysVector & LinSysSol, CGeometry *geometry, CConfig *config);


   void SetExternalSolve(CSysMatrix & Jacobian, CSysVector & LinSysRes, CSysVector & LinSysSol, CGeometry *geometry, CConfig *config);

 };

 #include "linear_solvers_structure.inl"
geometry_structure.hpp
Headers of the main subroutines for creating the geometrical structure. The subroutines and functions...

matrix_structure.hpp
Headers of the main subroutines for creating the sparse matrices-by-blocks. The subroutines and funct...

config_structure.hpp
All the information about the definition of the physical problem. The subroutines and functions are i...

CGeometry
Parent class for defining the geometry of the problem (complete geometry, multigrid agglomerated geom...
Definition: geometry_structure.hpp:72

option_structure.hpp
Defines classes for referencing options for easy input in CConfig.

std

mpi_structure.hpp
Headers of the mpi interface for generalized datatypes. The subroutines and functions are in the mpi_...

CMatrixVectorProduct
abstract base class for defining matrix-vector products
Definition: vector_structure.hpp:363

CPreconditioner
abstract base class for defining preconditioning operation
Definition: vector_structure.hpp:380

linear_solvers_structure.inl
inline subroutines of the linear_solvers_structure.hpp file.

CSysMatrix
Main class for defining sparse matrices-by-blocks with compressed row format.
Definition: matrix_structure.hpp:59

CSysSolve
Class for solving linear systems using classical and Krylov-subspace iterative methods.
Definition: linear_solvers_structure.hpp:67

vector_structure.hpp
Headers for the classes related to linear solvers (CG, FGMRES, etc) The subroutines and functions are...

CSysVector
Class for holding and manipulating vectors needed by linear solvers.
Definition: vector_structure.hpp:60

CConfig
Main class for defining the problem; basically this class reads the configuration file...
Definition: config_structure.hpp:68