DAMASK_EICMD/python/damask/solver/_marc.py

75 lines
2.3 KiB
Python
Raw Normal View History

2019-12-13 03:28:37 +05:30
import subprocess
import shlex
2020-10-31 17:29:46 +05:30
import re
import io
import os
2020-06-03 17:02:47 +05:30
from pathlib import Path
2019-12-13 03:28:37 +05:30
2020-03-22 00:48:46 +05:30
class Marc:
"""Wrapper to run DAMASK with MSCMarc."""
2019-12-13 03:28:37 +05:30
def __init__(self,version=os.environ['MSC_VERSION']):
"""
Create a Marc solver object.
2020-03-22 00:47:50 +05:30
Parameters
----------
version : float
Marc version
"""
self.solver = 'Marc'
self.version = version
2019-12-13 03:28:37 +05:30
2020-06-03 17:02:47 +05:30
@property
def library_path(self):
2019-12-13 03:28:37 +05:30
path_lib = Path(f'{os.environ["MSC_ROOT"]}/mentat{self.version}/shlib/linux64')
2020-10-31 17:29:46 +05:30
if not path_lib.is_dir():
raise FileNotFoundError(f'library path "{path_lib}" not found')
return path_lib
2019-12-13 03:28:37 +05:30
2020-06-03 17:02:47 +05:30
@property
def tools_path(self):
2019-12-13 03:28:37 +05:30
path_tools = Path(f'{os.environ["MSC_ROOT"]}/marc{self.version}/tools')
2020-10-31 17:29:46 +05:30
if not path_tools.is_dir():
raise FileNotFoundError(f'tools path "{path_tools}" not found')
return path_tools
2019-12-13 03:28:37 +05:30
def submit_job(self, model, job,
compile = False,
optimization = ''):
2019-12-13 03:28:37 +05:30
usersub = Path(os.environ['DAMASK_ROOT'])/'src/DAMASK_Marc'
2020-07-16 13:03:17 +05:30
usersub = usersub.parent/(usersub.name + ('.f90' if compile else '.marc'))
2020-06-25 01:04:51 +05:30
if not usersub.is_file():
2020-10-31 17:29:46 +05:30
raise FileNotFoundError(f'subroutine ({"source" if compile else "binary"}) "{usersub}" not found')
2019-12-13 03:28:37 +05:30
# Define options [see Marc Installation and Operation Guide, pp 23]
2020-06-25 01:04:51 +05:30
script = f'run_damask_{optimization}mp'
2020-03-22 00:47:50 +05:30
2020-10-31 17:29:46 +05:30
cmd = str(self.tools_path/script) + \
' -jid ' + model+'_'+job + \
' -nprocd 1 -autorst 0 -ci n -cr n -dcoup 0 -b no -v no'
cmd += ' -u ' + str(usersub) + ' -save y' if compile else \
' -prog ' + str(usersub.with_suffix(''))
print(cmd)
2020-10-31 17:29:46 +05:30
ret = subprocess.run(shlex.split(cmd),capture_output=True)
2020-10-31 17:29:46 +05:30
try:
if 3004 != int(re.search('Exit number ([0-9]+)',ret.stderr.decode()).group(1)):
print(ret.stderr.decode())
print(ret.stdout.decode())
raise RuntimeError(f'Marc simulation failed ({v})')
2020-10-31 17:29:46 +05:30
except (AttributeError,ValueError):
print(ret.stderr.decode())
print(ret.stdout.decode())
2020-10-31 17:29:46 +05:30
raise RuntimeError('Marc simulation failed (unknown return value)')