[petsc-users] matsetvaluesblocked4_

[Mark Adams]

>
>
>> 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/>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20200527/b58a178b/attachment.html>