It looks like my problem may be in using the older gcc 7.5.0 compiler.
The nans seem to go away when I compile petsc/slepc and my code on a
newer gcc 11.2 compiler.
Also, the petsc/slepc "make check" tests have errors on gcc 7.5 which
disappear on 11.2.
On 9/7/22 01:20, Jose E. Roman wrote:
DMDASetBlockFills() or DMDASetBlockFillsSparse() (they are just different ways
of providing the same information) will help you here enormously, your type of
case is exactly what they were designed for.
Barry
>
> On Sep 7, 2022, at 10:47 AM, Tu, Jiannan wrote:
>
> Barry and Hong,
>
>
Barry and Hong,
Thank you.
There are 26 components at each grid and there are not fully coupled in terms
of stiff functions. Mutual coupling is among about 6 components.
I would prefer not using matrix-free since the Jacobina is not difficult to
calculate and only up to 10 non-zeros at each
On Sun, Sep 4, 2022 at 5:50 PM Nicholas Arnold-Medabalimi <
narno...@umich.edu> wrote:
> Thank you. I'm able to access the label functions. Can you clarify what
> the Begin and End functions are doing?
>
> To clarify the behavior I am trying to achieve, I have identified the blue
> cell in the
You can run with -fp_trap to have the program stop as soon as the first Nan
or Inf appears, this can help track down why it is happening. In a debugger you
can also set the debugger to trap on floating point exceptions (syntax is
debugger dependent) to focus in on where it first happens.
> El 7 sept 2022, a las 6:18, Patrick Alken
> escribió:
>
> I sometimes get Nan output values in computed eigenvectors for the
> generalized symmetric eigenvalue problem produced by slepc. Is this a known
> issue, and is it related to the conditioning of the matrix pair (A,B)? Is
> there
