escript  Revision_
Namespaces | Classes | Typedefs | Enumerations | Functions | Variables
paso Namespace Reference

Namespaces

 util
 

Classes

struct  Connector
 
struct  Coupler
 
struct  DegreeAndIdx
 
struct  FCT_FluxLimiter
 
struct  FCT_Solver
 
struct  Function
 
struct  LinearSystem
 
struct  MergedSolver
 
struct  Options
 
class  PasoException
 PasoException exception class. More...
 
struct  Pattern
 
struct  Performance
 
struct  Preconditioner
 
struct  Preconditioner_AMG
 Local preconditioner. More...
 
struct  Preconditioner_AMG_Root
 
struct  Preconditioner_BoomerAMG
 
struct  Preconditioner_LocalAMG
 Local AMG preconditioner. More...
 
struct  Preconditioner_LocalSmoother
 
struct  Preconditioner_Smoother
 
struct  ReactiveSolver
 
struct  SharedComponents
 
struct  Solver_ILU
 ILU preconditioner. More...
 
struct  Solver_RILU
 RILU preconditioner. More...
 
struct  SparseMatrix
 
class  SystemMatrix
 this class holds a (distributed) stiffness matrix More...
 
struct  SystemMatrixPattern
 
class  TransportProblem
 
struct  UMFPACK_Handler
 

Typedefs

typedef boost::shared_ptr< ConnectorConnector_ptr
 
typedef boost::shared_ptr< const Connectorconst_Connector_ptr
 
typedef boost::shared_ptr< CouplerCoupler_ptr
 
typedef boost::shared_ptr< const Couplerconst_Coupler_ptr
 
typedef boost::shared_ptr< PatternPattern_ptr
 
typedef boost::shared_ptr< const Patternconst_Pattern_ptr
 
typedef boost::shared_ptr< PreconditionerPreconditioner_ptr
 
typedef boost::shared_ptr< const Preconditionerconst_Preconditioner_ptr
 
typedef boost::shared_ptr< SharedComponentsSharedComponents_ptr
 
typedef boost::shared_ptr< const SharedComponentsconst_SharedComponents_ptr
 
typedef boost::shared_ptr< SparseMatrixSparseMatrix_ptr
 
typedef boost::shared_ptr< const SparseMatrixconst_SparseMatrix_ptr
 
typedef int SparseMatrixType
 
typedef boost::shared_ptr< SystemMatrixSystemMatrix_ptr
 
typedef boost::shared_ptr< const SystemMatrixconst_SystemMatrix_ptr
 
typedef int SystemMatrixType
 
typedef boost::shared_ptr< SystemMatrixPatternSystemMatrixPattern_ptr
 
typedef boost::shared_ptr< const SystemMatrixPatternconst_SystemMatrixPattern_ptr
 
typedef boost::shared_ptr< TransportProblemTransportProblem_ptr
 
typedef boost::shared_ptr< const TransportProblemconst_TransportProblem_ptr
 

Enumerations

enum  SolverResult {
  NoError = 0, MaxIterReached, InputError, MemoryError,
  Breakdown, NegativeNormError, Divergence
}
 
enum  AMGBlockSelect { PASO_AMG_UNDECIDED =-1, PASO_AMG_IN_F =0, PASO_AMG_IN_C =1 }
 

Functions

double * createRandomVector (escript::const_Distribution_ptr dist)
 
void Preconditioner_AMG_free (Preconditioner_AMG *in)
 
int Preconditioner_AMG_getMaxLevel (const Preconditioner_AMG *in)
 
double Preconditioner_AMG_getCoarseLevelSparsity (const Preconditioner_AMG *in)
 
dim_t Preconditioner_AMG_getNumCoarseUnknowns (const Preconditioner_AMG *in)
 
Preconditioner_AMGPreconditioner_AMG_alloc (SystemMatrix_ptr A, int level, Options *options)
 
void Preconditioner_AMG_solve (SystemMatrix_ptr A, Preconditioner_AMG *amg, double *x, double *b)
 
void Preconditioner_AMG_setStrongConnections (SystemMatrix_ptr A, dim_t *degree_S, index_t *offset_S, index_t *S, double theta, double tau)
 
void Preconditioner_AMG_setStrongConnections_Block (SystemMatrix_ptr A, dim_t *degree_S, index_t *offset_S, index_t *S, double theta, double tau)
 
void Preconditioner_AMG_transposeStrongConnections (dim_t n, const dim_t *degree_S, const index_t *offset_S, const index_t *S, const dim_t nT, dim_t *degree_ST, index_t *offset_ST, index_t *ST)
 
void Preconditioner_AMG_CIJPCoarsening (dim_t n, dim_t my_n, AMGBlockSelect *split_marker, const dim_t *degree_S, const index_t *offset_S, const index_t *S, const dim_t *degree_ST, const index_t *offset_ST, const index_t *ST, const_Connector_ptr col_connector, escript::const_Distribution_ptr col_dist)
 
void Preconditioner_AMG_extendB (SystemMatrix_ptr A, SystemMatrix_ptr B)
 
void Preconditioner_AMG_CopyRemoteData (SystemMatrix_ptr P, index_t **p_ptr, index_t **p_idx, double **p_val, index_t *global_id, index_t block_size)
 
SystemMatrix_ptr Preconditioner_AMG_buildInterpolationOperator (SystemMatrix_ptr A, SystemMatrix_ptr P, SystemMatrix_ptr R)
 
SystemMatrix_ptr Preconditioner_AMG_buildInterpolationOperatorBlock (SystemMatrix_ptr A, SystemMatrix_ptr P, SystemMatrix_ptr R)
 
SystemMatrix_ptr Preconditioner_AMG_getProlongation (SystemMatrix_ptr A_p, const index_t *offset_S, const dim_t *degree_S, const index_t *S, const dim_t n_C, index_t *counter_C, const index_t interpolation_method)
 
void Preconditioner_AMG_setDirectProlongation (SystemMatrix_ptr P, SystemMatrix_ptr A, const index_t *offset_S, const dim_t *degree_S, const index_t *S, const index_t *counter_C)
 
void Preconditioner_AMG_setDirectProlongation_Block (SystemMatrix_ptr P, SystemMatrix_ptr A, const index_t *offset_S, const dim_t *degree_S, const index_t *S, const index_t *counter_C)
 
void Preconditioner_AMG_setClassicProlongation (SystemMatrix_ptr P, SystemMatrix_ptr A, const index_t *offset_S, const dim_t *degree_S, const index_t *S, const index_t *counter_C)
 
void Preconditioner_AMG_setClassicProlongation_Block (SystemMatrix_ptr P, SystemMatrix_ptr A, const index_t *offset_S, const dim_t *degree_S, const index_t *S, const index_t *counter_C)
 
SystemMatrix_ptr Preconditioner_AMG_getRestriction (SystemMatrix_ptr P)
 
void Preconditioner_AMG_Root_free (Preconditioner_AMG_Root *in)
 
