libMesh
Public Member Functions | Static Public Member Functions | Public Attributes | Protected Types | Protected Member Functions | Static Protected Attributes | Private Attributes | List of all members
SlepcSolverConfiguration Class Reference

A class that interfaces the SolverConfiguration to add the SLEPC option SetST. More...

Inheritance diagram for SlepcSolverConfiguration:
[legend]

Public Member Functions

 SlepcSolverConfiguration (libMesh::SlepcEigenSolver< libMesh::Number > &slepc_eigen_solver)
 Constructur: get a reference to the SlepcEigenSolver variable to be able to manipulate it. More...
 
 ~SlepcSolverConfiguration ()
 empty destructor More...
 
virtual void configure_solver () override
 Apply solver options to a particular solver. More...
 
void SetST (SpectralTransform st)
 This is the main functionality of this little class. More...
 
virtual void set_options_during_init ()
 Apply options during initialization of a solver. More...
 
virtual void respond_to_solve_failure (unsigned int)
 This method can be called after the solver has failed (e.g. More...
 

Static Public Member Functions

static std::string get_info ()
 Gets a string containing the reference information. More...
 
static void print_info (std::ostream &out=libMesh::out)
 Prints the reference information, by default to libMesh::out. More...
 
static unsigned int n_objects ()
 Prints the number of outstanding (created, but not yet destroyed) objects. More...
 
static void enable_print_counter_info ()
 Methods to enable/disable the reference counter output from print_info() More...
 
static void disable_print_counter_info ()
 

Public Attributes

std::map< std::string, Real > real_valued_data
 Store real-valued solver parameters in this map, e.g. More...
 
std::map< std::string, intint_valued_data
 Store integer solver parameters in this map, e.g. More...
 
std::map< std::string, std::string > string_data
 Store string data in this map, e.g. More...
 

Protected Types

typedef std::map< std::string, std::pair< unsigned int, unsigned int > > Counts
 Data structure to log the information. More...
 

Protected Member Functions

void increment_constructor_count (const std::string &name)
 Increments the construction counter. More...
 
void increment_destructor_count (const std::string &name)
 Increments the destruction counter. More...
 

Static Protected Attributes

static Counts _counts
 Actually holds the data. More...
 
static Threads::atomic< unsigned int_n_objects
 The number of objects. More...
 
static Threads::spin_mutex _mutex
 Mutual exclusion object to enable thread-safe reference counting. More...
 
static bool _enable_print_counter = true
 Flag to control whether reference count information is printed when print_info is called. More...
 

Private Attributes

libMesh::SlepcEigenSolver< libMesh::Number > & _slepc_solver
 The linear solver object that we are configuring. More...
 
SpectralTransform _st
 

Detailed Description

A class that interfaces the SolverConfiguration to add the SLEPC option SetST.

In case of dense eigenvalues, the application of SINVERT or CAYLEY is highly benefitial see the SLEPC-manual for details.

Definition at line 90 of file miscellaneous_ex14.C.

Member Typedef Documentation

◆ Counts

typedef std::map<std::string, std::pair<unsigned int, unsigned int> > libMesh::ReferenceCounter::Counts
protectedinherited

Data structure to log the information.

The log is identified by the class name.

Definition at line 117 of file reference_counter.h.

Constructor & Destructor Documentation

◆ SlepcSolverConfiguration()

SlepcSolverConfiguration::SlepcSolverConfiguration ( libMesh::SlepcEigenSolver< libMesh::Number > &  slepc_eigen_solver)
inline

Constructur: get a reference to the SlepcEigenSolver variable to be able to manipulate it.

Definition at line 97 of file miscellaneous_ex14.C.

97  :
98  _slepc_solver(slepc_eigen_solver),
100  {}

◆ ~SlepcSolverConfiguration()

SlepcSolverConfiguration::~SlepcSolverConfiguration ( )
inline

empty destructor

Definition at line 105 of file miscellaneous_ex14.C.

105 {}

Member Function Documentation

◆ configure_solver()

void SlepcSolverConfiguration::configure_solver ( )
overridevirtual

Apply solver options to a particular solver.

Override in subclasses to provide specific behavior.

Implements libMesh::SolverConfiguration.

Definition at line 680 of file miscellaneous_ex14.C.

681 {
682  // if a spectral transformation was requested
683  if (_st!=INVALID_ST)
684  {
685  // initialise the st with the default values and change the spectral transformation value.
686  ST st;
687  PetscErrorCode ierr = EPSGetST(_slepc_solver.eps(), &st);
688  if (ierr)
689  libmesh_error();
690 
691  // Set it to the desired type of spectral transformation.
692  // The value of the respective shift is chosen to be the target
693  // specified via \p set_position_of_spectrum().
694  switch (_st)
695  {
696  case SHIFT:
697  ierr = STSetType(st, STSHIFT);
698  break;
699  case SINVERT:
700  // this method has been renamed in 3.1
701 #if SLEPC_VERSION_LESS_THAN(3,1,0)
702  ierr = STSetType(st, STSINV);
703 #else
704  ierr = STSetType(st, STSINVERT);
705 #endif
706  break;
707  case CAYLEY:
708  ierr = STSetType(st, STCAYLEY);
709  break;
710  default:
711  // print a warning but do nothing more.
712  break;
713  }
714  if (ierr)
715  libmesh_error();
716 
717  // since st is a reference to the particular object used by \p _slepc_solver,
718  // we don't need to hand back the manipulated object. It will be applied before
719  // solving the system automatically.
720  }
721 }

References CAYLEY, libMesh::ierr, INVALID_ST, SHIFT, and SINVERT.

◆ disable_print_counter_info()

void libMesh::ReferenceCounter::disable_print_counter_info ( )
staticinherited

Definition at line 106 of file reference_counter.C.

107 {
108  _enable_print_counter = false;
109  return;
110 }

References libMesh::ReferenceCounter::_enable_print_counter.

Referenced by libMesh::LibMeshInit::LibMeshInit().

◆ enable_print_counter_info()

void libMesh::ReferenceCounter::enable_print_counter_info ( )
staticinherited

Methods to enable/disable the reference counter output from print_info()

Definition at line 100 of file reference_counter.C.

101 {
102  _enable_print_counter = true;
103  return;
104 }

References libMesh::ReferenceCounter::_enable_print_counter.

◆ get_info()

std::string libMesh::ReferenceCounter::get_info ( )
staticinherited

Gets a string containing the reference information.

Definition at line 47 of file reference_counter.C.

48 {
49 #if defined(LIBMESH_ENABLE_REFERENCE_COUNTING) && defined(DEBUG)
50 
51  std::ostringstream oss;
52 
53  oss << '\n'
54  << " ---------------------------------------------------------------------------- \n"
55  << "| Reference count information |\n"
56  << " ---------------------------------------------------------------------------- \n";
57 
58  for (const auto & pr : _counts)
59  {
60  const std::string name(pr.first);
61  const unsigned int creations = pr.second.first;
62  const unsigned int destructions = pr.second.second;
63 
64  oss << "| " << name << " reference count information:\n"
65  << "| Creations: " << creations << '\n'
66  << "| Destructions: " << destructions << '\n';
67  }
68 
69  oss << " ---------------------------------------------------------------------------- \n";
70 
71  return oss.str();
72 
73 #else
74 
75  return "";
76 
77 #endif
78 }

References libMesh::ReferenceCounter::_counts, and libMesh::Quality::name().

Referenced by libMesh::ReferenceCounter::print_info().

◆ increment_constructor_count()

void libMesh::ReferenceCounter::increment_constructor_count ( const std::string &  name)
inlineprotectedinherited

Increments the construction counter.

Should be called in the constructor of any derived class that will be reference counted.

Definition at line 181 of file reference_counter.h.

182 {
183  Threads::spin_mutex::scoped_lock lock(Threads::spin_mtx);
184  std::pair<unsigned int, unsigned int> & p = _counts[name];
185 
186  p.first++;
187 }

References libMesh::ReferenceCounter::_counts, libMesh::Quality::name(), and libMesh::Threads::spin_mtx.

Referenced by libMesh::ReferenceCountedObject< RBParametrized >::ReferenceCountedObject().

◆ increment_destructor_count()

void libMesh::ReferenceCounter::increment_destructor_count ( const std::string &  name)
inlineprotectedinherited

Increments the destruction counter.

Should be called in the destructor of any derived class that will be reference counted.

Definition at line 194 of file reference_counter.h.

195 {
196  Threads::spin_mutex::scoped_lock lock(Threads::spin_mtx);
197  std::pair<unsigned int, unsigned int> & p = _counts[name];
198 
199  p.second++;
200 }

References libMesh::ReferenceCounter::_counts, libMesh::Quality::name(), and libMesh::Threads::spin_mtx.

Referenced by libMesh::ReferenceCountedObject< RBParametrized >::~ReferenceCountedObject().

◆ n_objects()

static unsigned int libMesh::ReferenceCounter::n_objects ( )
inlinestaticinherited

Prints the number of outstanding (created, but not yet destroyed) objects.

Definition at line 83 of file reference_counter.h.

84  { return _n_objects; }

References libMesh::ReferenceCounter::_n_objects.

◆ print_info()

void libMesh::ReferenceCounter::print_info ( std::ostream &  out = libMesh::out)
staticinherited

Prints the reference information, by default to libMesh::out.

Definition at line 87 of file reference_counter.C.

88 {
90  out_stream << ReferenceCounter::get_info();
91 }

References libMesh::ReferenceCounter::_enable_print_counter, and libMesh::ReferenceCounter::get_info().

◆ respond_to_solve_failure()

virtual void libMesh::SolverConfiguration::respond_to_solve_failure ( unsigned int  )
inlinevirtualinherited

This method can be called after the solver has failed (e.g.

on catching an exception thrown by a solver). It allows an appropriate response to the failure, e.g. terminate the program, or change configuration options and try again. solve_failure_count specifies the number of times we've tried to recover from a failure.

Definition at line 74 of file solver_configuration.h.

75  {
76  libmesh_error_msg("Solver error"); // by default just throw an error
77  }

◆ set_options_during_init()

virtual void libMesh::SolverConfiguration::set_options_during_init ( )
inlinevirtualinherited

Apply options during initialization of a solver.

Default is a no-op. Override in subclasses to provide specific behavior.

Definition at line 58 of file solver_configuration.h.

58 {}

◆ SetST()

void SlepcSolverConfiguration::SetST ( SpectralTransform  st)
inline

This is the main functionality of this little class.

Definition at line 112 of file miscellaneous_ex14.C.

113  { _st=st;}

Referenced by main().

Member Data Documentation

◆ _counts

ReferenceCounter::Counts libMesh::ReferenceCounter::_counts
staticprotectedinherited

◆ _enable_print_counter

bool libMesh::ReferenceCounter::_enable_print_counter = true
staticprotectedinherited

Flag to control whether reference count information is printed when print_info is called.

Definition at line 141 of file reference_counter.h.

Referenced by libMesh::ReferenceCounter::disable_print_counter_info(), libMesh::ReferenceCounter::enable_print_counter_info(), and libMesh::ReferenceCounter::print_info().

◆ _mutex

Threads::spin_mutex libMesh::ReferenceCounter::_mutex
staticprotectedinherited

Mutual exclusion object to enable thread-safe reference counting.

Definition at line 135 of file reference_counter.h.

◆ _n_objects

Threads::atomic< unsigned int > libMesh::ReferenceCounter::_n_objects
staticprotectedinherited

The number of objects.

Print the reference count information when the number returns to 0.

Definition at line 130 of file reference_counter.h.

Referenced by libMesh::ReferenceCounter::n_objects(), libMesh::ReferenceCounter::ReferenceCounter(), and libMesh::ReferenceCounter::~ReferenceCounter().

◆ _slepc_solver

libMesh::SlepcEigenSolver<libMesh::Number>& SlepcSolverConfiguration::_slepc_solver
private

The linear solver object that we are configuring.

Definition at line 120 of file miscellaneous_ex14.C.

◆ _st

SpectralTransform SlepcSolverConfiguration::_st
private

Definition at line 121 of file miscellaneous_ex14.C.

◆ int_valued_data

std::map<std::string, int> libMesh::SolverConfiguration::int_valued_data
inherited

Store integer solver parameters in this map, e.g.

iteration limits.

Definition at line 87 of file solver_configuration.h.

◆ real_valued_data

std::map<std::string, Real> libMesh::SolverConfiguration::real_valued_data
inherited

Store real-valued solver parameters in this map, e.g.

solver tolerances.

Definition at line 82 of file solver_configuration.h.

◆ string_data

std::map<std::string, std::string> libMesh::SolverConfiguration::string_data
inherited

Store string data in this map, e.g.

"solver_type" : "gmres".

Definition at line 92 of file solver_configuration.h.


The documentation for this class was generated from the following file:
CAYLEY
Definition: miscellaneous_ex14.C:80
libMesh::ReferenceCounter::_counts
static Counts _counts
Actually holds the data.
Definition: reference_counter.h:122
libMesh::ReferenceCounter::_n_objects
static Threads::atomic< unsigned int > _n_objects
The number of objects.
Definition: reference_counter.h:130
libMesh::ReferenceCounter::get_info
static std::string get_info()
Gets a string containing the reference information.
Definition: reference_counter.C:47
SlepcSolverConfiguration::_st
SpectralTransform _st
Definition: miscellaneous_ex14.C:121
libMesh::SlepcEigenSolver::eps
EPS eps()
Definition: slepc_eigen_solver.h:235
libMesh::Threads::spin_mtx
spin_mutex spin_mtx
A convenient spin mutex object which can be used for obtaining locks.
Definition: threads.C:29
SHIFT
Definition: miscellaneous_ex14.C:78
SINVERT
Definition: miscellaneous_ex14.C:79
libMesh::ReferenceCounter::_enable_print_counter
static bool _enable_print_counter
Flag to control whether reference count information is printed when print_info is called.
Definition: reference_counter.h:141
ierr
PetscErrorCode ierr
Definition: petscdmlibmeshimpl.C:1028
SlepcSolverConfiguration::_slepc_solver
libMesh::SlepcEigenSolver< libMesh::Number > & _slepc_solver
The linear solver object that we are configuring.
Definition: miscellaneous_ex14.C:120
libMesh::Quality::name
std::string name(const ElemQuality q)
This function returns a string containing some name for q.
Definition: elem_quality.C:42
INVALID_ST
Definition: miscellaneous_ex14.C:82