Comment 3 for bug 1083720

On 11/27/2012 07:31 PM, Martin Sandve Alnæs wrote:
> Probably a caching issue somewhere in the PyDolfin layer.

Isn't it always somewhere in the PyDolfin layer :)

> The form sigature is unaffected by the mesh.

Could it be that it is the formdata that need to be recomputed?

Johan

> Den 27. nov. 2012 19:15 skrev "Simon Funke" <email address hidden>
> følgende:
>
>> Public bug reported:
>>
>> I submit this bug report following Marie's suggestions on the related
>> DOLFIN question (https://answers.launchpad.net/dolfin/+question/215120)
>>
>> DOLFIN crashes in the solve() routine after using replace() to replace
>> the function spaces of a form.
>>
>> While the solve succeeds if replaced function spaces have:
>> - a different polynomial degree or
>> - a different polynomial degree and a different mesh,
>> it does not work if the new function spaces is based on a new mesh only.
>>
>> This is the error message I get:
>> ***
>> -------------------------------------------------------------------------
>> *** Error: Unable to define linear variational problem a(u, v) = L(v) for
>> all v.
>> *** Reason: Expecting the solution variable u to be a member of the trial
>> space.
>> *** Where: This error was encountered inside LinearVariationalProblem.cpp.
>> ***
>> ------------------------------------------------------------------------
>>
>> The attached file demonstrates the problem.
>>
>> ** Affects: dolfin
>> Importance: Undecided
>> Status: New
>>
>> ** Attachment added: "poisson.py"
>>
>> https://bugs.launchpad.net/bugs/1083720/+attachment/3445290/+files/poisson.py
>>
>> --
>> You received this bug notification because you are a member of DOLFIN
>> Core Team, which is subscribed to DOLFIN.
>> https://bugs.launchpad.net/bugs/1083720
>>
>> Title:
>> replace() fails if the new function spaces have a different mesh
>>
>> Status in DOLFIN:
>> New
>>
>> Bug description:
>> I submit this bug report following Marie's suggestions on the related
>> DOLFIN question
>> (https://answers.launchpad.net/dolfin/+question/215120)
>>
>> DOLFIN crashes in the solve() routine after using replace() to replace
>> the function spaces of a form.
>>
>> While the solve succeeds if replaced function spaces have:
>> - a different polynomial degree or
>> - a different polynomial degree and a different mesh,
>> it does not work if the new function spaces is based on a new mesh only.
>>
>> This is the error message I get:
>> ***
>> -------------------------------------------------------------------------
>> *** Error: Unable to define linear variational problem a(u, v) = L(v)
>> for all v.
>> *** Reason: Expecting the solution variable u to be a member of the
>> trial space.
>> *** Where: This error was encountered inside
>> LinearVariationalProblem.cpp.
>> ***
>> ------------------------------------------------------------------------
>>
>> The attached file demonstrates the problem.
>>
>> To manage notifications about this bug go to:
>> https://bugs.launchpad.net/dolfin/+bug/1083720/+subscriptions
>>
>