Preconditioner_AMG_RootPreconditioner_AMG_Root_alloc (SystemMatrix_ptr A, Options *options)
 
void Preconditioner_AMG_Root_solve (SystemMatrix_ptr A, Preconditioner_AMG_Root *prec, double *x, double *b)
 
SolverResult Solver_BiCGStab (SystemMatrix_ptr A, double *r, double *x, dim_t *iter, double *tolerance, Performance *pp)
 
void BlockOps_Cpy_N (dim_t N, double *R, const double *V)
 
void BlockOps_SMV_2 (double *R, const double *mat, const double *V)
 performs operation R=R-mat*V (V and R are not overlapping) - 2x2 More...
 
void BlockOps_SMV_3 (double *R, const double *mat, const double *V)
 performs operation R=R-mat*V (V and R are not overlapping) - 3x3 More...
 
void BlockOps_SMV_N (dim_t N, double *R, const double *mat, const double *V)
 performs operation R=R-mat*V (V and R are not overlapping) - NxN More...
 
void BlockOps_MV_N (dim_t N, double *R, const double *mat, const double *V)
 
void BlockOps_invM_2 (double *invA, const double *A, int *failed)
 
void BlockOps_invM_3 (double *invA, const double *A, int *failed)
 
void BlockOps_invM_N (dim_t N, double *mat, index_t *pivot, int *failed)
 LU factorization of NxN matrix mat with partial pivoting. More...
 
void BlockOps_solve_N (dim_t N, double *X, double *mat, index_t *pivot, int *failed)
 solves system of linear equations A*X=B More...
 
void BlockOps_MViP_2 (const double *mat, double *V)
 inplace matrix vector product - order 2 More...
 
void BlockOps_MViP_3 (const double *mat, double *V)
 inplace matrix vector product - order 3 More...
 
void BlockOps_solveAll (dim_t n_block, dim_t n, double *D, index_t *pivot, double *x)
 
void Preconditioner_BoomerAMG_free (Preconditioner_BoomerAMG *in)
 
Preconditioner_BoomerAMGPreconditioner_BoomerAMG_alloc (SystemMatrix_ptr A, Options *options)
 
void Preconditioner_BoomerAMG_solve (SystemMatrix_ptr A, Preconditioner_BoomerAMG *amg, double *out, double *in)
 call the solver More...
 
SolverResult Solver_GMRES (SystemMatrix_ptr A, double *r, double *x, dim_t *iter, double *tolerance, dim_t Length_of_recursion, dim_t restart, Performance *pp)
 
SolverResult Solver_GMRES2 (Function *F, const double *f0, const double *x0, double *dx, dim_t *iter, double *tolerance, Performance *pp)
 
void Solver_ILU_free (Solver_ILU *in)
 
Solver_ILUSolver_getILU (SparseMatrix_ptr A, bool verbose)
 constructs the incomplete block factorization More...
 
void Solver_solveILU (SparseMatrix_ptr A, Solver_ILU *ilu, double *x, const double *b)
 
void Preconditioner_LocalAMG_free (Preconditioner_LocalAMG *in)
 
int Preconditioner_LocalAMG_getMaxLevel (const Preconditioner_LocalAMG *in)
 
double Preconditioner_LocalAMG_getCoarseLevelSparsity (const Preconditioner_LocalAMG *in)
 
dim_t Preconditioner_LocalAMG_getNumCoarseUnknowns (const Preconditioner_LocalAMG *in)
 
Preconditioner_LocalAMGPreconditioner_LocalAMG_alloc (SparseMatrix_ptr A_p, int level, Options *options)
 
void Preconditioner_LocalAMG_solve (SparseMatrix_ptr A, Preconditioner_LocalAMG *amg, double *x, double *b)
 
void Preconditioner_LocalAMG_setStrongConnections (SparseMatrix_ptr A, dim_t *degree_S, index_t *S, const double theta, const double tau)
 
void Preconditioner_LocalAMG_setStrongConnections_Block (SparseMatrix_ptr A, dim_t *degree_S, index_t *S, double theta, double tau)
 
void Preconditioner_LocalAMG_RungeStuebenSearch (dim_t n, const index_t *offset_S, const dim_t *degree_S, const index_t *S, AMGBlockSelect *split_marker, bool usePanel)
 
void Preconditioner_LocalAMG_enforceFFConnectivity (dim_t n, const index_t *offset_S, const dim_t *degree_S, const index_t *S, AMGBlockSelect *split_marker)
 ensures that two F nodes are connected via a C node More...
 
SparseMatrix_ptr Preconditioner_LocalAMG_getProlongation (SparseMatrix_ptr A_p, const index_t *offset_S, const dim_t *degree_S, const index_t *S, dim_t n_C, const index_t *counter_C, index_t interpolation_method)
 
void Preconditioner_LocalAMG_setDirectProlongation (SparseMatrix_ptr P_p, const_SparseMatrix_ptr A_p, const index_t *counter_C)
 
void Preconditioner_LocalAMG_setDirectProlongation_Block (SparseMatrix_ptr P_p, const_SparseMatrix_ptr A_p, const index_t *counter_C)
 
void Preconditioner_LocalAMG_setClassicProlongation (SparseMatrix_ptr P_p, SparseMatrix_ptr A_p, const index_t *offset_S, const dim_t *degree_S, const index_t *S, const index_t *counter_C)
 
void Preconditioner_LocalAMG_setClassicProlongation_Block (SparseMatrix_ptr P_p, SparseMatrix_ptr A_p, const index_t *offset_S, const dim_t *degree_S, const index_t *S, const index_t *counter_C)
 
SolverResult Solver_MINRES (SystemMatrix_ptr A, double *R, double *X, dim_t *iter, double *tolerance, Performance *pp)
 
void MKL_free (SparseMatrix *A)
 
void MKL_solve (SparseMatrix_ptr A, double *out, double *in, index_t reordering, dim_t numRefinements, bool verbose)
 
SolverResult Solver_NewtonGMRES (Function *F, double *x, Options *options, Performance *pp)
 
void checkPasoError ()
 
int comparDegreeAndIdx (const void *arg1, const void *arg2)
 
bool dropTree (index_t root, const Pattern *pattern, index_t *AssignedLevel, index_t *VerticesInTree, dim_t *numLevels, index_t *firstVertexInLevel, dim_t max_LevelWidth_abort, dim_t N)
 
SolverResult Solver_PCG (SystemMatrix_ptr A, double *r, double *x, dim_t *iter, double *tolerance, Performance *pp)
 
void Performance_open (Performance *pp, int verbose)
 sets up the monitoring process More...
 
int Performance_getEventIndex (Performance *pp, int event_id)
 find the index of an event in the list of monitored events More...
 
void Performance_close (Performance *pp, int verbose)
 shuts down the monitoring process More...
 
void Performance_startMonitor (Performance *pp, int monitor)
 switches on a monitor More...
 
void Performance_stopMonitor (Performance *pp, int monitor)
 switches off a monitor More...
 
