dolfin team mailing list archive

[Johan Hake <johan.hake@xxxxxxxxx>]

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. 

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
>