Ah I figured. Is it possible to one day add a feature where one could achieve BAIJ format if the fields all belong to same function space? I think there was a way to do this in DMPlex last I spoke with Matt Thanks, Justin On Thursday, January 14, 2016, Lawrence Mitchell < lawrence.mitchell@imperial.ac.uk> wrote:
On 14 Jan 2016, at 07:34, Justin Chang <jychang48@gmail.com <javascript:;>> 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 <javascript:;> https://mailman.ic.ac.uk/mailman/listinfo/firedrake
_______________________________________________ firedrake mailing list firedrake@imperial.ac.uk <javascript:;> https://mailman.ic.ac.uk/mailman/listinfo/firedrake