You may or may not need to do the whole chain: PETSc, petsc4py, FIAT, UFL, FFC, COFFEE, PyOP2, and Firedrake Good luck! On 4 March 2015 at 18:16, Anna Kalogirou <a.kalogirou@leeds.ac.uk> wrote:

I haven't updated after the new year, so that might be it. Thanks!

Anna.

On 04/03/15 18:14, Andrew McRae wrote:

Is your Firedrake completely up-to-date?

It sounds not-unrelated to https://github.com/firedrakeproject/firedrake/issues/472, which was fixed recently.

On 4 March 2015 at 18:09, Anna Kalogirou <a.kalogirou@leeds.ac.uk> wrote:

...

Dear all,

I am experiencing a problem with extruded meshes. When I try to solve the nonlinear VP:

m = IntervalMesh(N,L) mesh = ExtrudedMesh(m, layers=10)

F = ( v*(p1-p0)/dt + v*(h0-1) + sigma*inner(grad(v),grad(h0)) + 0.5*inner(grad(p0_5),grad(p0_5))*v )*dx + E_t*v*ds_v(1)

I get the error

File

"/opt/firedrake/2.0/lib/python2.7/site-packages/firedrake/variational_solver.py", line 97, in __init__ form_compiler_parameters=self._problem.form_compiler_parameters) File "/opt/firedrake/2.0/lib/python2.7/site-packages/firedrake/assemble.py", line 58, in assemble form_compiler_parameters=form_compiler_parameters) File "/opt/firedrake/2.0/lib/python2.7/site-packages/firedrake/assemble.py", line 79, in _assemble kernels = ffc_interface.compile_form(f, "form", parameters=form_compiler_parameters) File

"/opt/firedrake/2.0/lib/python2.7/site-packages/firedrake/ffc_interface.py", line 273, in compile_form ((kernel, needs_orientations), ) = ffc_kernel.kernels ValueError: need more than 0 values to unpack

However, this problem does not arise when I switch to a RectangleMesh ( and change the ds_v(1) to ds(1) ), OR when I solve the corresponding linear VP. Anyone has any ideas why this is happening?

Thank you.

Best, Anna.

_______________________________________________ firedrake mailing list firedrake@imperial.ac.uk https://mailman.ic.ac.uk/mailman/listinfo/firedrake

--

Dr Anna Kalogirou Research Fellow School of Mathematics University of Leeds