Changes between Version 8 and Version 9 of solvers


Ignore:
Timestamp:
10/29/20 11:28:06 (4 years ago)
Author:
tsantos
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • solvers

    v8 v9  
    66md.toolkits.DefaultAnalysis=bcgslbjacobioptions();
    77
    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.
    1313 
    1414md.toolkits=toolkits(); % to be sure that Mumps is the default solver for all analysis
     
    2020md.stressbalance.abstol=NaN; % no need
    2121
    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).
    2626}}}
    2727