[petsc-users] transform in a block matrix

2018-04-18 Thread Sonia Pozzi
Dear Petsc users and developers, I have a MPIAIJ matrix (already assembled) and I would like to transform it in a block matrix with of block dimension 3. Is there a easy way to do that? Thank you in advance, Sonia Ph.D. Student Group of Prof. Rolf Krause Institute of Computational Science Uni

Re: [petsc-users] PC FieldSplit - How to Count the Iterations

2018-03-15 Thread Sonia Pozzi
e: +41 58 666 4972 > On 14 Mar 2018, at 16:00, Matthew Knepley wrote: > > On Wed, Mar 14, 2018 at 11:21 PM, Sonia Pozzi <mailto:sonia.po...@usi.ch>> wrote: > Dear Barry, > > thank you for the answer. That helped a lot. Just a second curiosity. > I’m setting A00 to

Re: [petsc-users] PC FieldSplit - How to Count the Iterations

2018-03-14 Thread Sonia Pozzi
Dear Barry, thank you for the answer. That helped a lot. Just a second curiosity. I’m setting A00 to be solved with preonly+lu. I obtain the following ksp_0 KSPGetTotalIterations: 26 ksp_1 KSPGetTotalIterations: 22 Residual ksp_0: 0 Reason ksp_0: 4 Solution ksp_0 : Convergence in 1 iterations.

[petsc-users] PC FieldSplit - How to Count the Iterations

2018-03-14 Thread Sonia Pozzi
Schur complement. What is the number of its that I get from subksp_0? Is it the total number of calls or is it just related to the last solution for A00? In the second case, is there a strategy to have the total number? My best regards and thank you, Sonia Pozzi Ph.D. Student Group of Prof. Rolf

[petsc-users] ASM with matrix-free method

2017-02-02 Thread Sonia Pozzi
Dear Barry, At the moment I’m dealing with the following problem. Suppose I want to solve Sx=b, where is a product of matrices. In particular inside it contains also an inverse of a matrix. To compute the action of S I have created a shell matrix, but know I would like to apply the ASM as prec