Re: [petsc-users] Dropping single entries in matrix insertion for block matrices

2018-12-10 Thread Smith, Barry F. via petsc-users
> On Dec 10, 2018, at 8:30 PM, Jed Brown wrote: > > "Smith, Barry F." writes: > >>> On Dec 10, 2018, at 7:58 PM, Jed Brown wrote: >>> >>> "Smith, Barry F. via petsc-users" writes: >>> Well the thing is the indices you pass in are currently always "blocked base", they work aut

Re: [petsc-users] Dropping single entries in matrix insertion for block matrices

2018-12-10 Thread Jed Brown via petsc-users
"Smith, Barry F." writes: >> On Dec 10, 2018, at 7:58 PM, Jed Brown wrote: >> >> "Smith, Barry F. via petsc-users" writes: >> >>> Well the thing is the indices you pass in are currently always >>> "blocked base", they work automatically for both blocked and non >>> blocked versions of M

Re: [petsc-users] Dropping single entries in matrix insertion for block matrices

2018-12-10 Thread Smith, Barry F. via petsc-users
> On Dec 10, 2018, at 7:58 PM, Jed Brown wrote: > > "Smith, Barry F. via petsc-users" writes: > >> Well the thing is the indices you pass in are currently always >> "blocked base", they work automatically for both blocked and non >> blocked versions of MatSetValuesLocal(). > > What d

Re: [petsc-users] Dropping single entries in matrix insertion for block matrices

2018-12-10 Thread Jed Brown via petsc-users
"Smith, Barry F. via petsc-users" writes: >Well the thing is the indices you pass in are currently always >"blocked base", they work automatically for both blocked and non >blocked versions of MatSetValuesLocal(). What do you mean? You pass block indices to MatSetValuesBlocked and

Re: [petsc-users] Dropping single entries in matrix insertion for block matrices

2018-12-10 Thread Smith, Barry F. via petsc-users
> On Dec 10, 2018, at 5:57 AM, Lawrence Mitchell wrote: > > Hi Barry, > >> On 9 Dec 2018, at 19:25, Smith, Barry F. wrote: >> >> >> Lawrence, >> >> I understand what you want and it is a reasonable request. The problem >> is that currently ISLocalToGlobalMappingCreate() when used wi

Re: [petsc-users] Dropping single entries in matrix insertion for block matrices

2018-12-10 Thread Lawrence Mitchell via petsc-users
Hi Barry, > On 9 Dec 2018, at 19:25, Smith, Barry F. wrote: > > > Lawrence, > > I understand what you want and it is a reasonable request. The problem > is that currently ISLocalToGlobalMappingCreate() when used with block vectors > and matrices is always based on blocks, that is, fro

Re: [petsc-users] Dropping single entries in matrix insertion for block matrices

2018-12-09 Thread Smith, Barry F. via petsc-users
Lawrence, I understand what you want and it is a reasonable request. The problem is that currently ISLocalToGlobalMappingCreate() when used with block vectors and matrices is always based on blocks, that is, from the manual page, "There is one integer value in indices per block and i

[petsc-users] Dropping single entries in matrix insertion for block matrices

2018-12-09 Thread Lawrence Mitchell via petsc-users
Dear petsc-users, I have a matrix with a block size > 1. I would sometimes like to insert into it, applying Dirichlet conditions to one component of the block. Now, for normal block (or when the block size is 1) dirichlet conditions, I do this by swapping out the local to global map for one wit