void Preconditioner_free (Preconditioner *in)
 
PreconditionerPreconditioner_alloc (SystemMatrix_ptr A, Options *options)
 
void Preconditioner_solve (Preconditioner *prec, SystemMatrix_ptr A, double *x, double *b)
 
void Preconditioner_Smoother_free (Preconditioner_Smoother *in)
 
void Preconditioner_LocalSmoother_free (Preconditioner_LocalSmoother *in)
 
Preconditioner_SmootherPreconditioner_Smoother_alloc (SystemMatrix_ptr A, bool jacobi, bool is_local, bool verbose)
 constructs the symmetric Gauss-Seidel preconditioner More...
 
Preconditioner_LocalSmootherPreconditioner_LocalSmoother_alloc (SparseMatrix_ptr A, bool jacobi, bool verbose)
 
void Preconditioner_Smoother_solve (SystemMatrix_ptr A, Preconditioner_Smoother *gs, double *x, const double *b, dim_t sweeps, bool x_is_initial)
 
void Preconditioner_LocalSmoother_solve (SparseMatrix_ptr A, Preconditioner_LocalSmoother *gs, double *x, const double *b, dim_t sweeps, bool x_is_initial)
 
SolverResult Preconditioner_Smoother_solve_byTolerance (SystemMatrix_ptr A, Preconditioner_Smoother *gs, double *x, const double *b, double atol, dim_t *sweeps, bool x_is_initial)
 
void Preconditioner_LocalSmoother_Sweep (SparseMatrix_ptr A, Preconditioner_LocalSmoother *gs, double *x)
 
void Preconditioner_LocalSmoother_Sweep_sequential (SparseMatrix_ptr A, Preconditioner_LocalSmoother *smoother, double *x)
 inplace Gauss-Seidel sweep in sequential mode More...
 
void Preconditioner_LocalSmoother_Sweep_tiled (SparseMatrix_ptr A, Preconditioner_LocalSmoother *gs, double *x)
 
void Preconditioner_LocalSmoother_Sweep_colored (SparseMatrix_ptr A, Preconditioner_LocalSmoother *gs, double *x)
 
SparseMatrix_ptr Preconditioner_AMG_mergeSystemMatrix (SystemMatrix_ptr A)
 
void Preconditioner_AMG_mergeSolve (Preconditioner_AMG *amg)
 
void Solver_RILU_free (Solver_RILU *in)
 
Solver_RILUSolver_getRILU (SparseMatrix_ptr A, bool verbose)
 
void Solver_solveRILU (Solver_RILU *rilu, double *x, double *b)
 
void Solver_updateIncompleteSchurComplement (SparseMatrix_ptr A_CC, SparseMatrix_ptr A_CF, double *invA_FF, index_t *A_FF_pivot, SparseMatrix_ptr A_FC)
 
void solve_free (SystemMatrix *in)
 
void Solver_free (SystemMatrix *A)
 
SolverResult Solver (SystemMatrix_ptr A, double *x, double *b, Options *options, Performance *pp)
 calls the iterative solver More...
 
SolverResult Solver_TFQMR (SystemMatrix_ptr A, double *B, double *X, dim_t *iter, double *tolerance, Performance *pp)
 
void swap (index_t *r, index_t *c, double *v, int left, int right)
 
void q_sort (index_t *row, index_t *col, double *val, int begin, int end, int N)
 
void SparseMatrix_MatrixVector_CSC_OFFSET0 (const double alpha, const_SparseMatrix_ptr A, const double *in, const double beta, double *out)
 
void SparseMatrix_MatrixVector_CSC_OFFSET1 (const double alpha, const_SparseMatrix_ptr A, const double *in, const double beta, double *out)
 
void SparseMatrix_MatrixVector_CSR_OFFSET0 (const double alpha, const_SparseMatrix_ptr A, const double *in, const double beta, double *out)
 
void SparseMatrix_MatrixVector_CSR_OFFSET1 (const double alpha, const_SparseMatrix_ptr A, const double *in, const double beta, double *out)
 
void SparseMatrix_MatrixVector_CSR_OFFSET0_DIAG (const double alpha, const_SparseMatrix_ptr A, const double *in, const double beta, double *out)
 
SparseMatrix_ptr SparseMatrix_MatrixMatrix (const_SparseMatrix_ptr A, const_SparseMatrix_ptr B)
 
SparseMatrix_ptr SparseMatrix_MatrixMatrixTranspose (const_SparseMatrix_ptr A, const_SparseMatrix_ptr B, const_SparseMatrix_ptr T)
 
void SparseMatrix_MatrixMatrix_DD (SparseMatrix_ptr C, const_SparseMatrix_ptr A, const_SparseMatrix_ptr B)
 
void SparseMatrix_MatrixMatrix_DB (SparseMatrix_ptr C, const_SparseMatrix_ptr A, const_SparseMatrix_ptr B)
 
void SparseMatrix_MatrixMatrix_BD (SparseMatrix_ptr C, const_SparseMatrix_ptr A, const_SparseMatrix_ptr B)
 
void SparseMatrix_MatrixMatrix_BB (SparseMatrix_ptr C, const_SparseMatrix_ptr A, const_SparseMatrix_ptr B)
 
void SparseMatrix_MatrixMatrixTranspose_DD (SparseMatrix_ptr C, const_SparseMatrix_ptr A, const_SparseMatrix_ptr B, const_SparseMatrix_ptr T)
 
void SparseMatrix_MatrixMatrixTranspose_DB (SparseMatrix_ptr C, const_SparseMatrix_ptr A, const_SparseMatrix_ptr B, const_SparseMatrix_ptr T)
 
void SparseMatrix_MatrixMatrixTranspose_BD (SparseMatrix_ptr C, const_SparseMatrix_ptr A, const_SparseMatrix_ptr B, const_SparseMatrix_ptr T)
 
void SparseMatrix_MatrixMatrixTranspose_BB (SparseMatrix_ptr C, const_SparseMatrix_ptr A, const_SparseMatrix_ptr B, const_SparseMatrix_ptr T)
 
void SparseMatrix_MatrixVector_CSR_OFFSET0_stripe (double alpha, dim_t nRows, dim_t row_block_size, dim_t col_block_size, const index_t *ptr, const index_t *index, const double *val, const double *in, double beta, double *out)
 
static int calc_digits (int)
 
static void fmt_str (int, int, int *, int *, int *, char *, char *)
 
static void print_data (std::ofstream &, int, int, int, char *, const void *, int, int)
 
static void generate_HB (std::ofstream &, dim_t *, dim_t *, const double *)
 
void RHS_loadMM_toCSR (const char *filename, double *b, dim_t size)
 
static void q_sort (index_t *, index_t *, double *, int, int)
 
void UMFPACK_free (SparseMatrix *A)
 frees any UMFPACK related data from the matrix More...
 
void UMFPACK_solve (SparseMatrix_ptr A, double *out, double *in, dim_t numRefinements, bool verbose)
 calls the solver More...
 

