Template Numerical Library version\ main:be918e6f
Loading...
Searching...
No Matches
Public Types | Public Member Functions | Protected Member Functions | Protected Attributes | List of all members
TNL::Solvers::IterativeSolverMonitor< Real, Index > Class Template Reference

Object for monitoring convergence of iterative solvers. More...

#include <TNL/Solvers/IterativeSolverMonitor.h>

Inheritance diagram for TNL::Solvers::IterativeSolverMonitor< Real, Index >:
Inheritance graph
[legend]
Collaboration diagram for TNL::Solvers::IterativeSolverMonitor< Real, Index >:
Collaboration graph
[legend]

Public Types

using IndexType = Index
 A type for indexing.
 
using RealType = Real
 A type of the floating-point arithmetics.
 

Public Member Functions

 IterativeSolverMonitor ()=default
 Construct with no parameters.
 
void refresh () override
 Causes that the monitor prints out the status of the solver.
 
void setIterations (const IndexType &iterations)
 Set number of the current iteration.
 
void setNodesPerIteration (const IndexType &nodes)
 Set the number of nodes of the numerical mesh or lattice.
 
void setResidue (const RealType &residue)
 Set residue of the current approximation of the solution.
 
void setStage (const std::string &stage)
 This method can be used for naming a stage of the monitored solver.
 
void setTime (const RealType &time)
 Set the time of the simulated evolution if it is time dependent.
 
void setTimeStep (const RealType &timeStep)
 Set the time step for time dependent iterative solvers.
 
void setVerbose (const IndexType &verbose)
 Set up the verbosity of the monitor.
 
- Public Member Functions inherited from TNL::Solvers::SolverMonitor
 SolverMonitor ()=default
 Basic construct with no arguments.
 
bool isStopped () const
 Checks whether the main loop was stopped.
 
void runMainLoop ()
 Starts the main loop from which the method SolverMonitor::refresh is called in given time periods.
 
void setRefreshRate (const int &refreshRate)
 Set the time interval between two consecutive calls of SolverMonitor::refresh.
 
void setTimer (Timer &timer)
 Set a timer object for the solver monitor.
 
void stopMainLoop ()
 Stops the main loop of the monitor. See runMainLoop.
 

Protected Member Functions

int getLineWidth ()
 
- Protected Member Functions inherited from TNL::Solvers::SolverMonitor
double getElapsedTime ()
 

Protected Attributes

std::atomic_bool attributes_changed { false }
 
RealType elapsed_time_before_refresh = 0
 
IndexType iterations = 0
 
IndexType iterations_before_refresh = 0
 
RealType last_mlups = 0
 
IndexType nodesPerIteration = 0
 
RealType residue = 0
 
std::atomic_bool saved { false }
 
IndexType saved_iterations = 0
 
RealType saved_residue = 0
 
std::string saved_stage
 
RealType saved_time = 0
 
RealType saved_timeStep = 0
 
std::string stage
 
RealType time = 0
 
RealType timeStep = 0
 
IndexType verbose = 2
 
- Protected Attributes inherited from TNL::Solvers::SolverMonitor
std::atomic_bool started { false }
 
std::atomic_bool stopped { false }
 
std::atomic_int timeout_milliseconds { 500 }
 
Timer * timer = nullptr
 

Detailed Description

template<typename Real = double, typename Index = int>
class TNL::Solvers::IterativeSolverMonitor< Real, Index >

Object for monitoring convergence of iterative solvers.

Template Parameters
Realis a type of the floating-point arithmetics.
Indexis an indexing type.

The following example shows how to use the iterative solver monitor for monitoring convergence of linear iterative solver:

