|
Last change
on this file was 26744, checked in by Mathieu Morlighem, 4 years ago |
|
merged trunk-jpl and trunk for revision 26742
|
|
File size:
891 bytes
|
| Line | |
|---|
| 1 | function asm=asmstokesoptions(varargin)
|
|---|
| 2 | %ASMSTOKESOPTIONS - return Additive Schwartz Method Stokes PETSc options
|
|---|
| 3 | %
|
|---|
| 4 | % Usage:
|
|---|
| 5 | % options=asmstokesoptions;
|
|---|
| 6 |
|
|---|
| 7 | %retrieve options provided in varargin
|
|---|
| 8 | options=pairoptions(varargin{:});
|
|---|
| 9 | asm=struct();
|
|---|
| 10 |
|
|---|
| 11 | %default asm options
|
|---|
| 12 | asm.toolkit='petsc';
|
|---|
| 13 | asm.mat_type=getfieldvalue(options,'mat_type','mpiaij');
|
|---|
| 14 | asm.ksp_type=getfieldvalue(options,'ksp_type','gmres');
|
|---|
| 15 | asm.pc_type=getfieldvalue(options,'pc_type','asm');
|
|---|
| 16 | asm.sub_pc_type=getfieldvalue(options,'sub_pc_type','lu');
|
|---|
| 17 | asm.pc_asm_overlap=getfieldvalue(options,'pc_asm_overlap',1); % COMSOL's default
|
|---|
| 18 | asm.ksp_max_it=getfieldvalue(options,'ksp_max_it',100);
|
|---|
| 19 | asm.ksp_rtol=getfieldvalue(options,'ksp_rtol',1e-7); %tuned for best performance and to fit ISMIP-HOM-C 5km with MUMPS
|
|---|
| 20 | asm.ksp_atol=getfieldvalue(options,'ksp_atol',1e-10); %tuned for best performance and to fit ISMIP-HOM-C 5km with MUMPS
|
|---|
Note:
See
TracBrowser
for help on using the repository browser.