Variables

static const real_t LARGE_POSITIVE_FLOAT = escript::DataTypes::real_t_max()
 
static double Pattern_mis_seed =.4142135623730951
 
static const real_t EPSILON = escript::DataTypes::real_t_eps()
 
static const real_t PASO_RT_EXP_LIM_MIN = sqrt(EPSILON)
 
static const real_t PASO_RT_EXP_LIM_MAX = log(1./sqrt(EPSILON))
 
static dim_t M
 
static dim_t N
 
static dim_t nz
 
static int M
 
static int N
 
static int nz
 
static const real_t LARGE_POSITIVE_FLOAT = escript::DataTypes::real_t_max()
 

Typedef Documentation

◆ Connector_ptr

typedef boost::shared_ptr<Connector> paso::Connector_ptr

◆ const_Connector_ptr

typedef boost::shared_ptr<const Connector> paso::const_Connector_ptr

◆ const_Coupler_ptr

typedef boost::shared_ptr<const Coupler> paso::const_Coupler_ptr

◆ const_Pattern_ptr

typedef boost::shared_ptr<const Pattern> paso::const_Pattern_ptr

◆ const_Preconditioner_ptr

typedef boost::shared_ptr<const Preconditioner> paso::const_Preconditioner_ptr

◆ const_SharedComponents_ptr

typedef boost::shared_ptr<const SharedComponents> paso::const_SharedComponents_ptr

◆ const_SparseMatrix_ptr

typedef boost::shared_ptr<const SparseMatrix> paso::const_SparseMatrix_ptr

◆ const_SystemMatrix_ptr

typedef boost::shared_ptr<const SystemMatrix> paso::const_SystemMatrix_ptr

◆ const_SystemMatrixPattern_ptr

typedef boost::shared_ptr<const SystemMatrixPattern> paso::const_SystemMatrixPattern_ptr

◆ const_TransportProblem_ptr

typedef boost::shared_ptr<const TransportProblem> paso::const_TransportProblem_ptr

◆ Coupler_ptr

typedef boost::shared_ptr<Coupler> paso::Coupler_ptr

◆ Pattern_ptr

typedef boost::shared_ptr<Pattern> paso::Pattern_ptr

◆ Preconditioner_ptr

typedef boost::shared_ptr<Preconditioner> paso::Preconditioner_ptr

◆ SharedComponents_ptr

typedef boost::shared_ptr<SharedComponents> paso::SharedComponents_ptr

◆ SparseMatrix_ptr

typedef boost::shared_ptr<SparseMatrix> paso::SparseMatrix_ptr

◆ SparseMatrixType

◆ SystemMatrix_ptr

typedef boost::shared_ptr<SystemMatrix> paso::SystemMatrix_ptr

◆ SystemMatrixPattern_ptr

◆ SystemMatrixType

◆ TransportProblem_ptr

typedef boost::shared_ptr<TransportProblem> paso::TransportProblem_ptr

Enumeration Type Documentation

◆ AMGBlockSelect

Enumerator
PASO_AMG_UNDECIDED 
PASO_AMG_IN_F 
PASO_AMG_IN_C 

◆ SolverResult

Enumerator
NoError 
MaxIterReached 
InputError 
MemoryError 
Breakdown 
NegativeNormError 
Divergence 

Function Documentation

◆ BlockOps_Cpy_N()

void paso::BlockOps_Cpy_N ( dim_t  N,
double *  R,
const double *  V 
)
inline

◆ BlockOps_invM_2()

void paso::BlockOps_invM_2 ( double *  invA,
const double *  A,
int *  failed 
)
inline

◆ BlockOps_invM_3()

void paso::BlockOps_invM_3 ( double *  invA,
const double *  A,
int *  failed 
)
inline

◆ BlockOps_invM_N()

void paso::BlockOps_invM_N ( dim_t  N,
double *  mat,
index_t *  pivot,
int *  failed 
)
inline

LU factorization of NxN matrix mat with partial pivoting.

References PASO_MISSING_CLAPACK.

Referenced by paso::SparseMatrix::invMain().

◆ BlockOps_MV_N()

void paso::BlockOps_MV_N ( dim_t  N,
double *  R,
const double *  mat,
const double *  V 
)
inline

◆ BlockOps_MViP_2()

void paso::BlockOps_MViP_2 ( const double *  mat,
double *  V 
)
inline

◆ BlockOps_MViP_3()

void paso::BlockOps_MViP_3 ( const double *  mat,
double *  V 
)
inline

◆ BlockOps_SMV_2()

void paso::BlockOps_SMV_2 ( double *  R,
const double *  mat,
const double *  V 
)
inline

performs operation R=R-mat*V (V and R are not overlapping) - 2x2

Referenced by Preconditioner_LocalSmoother_Sweep_colored(), and Preconditioner_LocalSmoother_Sweep_sequential().

◆ BlockOps_SMV_3()

void paso::BlockOps_SMV_3 ( double *  R,
const double *  mat,
const double *  V 
)
inline

performs operation R=R-mat*V (V and R are not overlapping) - 3x3

Referenced by Preconditioner_LocalSmoother_Sweep_colored(), and Preconditioner_LocalSmoother_Sweep_sequential().

◆ BlockOps_SMV_N()

void paso::BlockOps_SMV_N ( dim_t  N,
double *  R,
const double *  mat,
const double *  V 
)
inline

performs operation R=R-mat*V (V and R are not overlapping) - NxN

References PASO_MISSING_CLAPACK.

Referenced by Preconditioner_LocalSmoother_Sweep_colored(), and Preconditioner_LocalSmoother_Sweep_sequential().

◆ BlockOps_solve_N()

void paso::BlockOps_solve_N ( dim_t  N,
double *  X,
double *  mat,
index_t *  pivot,
int *  failed 
)
inline

◆ BlockOps_solveAll()

void paso::BlockOps_solveAll ( dim_t  n_block,
dim_t  n,
double *  D,
index_t *  pivot,
double *  x 
)
inline

◆ calc_digits()

int paso::calc_digits ( int  var)
static

Referenced by generate_HB().

◆ checkPasoError()

void paso::checkPasoError ( )

◆ comparDegreeAndIdx()

int paso::comparDegreeAndIdx ( const void *  arg1,
const void *  arg2 
)

◆ createRandomVector()

double* paso::createRandomVector ( escript::const_Distribution_ptr  dist)
inline

◆ dropTree()

bool paso::dropTree ( index_t  root,
const Pattern pattern,
index_t *  AssignedLevel,
index_t *  VerticesInTree,
dim_t *  numLevels,
index_t *  firstVertexInLevel,
dim_t  max_LevelWidth_abort,
dim_t  N 
)

◆ fmt_str()

void paso::fmt_str ( int  nvalues,
int  integer,
int *  width,
int *  nlines,
int *  nperline,
char *  pfmt,
char *  fmt 
)
static

Referenced by generate_HB().

◆ generate_HB()

