Skip to content
Snippets Groups Projects
Commit c57673c9 authored by Timo Koch's avatar Timo Koch
Browse files

Merge branch 'cherry-pick-bec8d3f4' into 'releases/3.0'

Merge branch 'fix/facet-1p-convtest-parallel' into 'master'

See merge request !1447
parents e66c2a35 ed1b3872
No related branches found
No related tags found
1 merge request!1447Merge branch 'fix/facet-1p-convtest-parallel' into 'master'
...@@ -5,10 +5,10 @@ import os.path ...@@ -5,10 +5,10 @@ import os.path
import subprocess import subprocess
import sys import sys
if len(sys.argv) > 5 or len(sys.argv) < 3: if len(sys.argv) < 3:
sys.stderr.write("Invalid number of arguments given. Please provide the following arguments (in this order):\n\ sys.stderr.write("Invalid number of arguments given. Please provide the following arguments (in this order):\n\
- the name of the executable (either test_facetcoupling_tpfa_1p1p or test_facetcoupling_box_1p1p)\n\ - the name of the executable (either test_facetcoupling_tpfa_1p1p or test_facetcoupling_box_1p1p)\n\
- between one and three fracture permeabilities you want to consider") - the fracture permeabilities you want to perform the convergence test for")
sys.exit(1) sys.exit(1)
# name of the executable # name of the executable
...@@ -41,7 +41,7 @@ for permIndex in range(0, len(k)): ...@@ -41,7 +41,7 @@ for permIndex in range(0, len(k)):
sys.exit(1) sys.exit(1)
# write a temporary geo-file using num elements # write a temporary geo-file using num elements
tmpGeoFile = open("grids/tmp.geo", 'w') tmpGeoFile = open("grids/" + execName + ".geo", 'w')
lines[0] = "numElemsPerSide = " + str(int(cells)) + ";\n" lines[0] = "numElemsPerSide = " + str(int(cells)) + ";\n"
for line in lines: for line in lines:
tmpGeoFile.write(line) tmpGeoFile.write(line)
...@@ -49,17 +49,17 @@ for permIndex in range(0, len(k)): ...@@ -49,17 +49,17 @@ for permIndex in range(0, len(k)):
geoFile.close() geoFile.close()
tmpGeoFile.close() tmpGeoFile.close()
subprocess.call(['gmsh', '-2', 'grids/tmp.geo']) subprocess.call(['gmsh', '-2', "grids/" + execName + ".geo"])
subprocess.call(['./' + execName, 'params.input', subprocess.call(['./' + execName, 'params.input',
'-Vtk.OutputName', execName, '-Vtk.OutputName', execName,
'-Grid.File', 'grids/tmp.msh', '-Grid.File', "grids/" + execName + ".msh",
'-Grid.NumElemsPerSide', str(int(cells)), '-Grid.NumElemsPerSide', str(int(cells)),
'-LowDim.SpatialParams.Permeability', str(k[permIndex]), '-LowDim.SpatialParams.Permeability', str(k[permIndex]),
'-Problem.OutputFileName', execName + '.log']) '-Problem.OutputFileName', execName + '.log'])
# remove geo and msh file # remove geo and msh file
subprocess.call(['rm', 'grids/tmp.geo']) subprocess.call(['rm', "grids/" + execName + ".geo"])
subprocess.call(['rm', 'grids/tmp.msh']) subprocess.call(['rm', "grids/" + execName + ".msh"])
# check the rates and append them to the log file # check the rates and append them to the log file
logfile = open(execName + '.log', "r+") logfile = open(execName + '.log', "r+")
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment