dolfin team mailing list archive

[Marie Rognes <meg@xxxxxxxxx>]

Anders Logg wrote:
> On Sun, Jan 31, 2010 at 03:26:30PM -0800, Johan Hake wrote:
>   
> > Hello.
> >
> > The mesh registered in the forms (a and L) is the original mesh from the first
> > itteration, as it is cached in the form_data. ffc.jit does not recognize when
> > a mesh is changed.
> >     
>
> This shouldn't matter. The code should be the same when assembling on
> different meshes. Something else must be going wrong. Don't know what
> yet.
>
>   

A quick fix might be to not cache the form_data in jit_form?


--
Marie

> --
> Anders
>
>
>   
> > Not sure how this can be fixed. We might use the form together with the
> > id(mesh) as keys in the cache.
> >
> > Johan
> >
> > On Sunday 31 January 2010 13:34:00 Marie Rognes wrote:
> >     
> > > This gives me a terrible PETSc crash:
> > >
> > > from dolfin import *
> > >
> > > def run(mesh):
> > >
> > >     V = FunctionSpace(mesh, "CG", 1)
> > >     v = TestFunction(V)
> > >     u = TrialFunction(V)
> > >     f = Expression("1.0")
> > >
> > >     a = dot(v, u)*dx
> > >     L = dot(v, f)*dx
> > >
> > >     pde = VariationalProblem(a, L)
> > >     w = pde.solve()
> > >     plot(w)
> > >
> > > meshes = [UnitSquare(n, n) for n in range(1, 3)]
> > >
> > > for mesh in meshes:
> > >     run(mesh)
> > >
> > >
> > >
> > > _______________________________________________
> > > Mailing list: https://launchpad.net/~dolfin
> > > Post to     : dolfin@xxxxxxxxxxxxxxxxxxx
> > > Unsubscribe : https://launchpad.net/~dolfin
> > > More help   : https://help.launchpad.net/ListHelp
> > >
> > >       
> > _______________________________________________
> > Mailing list: https://launchpad.net/~dolfin
> > Post to     : dolfin@xxxxxxxxxxxxxxxxxxx
> > Unsubscribe : https://launchpad.net/~dolfin
> > More help   : https://help.launchpad.net/ListHelp
> >     
> > ------------------------------------------------------------------------
> >
> > _______________________________________________
> > Mailing list: https://launchpad.net/~dolfin
> > Post to     : dolfin@xxxxxxxxxxxxxxxxxxx
> > Unsubscribe : https://launchpad.net/~dolfin
> > More help   : https://help.launchpad.net/ListHelp
> >