void paso::generate_HB ( std::ofstream &  fp,
dim_t *  col_ptr,
dim_t *  row_ind,
const double *  val 
)
static

◆ MKL_free()

void paso::MKL_free ( SparseMatrix A)

◆ MKL_solve()

void paso::MKL_solve ( SparseMatrix_ptr  A,
double *  out,
double *  in,
index_t  reordering,
dim_t  numRefinements,
bool  verbose 
)

◆ Performance_close()

void paso::Performance_close ( Performance pp,
int  verbose 
)

◆ Performance_getEventIndex()

int paso::Performance_getEventIndex ( Performance pp,
int  event_id 
)

find the index of an event in the list of monitored events

References PERFORMANCE_UNMONITORED_EVENT.

Referenced by Performance_close().

◆ Performance_open()

void paso::Performance_open ( Performance pp,
int  verbose 
)

sets up the monitoring process

References PERFORMANCE_NUM_EVENTS, PERFORMANCE_NUM_MONITORS, and PERFORMANCE_UNUSED.

Referenced by paso::SystemMatrix::solve().

◆ Performance_startMonitor()

void paso::Performance_startMonitor ( Performance pp,
int  monitor 
)

◆ Performance_stopMonitor()

void paso::Performance_stopMonitor ( Performance pp,
int  monitor 
)

◆ Preconditioner_alloc()

Preconditioner * paso::Preconditioner_alloc ( SystemMatrix_ptr  A,
Options options 
)

◆ Preconditioner_AMG_alloc()

Preconditioner_AMG * paso::Preconditioner_AMG_alloc ( SystemMatrix_ptr  A,
int  level,
Options options 
)

◆ Preconditioner_AMG_buildInterpolationOperator()

SystemMatrix_ptr paso::Preconditioner_AMG_buildInterpolationOperator ( SystemMatrix_ptr  A,
SystemMatrix_ptr  P,
SystemMatrix_ptr  R 
)

◆ Preconditioner_AMG_buildInterpolationOperatorBlock()

SystemMatrix_ptr paso::Preconditioner_AMG_buildInterpolationOperatorBlock ( SystemMatrix_ptr  A,
SystemMatrix_ptr  P,
SystemMatrix_ptr  R 
)

◆ Preconditioner_AMG_CIJPCoarsening()

void paso::Preconditioner_AMG_CIJPCoarsening ( dim_t  n,
dim_t  my_n,
AMGBlockSelect split_marker,
const dim_t *  degree_S,
const index_t *  offset_S,
const index_t *  S,
const dim_t *  degree_ST,
const index_t *  offset_ST,
const index_t *  ST,
const_Connector_ptr  col_connector,
escript::const_Distribution_ptr  col_dist 
)

◆ Preconditioner_AMG_CopyRemoteData()

void paso::Preconditioner_AMG_CopyRemoteData ( SystemMatrix_ptr  P,
index_t **  p_ptr,
index_t **  p_idx,
double **  p_val,
index_t *  global_id,
index_t  block_size 
)

◆ Preconditioner_AMG_extendB()

void paso::Preconditioner_AMG_extendB ( SystemMatrix_ptr  A,
SystemMatrix_ptr  B 
)

◆ Preconditioner_AMG_free()

void paso::Preconditioner_AMG_free ( Preconditioner_AMG in)

◆ Preconditioner_AMG_getCoarseLevelSparsity()

double paso::Preconditioner_AMG_getCoarseLevelSparsity ( const Preconditioner_AMG in)

◆ Preconditioner_AMG_getMaxLevel()

int paso::Preconditioner_AMG_getMaxLevel ( const Preconditioner_AMG in)

◆ Preconditioner_AMG_getNumCoarseUnknowns()

dim_t paso::Preconditioner_AMG_getNumCoarseUnknowns ( const Preconditioner_AMG in)

◆ Preconditioner_AMG_getProlongation()

SystemMatrix_ptr paso::Preconditioner_AMG_getProlongation ( SystemMatrix_ptr  A_p,
const index_t *  offset_S,
const dim_t *  degree_S,
const index_t *  S,
const dim_t  n_C,
index_t *  counter_C,
const index_t  interpolation_method 
)

◆ Preconditioner_AMG_getRestriction()

SystemMatrix_ptr paso::Preconditioner_AMG_getRestriction ( SystemMatrix_ptr  P)

◆ Preconditioner_AMG_mergeSolve()

void paso::Preconditioner_AMG_mergeSolve ( Preconditioner_AMG amg)

◆ Preconditioner_AMG_mergeSystemMatrix()

SparseMatrix_ptr paso::Preconditioner_AMG_mergeSystemMatrix ( SystemMatrix_ptr  A)

◆ Preconditioner_AMG_Root_alloc()

Preconditioner_AMG_Root * paso::Preconditioner_AMG_Root_alloc ( SystemMatrix_ptr  A,
Options options 
)

◆ Preconditioner_AMG_Root_free()

void paso::Preconditioner_AMG_Root_free ( Preconditioner_AMG_Root in)

◆ Preconditioner_AMG_Root_solve()

void paso::Preconditioner_AMG_Root_solve ( SystemMatrix_ptr  A,
Preconditioner_AMG_Root prec,
double *  x,
double *  b 
)

◆ Preconditioner_AMG_setClassicProlongation()

void paso::Preconditioner_AMG_setClassicProlongation ( SystemMatrix_ptr  P,
SystemMatrix_ptr  A,
const index_t *  offset_S,
const dim_t *  degree_S,
const index_t *  S,
const index_t *  counter_C 
)

◆ Preconditioner_AMG_setClassicProlongation_Block()

void paso::Preconditioner_AMG_setClassicProlongation_Block ( SystemMatrix_ptr  P,
SystemMatrix_ptr  A,
const index_t *  offset_S,
const dim_t *  degree_S,
const index_t *  S,
const index_t *  counter_C 
)

◆ Preconditioner_AMG_setDirectProlongation()

void paso::Preconditioner_AMG_setDirectProlongation ( SystemMatrix_ptr  P,
SystemMatrix_ptr  A,
const index_t *  offset_S,
const dim_t *  degree_S,
const index_t *  S,
const index_t *  counter_C 
)

◆ Preconditioner_AMG_setDirectProlongation_Block()

void paso::Preconditioner_AMG_setDirectProlongation_Block ( SystemMatrix_ptr  P,
SystemMatrix_ptr  A,
const index_t *  offset_S,
const dim_t *  degree_S,
const index_t *  S,
const index_t *  counter_C 
)

◆ Preconditioner_AMG_setStrongConnections()

void paso::Preconditioner_AMG_setStrongConnections ( SystemMatrix_ptr  A,
dim_t *  degree_S,
index_t *  offset_S,
index_t *  S,
double  theta,
double  tau 
)

theta = threshold for strong connections tau = threshold for diagonal dominance

S_i={j N_i; i strongly coupled to j} in the sense that |A_{ij}| >= theta * max_k |A_{ik}|

