[Trilinos-Users] How can I print every iteration details in Belos solver?
Thornquist, Heidi K
hkthorn at sandia.gov
Thu Jun 26 11:01:01 MDT 2008
John,
That is true, Belos does not have AztecOO-like output because Belos uses an inherently different approach for convergence checking. In the future, we hope to have some more streamlined output features that look cleaner, but for now this is how you get the residual information for each iteration.
Thanks,
Heidi
On 6/26/08 10:49 AM, "john" <ax.equal.b at gmail.com> wrote:
Dear Heidi,
I have tried this option,but the output is not so clearly as in AztecOO.
Just as:
//---------------------------------------------------------
Belos::StatusTestOutput: Failed
(Num calls,Mod test,State test): (1, 1, Passed Failed Undefined)
Failed.......OR Combination ->
OK...........Number of Iterations = 0 < 100
Unconverged..(2-Norm Imp Res Vec) / (2-Norm Res0)
residual [ 0 ] = 1 > 1e-08
Belos::StatusTestOutput: Failed
(Num calls,Mod test,State test): (2, 1, Passed Failed Undefined)
Failed.......OR Combination ->
OK...........Number of Iterations = 1 < 100
Unconverged..(2-Norm Imp Res Vec) / (2-Norm Res0)
residual [ 0 ] = 0.287091 > 1e-08
Belos::StatusTestOutput: Failed
(Num calls,Mod test,State test): (3, 1, Passed Failed Undefined)
Failed.......OR Combination ->
OK...........Number of Iterations = 2 < 100
Unconverged..(2-Norm Imp Res Vec) / (2-Norm Res0)
residual [ 0 ] = 0.175683 > 1e-08
Belos::StatusTestOutput: Failed
(Num calls,Mod test,State test): (4, 1, Passed Failed Undefined)
Failed.......OR Combination ->
OK...........Number of Iterations = 3 < 100
Unconverged..(2-Norm Imp Res Vec) / (2-Norm Res0)
residual [ 0 ] = 0.0747333 > 1e-08
Belos::StatusTestOutput: Failed
(Num calls,Mod test,State test): (5, 1, Passed Failed Undefined)
Failed.......OR Combination ->
OK...........Number of Iterations = 4 < 100
Unconverged..(2-Norm Imp Res Vec) / (2-Norm Res0)
residual [ 0 ] = 0.0334115 > 1e-08
Belos::StatusTestOutput: Passed
(Num calls,Mod test,State test): (6, 1, Passed Failed Undefined)
Passed.......OR Combination ->
OK...........Number of Iterations = 5 < 100
Converged....(2-Norm Imp Res Vec) / (2-Norm Res0)
residual [ 0 ] = 1.55101e-17 < 1e-08
//------------------------------------
How can I remove the unrelated output?
"Belos::IterationDetails" is used for what?
Thanks
2008/6/26 Thornquist, Heidi K <hkthorn at sandia.gov>:
Hi,
To get the residual information for every iteration regardless of convergence, you need to add "+Belos::StatusTestDetails" to your verbosity.
Thanks,
Heidi
On 6/26/08 6:36 AM, "刘俊" <ax.equal.b at gmail.com> wrote:
Hi,
I am now testing ex1.cpp of Belos, but I can't print out every iteration details,
I have set the verbosity according to the manual,but it does not work. It print only the final results.
some code as follow:
..........
int verbosity = Belos::Errors+Belos::Warnings+Belos::Debug+Belos::IterationDetails;
Teuchos::RCP<Teuchos::ParameterList> myPL = Teuchos::rcp( new Teuchos::ParameterList() );
myPL->set( "Verbosity", verbosity );
myPL->set( "Block Size", blocksize );
myPL->set( "Maximum Iterations", 100 );
myPL->set( "Convergence Tolerance", 1.0e-8 );
myPL->set("Output Frequency",1);
Belos::BlockCGSolMgr<double,MV,OP> mySolver(myProblem, myPL);//CG方法
....................
the output:
--------------------------
Belos Example: Block CG
Problem info:
* Solving problem laplace_2d
* Number of global elements : -1
* Type of Map : linear
* Number of PDEs : 1
Belos::StatusTestOutput: Passed
(Num calls,Mod test,State test): (30, 1, Passed)
Passed.......OR Combination ->
OK...........Number of Iterations = 29 < 100
Converged....(2-Norm Imp Res Vec) / (2-Norm Res0)
residual [ 0 ] = 5.43159e-09 < 1e-08
residual [ 1 ] = 1.27115e-08 > 1e-08
Belos::StatusTestOutput: Passed
(Num calls,Mod test,State test): (32, 1, Passed)
Passed.......OR Combination ->
OK...........Number of Iterations = 30 < 100
Converged....(2-Norm Imp Res Vec) / (2-Norm Res0)
residual [ 0 ] = 5.43159e-09 < 1e-08
residual [ 1 ] = 4.938e-09 < 1e-08
Belos::BlockCGSolMgr::solve() converged!
////////////
I want to print the information about iteration as in AztecOO solver,such as:
*******************************************************
***** Problem: Epetra::CrsMatrix
***** Preconditioned CG solution
***** No preconditioning
***** No scaling
*******************************************************
iter: 0 residual = 1.000000e+00
iter: 1 residual = 2.990580e-01
iter: 2 residual = 1.203293e-01
iter: 3 residual = 7.955944e-02
iter: 4 residual = 4.113928e-02
iter: 5 residual = 2.821794e-02
iter: 6 residual = 3.141783e-02
iter: 7 residual = 3.228408e-02
iter: 8 residual = 2.419115e-02
iter: 9 residual = 1.694050e-02
iter: 10 residual = 1.324365e-02
iter: 11 residual = 7.659560e-03
iter: 12 residual = 4.284702e-03
iter: 13 residual = 2.637796e-03
iter: 14 residual = 1.648559e-03
iter: 15 residual = 1.207106e-03
iter: 16 residual = 9.482598e-04
iter: 17 residual = 6.114811e-04
iter: 18 residual = 3.355039e-04
iter: 19 residual = 2.229069e-04
iter: 20 residual = 9.825301e-05
iter: 21 residual = 6.336427e-05
iter: 22 residual = 3.948687e-05
iter: 23 residual = 1.821252e-05
iter: 24 residual = 1.239958e-05
iter: 25 residual = 7.315337e-06
iter: 26 residual = 3.092389e-06
iter: 27 residual = 1.322761e-06
iter: 28 residual = 5.867232e-07
iter: 29 residual = 2.594566e-07
iter: 30 residual = 8.810022e-08
iter: 31 residual = 3.242696e-08
iter: 32 residual = 9.748573e-09
Solution time: 0.000000 (sec.)
total iterations: 32
Solver performed 32iterations.
Norm of the true residual = 2.43177e-07
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://software.sandia.gov/mailman/private/trilinos-users/attachments/20080626/952460b9/attachment.html
More information about the Trilinos-Users
mailing list