> On Sep 19, 2016, at 10:43 AM, Cyrill Vonplanta <cyrill.von.pla...@usi.ch>
> Thanks a lot. I’d like to use this for a nonlinear variant of a
> block-gauss-seidel smoother. I would like to use MatInvertBlockDiagonal for
> speeding up my variant.
> I think I can work with this, however I also have the problem to turn my
> initial matrix into one with a blocksize of 3.When I call:
> MatConvert(A, MATBAIJ, MAT_INITIAL_MATRIX, Dinverse);
MatConvert() creates a new matrix so you should not create the Dinverse
beforehand; anything you put in Dinverse before is lost..m
block size > 1 really only makes sense if the block size is really greater
than one. So if A has blocks of size 3 you should create A as BAIJ and thus
never need to call the convert routine.
You can also set the block size for AIJ matrix to 3 and use
MatInvertBlockDiagonal() on that matrix and not use the BAIJ matrix.
Finally note that MatInvertBlockDiagonal() ends up calling (for block size
3) PetscKernel_A_gets_inverse_A_3() for each block. It sounds to me like that
is what you would want for nonlinear variant of a block-gauss-seidel smoother..
> Then the matrix Dinverse has blocksize 1 which comes from A. I checked the
> blocksize before the conversion and it was 3, so it seems to get lost.
> What is the correct (and elegant) way to turn a matrix into a block matrix?
> On 19/09/16 17:18, "Barry Smith" <bsm...@mcs.anl.gov> wrote:
>> This is very specialized for implementing point block Jacobi; I don't
>> think it is something you would want to use directly.
>> If you do want to use it, it simply returns the inverses of the block
>> diagonals in column major form. You can then call MatSetValues() with with
>> each of those blocks into another PETSc matrix. values[i*bs*bs] is the
>> starting point of each block in the array.
>>> On Sep 19, 2016, at 4:55 AM, Cyrill Vonplanta <cyrill.von.pla...@usi.ch>
>>> Dear PETSc-Users,
>>> I would like to use the inverted block diagonals of a a matrix. I have seen
>>> the function MatInvertBlockDiagonal() but I don’t know how to create a
>>> matrix out of them or an array of block matrizes.
>>> Does anyone have an example on how to use **values to create a PETSc