Referenced by Preconditioner_AMG_alloc().

◆ Preconditioner_AMG_setStrongConnections_Block()

void paso::Preconditioner_AMG_setStrongConnections_Block ( SystemMatrix_ptr  A,
dim_t *  degree_S,
index_t *  offset_S,
index_t *  S,
double  theta,
double  tau 
)

◆ Preconditioner_AMG_solve()

void paso::Preconditioner_AMG_solve ( SystemMatrix_ptr  A,
Preconditioner_AMG amg,
double *  x,
double *  b 
)

◆ Preconditioner_AMG_transposeStrongConnections()

void paso::Preconditioner_AMG_transposeStrongConnections ( dim_t  n,
const dim_t *  degree_S,
const index_t *  offset_S,
const index_t *  S,
const dim_t  nT,
dim_t *  degree_ST,
index_t *  offset_ST,
index_t *  ST 
)

◆ Preconditioner_BoomerAMG_alloc()

Preconditioner_BoomerAMG * paso::Preconditioner_BoomerAMG_alloc ( SystemMatrix_ptr  A,
Options options 
)

◆ Preconditioner_BoomerAMG_free()

void paso::Preconditioner_BoomerAMG_free ( Preconditioner_BoomerAMG in)

◆ Preconditioner_BoomerAMG_solve()

void paso::Preconditioner_BoomerAMG_solve ( SystemMatrix_ptr  A,
Preconditioner_BoomerAMG amg,
double *  out,
double *  in 
)

◆ Preconditioner_free()

void paso::Preconditioner_free ( Preconditioner in)

◆ Preconditioner_LocalAMG_alloc()

Preconditioner_LocalAMG * paso::Preconditioner_LocalAMG_alloc ( SparseMatrix_ptr  A_p,
int  level,
Options options 
)

References paso::Preconditioner_LocalAMG::A_C, paso::Preconditioner_LocalAMG::AMG_C, paso::Preconditioner_LocalAMG::b_C, paso::Options::coarse_matrix_refinements, paso::Options::coarsening_selection_time, paso::Options::coarsening_threshold, paso::util::cumsum_maskedFalse(), paso::util::cumsum_maskedTrue(), paso::Options::diagonal_dominance_threshold, escript::gettime(), paso::Options::interpolation_method, paso::Preconditioner_LocalAMG::level, paso::Options::level_max, MATRIX_FORMAT_BLK1, MATRIX_FORMAT_CSC, MATRIX_FORMAT_OFFSET1, paso::Options::min_coarse_matrix_size, paso::Options::min_coarse_sparsity, paso::Preconditioner_LocalAMG::P, PASO_AMG_IN_C, PASO_AMG_IN_F, PASO_CLASSIC_INTERPOLATION_WITH_FF_COUPLING, PASO_JACOBI, PASO_MKL, PASO_SMOOTHER, PASO_UMFPACK, paso::Options::post_sweeps, paso::Preconditioner_LocalAMG::post_sweeps, paso::Options::pre_sweeps, paso::Preconditioner_LocalAMG::pre_sweeps, Preconditioner_LocalAMG_enforceFFConnectivity(), Preconditioner_LocalAMG_getProlongation(), Preconditioner_LocalAMG_RungeStuebenSearch(), Preconditioner_LocalAMG_setStrongConnections(), Preconditioner_LocalAMG_setStrongConnections_Block(), Preconditioner_LocalSmoother_alloc(), paso::Preconditioner_LocalAMG::R, paso::Preconditioner_LocalAMG::r, paso::Preconditioner_LocalAMG::refinements, paso::Options::reordering, paso::Preconditioner_LocalAMG::reordering, S, SHOW_TIMING, paso::Options::smoother, paso::Preconditioner_LocalAMG::Smoother, SparseMatrix_MatrixMatrix(), SparseMatrix_MatrixMatrixTranspose(), USE_TRANSPOSE, paso::Options::usePanel, paso::Options::verbose, and paso::Preconditioner_LocalAMG::x_C.

Referenced by Preconditioner_AMG_Root_alloc().

◆ Preconditioner_LocalAMG_enforceFFConnectivity()

void paso::Preconditioner_LocalAMG_enforceFFConnectivity ( dim_t  n,
const index_t *  offset_S,
const dim_t *  degree_S,
const index_t *  S,
AMGBlockSelect split_marker 
)

ensures that two F nodes are connected via a C node

References paso::util::comparIndex(), PASO_AMG_IN_C, and PASO_AMG_IN_F.

Referenced by Preconditioner_LocalAMG_alloc().

◆ Preconditioner_LocalAMG_free()

void paso::Preconditioner_LocalAMG_free ( Preconditioner_LocalAMG in)

◆ Preconditioner_LocalAMG_getCoarseLevelSparsity()

double paso::Preconditioner_LocalAMG_getCoarseLevelSparsity ( const Preconditioner_LocalAMG in)

◆ Preconditioner_LocalAMG_getMaxLevel()

int paso::Preconditioner_LocalAMG_getMaxLevel ( const Preconditioner_LocalAMG in)

◆ Preconditioner_LocalAMG_getNumCoarseUnknowns()

dim_t paso::Preconditioner_LocalAMG_getNumCoarseUnknowns ( const Preconditioner_LocalAMG in)

◆ Preconditioner_LocalAMG_getProlongation()

SparseMatrix_ptr paso::Preconditioner_LocalAMG_getProlongation ( SparseMatrix_ptr  A_p,
const index_t *  offset_S,
const dim_t *  degree_S,
const index_t *  S,
dim_t  n_C,
const index_t *  counter_C,
index_t  interpolation_method 
)

◆ Preconditioner_LocalAMG_RungeStuebenSearch()

void paso::Preconditioner_LocalAMG_RungeStuebenSearch ( dim_t  n,
const index_t *  offset_S,
const dim_t *  degree_S,
const index_t *  S,
AMGBlockSelect split_marker,
bool  usePanel 
)

◆ Preconditioner_LocalAMG_setClassicProlongation()

void paso::Preconditioner_LocalAMG_setClassicProlongation ( SparseMatrix_ptr  P_p,
SparseMatrix_ptr  A_p,
const index_t *  offset_S,
const dim_t *  degree_S,
const index_t *  S,
const index_t *  counter_C 
)

◆ Preconditioner_LocalAMG_setClassicProlongation_Block()

void paso::Preconditioner_LocalAMG_setClassicProlongation_Block ( SparseMatrix_ptr  P_p,
SparseMatrix_ptr  A_p,
const index_t *  offset_S,
const dim_t *  degree_S,
const index_t *  S,
const index_t *  counter_C 
)

◆ Preconditioner_LocalAMG_setDirectProlongation()

void paso::Preconditioner_LocalAMG_setDirectProlongation ( SparseMatrix_ptr  P_p,
const_SparseMatrix_ptr  A_p,
const index_t *  counter_C 
)

◆ Preconditioner_LocalAMG_setDirectProlongation_Block()

