Re: [petsc-dev] Can MATAIJ support block size?

2018-09-20 Thread Fande Kong
Thanks Barry, I am going to take a try. Fande, On Thu, Sep 20, 2018 at 2:15 PM Smith, Barry F. wrote: > > Not a priority for us but if you implement clean code to do this we're > likely to accept it. > >Barry > > > > On Sep 20, 2018, at 12:27 PM, Fande Kong wrote: > > > > Hi

Re: [petsc-dev] Subsurface application and Algebraic Multigrid on GPUs

2018-09-20 Thread Smith, Barry F.
Brian, I have finished making the (relatively few) changes needed to get PETSc's GAMG to run on a combination of the CPU and GPU. Any of the AMG kernels that has a CUDA backed is run automatically on the GPU while the kernels without a CUDA backend are run on the CPU. In particular

Re: [petsc-dev] Can MATAIJ support block size?

2018-09-20 Thread Smith, Barry F.
Not a priority for us but if you implement clean code to do this we're likely to accept it. Barry > On Sep 20, 2018, at 12:27 PM, Fande Kong wrote: > > Hi Developers, > > MATBAIJ actually assumes that the point-block is dense. It is fine if the > block size is small for example less

Re: [petsc-dev] Can MATAIJ support block size?

2018-09-20 Thread Fande Kong
In order to make the discussion simple, we can just assume the sparse block is a diagonal matrix. How to take the advantage of the diagonal block matrix? Fande, On Thu, Sep 20, 2018 at 12:28 PM Fande Kong wrote: > > > On Thu, Sep 20, 2018 at 12:09 PM Matthew Knepley > wrote: > >> On Thu, Sep

Re: [petsc-dev] Can MATAIJ support block size?

2018-09-20 Thread Fande Kong
On Thu, Sep 20, 2018 at 12:09 PM Matthew Knepley wrote: > On Thu, Sep 20, 2018 at 1:26 PM Fande Kong wrote: > >> Hi Developers, >> >> MATBAIJ actually assumes that the point-block is dense. It is fine if the >> block size is small for example less than 50, but we may suffer from a >> memory

Re: [petsc-dev] Can MATAIJ support block size?

2018-09-20 Thread Matthew Knepley
On Thu, Sep 20, 2018 at 1:26 PM Fande Kong wrote: > Hi Developers, > > MATBAIJ actually assumes that the point-block is dense. It is fine if the > block size is small for example less than 50, but we may suffer from a > memory issue if the block is large such as 1000. The block is not necessary

[petsc-dev] Can MATAIJ support block size?

2018-09-20 Thread Fande Kong
Hi Developers, MATBAIJ actually assumes that the point-block is dense. It is fine if the block size is small for example less than 50, but we may suffer from a memory issue if the block is large such as 1000. The block is not necessary dense. In my application, neutron transport equations, the

Re: [petsc-dev] large bug in KSPSolve_Chebyshev() affects all multigrid solvers; all core developers including Mark please read this

2018-09-20 Thread Jed Brown
"Smith, Barry F." writes: > Problem 2) Less catastrophic > > When cheb->kspest is set the "regular" Chebyshev is also run (after the > eigenvalues are estimated). I am pretty sure this is wrong because there is > no reason to run the Chebyshev when estimating the eigenvalues? The >

Re: [petsc-dev] large bug in KSPSolve_Chebyshev() affects all multigrid solvers; all core developers including Mark please read this

2018-09-20 Thread Smith, Barry F.
Mark, Thanks for the explanation, I totally missed the > if (amatid != cheb->amatid || pmatid != cheb->pmatid || amatstate != > cheb->amatstate || pmatstate != cheb->pmatstate) { code that insures the computation of the eigenvalues is only once on each new matrix. Sorry for the

Re: [petsc-dev] large bug in KSPSolve_Chebyshev() affects all multigrid solvers; all core developers including Mark please read this

2018-09-20 Thread Patrick Sanan
2018-09-20 12:00 GMT+02:00 Mark Adams : > > > On Wed, Sep 19, 2018 at 7:44 PM Smith, Barry F. > wrote: > >> >> Look at the code in KSPSolve_Chebyshev(). >> >> Problem 1) VERY MAJOR >> >> Once you start running the eigenestimates it always runs them, this >> is because the routine begins

Re: [petsc-dev] large bug in KSPSolve_Chebyshev() affects all multigrid solvers; all core developers including Mark please read this

2018-09-20 Thread Mark Adams
On Wed, Sep 19, 2018 at 7:44 PM Smith, Barry F. wrote: > > Look at the code in KSPSolve_Chebyshev(). > > Problem 1) VERY MAJOR > > Once you start running the eigenestimates it always runs them, this is > because the routine begins with > > if (cheb->kspest) { > >but once

Re: [petsc-dev] Add dependency to configure

2018-09-20 Thread Marek Pecha
Hello, thank you a lot. Yesterday, I finished first version of OpenCV build script. Barry, thanks for you remark, however OpenCV should be compiled without dependency on OpenCL as well as CUDA etc. First version of OpenCV viewer, I am going to implement without OpenCL dependency. I’ve decided