Changes between Version 9 and Version 10 of solvers

Timestamp:

10/29/20 11:31:04 (4 years ago)

Author:

tsantos

Comment:

--

solvers

@@ -5 +5 @@
 % biconjugate gradient with block Jacobi preconditioner [Thiago's favorite]
 md.toolkits.DefaultAnalysis=bcgslbjacobioptions();
+}}}
 
-% Here is an example where bcg is used only for the stress balance analysis. 
-% The direct solver (MUMPS) is used for the rest (mass transport analysis, etc). 
-% Also, the recovery solver (MUMPS) is not called if bcg fails in a Picard iteration (linear solver). 
-% The only criterion is the mechanical equilibrium, which controls the non-linear solver as a whole. 
-% This scheme may compromise the monotonicity of the non-linear solver, but it is faster than calling MUMPS every time bgc fails.
- 
+Here is an example where bcg is used only for the stress balance analysis. The direct solver (MUMPS) is used for the rest (mass transport analysis, etc). Also, the recovery solver (MUMPS) is not called if bcg fails in a Picard iteration (linear solver). The only criterion is the mechanical equilibrium, which controls the non-linear solver as a whole. This scheme may compromise the monotonicity of the non-linear solver, but it is faster than calling MUMPS every time bgc fails.
+
+{{{
+#!m
 md.toolkits=toolkits(); % to be sure that Mumps is the default solver for all analysis
 md.toolkits=addoptions(md.toolkits,'StressbalanceAnalysis',bcgslbjacobioptions()); % add bcg for stress balance analysis only
@@ -19 +18 @@
 md.stressbalance.reltol=NaN; % no need
 md.stressbalance.abstol=NaN; % no need
-
-% Note: SSA is very stable, so bcg should convert fast. For HO(3D), bcg may not converge even after 50 iterations. 
-% In this case, look at the stress balance convergence (md.verbose=verbose('convergence',true,'solution',true);) 
-% and check if it is necessary to relax a little bit the restol value. 
-% The scheme above have being used for SSA (~250,000 elements) and for HO (~1,000,000 elements).
 }}}
+Note: SSA is very stable, so bcg should convert fast. For HO(3D), bcg may not converge even after 50 iterations. 
+In this case, look at the stress balance convergence (md.verbose=verbose('convergence',true,'solution',true);) 
+and check if it is necessary to relax a little bit the restol value. 
+The scheme above have being used for SSA (~250,000 elements) and for HO (~1,000,000 elements).
 
 {{{