[petsc-users] On the usage of MatSetValuesBlocked

Edoardo alinovi edoardo.alinovi at gmail.com
Thu Nov 3 11:30:52 CDT 2022


Hello  Matt,

I see, so from an operation point of view, I should pass a 1D array... How
should I unroll my nComp x nComp matrix? by row or by column?

e.g. say that my block is:

A=[
      1 2 3
      4 5 6
      7 8 9
     ]

than v = [1 2 3 4 5 6 7 8 9] or [1 4 7 2 5 8 3 6 9]?
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20221103/3b3313eb/attachment.html>


More information about the petsc-users mailing list