[firedrake] Problem with nonlinear solver on extruded mesh

Andrew McRae a.mcrae12 at imperial.ac.uk
Wed Mar 4 18:19:22 GMT 2015


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 at 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 at 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 at imperial.ac.uk
>> https://mailman.ic.ac.uk/mailman/listinfo/firedrake
>>
>
>
> --
>
>  Dr Anna Kalogirou
>  Research Fellow
>  School of Mathematics
>  University of Leeds
>
>
-------------- next part --------------
HTML attachment scrubbed and removed


More information about the firedrake mailing list