> On Nov 4, 2022, at 6:55 AM, Edoardo alinovi <edoardo.alin...@gmail.com> wrote: > > 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) Ah, likely the block size for the vector was not correct, leading to the memory corruption. MatCreateVecs() creates a vector compatible with the matrix, same block size and parallel layout so you don't need to worry about setting those values yourself. Barry > > 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
- Re: [petsc-users] On the usage of MatSetValuesBlocked Edoardo alinovi
- Re: [petsc-users] On the usage of MatSetValuesBlocked Edoardo alinovi
- Re: [petsc-users] On the usage of MatSetValuesBlocked Edoardo alinovi
- Re: [petsc-users] On the usage of MatSetValuesBlocked Barry Smith
- Re: [petsc-users] On the usage of MatSetValuesBlocked Edoardo alinovi
- Re: [petsc-users] On the usage of MatSetValuesBlocked Barry Smith
- Re: [petsc-users] On the usage of MatSetValuesBlocked Edoardo alinovi
- Re: [petsc-users] On the usage of MatSetValuesBlocked Edoardo alinovi
- Re: [petsc-users] On the usage of MatSetValuesBlocked Matthew Knepley
- Re: [petsc-users] On the usage of MatSetValuesBlocked Edoardo alinovi
- Re: [petsc-users] On the usage of MatSetValuesBlocked Barry Smith
- Re: [petsc-users] On the usage of MatSetValuesBlocked Edoardo alinovi
- Re: [petsc-users] On the usage of MatSetValuesBlocked Edoardo alinovi
- Re: [petsc-users] On the usage of MatSetValuesBlocked Pierre Jolivet
- Re: [petsc-users] On the usage of MatSetValuesBlocked Edoardo alinovi