void paso::Preconditioner_LocalAMG_setDirectProlongation_Block ( SparseMatrix_ptr  P_p,
const_SparseMatrix_ptr  A_p,
const index_t *  counter_C 
)

◆ Preconditioner_LocalAMG_setStrongConnections()

void paso::Preconditioner_LocalAMG_setStrongConnections ( SparseMatrix_ptr  A,
dim_t *  degree_S,
index_t *  S,
const double  theta,
const double  tau 
)

◆ Preconditioner_LocalAMG_setStrongConnections_Block()

void paso::Preconditioner_LocalAMG_setStrongConnections_Block ( SparseMatrix_ptr  A,
dim_t *  degree_S,
index_t *  S,
double  theta,
double  tau 
)

◆ Preconditioner_LocalAMG_solve()

void paso::Preconditioner_LocalAMG_solve ( SparseMatrix_ptr  A,
Preconditioner_LocalAMG amg,
double *  x,
double *  b 
)

◆ Preconditioner_LocalSmoother_alloc()

Preconditioner_LocalSmoother * paso::Preconditioner_LocalSmoother_alloc ( SparseMatrix_ptr  A,
bool  jacobi,
bool  verbose 
)

◆ Preconditioner_LocalSmoother_free()

void paso::Preconditioner_LocalSmoother_free ( Preconditioner_LocalSmoother in)

◆ Preconditioner_LocalSmoother_solve()

void paso::Preconditioner_LocalSmoother_solve ( SparseMatrix_ptr  A,
Preconditioner_LocalSmoother gs,
double *  x,
const double *  b,
dim_t  sweeps,
bool  x_is_initial 
)

◆ Preconditioner_LocalSmoother_Sweep()

void paso::Preconditioner_LocalSmoother_Sweep ( SparseMatrix_ptr  A,
Preconditioner_LocalSmoother gs,
double *  x 
)

◆ Preconditioner_LocalSmoother_Sweep_colored()

void paso::Preconditioner_LocalSmoother_Sweep_colored ( SparseMatrix_ptr  A,
Preconditioner_LocalSmoother gs,
double *  x 
)

◆ Preconditioner_LocalSmoother_Sweep_sequential()

void paso::Preconditioner_LocalSmoother_Sweep_sequential ( SparseMatrix_ptr  A,
Preconditioner_LocalSmoother gs,
double *  x 
)

◆ Preconditioner_LocalSmoother_Sweep_tiled()

void paso::Preconditioner_LocalSmoother_Sweep_tiled ( SparseMatrix_ptr  A,
Preconditioner_LocalSmoother gs,
double *  x 
)

◆ Preconditioner_Smoother_alloc()

Preconditioner_Smoother * paso::Preconditioner_Smoother_alloc ( SystemMatrix_ptr  A,
bool  jacobi,
bool  is_local,
bool  verbose 
)

◆ Preconditioner_Smoother_free()

void paso::Preconditioner_Smoother_free ( Preconditioner_Smoother in)

◆ Preconditioner_Smoother_solve()

void paso::Preconditioner_Smoother_solve ( SystemMatrix_ptr  A,
Preconditioner_Smoother gs,
double *  x,
const double *  b,
dim_t  sweeps,
bool  x_is_initial 
)

◆ Preconditioner_Smoother_solve_byTolerance()

SolverResult paso::Preconditioner_Smoother_solve_byTolerance ( SystemMatrix_ptr  A,
Preconditioner_Smoother gs,
double *  x,
const double *  b,
double  atol,
dim_t *  sweeps,
bool  x_is_initial 
)

◆ Preconditioner_solve()

void paso::Preconditioner_solve ( Preconditioner prec,
SystemMatrix_ptr  A,
double *  x,
double *  b 
)

◆ print_data()

void paso::print_data ( std::ofstream &  fp,
int  n_perline,
int  width,
int  nval,
char *  fmt,
const void *  ptr,
int  integer,
int  adjust 
)
static

Referenced by generate_HB().

◆ q_sort() [1/2]

void paso::q_sort ( index_t *  row,
index_t *  col,
double *  val,
int  begin,
int  end 
)
static

References q_sort(), and swap().

◆ q_sort() [2/2]

void paso::q_sort ( index_t *  row,
index_t *  col,
double *  val,
int  begin,
int  end,
int  N 
)

◆ RHS_loadMM_toCSR()

void paso::RHS_loadMM_toCSR ( const char *  filename,
double *  b,
dim_t  size 
)

◆ solve_free()

void paso::solve_free ( SystemMatrix in)

◆ Solver()

SolverResult paso::Solver ( SystemMatrix_ptr  A,
double *  x,
double *  b,
Options options,
Performance pp 
)

◆ Solver_BiCGStab()

SolverResult paso::Solver_BiCGStab ( SystemMatrix_ptr  A,
double *  r,
double *  x,
dim_t *  iter,
double *  tolerance,
Performance pp 
)

◆ Solver_free()

void paso::Solver_free ( SystemMatrix A)

◆ Solver_getILU()

Solver_ILU * paso::Solver_getILU ( SparseMatrix_ptr  A,
bool  verbose 
)

constructs the incomplete block factorization

References paso::Solver_ILU::factors, and escript::gettime().

Referenced by Preconditioner_alloc().

◆ Solver_getRILU()

Solver_RILU * paso::Solver_getRILU ( SparseMatrix_ptr  A,
bool  verbose 
)

◆ Solver_GMRES()

SolverResult paso::Solver_GMRES ( SystemMatrix_ptr  A,
double *  r,
double *  x,
dim_t *  iter,
double *  tolerance,
dim_t  Length_of_recursion,
dim_t  restart,
Performance pp 
)

◆ Solver_GMRES2()

SolverResult paso::Solver_GMRES2 ( Function F,
const double *  f0,
const double *  x0,
double *  dx,
dim_t *  iter,
double *  tolerance,
Performance pp 
)

◆ Solver_ILU_free()

void paso::Solver_ILU_free ( Solver_ILU in)

◆ Solver_MINRES()

SolverResult paso::Solver_MINRES ( SystemMatrix_ptr  A,
double *  R,
double *  X,
dim_t *  iter,
double *  tolerance,
Performance pp 
)

◆ Solver_NewtonGMRES()

SolverResult paso::Solver_NewtonGMRES ( Function F,
double *  x,
Options options,
Performance pp 
)

◆ Solver_PCG()

SolverResult paso::Solver_PCG ( SystemMatrix_ptr  A,
double *  r,
double *  x,
dim_t *  iter,
double *  tolerance,
Performance pp 
)

◆ Solver_RILU_free()

void paso::Solver_RILU_free ( Solver_RILU in)

◆ Solver_solveILU()

void paso::Solver_solveILU ( SparseMatrix_ptr  A,
Solver_ILU ilu,
double *  x,
const double *  b 
)

◆ Solver_solveRILU()

void paso::Solver_solveRILU ( Solver_RILU rilu,
double *  x,
double *  b 
)

◆ Solver_TFQMR()

