Changeset 15907
- Timestamp:
- 08/23/13 16:28:47 (12 years ago)
- Location:
- issm/trunk-jpl/src/m
- Files:
-
- 10 deleted
- 18 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/m/classes/clusters/generic.m
r15121 r15907 70 70 if ~isvalgrind, 71 71 if cluster.interactive 72 if ismpi,72 if IssmConfig('_HAVE_MPI_'), 73 73 fprintf(fid,'mpiexec -np %i %s/issm.exe %s %s %s ',cluster.np,cluster.codepath,EnumToString(solution),[cluster.executionpath '/' dirname],modelname); 74 74 else … … 76 76 end 77 77 else 78 if ismpi,78 if IssmConfig('_HAVE_MPI_'), 79 79 fprintf(fid,'mpiexec -np %i %s/issm.exe %s %s %s 2> %s.errlog >%s.outlog ',cluster.np,cluster.codepath,EnumToString(solution),[cluster.executionpath '/' dirname],modelname,modelname,modelname); 80 80 else … … 88 88 fprintf(fid,'LD_PRELOAD=%s \\\n',cluster.valgrindlib); 89 89 if ismac, 90 if ismpi,90 if IssmConfig('_HAVE_MPI_'), 91 91 fprintf(fid,'mpiexec -np %i %s --leak-check=full --dsymutil=yes --suppressions=%s %s/issm.exe %s %s %s 2> %s.errlog >%s.outlog ',... 92 92 cluster.np,cluster.valgrind,cluster.valgrindsup,cluster.codepath,EnumToString(solution),[cluster.executionpath '/' dirname], modelname,modelname,modelname); … … 96 96 end 97 97 else 98 if ismpi,98 if IssmConfig('_HAVE_MPI_'), 99 99 fprintf(fid,'mpiexec -np %i %s --leak-check=full --suppressions=%s %s/issm.exe %s %s %s 2> %s.errlog >%s.outlog ',... 100 100 cluster.np,cluster.valgrind,cluster.valgrindsup,cluster.codepath,EnumToString(solution),[cluster.executionpath '/' dirname],modelname,modelname,modelname); -
issm/trunk-jpl/src/m/classes/clusters/generic.py
r14477 r15907 1 1 import socket 2 import os3 2 import math 4 3 import subprocess 4 from IssmConfig import IssmConfig 5 5 from EnumToString import EnumToString 6 6 from issmdir import * … … 81 81 if not isvalgrind: 82 82 if self.interactive: 83 if ismpi():83 if IssmConfig('_HAVE_MPI_')[0]: 84 84 fid.write('mpiexec -np %i %s/issm.exe %s %s/%s %s ' % (self.np,self.codepath,EnumToString(solution)[0],self.executionpath,dirname,modelname)) 85 85 else: 86 86 fid.write('%s/issm.exe %s %s/%s %s ' % (self.codepath,EnumToString(solution)[0],self.executionpath,dirname,modelname)) 87 87 else: 88 if ismpi():88 if IssmConfig('_HAVE_MPI_')[0]: 89 89 fid.write('mpiexec -np %i %s/issm.exe %s %s/%s %s 2> %s.errlog >%s.outlog ' % (self.np,self.codepath,EnumToString(solution)[0],self.executionpath,dirname,modelname,modelname,modelname)) 90 90 else: … … 95 95 #Add --gen-suppressions=all to get suppression lines 96 96 fid.write('LD_PRELOAD=%s \\\n' % self.valgrindlib) 97 if ismpi():97 if IssmConfig('_HAVE_MPI_')[0]: 98 98 fid.write('mpiexec -np %i %s --leak-check=full --suppressions=%s %s/issm.exe %s %s/%s %s 2> %s.errlog >%s.outlog ' % \ 99 99 (self.np,self.valgrind,self.valgrindsup,self.codepath,EnumToString(solution)[0],self.executionpath,dirname,modelname,modelname,modelname)) -
issm/trunk-jpl/src/m/classes/model/model.py
r15771 r15907 38 38 from private import private 39 39 from EnumDefinitions import * 40 from ismumps import *41 40 from mumpsoptions import * 42 41 from iluasmoptions import * -
issm/trunk-jpl/src/m/classes/model/planet.m
r15771 r15907 38 38 md.settings = settings(); 39 39 md.solver = solver(); 40 if ismumps(),41 md.solver = addoptions(md.solver,StressbalanceVerticalAnalysisEnum(),mumpsoptions());42 else43 md.solver = addoptions(md.solver,StressbalanceVerticalAnalysisEnum(),iluasmoptions());44 end45 40 md.cluster = generic(); 46 41 md.balancethickness = balancethickness(); -
issm/trunk-jpl/src/m/classes/qmu/@dakota_method/dmeth_params_write.m
r15771 r15907 350 350 param_write(fid,sbeg,'seed',' = ','\n',dm.params); 351 351 param_write(fid,sbeg,'fixed_seed','','\n',dm.params); 352 dver=textscan( dakotaversion(),'%[0123456789].%[0123456789].%[0123456789]');352 dver=textscan(IssmConfig('_DAKOTA_VERSION_'),'%[0123456789].%[0123456789].%[0123456789]'); 353 353 if ((str2num(dver{1}{1})==4 && str2num(dver{2}{1})>2) || str2num(dver{1}{1})>4) 354 354 param_write(fid,sbeg,'rng',' ','\n',dm.params); -
issm/trunk-jpl/src/m/classes/toolkits.m
r15771 r15907 40 40 41 41 %default toolkits: 42 if ispetsc,42 if IssmConfig('_HAVE_PETSC_'), 43 43 %MUMPS is the default toolkits 44 if ismumps(),44 if IssmConfig('_HAVE_MUMPS_'), 45 45 obj.DefaultAnalysis = mumpsoptions(); 46 46 else -
issm/trunk-jpl/src/m/classes/toolkits.py
r15771 r15907 1 from ismumps import * 2 from ispetsc import * 1 from IssmConfig import IssmConfig 3 2 from mumpsoptions import * 4 3 from iluasmoptions import * … … 19 18 def __init__(self): # {{{ 20 19 #default toolkits 21 if ispetsc():20 if IssmConfig('_HAVE_PETSC_')[0]: 22 21 #MUMPS is the default toolkits 23 if ismumps():22 if IssmConfig('_HAVE_MUMPS_')[0]: 24 23 self.DefaultAnalysis = mumpsoptions() 25 24 else: -
issm/trunk-jpl/src/m/solve/solve.py
r15771 r15907 8 8 from ismodelselfconsistent import * 9 9 from marshall import * 10 from ispetsc import *11 10 from waitonlock import * 12 11 from loadresultsfromcluster import * -
issm/trunk-jpl/src/m/solvers/asmoptions.py
r14610 r15907 1 #module imports {{{2 1 from pairoptions import * 3 from petscversion import *4 #}}}5 2 def asmoptions(*args): 6 3 #ASMOPTIONS - return ASM petsc options -
issm/trunk-jpl/src/m/solvers/issmsolver.py
r14728 r15907 1 #module imports {{{2 1 from pairoptions import * 3 from petscversion import *4 #}}}5 2 def asmoptions(*args): 6 3 #ISSMSOLVE - return issm solver options -
issm/trunk-jpl/src/m/solvers/jacobiasmoptions.py
r14610 r15907 1 #module imports {{{2 1 from pairoptions import * 3 from petscversion import *4 #}}}5 2 def jacobiasmoptions(*args): 6 3 #ASMOPTIONS - return Additive Shwartz Method with Jacobi preconditioner petsc options -
issm/trunk-jpl/src/m/solvers/jacobicgoptions.py
r14610 r15907 1 #module imports {{{2 1 from pairoptions import * 3 from petscversion import *4 #}}}5 2 def jacobicgoptions(*args): 6 3 #ASMOPTIONS - return Additive Shwartz Method with Jacobi preconditioner petsc options -
issm/trunk-jpl/src/m/solvers/matlaboptions.py
r14610 r15907 1 #module imports {{{2 1 from pairoptions import * 3 from petscversion import *4 #}}}5 2 def matlaboptions(*args): 6 3 #MATLABOPTIONS - return Matlab petsc options -
issm/trunk-jpl/src/m/solvers/mumpsoptions.m
r14610 r15907 10 10 11 11 %default mumps options 12 PETSC_VERSION= petscversion();13 if PETSC_VERSION==2 ,12 PETSC_VERSION=IssmConfig('_PETSC_MAJOR_')[0]; 13 if PETSC_VERSION==2., 14 14 mumps.toolkit='petsc'; 15 15 mumps.mat_type=getfieldvalue(options,'mat_type','aijmumps'); … … 19 19 mumps.pc_factor_shift_positive_definite=getfieldvalue(options,'pc_factor_shift_positive_definite','true'); 20 20 end 21 if PETSC_VERSION==3 ,21 if PETSC_VERSION==3., 22 22 mumps.toolkit='petsc'; 23 23 mumps.mat_type=getfieldvalue(options,'mat_type','mpiaij'); -
issm/trunk-jpl/src/m/solvers/mumpsoptions.py
r14610 r15907 1 1 from collections import OrderedDict 2 2 from pairoptions import * 3 from petscversionimport *3 from IssmConfig import * 4 4 5 5 def mumpsoptions(*args): … … 16 16 17 17 #default mumps options 18 PETSC_VERSION= petscversion()19 if PETSC_VERSION==2 :18 PETSC_VERSION=IssmConfig('_PETSC_MAJOR_')[0] 19 if PETSC_VERSION==2.: 20 20 mumps['toolkit']='petsc' 21 21 mumps['mat_type']=options.getfieldvalue('mat_type','aijmumps') … … 24 24 mumps['mat_mumps_icntl_14']=options.getfieldvalue('mat_mumps_icntl_14',120) 25 25 mumps['pc_factor_shift_positive_definite']=options.getfieldvalue('pc_factor_shift_positive_definite','true') 26 if PETSC_VERSION==3 :26 if PETSC_VERSION==3.: 27 27 mumps['toolkit']='petsc' 28 28 mumps['mat_type']=options.getfieldvalue('mat_type','mpiaij') -
issm/trunk-jpl/src/m/solvers/soroptions.py
r14610 r15907 1 #module imports {{{2 1 from pairoptions import * 3 from petscversion import *4 #}}}5 2 def soroptions(*args): 6 3 #SOROPTIONS - return Relaxation Solver petsc options -
issm/trunk-jpl/src/m/solvers/stokesoptions.m
r15568 r15907 10 10 11 11 %default stokes options 12 PETSC_VERSION= petscversion();13 if PETSC_VERSION==2 ,12 PETSC_VERSION=IssmConfig('_PETSC_MAJOR_')[0]; 13 if PETSC_VERSION==2., 14 14 error('stokesoptions error message: multi-physics options not supported in Petsc 2'); 15 15 end 16 if PETSC_VERSION==3 ,16 if PETSC_VERSION==3., 17 17 stokes.toolkit='petsc'; 18 18 stokes.mat_type=getfieldvalue(options,'mat_type','mpiaij'); -
issm/trunk-jpl/src/m/solvers/stokesoptions.py
r15568 r15907 1 #module imports {{{2 1 from pairoptions import * 3 from petscversion import * 4 #}}} 2 from IssmConfig import * 5 3 def stokesoptions(*args): 6 4 #STOKESOPTIONS - return STOKES multi-physics solver petsc options … … 14 12 15 13 #default stokes options 16 PETSC_VERSION= petscversion()14 PETSC_VERSION=IssmConfig('_PETSC_MAJOR_')[0] 17 15 18 if PETSC_VERSION==2 :16 if PETSC_VERSION==2.: 19 17 raise RuntimeError('stokesoptions error message: multi-physics options not supported in Petsc 2') 20 if PETSC_VERSION==3 :18 if PETSC_VERSION==3.: 21 19 options=[['toolkit','petsc'],['mat_type','mpiaij'],['ksp_max_it',1000],['ksp_type','gmres'],['pc_type','fieldsplit'],['pc_field_split_type','schur'],\ 22 20 ['fieldsplit_0_pc_type','hypre'],['fieldsplit_0_ksp_type','gmres'],['fieldsplit_0_pc_hypre_type','boomerang'],\
Note:
See TracChangeset
for help on using the changeset viewer.