1#include <iostream>
2#include <memory>
3#include <chrono>
4#include <thread>
5#include <TNL/Matrices/SparseMatrix.h>
6#include <TNL/Devices/Sequential.h>
7#include <TNL/Devices/Cuda.h>
8#include <TNL/Solvers/Linear/Jacobi.h>
9
10template< typename Device >
11void
12iterativeLinearSolverExample()
13{
14 /***
15 * Set the following matrix (dots represent zero matrix elements):
16 *
17 * / 2.5 -1 . . . \
18 * | -1 2.5 -1 . . |
19 * | . -1 2.5 -1. . |
20 * | . . -1 2.5 -1 |
21 * \ . . . -1 2.5 /
22 */
23 using MatrixType = TNL::Matrices::SparseMatrix< double, Device >;
24 using Vector = TNL::Containers::Vector< double, Device >;
25 const int size( 5 );
26 auto matrix_ptr = std::make_shared< MatrixType >();
27 matrix_ptr->setDimensions( size, size );
28 matrix_ptr->setRowCapacities( Vector( { 2, 3, 3, 3, 2 } ) );
29
30 auto f = [ = ] __cuda_callable__( typename MatrixType::RowView & row ) mutable
31 {
32 const int rowIdx = row.getRowIndex();
33 if( rowIdx == 0 ) {
34 row.setElement( 0, rowIdx, 2.5 ); // diagonal element
35 row.setElement( 1, rowIdx + 1, -1 ); // element above the diagonal
36 }
37 else if( rowIdx == size - 1 ) {
38 row.setElement( 0, rowIdx - 1, -1.0 ); // element below the diagonal
39 row.setElement( 1, rowIdx, 2.5 ); // diagonal element
40 }
41 else {
42 row.setElement( 0, rowIdx - 1, -1.0 ); // element below the diagonal
43 row.setElement( 1, rowIdx, 2.5 ); // diagonal element
44 row.setElement( 2, rowIdx + 1, -1.0 ); // element above the diagonal
45 }
46 };
47
48 /***
49 * Set the matrix elements.
50 */
51 matrix_ptr->forAllRows( f );
52 std::cout << *matrix_ptr << std::endl;
53
54 /***
55 * Set the right-hand side vector.
56 */
57 Vector x( size, 1.0 );
58 Vector b( size );
59 matrix_ptr->vectorProduct( x, b );
60 x = 0.0;
61 std::cout << "Vector b = " << b << std::endl;
62
63 /***
64 * Setup solver of the linear system.
65 */
66 using LinearSolver = TNL::Solvers::Linear::Jacobi< MatrixType >;
67 LinearSolver solver;
68 solver.setMatrix( matrix_ptr );
69 solver.setOmega( 0.0005 );
70
71 /***
72 * Setup monitor of the iterative solver.
73 */
74 using IterativeSolverMonitorType = TNL::Solvers::IterativeSolverMonitor< double, int >;
75 IterativeSolverMonitorType monitor;
76 TNL::Solvers::SolverMonitorThread monitorThread( monitor );
77 monitor.setRefreshRate( 10 ); // refresh rate in milliseconds
78 monitor.setVerbose( 1 );
79 monitor.setStage( "Jacobi stage:" );
80 solver.setSolverMonitor( monitor );
81 solver.setConvergenceResidue( 1.0e-6 );
82 solver.solve( b, x );
83 monitor.stopMainLoop();
84 std::cout << "Vector x = " << x << std::endl;
85}
86
87int
88main( int argc, char* argv[] )
89{
90 std::cout << "Solving linear system on host: " << std::endl;
91 iterativeLinearSolverExample< TNL::Devices::Sequential >();
92
93#ifdef __CUDACC__
94 std::cout << "Solving linear system on CUDA device: " << std::endl;
95 iterativeLinearSolverExample< TNL::Devices::Cuda >();
96#endif
97}
__cuda_callable__
#define __cuda_callable__
Definition Macros.h:49
TNL::Containers::Vector
Vector extends Array with algebraic operations.
Definition Vector.h:36
TNL::Matrices::SparseMatrix
Implementation of sparse matrix, i.e. matrix storing only non-zero elements.
Definition SparseMatrix.h:57
TNL::Solvers::Linear::Jacobi
Iterative solver of linear systems based on the Jacobi method.
Definition Jacobi.h:21
TNL::Solvers::Linear::LinearSolver::setMatrix
void setMatrix(const MatrixPointer &matrix)
Set the matrix of the linear system.
Definition LinearSolver.h:120
TNL::Solvers::SolverMonitorThread
A RAII wrapper for launching the SolverMonitor's main loop in a separate thread.
Definition SolverMonitor.h:133
std::endl
T endl(T... args)
std::make_shared
T make_shared(T... args)

The result looks as follows:

