Comment 7 for bug 1083720
Revision history for this message
| Martin Sandve Alnæs (martinal) wrote Re: [Bug 1083720] [NEW] replace() fails if the new function spaces have a different mesh | #7 |
© 2004
Canonical Ltd.
•
Terms of use
•
Data privacy
•
Contact Launchpad Support
•
Blog
•
Careers
•
System status
•
0e1f616
(Get the code!)
I think I know... Create an Argument instead of TestFunction and it should
work fine.
Den 27. nov. 2012 21:55 skrev "Johan Hake" <email address hidden>
følgende:
> On 11/27/2012 09:04 PM, Martin Sandve Alnæs wrote:
> > The form data is not carried over to the new form returned by replace, so
> > no...
>
> :)
>
> It actually looks like it is replace or compute_form_data that is not
> doing its job. If you add the following lines:
>
> form_data = F_fine.
>
> print [V_fine == func.function_
> form_data.
> print [V == func.function_
>
> after:
> F_fine = replace(F, {u:u_fine, v:v_fine})
>
> you see that it is still the original functionspace that is used.
>
> Johan
>
> > Den 27. nov. 2012 20:00 skrev "Johan Hake" <email address hidden>
> > følgende:
> >
> >> 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:/
> >>>>
> >>>> 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
> >> LinearVariation
> >>>> ***
> >>>>
> -------
> >>>>
> >>>> The attached file demonstrates the problem.
> >>>>
> >>>> ** Affects: dolfin
> >>>> Importance: Undecided
> >>>> Status: New
> >>>>
> >>>> ** Attachment added: "poisson.py"
> >>>>
> >>>>
> >>
> https:/
> >>>>
> >>>> --
> >>>> You received this bug notification because you are a member of DOLFIN
> >>>> Core Team, which is subscribed to DOLFIN.
> >>>> https:/
> >>>>
> >>>> 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:/
> >>>>
> >>>> 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
> >>>> LinearVariation
> >>>> ***
> >>>>
> -------
> >>>>
> >>>> The attached file demonstrates the problem.
> >>>>
> >>>> To manage notifications about this bug go to:
> >>>> https:/
> >>>>
> >>>
> >>
> >> --
> >> You received this bug notification because you are a member of DOLFIN
> >> Core Team, which is subscribed to DOLFIN.
> >> https:/
> >>
> >> 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:/
> >>
> >> 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
> >> LinearVariation
> >> ***
> >> -------
> >>
> >> The attached file demonstrates the problem.
> >>
> >> To manage notifications about this bug go to:
> >> https:/
> >>
> >
>
> --
> You received this bug notification because you are a member of DOLFIN
> Core Team, which is subscribed to DOLFIN.
> https:/
>
> 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:/
>
> 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
> LinearVariation
> ***
> -------
>
> The attached file demonstrates the problem.
>
> To manage notifications about this bug go to:
> https:/
>