Mixed Helmholtz attempt.
Using branches hdivcurl of FIAT/UFL/FFC, extrusion-vfs of PyOP2, extrusion-hdivcurl-proper of firedrake
from firedrake import *
m = UnitSquareMesh(4, 4)
mesh = ExtrudedMesh(m, 5, layer_height=0.25)
horiz_elt = FiniteElement("RT", "triangle", 1)
vert_elt = FiniteElement("DG", "interval", 0)
product_elt = HDiv(OuterProductElement(horiz_elt, vert_elt))
V1 = FunctionSpace(mesh, product_elt)
horiz_elt = FiniteElement("DG", "triangle", 0)
vert_elt = FiniteElement("DG", "interval", 0)
product_elt = OuterProductElement(horiz_elt, vert_elt)
V2 = FunctionSpace(mesh, product_elt)
f = Function(V2)
exact = Function(V2)
f.interpolate(Expression("(1+8*pi*pi)*sin(x[0]*pi*2)*sin(x[1]*pi*2)"))
exact.interpolate(Expression("sin(x[0]*pi*2)*sin(x[1]*pi*2)"))
W = V1 * V2
u, p = TrialFunctions(W)
v, q = TestFunctions(W)
a = (p*q - q*div(u) + dot(v, u) + div(v)*p)*dx
L = f*q*dx
out = Function(W)
solve(a == L, out)
print out.dat.data
# print sqrt(assemble((out[2]-exact)*(out[2]-exact)*dx))
atm112@ubuntu:~$ python mixedbork.py
In instant.recompile: The module did not compile with command 'python setup.py build_ext install --install-platlib=.', see '/home/atm112/.instant/error/instant_module_c6462c135a902b9c4679ed9d63d5d267d43b9876/compile.log'
Traceback (most recent call last):
File "mixedbork.py", line 28, in <module>
solve(a == L, out)
File "/home/atm112/local/firedrake/python/firedrake/solving.py", line 644, in solve
_solve_varproblem(*args, **kwargs)
File "/home/atm112/local/firedrake/python/firedrake/solving.py", line 672, in _solve_varproblem
solver.solve()
File "/home/atm112/local/firedrake/python/firedrake/solving.py", line 229, in solve
self.snes.solve(None, v)
File "SNES.pyx", line 404, in petsc4py.PETSc.SNES.solve (src/petsc4py.PETSc.c:124245)
File "petscsnes.pxi", line 223, in petsc4py.PETSc.SNES_Function (src/petsc4py.PETSc.c:27067)
File "/home/atm112/local/firedrake/python/firedrake/solving.py", line 156, in form_function
with self._F_tensor.dat.vec_ro as v:
File "/usr/lib/python2.7/contextlib.py", line 17, in __enter__
return self.gen.next()
File "/home/atm112/local/PyOP2/pyop2/petsc_base.py", line 149, in vecscatter
with acc(d) as v:
File "/usr/lib/python2.7/contextlib.py", line 17, in __enter__
return self.gen.next()
File "/home/atm112/local/PyOP2/pyop2/petsc_base.py", line 84, in vec_context
self._force_evaluation()
File "/home/atm112/local/PyOP2/pyop2/base.py", line 1335, in _force_evaluation
_trace.evaluate(set([self]), set([self]))
File "/home/atm112/local/PyOP2/pyop2/base.py", line 145, in evaluate
comp._run()
File "/home/atm112/local/PyOP2/pyop2/base.py", line 2881, in _run
return self.compute()
File "/home/atm112/local/PyOP2/pyop2/base.py", line 2888, in compute
self._compute_if_not_empty(self.it_space.iterset.core_part)
File "/home/atm112/local/PyOP2/pyop2/base.py", line 2900, in _compute_if_not_empty
self._compute(part)
File "/home/atm112/local/PyOP2/pyop2/sequential.py", line 101, in _compute
fun(*self._jit_args)
File "/home/atm112/local/PyOP2/pyop2/host.py", line 411, in __call__
return self.compile()(*args)
File "/home/atm112/local/PyOP2/pyop2/host.py", line 452, in compile
modulename=self._kernel.name if configuration["debug"] else None)
File "/usr/local/lib/python2.7/dist-packages/instant/inlining.py", line 95, in inline_with_numpy
module = build_module(**kwargs)
File "/usr/local/lib/python2.7/dist-packages/instant/build.py", line 542, in build_module
recompile(modulename, module_path, new_compilation_checksum, build_system)
File "/usr/local/lib/python2.7/dist-packages/instant/build.py", line 123, in recompile
instant_error(msg % (cmd, compile_log_filename_dest))
File "/usr/local/lib/python2.7/dist-packages/instant/output.py", line 57, in instant_error
raise RuntimeError(text)
RuntimeError: In instant.recompile: The module did not compile with command 'python setup.py build_ext install --install-platlib=.', see '/home/atm112/.instant/error/instant_module_c6462c135a902b9c4679ed9d63d5d267d43b9876/compile.log'
running build_ext
building '_instant_module_c6462c135a902b9c4679ed9d63d5d267d43b9876' extension
creating build
creating build/temp.linux-x86_64-2.7
mpicc -fno-strict-aliasing -DNDEBUG -g -fwrapv -O2 -Wall -Wstrict-prototypes -fPIC -I/usr/local/lib/python2.7/dist-packages/petsc/include -I/usr/lib/python2.7/dist-packages/numpy/core/include -I/usr/include/python2.7 -c mat_utils.cxx -o build/temp.linux-x86_64-2.7/mat_utils.o
cc1plus: warning: command line option ‘-Wstrict-prototypes’ is valid for C/ObjC but not for C++ [enabled by default]
mpicc -fno-strict-aliasing -DNDEBUG -g -fwrapv -O2 -Wall -Wstrict-prototypes -fPIC -I/usr/local/lib/python2.7/dist-packages/petsc/include -I/usr/lib/python2.7/dist-packages/numpy/core/include -I/usr/include/python2.7 -c instant_module_c6462c135a902b9c4679ed9d63d5d267d43b9876_wrap.cxx -o build/temp.linux-x86_64-2.7/instant_module_c6462c135a902b9c4679ed9d63d5d267d43b9876_wrap.o
cc1plus: warning: command line option ‘-Wstrict-prototypes’ is valid for C/ObjC but not for C++ [enabled by default]
In file included from /usr/lib/python2.7/dist-packages/numpy/core/include/numpy/ndarraytypes.h:1728:0,
from /usr/lib/python2.7/dist-packages/numpy/core/include/numpy/ndarrayobject.h:17,
from /usr/lib/python2.7/dist-packages/numpy/core/include/numpy/arrayobject.h:15,
from instant_module_c6462c135a902b9c4679ed9d63d5d267d43b9876_wrap.cxx:3197:
/usr/lib/python2.7/dist-packages/numpy/core/include/numpy/npy_deprecated_api.h:11:2: warning: #warning "Using deprecated NumPy API, disable it by #defining NPY_NO_DEPRECATED_API NPY_1_7_API_VERSION" [-Wcpp]
#warning "Using deprecated NumPy API, disable it by #defining NPY_NO_DEPRECATED_API NPY_1_7_API_VERSION"
^
instant_module_c6462c135a902b9c4679ed9d63d5d267d43b9876_wrap.cxx: In function ‘void wrap_form_cell_integral_0_otherwise__(PyObject*, PyObject*, PyObject*, PyObject*, PyObject*, PyObject*, PyObject*, PyObject*, PyObject*, PyObject*, PyObject*, PyObject*, PyObject*, PyObject*, PyObject*, PyObject*, PyObject*, PyObject*, PyObject*, PyObject*)’:
instant_module_c6462c135a902b9c4679ed9d63d5d267d43b9876_wrap.cxx:3254:11: error: ‘xtr_arg0_0_map1_0’ was not declared in this scope
xtr_arg0_0_map1_0[0] = *(arg0_0_map1_0 + i * 1 + 0);;
^
instant_module_c6462c135a902b9c4679ed9d63d5d267d43b9876_wrap.cxx:3257:74: warning: left operand of comma operator has no effect [-Wunused-value]
form_cell_integral_0_otherwise(arg0_0 + xtr_arg0_0_map0_0[i_0]*(1, 1), arg1_0_vec, arg2_0_vec, arg3_0_vec, i_0 + 0);
^
instant_module_c6462c135a902b9c4679ed9d63d5d267d43b9876_wrap.cxx:3280:40: warning: left operand of comma operator has no effect [-Wunused-value]
arg3_0_vec[0] += _off30[0] * (1, 1);
^
instant_module_c6462c135a902b9c4679ed9d63d5d267d43b9876_wrap.cxx:3281:40: warning: left operand of comma operator has no effect [-Wunused-value]
arg3_0_vec[1] += _off30[1] * (1, 1);
^
instant_module_c6462c135a902b9c4679ed9d63d5d267d43b9876_wrap.cxx:3282:40: warning: left operand of comma operator has no effect [-Wunused-value]
arg3_0_vec[2] += _off30[2] * (1, 1);
^
instant_module_c6462c135a902b9c4679ed9d63d5d267d43b9876_wrap.cxx:3283:40: warning: left operand of comma operator has no effect [-Wunused-value]
arg3_0_vec[3] += _off30[3] * (1, 1);
^
instant_module_c6462c135a902b9c4679ed9d63d5d267d43b9876_wrap.cxx:3291:47: error: ‘xtr_arg0_0_map0_1’ was not declared in this scope
form_cell_integral_0_otherwise(arg0_0 + xtr_arg0_0_map0_1[i_0]*(1, 1), arg1_0_vec, arg2_0_vec, arg3_0_vec, i_0 + 3);
^
instant_module_c6462c135a902b9c4679ed9d63d5d267d43b9876_wrap.cxx:3291:74: warning: left operand of comma operator has no effect [-Wunused-value]
form_cell_integral_0_otherwise(arg0_0 + xtr_arg0_0_map0_1[i_0]*(1, 1), arg1_0_vec, arg2_0_vec, arg3_0_vec, i_0 + 3);
^
instant_module_c6462c135a902b9c4679ed9d63d5d267d43b9876_wrap.cxx:3314:40: warning: left operand of comma operator has no effect [-Wunused-value]
arg3_0_vec[0] += _off30[0] * (1, 1);
^
instant_module_c6462c135a902b9c4679ed9d63d5d267d43b9876_wrap.cxx:3315:40: warning: left operand of comma operator has no effect [-Wunused-value]
arg3_0_vec[1] += _off30[1] * (1, 1);
^
instant_module_c6462c135a902b9c4679ed9d63d5d267d43b9876_wrap.cxx:3316:40: warning: left operand of comma operator has no effect [-Wunused-value]
arg3_0_vec[2] += _off30[2] * (1, 1);
^
instant_module_c6462c135a902b9c4679ed9d63d5d267d43b9876_wrap.cxx:3317:40: warning: left operand of comma operator has no effect [-Wunused-value]
arg3_0_vec[3] += _off30[3] * (1, 1);
^
instant_module_c6462c135a902b9c4679ed9d63d5d267d43b9876_wrap.cxx:3205:11: warning: unused variable ‘arg0_1’ [-Wunused-variable]
double *arg0_1 = (double *)(((PyArrayObject *)_arg0_1)->data);
^
instant_module_c6462c135a902b9c4679ed9d63d5d267d43b9876_wrap.cxx:3224:9: warning: unused variable ‘_off31’ [-Wunused-variable]
int * _off31 = (int *)(((PyArrayObject *)off31)->data);
^
error: command 'mpicc' failed with exit status 1