Solving linear system on host:
Row: 0 -> 0:2.5 1:-1
Row: 1 -> 0:-1 1:2.5 2:-1
Row: 2 -> 1:-1 2:2.5 3:-1
Row: 3 -> 2:-1 3:2.5 4:-1
Row: 4 -> 3:-1 4:2.5
Vector b = [ 1.5, 0.5, 0.5, 0.5, 1.5 ]
Jacobi stage: ITER: 40899 RES: 0.0013506
Jacobi stage: ITER: 82863 RES: 2.1425e-06
Vector x = [ 0.999999, 0.999999, 0.999998, 0.999999, 0.999999 ]
Solving linear system on CUDA device:
Row: 0 -> 0:2.5 1:-1
Row: 1 -> 0:-1 1:2.5 2:-1
Row: 2 -> 1:-1 2:2.5 3:-1
Row: 3 -> 2:-1 3:2.5 4:-1
Row: 4 -> 3:-1 4:2.5
Vector b = [ 1.5, 0.5, 0.5, 0.5, 1.5 ]
Jacobi stage: ITER: 5675 RES: 0.30444
Jacobi stage: ITER: 10887 RES: 0.13573
Jacobi stage: ITER: 11675 RES: 0.12025
Jacobi stage: ITER: 12349 RES: 0.10838
Jacobi stage: ITER: 12987 RES: 0.098292
Jacobi stage: ITER: 13591 RES: 0.089582
Jacobi stage: ITER: 14131 RES: 0.082451
Jacobi stage: ITER: 14747 RES: 0.075006
Jacobi stage: ITER: 15355 RES: 0.068318
Jacobi stage: ITER: 15967 RES: 0.062189
Jacobi stage: ITER: 16583 RES: 0.056574
Jacobi stage: ITER: 17195 RES: 0.051498
Jacobi stage: ITER: 17811 RES: 0.046849
Jacobi stage: ITER: 18429 RES: 0.042593
Jacobi stage: ITER: 19051 RES: 0.038724
Jacobi stage: ITER: 19589 RES: 0.035642
Jacobi stage: ITER: 20225 RES: 0.032324
Jacobi stage: ITER: 20853 RES: 0.029352
Jacobi stage: ITER: 21445 RES: 0.026801
Jacobi stage: ITER: 22049 RES: 0.024426
Jacobi stage: ITER: 22663 RES: 0.022235
Jacobi stage: ITER: 23303 RES: 0.020153
Jacobi stage: ITER: 23939 RES: 0.018277
Jacobi stage: ITER: 24575 RES: 0.016576
Jacobi stage: ITER: 25211 RES: 0.015033
Jacobi stage: ITER: 25849 RES: 0.013626
Jacobi stage: ITER: 26487 RES: 0.012358
Jacobi stage: ITER: 27123 RES: 0.011208
Jacobi stage: ITER: 27759 RES: 0.010164
Jacobi stage: ITER: 28395 RES: 0.0092184
Jacobi stage: ITER: 29031 RES: 0.0083605
Jacobi stage: ITER: 29667 RES: 0.0075823
Jacobi stage: ITER: 30305 RES: 0.0068724
Jacobi stage: ITER: 30941 RES: 0.0062328
Jacobi stage: ITER: 31577 RES: 0.0056527
Jacobi stage: ITER: 32213 RES: 0.0051266
Jacobi stage: ITER: 32845 RES: 0.0046523
Jacobi stage: ITER: 33475 RES: 0.0042245
Jacobi stage: ITER: 34111 RES: 0.0038313
Jacobi stage: ITER: 34743 RES: 0.0034769
Jacobi stage: ITER: 35375 RES: 0.0031552
Jacobi stage: ITER: 36007 RES: 0.0028633
Jacobi stage: ITER: 36639 RES: 0.0025984
Jacobi stage: ITER: 37271 RES: 0.0023581
Jacobi stage: ITER: 37903 RES: 0.0021399
Jacobi stage: ITER: 38533 RES: 0.0019419
Jacobi stage: ITER: 39167 RES: 0.0017623
Jacobi stage: ITER: 39799 RES: 0.0015992
Jacobi stage: ITER: 40429 RES: 0.0014513
Jacobi stage: ITER: 41061 RES: 0.001317
Jacobi stage: ITER: 41695 RES: 0.0011952
Jacobi stage: ITER: 42327 RES: 0.0010846
Jacobi stage: ITER: 42959 RES: 0.00098427
Jacobi stage: ITER: 43591 RES: 0.00089322
Jacobi stage: ITER: 44223 RES: 0.00081058
Jacobi stage: ITER: 44855 RES: 0.00073559
Jacobi stage: ITER: 45487 RES: 0.00066754
Jacobi stage: ITER: 46119 RES: 0.00060578
Jacobi stage: ITER: 46751 RES: 0.00054974
Jacobi stage: ITER: 47383 RES: 0.00049888
Jacobi stage: ITER: 48015 RES: 0.00045273
Jacobi stage: ITER: 48649 RES: 0.00041059
Jacobi stage: ITER: 49281 RES: 0.00037261
Jacobi stage: ITER: 49913 RES: 0.00033814
Jacobi stage: ITER: 50547 RES: 0.00030686
Jacobi stage: ITER: 51179 RES: 0.00027847
Jacobi stage: ITER: 51811 RES: 0.00025271
Jacobi stage: ITER: 52443 RES: 0.00022933
Jacobi stage: ITER: 53075 RES: 0.00020811
Jacobi stage: ITER: 53705 RES: 0.00018886
Jacobi stage: ITER: 54337 RES: 0.00017139
Jacobi stage: ITER: 54977 RES: 0.00015534
Jacobi stage: ITER: 55611 RES: 0.00014097
Jacobi stage: ITER: 56243 RES: 0.00012793
Jacobi stage: ITER: 56879 RES: 0.00011602
Jacobi stage: ITER: 57517 RES: 0.00010516
Jacobi stage: ITER: 58153 RES: 9.5371e-05
Jacobi stage: ITER: 58789 RES: 8.6495e-05
Jacobi stage: ITER: 59427 RES: 7.8445e-05
Jacobi stage: ITER: 60063 RES: 7.1144e-05
Jacobi stage: ITER: 60699 RES: 6.4522e-05
Jacobi stage: ITER: 61337 RES: 5.8481e-05
Jacobi stage: ITER: 61975 RES: 5.3038e-05
Jacobi stage: ITER: 62611 RES: 4.8102e-05
Jacobi stage: ITER: 63247 RES: 4.3625e-05
Jacobi stage: ITER: 63883 RES: 3.9565e-05
Jacobi stage: ITER: 64521 RES: 3.5861e-05
Jacobi stage: ITER: 65157 RES: 3.2523e-05
Jacobi stage: ITER: 65795 RES: 2.9496e-05
Jacobi stage: ITER: 66429 RES: 2.6751e-05
Jacobi stage: ITER: 67059 RES: 2.4291e-05
Jacobi stage: ITER: 67691 RES: 2.2044e-05
Jacobi stage: ITER: 68323 RES: 2.0004e-05
Jacobi stage: ITER: 68957 RES: 1.8143e-05
Jacobi stage: ITER: 69589 RES: 1.6464e-05
Jacobi stage: ITER: 70221 RES: 1.4941e-05
Jacobi stage: ITER: 70853 RES: 1.3559e-05
Jacobi stage: ITER: 71483 RES: 1.2312e-05
Jacobi stage: ITER: 72119 RES: 1.1166e-05
Jacobi stage: ITER: 72755 RES: 1.0127e-05
Jacobi stage: ITER: 73391 RES: 9.1843e-06
Jacobi stage: ITER: 74029 RES: 8.3244e-06
Jacobi stage: ITER: 74665 RES: 7.5497e-06
Jacobi stage: ITER: 75303 RES: 6.847e-06
Jacobi stage: ITER: 75939 RES: 6.2098e-06
Jacobi stage: ITER: 76577 RES: 5.6284e-06
Jacobi stage: ITER: 77213 RES: 5.1045e-06
Jacobi stage: ITER: 77849 RES: 4.6294e-06
Jacobi stage: ITER: 78487 RES: 4.1986e-06
Jacobi stage: ITER: 79123 RES: 3.8078e-06
Jacobi stage: ITER: 79759 RES: 3.4534e-06
Jacobi stage: ITER: 80399 RES: 3.1301e-06
Jacobi stage: ITER: 81035 RES: 2.8388e-06
Jacobi stage: ITER: 81673 RES: 2.573e-06
Jacobi stage: ITER: 82309 RES: 2.3335e-06
Jacobi stage: ITER: 82943 RES: 2.1176e-06
Jacobi stage: ITER: 83575 RES: 1.9217e-06
Jacobi stage: ITER: 84207 RES: 1.7439e-06
Jacobi stage: ITER: 84839 RES: 1.5826e-06
Jacobi stage: ITER: 85471 RES: 1.4362e-06
Jacobi stage: ITER: 86103 RES: 1.3033e-06
Jacobi stage: ITER: 86735 RES: 1.1827e-06
Jacobi stage: ITER: 87371 RES: 1.0727e-06
Vector x = [ 0.999999, 0.999999, 0.999998, 0.999999, 0.999999 ]

