This shows you the differences between two versions of the page.
| Both sides previous revision Previous revision Next revision | Previous revision | ||
|
appendices:a19 [2019/06/27 17:45] helene |
appendices:a19 [2020/08/25 15:46] (current) |
||
|---|---|---|---|
| Line 31: | Line 31: | ||
| ===== Guideline to using GMRES method ===== | ===== Guideline to using GMRES method ===== | ||
| - | KNSYM must be equal to 6 at the 1st line and the 15th column in the loading file. Then the parameters for the iterative method must be entered on the second line with the following format: | + | KNSYM must be equal to 6 at the 1st line and the 15th column in the [[:lagamex|loading file]]. Then the parameters for the iterative method must be entered on the second line (this line is added in the loading file) with the following format: |
| ^Format^Name^Description^Use^Default value^ | ^Format^Name^Description^Use^Default value^ | ||
| | I5 | IM |size of the Krylov subspace (must be inferior to 1000)|GMRES parameter|20| | | I5 | IM |size of the Krylov subspace (must be inferior to 1000)|GMRES parameter|20| | ||
| | I5 | MAXITS |Maximum of iterations allowed |GMRES parameters|5000| | | I5 | MAXITS |Maximum of iterations allowed |GMRES parameters|5000| | ||
| + | | I5 | LFIL | Fill-in tolerance parameter |Preconditioning parameters|50| | ||
| + | | I5 | IOUT | Output unit number for printing intermediate results (=1 printing, =0 no printing)|GMRES parameters| 0| | ||
| + | | G10.0 | EPS | Tolerance for stopping criterion |GMRES parameter|10<sup>-4</sup>| | ||
| + | | G10.0 | DROPTOL |Sets the threshold for dropping small terms in the factorization|Preconditioning parameters| 10<sup>-8</sup>| | ||
| + | | G10.0 | RESTOL |Sets the threshold value of the residual which consider the convergence|GMRES parameters|10<sup>-8</sup>| | ||
| + | ===== Example of parameters as function of number of DOF ===== | ||
| + | === Example 1 === | ||
| + | ^DOF^IM^EPS^LFIL^DROPTOL^ | ||
| + | |26392 |2000 |1.10<sup>-4</sup>| 50| 1.10<sup>-6</sup>| | ||
| + | |14637| 500 |1.10<sup>-4</sup> |100 |5.10<sup>-3</sup>| | ||
| + | |||
| + | === Example 2 === | ||
| + | ^DOF^IM^MAXITS^LFIL^EPS^DROPTOL^ | ||
| + | |42700|20|5000|15|10<sup>-4</sup>|10<sup>-2</sup>| | ||
| + | |58483|175|5000|200|10<sup>-5</sup>|10<sup>-6</sup>| | ||
| + | |||
| + | ===== Effect of the parameters on the convergence ===== | ||
| + | |||
| + | * Higher are IM and LFIL, lower the number of iterations before convergence is. However the memory cost is more important. | ||
| + | * Higher EPS is, lower the number of iterations before convergence for the iterative solver is. However if the residual is not enough low, the out-of-balanced forces diminish slower and the Newton-Raphson can converge slower. | ||
| + | * Higher DROTOL is, lower the memory cost of the preconditioning matrix is. However the number of iterations before convergence for the iterative solver increases. | ||