Say if I had this mixed function space:

parameters["matnest"] = False

Q = FunctionSpace(mesh,"CG",1)

W = Q*Q*Q

u1,u2,u3 = TrialFunctions(W)

v1,v2,v3 = TestFunctions(W)

What will the corresponding global matrix format be? Will it be of type MATMPIBAIJ of block size 3? Or will it simply be a monolithic MATMPIAIJ matrix?

Thanks,

Justin