Barry,

Can you please provide me with an example on how to use MatSetValuesBlocked? 

To play it easy, let's say that I want to insert a 3x3 block matrix b into the matrix A, rows 0-2, columns 0-2. Up to what I've understood (very few apparently XD ), I would do like this:

b(3,3) = 11.0 

call MatSetValuesBlocked(A, 3, 0, 3, 0, b, INSERT_VALUES, ierr).  

This does not work at all, I get this result that does not make any sense 😫

It places 6 values instead of 9 and it they are in odd locations (0 1 2 9 10 11).

Also I noted that I am getting  different results if in place of  the zero in red I use a fortran integer 😥

Super thanks for the help