The following example shows how to employ timer (TNL::Timer) to the monitor of iterative solvers:

1#include <iostream>
2#include <memory>
3#include <chrono>
4#include <thread>
5#include <TNL/Timer.h>
6#include <TNL/Matrices/SparseMatrix.h>
7#include <TNL/Devices/Sequential.h>
8#include <TNL/Devices/Cuda.h>
9#include <TNL/Solvers/Linear/Jacobi.h>
10
11template< typename Device >
12void
13iterativeLinearSolverExample()
14{
15 /***
16 * Set the following matrix (dots represent zero matrix elements):
17 *
18 * / 2.5 -1 . . . \
19 * | -1 2.5 -1 . . |
20 * | . -1 2.5 -1. . |
21 * | . . -1 2.5 -1 |
22 * \ . . . -1 2.5 /
23 */
24 using MatrixType = TNL::Matrices::SparseMatrix< double, Device >;
25 using Vector = TNL::Containers::Vector< double, Device >;
26 const int size( 5 );
27 auto matrix_ptr = std::make_shared< MatrixType >();
28 matrix_ptr->setDimensions( size, size );
29 matrix_ptr->setRowCapacities( Vector( { 2, 3, 3, 3, 2 } ) );
30
31 auto f = [ = ] __cuda_callable__( typename MatrixType::RowView & row ) mutable
32 {
33 const int rowIdx = row.getRowIndex();
34 if( rowIdx == 0 ) {
35 row.setElement( 0, rowIdx, 2.5 ); // diagonal element
36 row.setElement( 1, rowIdx + 1, -1 ); // element above the diagonal
37 }
38 else if( rowIdx == size - 1 ) {
39 row.setElement( 0, rowIdx - 1, -1.0 ); // element below the diagonal
40 row.setElement( 1, rowIdx, 2.5 ); // diagonal element
41 }
42 else {
43 row.setElement( 0, rowIdx - 1, -1.0 ); // element below the diagonal
44 row.setElement( 1, rowIdx, 2.5 ); // diagonal element
45 row.setElement( 2, rowIdx + 1, -1.0 ); // element above the diagonal
46 }
47 };
48
49 /***
50 * Set the matrix elements.
51 */
52 matrix_ptr->forAllRows( f );
53 std::cout << *matrix_ptr << std::endl;
54
55 /***
56 * Set the right-hand side vector.
57 */
58 Vector x( size, 1.0 );
59 Vector b( size );
60 matrix_ptr->vectorProduct( x, b );
61 x = 0.0;
62 std::cout << "Vector b = " << b << std::endl;
63
64 /***
65 * Setup solver of the linear system.
66 */
67 using LinearSolver = TNL::Solvers::Linear::Jacobi< MatrixType >;
68 LinearSolver solver;
69 solver.setMatrix( matrix_ptr );
70 solver.setOmega( 0.0005 );
71
72 /***
73 * Setup monitor of the iterative solver.
74 */
75 using IterativeSolverMonitorType = TNL::Solvers::IterativeSolverMonitor< double, int >;
76 IterativeSolverMonitorType monitor;
77 TNL::Solvers::SolverMonitorThread mmonitorThread( monitor );
78 monitor.setRefreshRate( 10 ); // refresh rate in milliseconds
79 monitor.setVerbose( 1 );
80 monitor.setStage( "Jacobi stage:" );
81 TNL::Timer timer;
82 monitor.setTimer( timer );
83 timer.start();
84 solver.setSolverMonitor( monitor );
85 solver.setConvergenceResidue( 1.0e-6 );
86 solver.solve( b, x );
87 monitor.stopMainLoop();
88 std::cout << "Vector x = " << x << std::endl;
89}
90
91int
92main( int argc, char* argv[] )
93{
94 std::cout << "Solving linear system on host: " << std::endl;
95 iterativeLinearSolverExample< TNL::Devices::Sequential >();
96
97#ifdef __CUDACC__
98 std::cout << "Solving linear system on CUDA device: " << std::endl;
99 iterativeLinearSolverExample< TNL::Devices::Cuda >();
100#endif
101}
TNL::Timer
Class for real time, CPU time and CPU cycles measuring.
Definition Timer.h:25
TNL::Timer::start
void start()
Starts timer.
Definition Timer.hpp:42

