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
