Context Navigation

asmstokesoptions.m

Last change on this file was 26744, checked in by Mathieu Morlighem, 3 years ago
merged trunk-jpl and trunk for revision 26742
File size: 891 bytes

Rev	Line
[26661]	1	function asm=asmstokesoptions(varargin)
	2	%ASMSTOKESOPTIONS - return Additive Schwartz Method Stokes PETSc options
[25935]	3	%
	4	% Usage:
[26661]	5	% options=asmstokesoptions;
[25935]	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);
[26607]	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.