Assembling of Navier-Stokes gradient form segfaults

I came across this curious issue while playing with NavierStokes, and I can't explain it to myself.

The following code segfaults at the last instruction on my debian testing, dolfin 1.4.0 from repos

V = VectorElement("CG", triangle, 2)
Q = FiniteElement("CG", triangle, 1)

u = TrialFunction(V)
v = TestFunction(V)
p = TrialFunction(Q)
q = TestFunction(Q)

b = -q*div(u)*dx
g = -p*div(v)*dx

forms = [b,g]

#include <dolfin.h>

#ifdef USING_MSHR
#include <mshr.h>
#endif

#include "NavierStokesYosida.h"

using namespace dolfin;

int main(int argc, char *argv[])
{
    int resolution = 50;
#ifndef USING_MSHR
    /*
     * create mesh legacy version
     */
    Rectangle cavity(0,0,1,1);
    Mesh mesh(cavity,resolution);
#else
    /*
     * create mesh mshr version
     */
    mshr::Rectangle cavity(Point(0,0,0),Point(1,1,0));
    Mesh mesh;
    mshr::generate(mesh,cavity,resolution);
#endif

    NavierStokesYosida::Form_g::TestSpace V(mesh);
    NavierStokesYosida::Form_b::TestSpace Q(mesh);
    NavierStokesYosida::Form_b b(V,Q);
    NavierStokesYosida::Form_g g(Q,V);
    PETScMatrix B;
    PETScMatrix G;
    assemble(B,b);
    assemble(G,g);
}

What's most puzzling, it crashes when assembling G, but not B, which is very similar.

I attach main.cpp, ufl file and CMakeLists.txt. There's code for both meshing with mshr and legacy meshing.

Comments (9)