8 | | Here is an example where bcg is used only for the stress balance analysis. |
9 | | The direct solver (MUMPS) is used for the rest (mass transport analysis, etc). |
10 | | Also, the recovery solver (MUMPS) is not called if bcg fails in a Picard iteration (linear solver). |
11 | | The only criterion is the mechanical equilibrium, which controls the non-linear solver as a whole. |
12 | | This scheme may compromise the monotonicity of the non-linear solver, but it is faster than calling MUMPS every time bgc fails. |
| 8 | % Here is an example where bcg is used only for the stress balance analysis. |
| 9 | % The direct solver (MUMPS) is used for the rest (mass transport analysis, etc). |
| 10 | % Also, the recovery solver (MUMPS) is not called if bcg fails in a Picard iteration (linear solver). |
| 11 | % The only criterion is the mechanical equilibrium, which controls the non-linear solver as a whole. |
| 12 | % This scheme may compromise the monotonicity of the non-linear solver, but it is faster than calling MUMPS every time bgc fails. |
22 | | Note: SSA is very stable, so bcg should convert fast. For HO(3D), bcg may not converge even after 50 iterations. |
23 | | In this case, look at the stress balance convergence (md.verbose=verbose('convergence',true,'solution',true);) |
24 | | and check if it is necessary to relax a little bit the restol value. |
25 | | The scheme above have being used for SSA (~250,000 elements) and for HO (~1,000,000 elements). |
| 22 | % Note: SSA is very stable, so bcg should convert fast. For HO(3D), bcg may not converge even after 50 iterations. |
| 23 | % In this case, look at the stress balance convergence (md.verbose=verbose('convergence',true,'solution',true);) |
| 24 | % and check if it is necessary to relax a little bit the restol value. |
| 25 | % The scheme above have being used for SSA (~250,000 elements) and for HO (~1,000,000 elements). |