dolfin team mailing list archive

[Marie Rognes <meg@xxxxxxxxx>]

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)