Commit 51004f3f authored by Matthieu Constant's avatar Matthieu Constant
Browse files

inject

parent db7229d8
Pipeline #4617 passed with stage
in 22 seconds
......@@ -36,6 +36,7 @@ if not os.path.isdir(outputdir) :
t = 0
ii = 0
use_lmgc = True
def genInitialPosition(filename, N, r, lx, ly, rhop) :
......@@ -70,13 +71,15 @@ dt = .001 # time step
shutil.copy("mesh.msh", outputdir +"/mesh.msh")
genInitialPosition(outputdir, N, r, lx, ly, rhop)
friction=0.3
lmgc90Interface.scontactTolmgc90(outputdir, 2, 0, friction)
p = lmgc90Interface.ParticleProblem(2)
p.write_vtk(outputdir,0,0)
if use_lmgc:
friction=0.3
lmgc90Interface.scontactTolmgc90(outputdir, 2, 0, friction)
p = lmgc90Interface.ParticleProblem(2)
p.write_vtk(outputdir,0,0)
else:
p = scontact.ParticleProblem(2)
p.read_vtk(outputdir,0)
p = scontact.ParticleProblem(2)
p.read_vtk(outputdir,0)
# number of iterations between output files
outf = 25
ii = 0
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment