Re: [petsc-users] MatGetDiagonalBlock and shell matrices

2016-08-29 Thread Steven Dargaville
HI all Thanks for this Barry, the branch you've made seems to be working perfectly. Do you think this change will make it into a future PETSc release? I can imagine this would be very useful for a number of people. Thanks Steven On Sat, Aug 27, 2016 at 11:02 PM, Jed Brown

Re: [petsc-users] MatGetDiagonalBlock and shell matrices

2016-08-26 Thread Barry Smith
I have added support for this operation in the branch barry/add-matshell-matgetdiagonalblock-fortran Barry Because this required changing the handling of the method MatGetDiagonalBlock() the addition cannot be back ported to PETSc 3.7.x > On Aug 26, 2016, at 8:22 AM, Dave May

Re: [petsc-users] MatGetDiagonalBlock and shell matrices

2016-08-26 Thread Dave May
On Friday, 26 August 2016, Steven Dargaville wrote: > Hi Dave > > Thanks for the response. I'm actually using fortran and I wasn't sure that > PetscObjectComposeFunction would be accessible, and if so, what sort of > fortran magic I might need to call this function

Re: [petsc-users] MatGetDiagonalBlock and shell matrices

2016-08-26 Thread Steven Dargaville
Hi Dave Thanks for the response. I'm actually using fortran and I wasn't sure that PetscObjectComposeFunction would be accessible, and if so, what sort of fortran magic I might need to call this function (possibly an interface, with or without c_opt). Do you know if it is possible to call that

Re: [petsc-users] MatGetDiagonalBlock and shell matrices

2016-08-26 Thread Dave May
On 26 August 2016 at 14:34, Dave May wrote: > > > On 26 August 2016 at 14:14, Steven Dargaville > wrote: > >> Hi all >> >> I'm just wondering if there is any plans in the future for >> MatGetDiagonalBlock to support shell matrices by

Re: [petsc-users] MatGetDiagonalBlock and shell matrices

2016-08-26 Thread Dave May
On 26 August 2016 at 14:14, Steven Dargaville wrote: > Hi all > > I'm just wondering if there is any plans in the future for > MatGetDiagonalBlock to support shell matrices by registering a > user-implemented MATOP? MatGetDiagonal supports MATOP, but the block >

[petsc-users] MatGetDiagonalBlock and shell matrices

2016-08-26 Thread Steven Dargaville
Hi all I'm just wondering if there is any plans in the future for MatGetDiagonalBlock to support shell matrices by registering a user-implemented MATOP? MatGetDiagonal supports MATOP, but the block version of this does not. I found a previous query on the user list which touched on this and