On 14 Jan 2016, at 07:34, Justin Chang <jychang48@gmail.com> wrote:
Hi all,
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?
It will be an mpiaij matrix. Due to how we build the dof maps, baij does not work here because successive dofs from the different fields in the mixed space do not live next to one another. Lawrence
Thanks, Justin _______________________________________________ firedrake mailing list firedrake@imperial.ac.uk https://mailman.ic.ac.uk/mailman/listinfo/firedrake