Re: [petsc-users] definition of level in PCILU

2019-01-03 Thread Smith, Barry F. via petsc-users
There is a hint in a comment in matimpl.h Note: the level of factor(i,j) is set as lvl(i,j) = min{ lvl(i,j), lvl(i,prow)+lvl(prow,j)+1) So it uses the "sum" rule. Barry > On Jan 3, 2019, at 9:38 AM, Tim Steinhoff via petsc-users > wrote: > > Dear PETSc Team > > I checked

[petsc-users] Set diagonals other than main diagonal

2019-01-03 Thread Sajid Ali via petsc-users
Hi, Is there any simple way of setting minor diagonals. The main diagonal can be set by MatDiagonalSet but there's no equivalent way of doing it for for minor diagonals. Does the preferred way to do this involve using MatSetValues or there a simpler way? Thank You, Sajid Ali Applied Physics

Re: [petsc-users] MUMPS Error

2019-01-03 Thread Zhang, Hong via petsc-users
You may try different matrix orderings, or try superlu_dist. Hong On Thu, Jan 3, 2019 at 7:39 AM Matthew Knepley via petsc-users mailto:petsc-users@mcs.anl.gov>> wrote: On Thu, Jan 3, 2019 at 5:58 AM Sal Am via petsc-users mailto:petsc-users@mcs.anl.gov>> wrote: [0]PETSC ERROR:

[petsc-users] Check vector for zeros

2019-01-03 Thread Florian Lindner via petsc-users
Hello, Happy New Year Everybody! I get a vector from a linear solve and is used a divisor in VecPointwiseDivide. Clearly, I must check for zero entries before doing the division. What is the best way to do so, especially performance wise? The only way I come up with so far, is to VecCopy,

Re: [petsc-users] Problems about Picard and NolinearGS

2019-01-03 Thread Matthew Knepley via petsc-users
On Thu, Jan 3, 2019 at 7:36 AM Yingjie Wu via petsc-users < petsc-users@mcs.anl.gov> wrote: > Thanks for your reply. > I read the article you provided. This is my first contact with the > quasi-Newton method. > I have some problem: > 1. From the point of view of algorithm, the quasi-Newton method

Re: [petsc-users] definition of level in PCILU

2019-01-03 Thread Smith, Barry F. via petsc-users
There is a hint in a comment in matimpl.h Note: the level of factor(i,j) is set as lvl(i,j) = min{ lvl(i,j), lvl(i,prow)+lvl(prow,j)+1) So it is the sum approach. Barry > On Jan 3, 2019, at 9:38 AM, Tim Steinhoff via petsc-users > wrote: > > Dear PETSc Team > > I checked

[petsc-users] Question about correctly catching fp_trap

2019-01-03 Thread Sajid Ali via petsc-users
Hi, I've compiled a program using petsc with debugging enabled. Confirmation of the same : https://pastebin.com/aa0XDheD (check if the loaded module is indeed petsc+debug, compile) When I run it in a debugger, I get the error as shown below: https://pastebin.com/GJtB2Ghz What am I missing?

Re: [petsc-users] Set diagonals other than main diagonal

2019-01-03 Thread Sajid Ali via petsc-users
Got it. Thank you! On Thu, Jan 3, 2019 at 11:46 AM Matthew Knepley wrote: > On Thu, Jan 3, 2019 at 12:24 PM Sajid Ali via petsc-users < > petsc-users@mcs.anl.gov> wrote: > >> Hi, >> >> Is there any simple way of setting minor diagonals. The main diagonal can >> be set by MatDiagonalSet but

Re: [petsc-users] Check vector for zeros

2019-01-03 Thread Jed Brown via petsc-users
Florian Lindner via petsc-users writes: > I get a vector from a linear solve and is used a divisor in > VecPointwiseDivide. Clearly, I must check for zero entries before doing the > division. Note that VecPointwiseDivide does not divide by values that are identically zero. for (i=0; i What

Re: [petsc-users] Installation error on macOS Mojave using GNU compiler

2019-01-03 Thread Balay, Satish via petsc-users
On Thu, 3 Jan 2019, Matthew Knepley via petsc-users wrote: > On Thu, Jan 3, 2019 at 7:02 PM Danyang Su via petsc-users < > petsc-users@mcs.anl.gov> wrote: > > > Hi All, > > > > I am trying to install PETSc on macOS Mojave using GNU compiler. > > First, I tried the debug version using the

Re: [petsc-users] Question about correctly catching fp_trap

2019-01-03 Thread Smith, Barry F. via petsc-users
This appears to be more a spack question than a PETSc question. That make (that doesn't have the -g) is controlled by spack, not PETSc. Barry > On Jan 3, 2019, at 6:15 PM, Sajid Ali via petsc-users > wrote: > > Hi, > > I've compiled a program using petsc with debugging enabled. >

Re: [petsc-users] Question about correctly catching fp_trap

2019-01-03 Thread Balay, Satish via petsc-users
I guess the primary bug here is 'petsc+debug' option is broken in spack. Its not clear to me if a package in spack should have 'debug' option. [so I should remove this?] since spack way is perhaps: spack install CFLAGS=-g FFLAGS=-g CXXFLAGS=-g petsc Satish On Fri, 4 Jan 2019, Smith, Barry F.

Re: [petsc-users] Question about correctly catching fp_trap

2019-01-03 Thread Smith, Barry F. via petsc-users
> On Jan 3, 2019, at 8:44 PM, Balay, Satish wrote: > > I guess the primary bug here is 'petsc+debug' option is broken in spack. > > Its not clear to me if a package in spack should have 'debug' option. [so I > should remove this?] Question you should ask the spack folks. Many systems