Thanks Matt, I have found out that setValuesblocked will work if I do:
call MatCreateVecs(A, x, y, ierr) call setValuesBlocked(x, nblocks, varray, ierr) However, there is nogetValuesBlocked. Not the end of the world, it is handy to set and get stuff by block and not by single entry :) Cheers
