2010-10-26 20:32:20 +05:30
|
|
|
#!/usr/bin/env python
|
|
|
|
|
2011-11-09 23:34:51 +05:30
|
|
|
# Makes postprocessing routines acessible from everywhere.
|
2010-10-26 20:32:20 +05:30
|
|
|
|
2011-12-15 14:23:18 +05:30
|
|
|
import os,sys,glob,string,damask
|
2011-12-06 22:28:17 +05:30
|
|
|
from optparse import OptionParser, Option
|
2010-10-26 20:32:20 +05:30
|
|
|
|
2011-11-09 21:50:52 +05:30
|
|
|
|
2011-12-06 22:28:17 +05:30
|
|
|
# -----------------------------
|
|
|
|
class extendableOption(Option):
|
|
|
|
# -----------------------------
|
|
|
|
# used for definition of new option parser action 'extend', which enables to take multiple option arguments
|
|
|
|
# taken from online tutorial http://docs.python.org/library/optparse.html
|
|
|
|
|
|
|
|
ACTIONS = Option.ACTIONS + ("extend",)
|
|
|
|
STORE_ACTIONS = Option.STORE_ACTIONS + ("extend",)
|
|
|
|
TYPED_ACTIONS = Option.TYPED_ACTIONS + ("extend",)
|
|
|
|
ALWAYS_TYPED_ACTIONS = Option.ALWAYS_TYPED_ACTIONS + ("extend",)
|
|
|
|
|
|
|
|
def take_action(self, action, dest, opt, value, values, parser):
|
|
|
|
if action == "extend":
|
|
|
|
lvalue = value.split(",")
|
|
|
|
values.ensure_value(dest, []).extend(lvalue)
|
|
|
|
else:
|
|
|
|
Option.take_action(self, action, dest, opt, value, values, parser)
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
########################################################
|
|
|
|
# MAIN
|
|
|
|
########################################################
|
|
|
|
|
|
|
|
parser = OptionParser(option_class=extendableOption, usage='%prog options', description = """
|
|
|
|
Sets up the pre and post processing tools of DAMASK
|
|
|
|
|
|
|
|
""" + string.replace('$Id: addDivergence.py 1129 2011-12-01 12:01:13Z MPIE\m.diehl $','\n','\\n')
|
|
|
|
)
|
|
|
|
|
|
|
|
parser.add_option('--F90', dest='compiler', type='string', \
|
|
|
|
help='name of F90 compiler [%default]')
|
|
|
|
|
|
|
|
parser.set_defaults(compiler = 'ifort')
|
|
|
|
(options,filenames) = parser.parse_args()
|
|
|
|
|
|
|
|
#translating name of compiler for use with f2py and setting subdirname of acml
|
|
|
|
if options.compiler == 'gfortran':
|
|
|
|
f2py_compiler='gnu95 --f90flags="-fno-range-check"'
|
|
|
|
else:
|
|
|
|
f2py_compiler='intelem'
|
2011-12-07 00:08:03 +05:30
|
|
|
|
|
|
|
acml_subdir='%s64/lib'%options.compiler
|
2011-12-06 22:28:17 +05:30
|
|
|
|
|
|
|
|
2011-12-15 14:23:18 +05:30
|
|
|
damaskEnv = damask.Environment()
|
|
|
|
baseDir = damaskEnv.relPath('processing/')
|
|
|
|
codeDir = damaskEnv.relPath('code/')
|
2011-12-06 22:28:17 +05:30
|
|
|
|
|
|
|
#define ToDo list
|
2011-01-26 20:47:03 +05:30
|
|
|
bin_link = { \
|
2011-11-09 23:21:00 +05:30
|
|
|
'pre' : [
|
|
|
|
'marc_addUserOutput.py',
|
|
|
|
'mentat_pbcOnBoxMesh.py',
|
|
|
|
'mentat_spectralBox.py',
|
|
|
|
'patchFromReconstructedBoundaries.py',
|
|
|
|
'spectral_geomCheck.py',
|
|
|
|
'spectral_minimalSurface.py',
|
|
|
|
'spectral_vicinityOffset.py',
|
2011-12-06 22:28:17 +05:30
|
|
|
'voronoi_randomSeeding.exe',
|
|
|
|
'voronoi_tessellation.exe',
|
2011-11-09 23:21:00 +05:30
|
|
|
],
|
|
|
|
'post' : [
|
|
|
|
'3Dvisualize.py',
|
|
|
|
'addCauchy.py',
|
2011-12-02 20:45:36 +05:30
|
|
|
'addCalculation.py',
|
2011-11-09 23:21:00 +05:30
|
|
|
'addDeterminant.py',
|
|
|
|
'addDivergence.py',
|
|
|
|
'addMises.py',
|
|
|
|
'addNorm.py',
|
|
|
|
'addStrainTensors.py',
|
|
|
|
'addCompatibilityMismatch.py',
|
|
|
|
'averageDown.py',
|
|
|
|
'mentat_colorMap.py',
|
|
|
|
'postResults.py',
|
|
|
|
'spectral_iterationCount.py',
|
|
|
|
'spectral_convergence.py',
|
|
|
|
'spectral_scaleGeom.py',
|
|
|
|
],
|
|
|
|
}
|
2010-10-26 20:32:20 +05:30
|
|
|
|
2011-01-26 20:47:03 +05:30
|
|
|
compile = { \
|
2011-11-09 23:21:00 +05:30
|
|
|
'pre' : [
|
|
|
|
'voronoi_randomSeeding.f90',
|
|
|
|
'voronoi_tessellation.f90',
|
|
|
|
],
|
|
|
|
'post' : [
|
|
|
|
],
|
|
|
|
}
|
2011-12-06 22:28:17 +05:30
|
|
|
|
2011-01-26 20:47:03 +05:30
|
|
|
|
2011-03-11 14:50:37 +05:30
|
|
|
execute = { \
|
2011-12-06 22:28:17 +05:30
|
|
|
'postMath' : [
|
2011-12-15 14:23:18 +05:30
|
|
|
'rm %s'%(os.path.join(damaskEnv.relPath('lib/'),'DAMASK.so')),
|
2011-12-06 22:28:17 +05:30
|
|
|
# The following command is used to compile math.f90 and make the functions defined in DAMASK_math.pyf
|
2011-12-07 00:08:03 +05:30
|
|
|
# available for python in the module DAMASK_math.so
|
2011-12-06 22:28:17 +05:30
|
|
|
# It uses the fortran wrapper f2py that is included in the numpy package to construct the
|
|
|
|
# module postprocessingMath.so out of the fortran code postprocessingMath.f90
|
|
|
|
# for the generation of the pyf file:
|
|
|
|
#f2py -m DAMASK -h DAMASK.pyf --overwrite-signature ../../code/math.f90 \
|
2011-12-07 00:08:03 +05:30
|
|
|
'f2py %s '%(os.path.join(codeDir,'DAMASK.pyf')) +\
|
2011-12-06 22:28:17 +05:30
|
|
|
'-c --fcompiler=%s '%(f2py_compiler) +\
|
2011-12-07 00:08:03 +05:30
|
|
|
'%s ' %(os.path.join(codeDir,'DAMASK2Python_helper.f90'))+\
|
|
|
|
'%s ' %(os.path.join(codeDir,'math.f90'))+\
|
2011-12-15 14:23:18 +05:30
|
|
|
'%s ' %(os.path.join(damaskEnv.pathInfo['fftw'],'libfftw3.a'))+\
|
|
|
|
'%s' %(os.path.join(damaskEnv.pathInfo['acml'],acml_subdir,'libacml.a')),
|
|
|
|
'mv %s %s' %(os.path.join(codeDir,'DAMASK.so'),damaskEnv.relPath('lib/')),
|
2011-12-06 22:28:17 +05:30
|
|
|
]
|
2011-11-09 23:21:00 +05:30
|
|
|
}
|
2011-03-11 14:50:37 +05:30
|
|
|
|
2011-11-09 21:07:45 +05:30
|
|
|
|
2011-12-06 22:28:17 +05:30
|
|
|
for dir in compile:
|
|
|
|
for file in compile[dir]:
|
|
|
|
src = os.path.abspath(os.path.join(baseDir,dir,file))
|
|
|
|
if os.path.isfile(src):
|
|
|
|
try:
|
|
|
|
os.system('rm %s.exe'%(os.path.splitext(src)[0]))
|
|
|
|
print 'removing %s.exe '%(os.path.splitext(src)[0])
|
|
|
|
except:
|
|
|
|
pass
|
|
|
|
print 'compiling ',src,'using',options.compiler
|
|
|
|
os.system('%s -O2 -o%s.exe %s'%(options.compiler,os.path.splitext(src)[0],src))
|
|
|
|
|
|
|
|
os.chdir(codeDir) # needed for compilation with gfortran and f2py
|
|
|
|
for tasks in execute:
|
|
|
|
for cmd in execute[tasks]:
|
|
|
|
os.system(cmd)
|
2011-12-15 14:23:18 +05:30
|
|
|
os.chdir(damaskEnv.relPath('processing/setup/'))
|
2011-12-06 22:28:17 +05:30
|
|
|
|
|
|
|
modules = glob.glob('*.mod')
|
|
|
|
for module in modules:
|
|
|
|
print 'removing', module
|
|
|
|
os.remove(module)
|
|
|
|
|
2011-01-26 20:47:03 +05:30
|
|
|
for dir in bin_link:
|
|
|
|
for file in bin_link[dir]:
|
2011-12-06 22:28:17 +05:30
|
|
|
src = os.path.abspath(os.path.join(baseDir,dir,file))
|
2010-10-26 20:32:20 +05:30
|
|
|
if (file == ''):
|
2011-12-15 14:23:18 +05:30
|
|
|
sym_link = os.path.abspath(os.path.join(damaskEnv.binDir(),dir))
|
2010-10-26 20:32:20 +05:30
|
|
|
else:
|
2011-12-15 14:23:18 +05:30
|
|
|
sym_link = os.path.abspath(os.path.join(damaskEnv.binDir(),os.path.splitext(file)[0]))
|
2011-11-09 21:50:52 +05:30
|
|
|
print sym_link,'-->',src
|
|
|
|
if os.path.lexists(sym_link):
|
|
|
|
os.remove(sym_link)
|
|
|
|
os.symlink(src,sym_link)
|
|
|
|
|