Re: [petsc-dev] ksp_error_if_not_converged in multilevel solvers

2019-10-20 Thread Pierre Jolivet via petsc-dev
On Oct 20, 2019, at 6:07 PM, "Smith, Barry F." wrote: > > The reason the code works this way is that normally > -ksp_error_if_not_converged is propagated into the inner (and innerer) solves > and normally it is desirable that these inner solves do not error simply > because they reach

Re: [petsc-dev] BlockGetIndices and GetBlockIndices

2019-10-20 Thread Smith, Barry F. via petsc-dev
> On Oct 21, 2019, at 12:23 AM, Pierre Jolivet > wrote: > > > >> On 21 Oct 2019, at 7:11 AM, Smith, Barry F. wrote: >> >> >> >>> On Oct 20, 2019, at 11:52 PM, Pierre Jolivet >>> wrote: >>> >>> >>> On 21 Oct 2019, at 6:42 AM, Smith, Barry F. wrote: Could you provide

Re: [petsc-dev] BlockGetIndices and GetBlockIndices

2019-10-20 Thread Pierre Jolivet via petsc-dev
> On 21 Oct 2019, at 7:11 AM, Smith, Barry F. wrote: > > > >> On Oct 20, 2019, at 11:52 PM, Pierre Jolivet >> wrote: >> >> >> >>> On 21 Oct 2019, at 6:42 AM, Smith, Barry F. wrote: >>> >>> Could you provide a use case where you want to access/have a block size of >>> a IS that is

Re: [petsc-dev] BlockGetIndices and GetBlockIndices

2019-10-20 Thread Smith, Barry F. via petsc-dev
> On Oct 20, 2019, at 11:52 PM, Pierre Jolivet > wrote: > > > >> On 21 Oct 2019, at 6:42 AM, Smith, Barry F. wrote: >> >> Could you provide a use case where you want to access/have a block size of >> a IS that is not an ISBlock? > > In the end, all I really want is get access to the

Re: [petsc-dev] BlockGetIndices and GetBlockIndices

2019-10-20 Thread Pierre Jolivet via petsc-dev
> On 21 Oct 2019, at 6:42 AM, Smith, Barry F. wrote: > > Could you provide a use case where you want to access/have a block size of a > IS that is not an ISBlock? In the end, all I really want is get access to the underlying is->data->idx without having to worry about the subclass of is.

Re: [petsc-dev] BlockGetIndices and GetBlockIndices

2019-10-20 Thread Smith, Barry F. via petsc-dev
> On Oct 16, 2019, at 9:41 AM, Stefano Zampini via petsc-dev > wrote: > > I just took a look at the ISGENERAL code. ISSetBlockSize_General just sets > the block size of the layout (??) > ISGetIndices always return the data->idx memory. > So, a more profound question is: what is the model

Re: [petsc-dev] BlockGetIndices and GetBlockIndices

2019-10-20 Thread Smith, Barry F. via petsc-dev
> On Oct 16, 2019, at 2:50 AM, Pierre Jolivet via petsc-dev > wrote: > > Hello, > I’m trying to understand what is the rationale for naming a function > ISBlockGetIndices and another ISLocalToGlobalMappingGetBlockIndices (BlockGet > vs. GetBlock). ISBlockGetIndices returns the indices

Re: [petsc-dev] ksp_error_if_not_converged in multilevel solvers

2019-10-20 Thread Smith, Barry F. via petsc-dev
The reason the code works this way is that normally -ksp_error_if_not_converged is propagated into the inner (and innerer) solves and normally it is desirable that these inner solves do not error simply because they reach the maximum number of iterations since for nested iterative methods

Re: [petsc-dev] ksp_error_if_not_converged in multilevel solvers

2019-10-20 Thread Mark Adams via petsc-dev
> If one just wants to run a fixed number of iterations, not checking for > convergence, why would one set ksp->errorifnotconverged to true? > > Good question. I can see not worrying too much about convergence on the coarse grids, but to not allow it ... and now that I think about it, it seems

[petsc-dev] ksp_error_if_not_converged in multilevel solvers

2019-10-20 Thread Pierre Jolivet via petsc-dev
Hello, I’m trying to get multilevel solvers to error out when coarse levels are not converging, but I’m failing… Could someone either tell me if this is not possible to do so, or help me find what the problem in my options is, please? (in src/ksp/ksp/examples/tutorials) $ mpirun -n 8 ./ex71