MatSolve changed authored by Dmitry Bagaev's avatar Dmitry Bagaev
# Solve the Matrix stored in mtx format
The code for this example is located in `examples/MatSolve`
The code for this example is located in `Examples/MatSolve`
## Brief
......@@ -13,26 +13,24 @@ The code loads the matrix and right-hand side from files. If right-hand side is
## Arguments
```
Usage: ./MatSolve method_number<0:INNER_ILU2,1:INNER_MLILUC,2:PETSc,3:Trilinos_Aztec,4:Trilinos_Belos,5:Trilinos_Ifpack,6:Trilinos_ML,7:ANI> matrix.mtx [right_hand_side.rhs] [solver_options.txt]
Usage: ./MatSolve -m "matrix.mtx" [arguments]
```
- First parameter is the Solver type:
+ 0 – `INNER_ILU2`, inner Solver based on BiCGStab(L) solver with second order ILU factorization as preconditioner;
+ 1 – `INNER_MLILUC`, inner Solver based on BiCGStab(L) solver with second order Crout-ILU with inversed-based condition estimation and unsymmetric reordering for diagonal dominance as preconditioner;
+ 2 – `PETSc`, external Solver AztecOO from Trilinos package;
+ 3 – `Trilinos_Aztec`, external Solver Belos from Trilinos package, currently without preconditioner;
+ 4 – `Trilinos_Belos`, external Solver AztecOO with ML preconditioner;
+ 5 – `Trilinos_Ifpack`, external Solver AztecOO with Ifpack preconditioner;
+ 6 – `Trilinos_ML`, external Solver PETSc;
+ 7 – `ANI`, external Solver from ANI3D based on ILU2 (sequential Fortran version).
- Second parameter is the file with matrix in mtx format.
- Third optional parameter is the file with right-hand side in mtx format.
- Fourth optional parameter is the file with solver parameters, see `examples/MatSolve/database.txt` as example.
- Available arguments
* Required
+ -m, --matrix, specifies path to the matrix file, e.g. "/tmp/A.mtx"
* Optional
+ -b, --bvector, specifies path to the rhs file, e.g. "/tmp/b.rhs"
+ -x, --xvector, specifies path to the solution
+ -p, --parameters, specifies path to the database file, see `Examples/MatSolve/database.xml` as example;
+ -t, --type, specifies solver type, "inner_ilu2" by default
+ -h, --help, prints help message
## Running example
You can generate matrices using [[FVDiscr|1710-FVDiscr-Example]] example. The following line uses matrix `/tmp/A.mtx` and right-hand side `/tmp/b.rhs` from [[FVDiscr example|1710-FVDiscr-Example#running-example]].
```
$ cd examples/MatSolve
$ mpirun -np 4 ./MatSolve 0 /tmp/A.mtx /tmp/b.rhs
$ cd Examples/MatSolve
$ mpirun -np 4 ./MatSolve -m /tmp/A.mtx -b /tmp/b.rhs
load matrix: 0.346156
load vector: 0.0542159
preconditioner: 0.516358
......@@ -43,4 +41,3 @@ norms: 0.00276518
reason: converged due to relative tolerance
```