There is a reason PETSc does not have variable block matrices. It is messy
and rarely a big win. It is doable, ML/Aztec does it, but it is messy.
On Sat, Dec 1, 2018 at 6:20 PM Smith, Barry F. via petsc-users <
petsc-users@mcs.anl.gov> wrote:
>
>Well, you need to start somewhere. It is going
Well, you need to start somewhere. It is going to take a combination of
this, MatLUFactorSymbolic_SeqAIJ(), MatLUFactorSymbolic_SeqAIJ() and
MatLUFactorNumerical_SeqAIJ() to get what you want. Note that one of the first
things you need to do is given the MatGetVariableBlockSizes() determine
This would use the blocks as indicated by MatSetVariableBlockSizes, not
by using a matrix format that explicitly blocks? The information it
needs is too big/complex to yield to the preconditioner via a call like
MatInvertVariableBlockDiagonal. Instead, I would model it off MatFactor
and the ILU c
Hi,
I want to implement a variable block size BILU preconditioner in PETSc.
I am wondering if modifying the PCVPBJACOBI preconditioner
(vpbjacobi.c), to get what I want, is a good idea.
Any help is much appreciated.
--
Ali Reza Khaz’ali
Assistant Professor of Petroleum Engineering,
Departmen
On Fri, Nov 30, 2018 at 10:19 AM NENNIG Benoit via petsc-users <
petsc-users@mcs.anl.gov> wrote:
> Dear petsc users,
>
> I have parallel matrix A (mpiaij) and I would like to create a matrix B
> like
> B = [A v
> wT 0 ]
> where wT and v are vectors.
> A is involved in eigenvalue computation (