Construct and properly attach subfunctions
Default to 0 for incomplete CompositeFunctions

I'm not sure this is the right API.  Return a special code which
libMesh can treat as "leave existing values unchanged"?  Throw an
error, but trust to the MultiphysicsSystem to test for and fill in
missing indices with a ZeroFunction?