>
>
>> I think I may know what your problem is. Plex evaluates the blocksize by
> looking for an equal number of dofs
> on each point. This is sufficient, but not necessary. If you are using
> higher order methods, there is block structure
> there that I will not see.
>
I don't understand what the order has to do with it.
I use code like this to setup the dofs:
for (ii=0;ii<ctx->num_species;ii++) {
ierr = PetscFECreateDefault(PetscObjectComm((PetscObject) dm), dim, 1,
ctx->simplex, NULL, PETSC_DECIDE, &ctx->fe[ii]);CHKERRQ(ierr);
ierr = DMSetField(dm, ii, NULL, (PetscObject)
ctx->fe[ii]);CHKERRQ(ierr);
}
Everything is constant, elements (eg, Q3) and dofs/vertex.
>
> Jed, is there an obvious way to see that structure that I am missing?
>
> Thanks,
>
> Matt
>
> --
> What most experimenters take for granted before they begin their
> experiments is infinitely more interesting than any results to which their
> experiments lead.
> -- Norbert Wiener
>
> https://www.cse.buffalo.edu/~knepley/
> <http://www.cse.buffalo.edu/~knepley/>
>
