Changeset 14610
- Timestamp:
- 04/16/13 16:12:55 (12 years ago)
- Location:
- issm/trunk-jpl/src/m/solvers
- Files:
-
- 16 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/m/solvers/asmoptions.m
r13646 r14610 10 10 11 11 %default asm options 12 asm.toolkit='petsc'; 12 13 asm.mat_type=getfieldvalue(options,'mat_type','aij'); 13 14 asm.ksp_type=getfieldvalue(options,'ksp_type','gmres'); -
issm/trunk-jpl/src/m/solvers/asmoptions.py
r12929 r14610 12 12 arguments=pairoptions(*args) 13 13 14 options=[[' mat_type','aij'],['ksp_type','gmres'],['pc_type','asm'],['sub_pc_type','lu'],['pc_asm_overlap',3],['ksp_max_it',100],['ksp_rtol',1e-30]];14 options=[['toolkit','petsc'],['mat_type','aij'],['ksp_type','gmres'],['pc_type','asm'],['sub_pc_type','lu'],['pc_asm_overlap',3],['ksp_max_it',100],['ksp_rtol',1e-30]]; 15 15 16 16 #now, go through our arguments, and write over default options. -
issm/trunk-jpl/src/m/solvers/iluasmoptions.m
r13646 r14610 10 10 11 11 %default iluasm options 12 iluasm.toolkit='petsc'; 12 13 iluasm.mat_type=getfieldvalue(options,'mat_type','aij'); 13 14 iluasm.ksp_type=getfieldvalue(options,'ksp_type','gmres'); -
issm/trunk-jpl/src/m/solvers/iluasmoptions.py
r13107 r14610 15 15 16 16 #default iluasm options 17 iluasm['toolkit']='petsc' 17 18 iluasm['mat_type']=options.getfieldvalue('mat_type','aij') 18 19 iluasm['ksp_type']=options.getfieldvalue('ksp_type','gmres') -
issm/trunk-jpl/src/m/solvers/jacobiasmoptions.m
r13646 r14610 10 10 11 11 %default jacobiasm options 12 jacobiasm.toolkit='petsc'; 12 13 jacobiasm.mat_type=getfieldvalue(options,'mat_type','aij'); 13 14 jacobiasm.ksp_type=getfieldvalue(options,'ksp_type','gmres'); -
issm/trunk-jpl/src/m/solvers/jacobiasmoptions.py
r12929 r14610 12 12 arguments=pairoptions(*args) 13 13 14 options=[[' mat_type','aij'],['ksp_type','gmres'],['pc_type','asm'],['sub_pc_type','jacobi'],['pc_asm_overlap',3],['ksp_max_it',100],['ksp_rtol',1e-15]];14 options=[['toolkit','petsc'],['mat_type','aij'],['ksp_type','gmres'],['pc_type','asm'],['sub_pc_type','jacobi'],['pc_asm_overlap',3],['ksp_max_it',100],['ksp_rtol',1e-15]]; 15 15 16 16 #now, go through our arguments, and write over default options. -
issm/trunk-jpl/src/m/solvers/jacobicgoptions.m
r13646 r14610 10 10 11 11 %default jacobiasm options 12 jacobicg.toolkit='petsc'; 12 13 jacobicg.mat_type=getfieldvalue(options,'mat_type','aij'); 13 14 jacobicg.ksp_type=getfieldvalue(options,'ksp_type','cg'); -
issm/trunk-jpl/src/m/solvers/jacobicgoptions.py
r12929 r14610 12 12 arguments=pairoptions(*args) 13 13 14 options=[[' mat_type','aij'],['ksp_type','cg'],['ksp_max_it',100],['ksp_rtol',1e-15]];14 options=[['toolkit','petsc'],['mat_type','aij'],['ksp_type','cg'],['ksp_max_it',100],['ksp_rtol',1e-15]]; 15 15 16 16 #now, go through our arguments, and write over default options. -
issm/trunk-jpl/src/m/solvers/matlaboptions.m
r13646 r14610 10 10 11 11 %default matlab options 12 maltab.toolkit='petsc'; 12 13 maltab.ksp_type='matlab'; -
issm/trunk-jpl/src/m/solvers/matlaboptions.py
r12929 r14610 12 12 arguments=pairoptions(*args) 13 13 14 options=[[' ksp_type','matlab']];14 options=[['toolkit','petsc'],['ksp_type','matlab']]; 15 15 16 16 #now, go through our arguments, and write over default options. -
issm/trunk-jpl/src/m/solvers/mumpsoptions.m
r11871 r14610 12 12 PETSC_VERSION=petscversion(); 13 13 if PETSC_VERSION==2, 14 mumps.toolkit='petsc'; 14 15 mumps.mat_type=getfieldvalue(options,'mat_type','aijmumps'); 15 16 mumps.ksp_type=getfieldvalue(options,'ksp_type','preonly'); … … 19 20 end 20 21 if PETSC_VERSION==3, 22 mumps.toolkit='petsc'; 21 23 mumps.mat_type=getfieldvalue(options,'mat_type','mpiaij'); 22 24 mumps.ksp_type=getfieldvalue(options,'ksp_type','preonly'); -
issm/trunk-jpl/src/m/solvers/mumpsoptions.py
r13107 r14610 18 18 PETSC_VERSION=petscversion() 19 19 if PETSC_VERSION==2: 20 mumps['toolkit']='petsc' 20 21 mumps['mat_type']=options.getfieldvalue('mat_type','aijmumps') 21 22 mumps['ksp_type']=options.getfieldvalue('ksp_type','preonly') … … 24 25 mumps['pc_factor_shift_positive_definite']=options.getfieldvalue('pc_factor_shift_positive_definite','true') 25 26 if PETSC_VERSION==3: 27 mumps['toolkit']='petsc' 26 28 mumps['mat_type']=options.getfieldvalue('mat_type','mpiaij') 27 29 mumps['ksp_type']=options.getfieldvalue('ksp_type','preonly') -
issm/trunk-jpl/src/m/solvers/soroptions.m
r13646 r14610 10 10 11 11 %default sor options 12 sor.toolkit='petsc'; 12 13 sor.mat_type=getfieldvalue(options,'mat_type','aij'); 13 14 sor.ksp_type=getfieldvalue(options,'ksp_type','cg'); -
issm/trunk-jpl/src/m/solvers/soroptions.py
r12929 r14610 12 12 arguments=pairoptions(*args) 13 13 14 options=[[' mat_type','aij'],['ksp_type','cg'],['pc_type','sor'],['pc_sor_omega',1.1],['pc_sor_its',2]];14 options=[['toolkit','petsc'],['mat_type','aij'],['ksp_type','cg'],['pc_type','sor'],['pc_sor_omega',1.1],['pc_sor_its',2]]; 15 15 16 16 #now, go through our arguments, and write over default options. -
issm/trunk-jpl/src/m/solvers/stokesoptions.m
r12645 r14610 15 15 end 16 16 if PETSC_VERSION==3, 17 stokes.toolkit='petsc'; 17 18 stokes.mat_type=getfieldvalue(options,'mat_type','mpiaij'); 18 19 stokes.ksp_max_it=getfieldvalue(options,'ksp_max_it',1000); -
issm/trunk-jpl/src/m/solvers/stokesoptions.py
r12929 r14610 19 19 raise RuntimeError('stokesoptions error message: multi-physics options not supported in Petsc 2') 20 20 if PETSC_VERSION==3: 21 options=[[' mat_type','mpiaij'],['ksp_max_it',1000],['ksp_type','gmres'],['pc_type','fieldsplit'],['pc_field_split_type','schur'],\21 options=[['toolkit','petsc'],['mat_type','mpiaij'],['ksp_max_it',1000],['ksp_type','gmres'],['pc_type','fieldsplit'],['pc_field_split_type','schur'],\ 22 22 ['fieldsplit_0_pc_type','hypre'],['fieldsplit_0_ksp_type','gmres'],['fieldsplit_0_pc_hypre_type','boomerang'],\ 23 23 ['fieldsplit_1_pc_type','jacobi'],['fieldsplit_1_ksp_type','preonly'],['issm_option_solver','stokes']]
Note:
See TracChangeset
for help on using the changeset viewer.