[petsc-users] On the usage of MatSetValuesBlocked

[Edoardo alinovi]

Hi Barry,

Thanks for popping in. This is my code:


*                blockValues =
this%getDiagonalBlockValues(iElement=iElement)                call
MatSetValuesBlocked(this%A, 4-bdim, mesh%cellGlobalAddr(iElement)-1,
4-bdim, mesh%cellGlobalAddr(iElement)-1, blockValues, INSERT_VALUES, ierr)*

*blockValues *is a 3x3 or 4x4 matrix and I am passing it straight away.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20221103/58d63cc8/attachment-0001.html>