The result looks as follows:

Solving linear system on host:
Row: 0 -> 0:2.5 1:-1
Row: 1 -> 0:-1 1:2.5 2:-1
Row: 2 -> 1:-1 2:2.5 3:-1
Row: 3 -> 2:-1 3:2.5 4:-1
Row: 4 -> 3:-1 4:2.5
Vector b = [ 1.5, 0.5, 0.5, 0.5, 1.5 ]
ELA:1.451e-0
ELA: 0.10013 Jacobi stage: ITER: 41531 RES: 0.0012249
ELA: 0.20026 Jacobi stage: ITER: 82731 RES: 2.1877e-06
Vector x = [ 0.999999, 0.999999, 0.999998, 0.999999, 0.999999 ]
Solving linear system on CUDA device:
Row: 0 -> 0:2.5 1:-1
Row: 1 -> 0:-1 1:2.5 2:-1
Row: 2 -> 1:-1 2:2.5 3:-1
Row: 3 -> 2:-1 3:2.5 4:-1
Row: 4 -> 3:-1 4:2.5
Vector b = [ 1.5, 0.5, 0.5, 0.5, 1.5 ]
ELA:1.5231e-
ELA: 0.10012 Jacobi stage: ITER: 691 RES: 0.80439
ELA: 0.20023 Jacobi stage: ITER: 1451 RES: 0.65919
ELA: 0.30036 Jacobi stage: ITER: 2039 RES: 0.57741
ELA: 0.40047 Jacobi stage: ITER: 2641 RES: 0.51089
ELA: 0.50058 Jacobi stage: ITER: 3175 RES: 0.46241
ELA: 0.60072 Jacobi stage: ITER: 3783 RES: 0.4153
ELA: 0.70084 Jacobi stage: ITER: 4391 RES: 0.37476
ELA: 0.80098 Jacobi stage: ITER: 5003 RES: 0.33903
ELA: 0.90114 Jacobi stage: ITER: 5611 RES: 0.30755
ELA: 1.0013 Jacobi stage: ITER: 6219 RES: 0.27938
ELA: 1.1014 Jacobi stage: ITER: 6829 RES: 0.25387
ELA: 1.2015 Jacobi stage: ITER: 7443 RES: 0.23082
ELA: 1.3016 Jacobi stage: ITER: 8055 RES: 0.20996
ELA: 1.4017 Jacobi stage: ITER: 8599 RES: 0.19304
ELA: 1.5018 Jacobi stage: ITER: 9193 RES: 0.17609
ELA: 1.6019 Jacobi stage: ITER: 9813 RES: 0.16006
ELA: 1.702 Jacobi stage: ITER: 10407 RES: 0.14613
ELA: 1.8021 Jacobi stage: ITER: 10997 RES: 0.13341
ELA: 1.9022 Jacobi stage: ITER: 11605 RES: 0.12151
ELA: 2.0023 Jacobi stage: ITER: 12239 RES: 0.11026
ELA: 2.1025 Jacobi stage: ITER: 12875 RES: 0.099998
ELA: 2.2026 Jacobi stage: ITER: 13511 RES: 0.09069
ELA: 2.3026 Jacobi stage: ITER: 14147 RES: 0.082248
ELA: 2.4027 Jacobi stage: ITER: 14785 RES: 0.074547
ELA: 2.5028 Jacobi stage: ITER: 15421 RES: 0.067608
ELA: 2.6029 Jacobi stage: ITER: 16059 RES: 0.061316
ELA: 2.703 Jacobi stage: ITER: 16695 RES: 0.055609
ELA: 2.8031 Jacobi stage: ITER: 17331 RES: 0.050433
ELA: 2.9032 Jacobi stage: ITER: 17967 RES: 0.04574
ELA: 3.0033 Jacobi stage: ITER: 18603 RES: 0.041483
ELA: 3.1034 Jacobi stage: ITER: 19239 RES: 0.037622
ELA: 3.2035 Jacobi stage: ITER: 19875 RES: 0.03412
ELA: 3.3036 Jacobi stage: ITER: 20513 RES: 0.030926
ELA: 3.4037 Jacobi stage: ITER: 21149 RES: 0.028047
ELA: 3.5038 Jacobi stage: ITER: 21781 RES: 0.025453
ELA: 3.6039 Jacobi stage: ITER: 22411 RES: 0.023112
ELA: 3.704 Jacobi stage: ITER: 23047 RES: 0.020961
ELA: 3.8041 Jacobi stage: ITER: 23679 RES: 0.019022
ELA: 3.9042 Jacobi stage: ITER: 24311 RES: 0.017262
ELA: 4.0043 Jacobi stage: ITER: 24941 RES: 0.015665
ELA: 4.1044 Jacobi stage: ITER: 25573 RES: 0.014216
ELA: 4.2045 Jacobi stage: ITER: 26205 RES: 0.012901
ELA: 4.3046 Jacobi stage: ITER: 26837 RES: 0.011707
ELA: 4.4048 Jacobi stage: ITER: 27467 RES: 0.010631
ELA: 4.5049 Jacobi stage: ITER: 28101 RES: 0.0096413
ELA: 4.605 Jacobi stage: ITER: 28733 RES: 0.0087493
ELA: 4.7051 Jacobi stage: ITER: 29365 RES: 0.0079399
ELA: 4.8052 Jacobi stage: ITER: 29995 RES: 0.0072098
ELA: 4.9053 Jacobi stage: ITER: 30627 RES: 0.0065428
ELA: 5.0054 Jacobi stage: ITER: 31259 RES: 0.0059375
ELA: 5.1055 Jacobi stage: ITER: 31891 RES: 0.0053882
ELA: 5.2056 Jacobi stage: ITER: 32523 RES: 0.0048897
ELA: 5.3057 Jacobi stage: ITER: 33155 RES: 0.0044374
ELA: 5.4058 Jacobi stage: ITER: 33787 RES: 0.0040268
ELA: 5.5059 Jacobi stage: ITER: 34419 RES: 0.0036543
ELA: 5.606 Jacobi stage: ITER: 35051 RES: 0.0033162
ELA: 5.7061 Jacobi stage: ITER: 35683 RES: 0.0030094
ELA: 5.8062 Jacobi stage: ITER: 36315 RES: 0.002731
ELA: 5.9063 Jacobi stage: ITER: 36947 RES: 0.0024784
ELA: 6.0064 Jacobi stage: ITER: 37579 RES: 0.0022491
ELA: 6.1065 Jacobi stage: ITER: 38211 RES: 0.002041
ELA: 6.2066 Jacobi stage: ITER: 38843 RES: 0.0018522
ELA: 6.3067 Jacobi stage: ITER: 39475 RES: 0.0016808
ELA: 6.4068 Jacobi stage: ITER: 40107 RES: 0.0015253
ELA: 6.5069 Jacobi stage: ITER: 40739 RES: 0.0013842
ELA: 6.6069 Jacobi stage: ITER: 41371 RES: 0.0012562
ELA: 6.707 Jacobi stage: ITER: 42003 RES: 0.00114
ELA: 6.8071 Jacobi stage: ITER: 42635 RES: 0.0010345
ELA: 6.9072 Jacobi stage: ITER: 43265 RES: 0.00093879
ELA: 7.0073 Jacobi stage: ITER: 43895 RES: 0.00085247
ELA: 7.1074 Jacobi stage: ITER: 44529 RES: 0.00077313
ELA: 7.2075 Jacobi stage: ITER: 45161 RES: 0.0007016
ELA: 7.3076 Jacobi stage: ITER: 45797 RES: 0.0006363
ELA: 7.4077 Jacobi stage: ITER: 46433 RES: 0.00057708
ELA: 7.5078 Jacobi stage: ITER: 47071 RES: 0.00052337
ELA: 7.6079 Jacobi stage: ITER: 47707 RES: 0.00047466
ELA: 7.708 Jacobi stage: ITER: 48343 RES: 0.00043049
ELA: 7.8081 Jacobi stage: ITER: 48979 RES: 0.00039042
ELA: 7.9082 Jacobi stage: ITER: 49617 RES: 0.00035387
ELA: 8.0083 Jacobi stage: ITER: 50253 RES: 0.00032093
ELA: 8.1084 Jacobi stage: ITER: 50891 RES: 0.00029106
ELA: 8.2085 Jacobi stage: ITER: 51527 RES: 0.00026397
ELA: 8.3086 Jacobi stage: ITER: 52163 RES: 0.0002394
ELA: 8.4087 Jacobi stage: ITER: 52799 RES: 0.00021712
ELA: 8.5087 Jacobi stage: ITER: 53437 RES: 0.00019679
ELA: 8.6088 Jacobi stage: ITER: 54073 RES: 0.00017848
ELA: 8.7089 Jacobi stage: ITER: 54709 RES: 0.00016187
ELA: 8.809 Jacobi stage: ITER: 55345 RES: 0.0001468
ELA: 8.9091 Jacobi stage: ITER: 55977 RES: 0.00013322
ELA: 9.0092 Jacobi stage: ITER: 56607 RES: 0.00012097
ELA: 9.1093 Jacobi stage: ITER: 57239 RES: 0.00010978
ELA: 9.2094 Jacobi stage: ITER: 57871 RES: 9.9623e-05
ELA: 9.3095 Jacobi stage: ITER: 58503 RES: 9.0407e-05
ELA: 9.4096 Jacobi stage: ITER: 59135 RES: 8.2043e-05
ELA: 9.5097 Jacobi stage: ITER: 59767 RES: 7.4453e-05
ELA: 9.6098 Jacobi stage: ITER: 60399 RES: 6.7565e-05
ELA: 9.7099 Jacobi stage: ITER: 61033 RES: 6.1277e-05
ELA: 9.81 Jacobi stage: ITER: 61669 RES: 5.5574e-05
ELA: 9.9101 Jacobi stage: ITER: 62305 RES: 5.0401e-05
ELA: 10.01 Jacobi stage: ITER: 62943 RES: 4.5711e-05
ELA: 10.11 Jacobi stage: ITER: 63579 RES: 4.1456e-05
ELA: 10.21 Jacobi stage: ITER: 64215 RES: 3.7598e-05
ELA: 10.31 Jacobi stage: ITER: 64851 RES: 3.4099e-05
ELA: 10.411 Jacobi stage: ITER: 65487 RES: 3.0925e-05
ELA: 10.511 Jacobi stage: ITER: 66127 RES: 2.803e-05
ELA: 10.611 Jacobi stage: ITER: 66763 RES: 2.5421e-05
ELA: 10.711 Jacobi stage: ITER: 67399 RES: 2.3055e-05
ELA: 10.811 Jacobi stage: ITER: 68035 RES: 2.0909e-05
ELA: 10.911 Jacobi stage: ITER: 68673 RES: 1.8952e-05
ELA: 11.011 Jacobi stage: ITER: 69307 RES: 1.7198e-05
ELA: 11.111 Jacobi stage: ITER: 69943 RES: 1.5598e-05
ELA: 11.211 Jacobi stage: ITER: 70581 RES: 1.4137e-05
ELA: 11.311 Jacobi stage: ITER: 71217 RES: 1.2821e-05
ELA: 11.411 Jacobi stage: ITER: 71851 RES: 1.1635e-05
ELA: 11.512 Jacobi stage: ITER: 72483 RES: 1.0559e-05
ELA: 11.612 Jacobi stage: ITER: 73115 RES: 9.5821e-06
ELA: 11.712 Jacobi stage: ITER: 73747 RES: 8.6956e-06
ELA: 11.812 Jacobi stage: ITER: 74379 RES: 7.8911e-06
ELA: 11.912 Jacobi stage: ITER: 75011 RES: 7.1611e-06
ELA: 12.012 Jacobi stage: ITER: 75643 RES: 6.4986e-06
ELA: 12.112 Jacobi stage: ITER: 76279 RES: 5.8938e-06
ELA: 12.212 Jacobi stage: ITER: 76935 RES: 5.3289e-06
ELA: 12.312 Jacobi stage: ITER: 77893 RES: 4.5983e-06
ELA: 12.412 Jacobi stage: ITER: 78855 RES: 3.9678e-06
ELA: 12.513 Jacobi stage: ITER: 79815 RES: 3.4238e-06
ELA: 12.613 Jacobi stage: ITER: 80723 RES: 2.9781e-06
ELA: 12.713 Jacobi stage: ITER: 81683 RES: 2.5698e-06
ELA: 12.813 Jacobi stage: ITER: 82643 RES: 2.2175e-06
ELA: 12.913 Jacobi stage: ITER: 83603 RES: 1.9135e-06
ELA: 13.013 Jacobi stage: ITER: 84563 RES: 1.6511e-06
ELA: 13.113 Jacobi stage: ITER: 85523 RES: 1.4248e-06
ELA: 13.213 Jacobi stage: ITER: 86485 RES: 1.2287e-06
ELA: 13.313 Jacobi stage: ITER: 87467 RES: 1.057e-06
Vector x = [ 0.999999, 0.999999, 0.999998, 0.999999, 0.999999 ]

