Skip to content
Snippets Groups Projects
Commit 7692965d authored by Andreas Dedner's avatar Andreas Dedner
Browse files

set up things for multithreading

parent b1e31625
No related branches found
No related tags found
1 merge request!4Latest features added to dune-fem-dg.
import mpi4py.rc
mpi4py.rc.threaded = True
from dune.fem import parameter
from dune.femdg.testing import run
from euler import sod as problem
# from euler import vortex as problem
# from euler import sod as problem
from euler import vortex as problem
# from euler import leVeque as problem
# from euler import radialSod3 as problem
......@@ -17,4 +19,4 @@ primitive=lambda Model,uh: {"pressure":Model.toPrim(uh)[2]}
run(*problem(),
startLevel=0, polOrder=2, limiter="default",
primitive=primitive, saveStep=0.1, subsamp=2,
dt=None)
dt=None,nThreads=4,grid="yasp")
......@@ -6,7 +6,7 @@ fem.parallel.numberofthreads: 1
# write diagnostics file (
# 0 don't, 1 only speedup file, 2 write all runfiles
# 3 only write 0, others at end, 4 all files at end for scaling)
fem.parallel.diagnostics: 1
fem.parallel.diagnostics: 4
# if true non-blocking communication is enabled
femdg.nonblockingcomm: true
fem.threads.verbose: true
......
......@@ -9,9 +9,9 @@ from ufl import dot, SpatialCoordinate
def run(Model, initial, x0,x1,N, endTime, name, exact,
polOrder, limiter="default", startLevel=0,
primitive=None, saveStep=None, subsamp=0,
dt=None):
dt=None,grid="yasp",nThreads=1):
domain = cartesianDomain(x0,x1,N,periodic=[False,False])
grid = create.grid("alucube",domain)
grid = create.grid(grid,domain)
grid.hierarchicalGrid.globalRefine(startLevel)
dimR = Model.dimension
t = 0
......@@ -20,7 +20,7 @@ def run(Model, initial, x0,x1,N, endTime, name, exact,
space = create.space("dgonb", grid, order=polOrder, dimrange=dimR)
u_h = space.interpolate(initial, name='u_h')
operator = createFemDGSolver( Model, space, limiter=limiter )
operator = createFemDGSolver( Model, space, limiter=limiter, nThreads=nThreads )
operator.applyLimiter( u_h );
print("number of elements: ",grid.size(0),flush=True)
if saveStep is not None:
......
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