diff --git a/firedrake/randomfields.py b/firedrake/randomfields.py index d4a8af8392..34ee5bb3a2 100644 --- a/firedrake/randomfields.py +++ b/firedrake/randomfields.py @@ -27,7 +27,6 @@ from pyop2 import op2 from pyop2.mpi import COMM_WORLD -from pyop2.parloop import DatParloopArg _default_pcg = PCG64() @@ -138,22 +137,9 @@ def WhiteNoise(V, rng=None, scale=1.0): include_dirs=BLASLAPACK_INCLUDE.split(), ldargs=BLASLAPACK_LIB.split()) - # Construct arguments for par loop - # ~ def get_map(x): - # ~ return x.cell_node_map() i, _ = mass_ker.indices - # ~ (access, get_map, i, *, function, V) - # ~ if i is None: - # ~ map_ = get_map(V) - # ~ return function.dat(access, map_) - # ~ else: - # ~ map_ = get_map(V[i]) - # ~ return function.dat[i](access, map_) - - # ~ z_arg = _vector_arg(op2.READ, get_map, i, function=iid_normal, V=Vbrok) z_arg = iid_normal.dat(op2.READ, Vbrok.cell_node_map()) - # ~ b_arg = _vector_arg(op2.INC, get_map, i, function=wnoise, V=V) b_arg = wnoise.dat(op2.INC, V.cell_node_map()) coords = mesh.coordinates