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 : [email protected]
Unsubscribe : https://launchpad.net/~dolfin
More help : https://help.launchpad.net/ListHelp
_______________________________________________
Mailing list: https://launchpad.net/~dolfin
Post to : [email protected]
Unsubscribe : https://launchpad.net/~dolfin
More help : https://help.launchpad.net/ListHelp
------------------------------------------------------------------------
_______________________________________________
Mailing list: https://launchpad.net/~dolfin
Post to : [email protected]
Unsubscribe : https://launchpad.net/~dolfin
More help : https://help.launchpad.net/ListHelp
_______________________________________________
Mailing list: https://launchpad.net/~dolfin
Post to : [email protected]
Unsubscribe : https://launchpad.net/~dolfin
More help : https://help.launchpad.net/ListHelp