Member Function Documentation

◆ refresh()

template<typename Real, typename Index>
void TNL::Solvers::IterativeSolverMonitor< Real, Index >::refresh ( )
overridevirtual

Causes that the monitor prints out the status of the solver.

Implements TNL::Solvers::SolverMonitor.

◆ setIterations()

template<typename Real, typename Index>
void TNL::Solvers::IterativeSolverMonitor< Real, Index >::setIterations ( const IndexType & iterations)

Set number of the current iteration.

Parameters
iterationsis number of the current iteration.

◆ setNodesPerIteration()

template<typename Real, typename Index>
void TNL::Solvers::IterativeSolverMonitor< Real, Index >::setNodesPerIteration ( const IndexType & nodes)

Set the number of nodes of the numerical mesh or lattice.

This can be used to compute the number of nodes processed per one second.

Parameters
nodesis number of nodes of the numerical mesh or lattice.

◆ setResidue()

template<typename Real = double, typename Index = int>
void TNL::Solvers::IterativeSolverMonitor< Real, Index >::setResidue ( const RealType & residue)

Set residue of the current approximation of the solution.

Parameters
residueis a residue of the current approximation of the solution.

◆ setStage()

template<typename Real, typename Index>
void TNL::Solvers::IterativeSolverMonitor< Real, Index >::setStage ( const std::string & stage)

This method can be used for naming a stage of the monitored solver.

The stage name can be used to differ between various stages of iterative solvers.

Parameters
stageis name of the solver stage.

◆ setTime()

template<typename Real, typename Index>
void TNL::Solvers::IterativeSolverMonitor< Real, Index >::setTime ( const RealType & time)

Set the time of the simulated evolution if it is time dependent.

This can be used for example when solving parabolic or hyperbolic PDEs.

Parameters
timetime of the simulated evolution.

◆ setTimeStep()

template<typename Real, typename Index>
void TNL::Solvers::IterativeSolverMonitor< Real, Index >::setTimeStep ( const RealType & timeStep)

Set the time step for time dependent iterative solvers.

Parameters
timeSteptime step of the time dependent iterative solver.

◆ setVerbose()

template<typename Real, typename Index>
void TNL::Solvers::IterativeSolverMonitor< Real, Index >::setVerbose ( const IndexType & verbose)

Set up the verbosity of the monitor.

Parameters
verboseis the new value of the verbosity of the monitor.
The documentation for this class was generated from the following files:
Generated on Sun Mar 9 2025 21:30:32 for Template Numerical Library by doxygen 1.13.2