SolverResult paso::Solver_TFQMR ( SystemMatrix_ptr  A,
double *  B,
double *  X,
dim_t *  iter,
double *  tolerance,
Performance pp 
)

◆ Solver_updateIncompleteSchurComplement()

void paso::Solver_updateIncompleteSchurComplement ( SparseMatrix_ptr  A_CC,
SparseMatrix_ptr  A_CF,
double *  invA_FF,
index_t *  A_FF_pivot,
SparseMatrix_ptr  A_FC 
)

References paso::util::comparIndex().

Referenced by Solver_getRILU().

◆ SparseMatrix_MatrixMatrix()

SparseMatrix_ptr paso::SparseMatrix_MatrixMatrix ( const_SparseMatrix_ptr  A,
const_SparseMatrix_ptr  B 
)

◆ SparseMatrix_MatrixMatrix_BB()

void paso::SparseMatrix_MatrixMatrix_BB ( SparseMatrix_ptr  C,
const_SparseMatrix_ptr  A,
const_SparseMatrix_ptr  B 
)

◆ SparseMatrix_MatrixMatrix_BD()

void paso::SparseMatrix_MatrixMatrix_BD ( SparseMatrix_ptr  C,
const_SparseMatrix_ptr  A,
const_SparseMatrix_ptr  B 
)

◆ SparseMatrix_MatrixMatrix_DB()

void paso::SparseMatrix_MatrixMatrix_DB ( SparseMatrix_ptr  C,
const_SparseMatrix_ptr  A,
const_SparseMatrix_ptr  B 
)

◆ SparseMatrix_MatrixMatrix_DD()

void paso::SparseMatrix_MatrixMatrix_DD ( SparseMatrix_ptr  C,
const_SparseMatrix_ptr  A,
const_SparseMatrix_ptr  B 
)

◆ SparseMatrix_MatrixMatrixTranspose()

SparseMatrix_ptr paso::SparseMatrix_MatrixMatrixTranspose ( const_SparseMatrix_ptr  A,
const_SparseMatrix_ptr  B,
const_SparseMatrix_ptr  T 
)

◆ SparseMatrix_MatrixMatrixTranspose_BB()

void paso::SparseMatrix_MatrixMatrixTranspose_BB ( SparseMatrix_ptr  C,
const_SparseMatrix_ptr  A,
const_SparseMatrix_ptr  B,
const_SparseMatrix_ptr  T 
)

◆ SparseMatrix_MatrixMatrixTranspose_BD()

void paso::SparseMatrix_MatrixMatrixTranspose_BD ( SparseMatrix_ptr  C,
const_SparseMatrix_ptr  A,
const_SparseMatrix_ptr  B,
const_SparseMatrix_ptr  T 
)

◆ SparseMatrix_MatrixMatrixTranspose_DB()

void paso::SparseMatrix_MatrixMatrixTranspose_DB ( SparseMatrix_ptr  C,
const_SparseMatrix_ptr  A,
const_SparseMatrix_ptr  B,
const_SparseMatrix_ptr  T 
)

◆ SparseMatrix_MatrixMatrixTranspose_DD()

void paso::SparseMatrix_MatrixMatrixTranspose_DD ( SparseMatrix_ptr  C,
const_SparseMatrix_ptr  A,
const_SparseMatrix_ptr  B,
const_SparseMatrix_ptr  T 
)

◆ SparseMatrix_MatrixVector_CSC_OFFSET0()

void paso::SparseMatrix_MatrixVector_CSC_OFFSET0 ( const double  alpha,
const_SparseMatrix_ptr  A,
const double *  in,
const double  beta,
double *  out 
)

◆ SparseMatrix_MatrixVector_CSC_OFFSET1()

void paso::SparseMatrix_MatrixVector_CSC_OFFSET1 ( const double  alpha,
const_SparseMatrix_ptr  A,
const double *  in,
const double  beta,
double *  out 
)

◆ SparseMatrix_MatrixVector_CSR_OFFSET0()

void paso::SparseMatrix_MatrixVector_CSR_OFFSET0 ( const double  alpha,
const_SparseMatrix_ptr  A,
const double *  in,
const double  beta,
double *  out 
)

◆ SparseMatrix_MatrixVector_CSR_OFFSET0_DIAG()

void paso::SparseMatrix_MatrixVector_CSR_OFFSET0_DIAG ( const double  alpha,
const_SparseMatrix_ptr  A,
const double *  in,
const double  beta,
double *  out 
)

◆ SparseMatrix_MatrixVector_CSR_OFFSET0_stripe()

void paso::SparseMatrix_MatrixVector_CSR_OFFSET0_stripe ( double  alpha,
dim_t  nRows,
dim_t  row_block_size,
dim_t  col_block_size,
const index_t *  ptr,
const index_t *  index,
const double *  val,
const double *  in,
double  beta,
double *  out 
)

◆ SparseMatrix_MatrixVector_CSR_OFFSET1()

void paso::SparseMatrix_MatrixVector_CSR_OFFSET1 ( const double  alpha,
const_SparseMatrix_ptr  A,
const double *  in,
const double  beta,
double *  out 
)

◆ swap()

void paso::swap ( index_t *  r,
index_t *  c,
double *  v,
int  left,
int  right 
)

◆ UMFPACK_free()

void paso::UMFPACK_free ( SparseMatrix A)

frees any UMFPACK related data from the matrix

References paso::UMFPACK_Handler::numeric, paso::SparseMatrix::solver_p, and paso::UMFPACK_Handler::symbolic.

Referenced by solve_free(), and paso::SparseMatrix::~SparseMatrix().

◆ UMFPACK_solve()

void paso::UMFPACK_solve ( SparseMatrix_ptr  A,
double *  out,
double *  in,
dim_t  numRefinements,
bool  verbose 
)

Variable Documentation

◆ EPSILON

const real_t paso::EPSILON = escript::DataTypes::real_t_eps()
static

◆ LARGE_POSITIVE_FLOAT [1/2]

const real_t paso::LARGE_POSITIVE_FLOAT = escript::DataTypes::real_t_max()
static

◆ LARGE_POSITIVE_FLOAT [2/2]

const real_t paso::LARGE_POSITIVE_FLOAT = escript::DataTypes::real_t_max()
static

◆ M [1/2]

int paso::M
static

◆ M [2/2]

dim_t paso::M
static

◆ N [1/2]

int paso::N
static

◆ N [2/2]

dim_t paso::N
static

◆ nz [1/2]

dim_t paso::nz
static

◆ nz [2/2]

int paso::nz
static

◆ PASO_RT_EXP_LIM_MAX

const real_t paso::PASO_RT_EXP_LIM_MAX = log(1./sqrt(EPSILON))
static

◆ PASO_RT_EXP_LIM_MIN

const real_t paso::PASO_RT_EXP_LIM_MIN = sqrt(EPSILON)
static

◆ Pattern_mis_seed

double paso::Pattern_mis_seed =.4142135623730951
static