search

ForestClaw / forestclaw

Quadtree/octree adaptive PDE solver based based on p4est.
http://www.forestclaw.org
BSD 2-Clause "Simplified" License
58 stars 21 forks source link

Floating point overflow in one test #161

Closed tim-griesbach closed 1 year ago

tim-griesbach commented 2 years ago

Description

On my system, a test fails when I call make check for the current develop branch (d618b0dae87396fcf261b26242ee67f192aff0ab). applications/clawpack/shallow/2d/radialdam/regressions.sh fails with a floating point overflow. The details are in the test suite output below. I used gcc (Ubuntu 7.5.0-3ubuntu1~18.04) 7.5.0 as compiler.

Configuration

../forestclaw/configure \
    --enable-mpi \
    --enable-clawpack \
    --without-blas \
    --enable-debug \
    CC=mpicc \
    CFLAGS="-O2 -Wall" \
    CXX=mpicxx \
    CXXFLAGS="-O2 -Wall" \
    F77=gfortran \
    FC=gfortran \
    FFLAGS="-O2 -cpp" \
    --disable-shared \
    LIBS="-lgfortran"

Test Suite Output

=======================================================
   ForestClaw 0.1.6358-d618b-dirty: ./test-suite.log
=======================================================

# TOTAL: 24
# PASS:  23
# SKIP:  0
# XFAIL: 0
# FAIL:  1
# XPASS: 0
# ERROR: 0

.. contents:: :depth: 2

FAIL: applications/clawpack/shallow/2d/radialdam/regressions.sh
===============================================================

[libsc] This is libsc 2.8.2.85-1ae8
[libsc] CPP                      mpicc -E
[libsc] CPPFLAGS                 
[libsc] CC                       mpicc
[libsc] CFLAGS                   -O2 -Wall
[libsc] LDFLAGS                  
[libsc] LIBS                     -lz -lm -lgfortran
[p4est] This is p4est 2.8.7-6c4d-dirty
[p4est] CPP                      mpicc -E
[p4est] CPPFLAGS                 
[p4est] CC                       mpicc
[p4est] CFLAGS                   -O2 -Wall
[p4est] LDFLAGS                  
[p4est] LIBS                     -lz -lm -lgfortran
[fclaw] This is ForestClaw 0.1.6358-d618b-dirty
[fclaw] CPP                      mpicc -E
[fclaw] CPPFLAGS                 
[fclaw] F77                      gfortran
[fclaw] FFLAGS                   -O2 -cpp
[fclaw] CC                       mpicc
[fclaw] CFLAGS                   -O2 -Wall
[fclaw] CXX                      mpicxx
[fclaw] CXXFLAGS                 -O2 -Wall
[fclaw] LDFLAGS                  
[fclaw] FLIBS                     -L/usr/lib/gcc/x86_64-linux-gnu/7 -L/usr/lib/gcc/x86_64-linux-gnu/7/../../../x86_64-linux-gnu -L/usr/lib/gcc/x86_64-linux-gnu/7/../../../../lib -L/lib/x86_64-linux-gnu -L/lib/../lib -L/usr/lib/x86_64-linux-gnu -L/usr/lib/../lib -L/usr/lib/gcc/x86_64-linux-gnu/7/../../.. -lgfortran -lm -lquadmath
[fclaw] LIBS                     -lz -lm -lgfortran
[fclaw] Reading file fclaw_options.ini.
[fclaw] Options:
[fclaw]    help                                 false
[fclaw]    version                              false
[fclaw]    print-options                        false
[fclaw]    verbosity                            production
[fclaw]    lib-verbosity                        essential
[fclaw]    initial_dt                           0.0625
[fclaw]    max_cfl                              1
[fclaw]    desired_cfl                          0.9
[fclaw]    reduce-cfl                           true
[fclaw]    use_fixed_dt                         false
[fclaw]    outstyle                             1
[fclaw]    tfinal                               0.4
[fclaw]    nout                                 1
[fclaw]    nstep                                1
[fclaw]    advance-one-step                     false
[fclaw]    outstyle-uses-maxlevel               false
[fclaw]    subcycle                             true
[fclaw]    weighted_partition                   true
[fclaw]    time-sync                            true
[fclaw]    flux-correction                      true
[fclaw]    fluctuation-correction               true
[fclaw]    output                               false
[fclaw]    output-gauges                        false
[fclaw]    gauge-buffer-length                  1
[fclaw]    tikz-out                             false
[fclaw]    tikz-figsize                         8 6
[fclaw]    tikz-plot-prefix                     plot
[fclaw]    tikz-plot-suffix                     png
[fclaw]    tikz-mesh-only                       false
[fclaw]    tikz-plot-fig                        true
[fclaw]    prefix                               fort
[fclaw]    vtkspace                             0
[fclaw]    init_ghostcell                       false
[fclaw]    minlevel                             1
[fclaw]    maxlevel                             5
[fclaw]    regrid_interval                      1
[fclaw]    refratio                             2
[fclaw]    smooth-refine                        true
[fclaw]    smooth-level                         4
[fclaw]    coarsen-delay                        0
[fclaw]    refine_threshold                     0.5
[fclaw]    coarsen_threshold                    0.25
[fclaw]    run-user-diagnostics                 false
[fclaw]    compute-error                        false
[fclaw]    conservation-check                   true
[fclaw]    report-timing                        true
[fclaw]    report-timing-verbosity              wall
[fclaw]    ghost_patch_pack_area                true
[fclaw]    ghost_patch_pack_extra               false
[fclaw]    ghost_patch_pack_numextrafields      0
[fclaw]    trapfpe                              true
[fclaw]    mpi_debug                            false
[fclaw]    ax                                   -2.5
[fclaw]    bx                                   2.5
[fclaw]    ay                                   -2.5
[fclaw]    by                                   2.5
[fclaw]    az                                   0
[fclaw]    bz                                   1
[fclaw]    manifold                             false
[fclaw]    mi                                   1
[fclaw]    mj                                   1
[fclaw]    periodic_x                           false
[fclaw]    periodic_y                           false
[fclaw]    scale                                1 1 1
[fclaw]    shift                                0 0 0
[fclaw]    phi                                  0
[fclaw]    theta                                0
[fclaw]    clawpatch:mx                         16
[fclaw]    clawpatch:my                         16
[fclaw]    clawpatch:maux                       0
[fclaw]    clawpatch:mbc                        2
[fclaw]    clawpatch:meqn                       3
[fclaw]    clawpatch:rhs-fields                 1
[fclaw]    clawpatch:interp_stencil_width       3
[fclaw]    clawpatch:ghost_patch_pack_aux       true
[fclaw]    clawpatch:refinement-criteria        gradient
[fclaw]    clawpack46:order                     2 2
[fclaw]    clawpack46:mcapa                     0
[fclaw]    clawpack46:src_term                  false
[fclaw]    clawpack46:use-fwaves                false
[fclaw]    clawpack46:mwaves                    3
[fclaw]    clawpack46:mthlim                    4 4 4
[fclaw]    clawpack46:mthbc                     1 1 1 1
[fclaw]    clawpack46:ascii-out                 true
[fclaw]    clawpack46:vtk-out                   false
[fclaw]    clawpack5:order                      2 2
[fclaw]    clawpack5:mcapa                      0
[fclaw]    clawpack5:src_term                   false
[fclaw]    clawpack5:use-fwaves                 false
[fclaw]    clawpack5:mwaves                     3
[fclaw]    clawpack5:mthlim                     4 4 4
[fclaw]    clawpack5:mthbc                      1 1 1 1
[fclaw]    clawpack5:ascii-out                  true
[fclaw]    clawpack5:vtk-out                    false
[fclaw]    user:example                         0
[fclaw]    user:g                               1
[fclaw]    user:x0                              0
[fclaw]    user:y0                              0
[fclaw]    user:r0                              0.5
[fclaw]    user:hin                             2
[fclaw]    user:hout                            1
[fclaw]    user:alpha                           0.4
[fclaw]    user:claw-version                    4
[fclaw] Arguments: none
[p4est 1] Local minimum/maximum levels:  1  1
[p4est 0] Local minimum/maximum levels:  1  1
[p4est] Global minimum/maximum levels:  1  1
[p4est 0] Patches on level  1:         2
[fclaw] Max threads set to 0
[p4est 1] Patches on level  1:         2
[fclaw] sum[0] =    2.5785351818898182e+01    0.0000000000000000e+00
[fclaw] sum[1] =    0.0000000000000000e+00    0.0000000000000000e+00
[fclaw] sum[2] =    0.0000000000000000e+00    0.0000000000000000e+00
[fclaw] Time step    0  at time   0.00000000e+00

[fclaw] Level 1 (3-5) step     1 : dt =    6.250e-02; maxcfl (step) =    0.648; Final time =       0.0625
[fclaw] Level 1 (3-5) step     2 : dt =    8.681e-02; maxcfl (step) =    0.905; Final time =       0.1493
[fclaw] Level 1 (3-5) step     3 : dt =    8.628e-02; maxcfl (step) =    0.903; Final time =       0.2356
[fclaw] Level 1 (3-5) step     4 : dt =    8.604e-02; maxcfl (step) =    0.900; Final time =       0.3216
[fclaw] Level 1 (3-5) step     5 : dt =    7.837e-02; maxcfl (step) =    0.814; Final time =       0.4000
[fclaw] sum[0] =    2.5785351818898182e+01    0.0000000000000000e+00
[fclaw] sum[1] =    7.2123498963688636e-19    7.2123498963688636e-19
[fclaw] sum[2] =    0.0000000000000000e+00    0.0000000000000000e+00
[fclaw] Time step    1  at time   4.00000000e-01

[fclaw] Finalizing run
[libsc] Statistics for   WALLTIME
[libsc]    Global number of values:       2
[libsc]    Mean value (std. dev.):           5.77838 (1.46e-07 = 2.53e-06%)
[libsc]    Minimum attained at rank       0: 5.77838
[libsc]    Maximum attained at rank       1: 5.77838
[libsc] Procs 2 advance 5 3.62285 exchange 83 0.623892 regrid 5 3 0.016771
[libsc] Max/P 2 advance 5 3.64783 exchange 83 0.72021 regrid 5 3 0.0171044
[libsc] This is libsc 2.8.2.85-1ae8
[libsc] CPP                      mpicc -E
[libsc] CPPFLAGS                 
[libsc] CC                       mpicc
[libsc] CFLAGS                   -O2 -Wall
[libsc] LDFLAGS                  
[libsc] LIBS                     -lz -lm -lgfortran
[p4est] This is p4est 2.8.7-6c4d-dirty
[p4est] CPP                      mpicc -E
[p4est] CPPFLAGS                 
[p4est] CC                       mpicc
[p4est] CFLAGS                   -O2 -Wall
[p4est] LDFLAGS                  
[p4est] LIBS                     -lz -lm -lgfortran
[fclaw] This is ForestClaw 0.1.6358-d618b-dirty
[fclaw] CPP                      mpicc -E
[fclaw] CPPFLAGS                 
[fclaw] F77                      gfortran
[fclaw] FFLAGS                   -O2 -cpp
[fclaw] CC                       mpicc
[fclaw] CFLAGS                   -O2 -Wall
[fclaw] CXX                      mpicxx
[fclaw] CXXFLAGS                 -O2 -Wall
[fclaw] LDFLAGS                  
[fclaw] FLIBS                     -L/usr/lib/gcc/x86_64-linux-gnu/7 -L/usr/lib/gcc/x86_64-linux-gnu/7/../../../x86_64-linux-gnu -L/usr/lib/gcc/x86_64-linux-gnu/7/../../../../lib -L/lib/x86_64-linux-gnu -L/lib/../lib -L/usr/lib/x86_64-linux-gnu -L/usr/lib/../lib -L/usr/lib/gcc/x86_64-linux-gnu/7/../../.. -lgfortran -lm -lquadmath
[fclaw] LIBS                     -lz -lm -lgfortran
[fclaw] Reading file fclaw_options.ini.
[fclaw] Options:
[fclaw]    help                                 false
[fclaw]    version                              false
[fclaw]    print-options                        false
[fclaw]    verbosity                            production
[fclaw]    lib-verbosity                        essential
[fclaw]    initial_dt                           0.0625
[fclaw]    max_cfl                              1
[fclaw]    desired_cfl                          0.9
[fclaw]    reduce-cfl                           true
[fclaw]    use_fixed_dt                         false
[fclaw]    outstyle                             1
[fclaw]    tfinal                               0.4
[fclaw]    nout                                 1
[fclaw]    nstep                                1
[fclaw]    advance-one-step                     false
[fclaw]    outstyle-uses-maxlevel               false
[fclaw]    subcycle                             true
[fclaw]    weighted_partition                   true
[fclaw]    time-sync                            true
[fclaw]    flux-correction                      true
[fclaw]    fluctuation-correction               true
[fclaw]    output                               false
[fclaw]    output-gauges                        false
[fclaw]    gauge-buffer-length                  1
[fclaw]    tikz-out                             false
[fclaw]    tikz-figsize                         8 6
[fclaw]    tikz-plot-prefix                     plot
[fclaw]    tikz-plot-suffix                     png
[fclaw]    tikz-mesh-only                       false
[fclaw]    tikz-plot-fig                        true
[fclaw]    prefix                               fort
[fclaw]    vtkspace                             0
[fclaw]    init_ghostcell                       false
[fclaw]    minlevel                             1
[fclaw]    maxlevel                             5
[fclaw]    regrid_interval                      1
[fclaw]    refratio                             2
[fclaw]    smooth-refine                        true
[fclaw]    smooth-level                         4
[fclaw]    coarsen-delay                        0
[fclaw]    refine_threshold                     0.5
[fclaw]    coarsen_threshold                    0.25
[fclaw]    run-user-diagnostics                 false
[fclaw]    compute-error                        false
[fclaw]    conservation-check                   true
[fclaw]    report-timing                        true
[fclaw]    report-timing-verbosity              wall
[fclaw]    ghost_patch_pack_area                true
[fclaw]    ghost_patch_pack_extra               false
[fclaw]    ghost_patch_pack_numextrafields      0
[fclaw]    trapfpe                              true
[fclaw]    mpi_debug                            false
[fclaw]    ax                                   -2.5
[fclaw]    bx                                   2.5
[fclaw]    ay                                   -2.5
[fclaw]    by                                   2.5
[fclaw]    az                                   0
[fclaw]    bz                                   1
[fclaw]    manifold                             false
[fclaw]    mi                                   1
[fclaw]    mj                                   1
[fclaw]    periodic_x                           false
[fclaw]    periodic_y                           false
[fclaw]    scale                                1 1 1
[fclaw]    shift                                0 0 0
[fclaw]    phi                                  0
[fclaw]    theta                                0
[fclaw]    clawpatch:mx                         16
[fclaw]    clawpatch:my                         16
[fclaw]    clawpatch:maux                       0
[fclaw]    clawpatch:mbc                        2
[fclaw]    clawpatch:meqn                       3
[fclaw]    clawpatch:rhs-fields                 1
[fclaw]    clawpatch:interp_stencil_width       3
[fclaw]    clawpatch:ghost_patch_pack_aux       true
[fclaw]    clawpatch:refinement-criteria        gradient
[fclaw]    clawpack46:order                     2 2
[fclaw]    clawpack46:mcapa                     0
[fclaw]    clawpack46:src_term                  false
[fclaw]    clawpack46:use-fwaves                false
[fclaw]    clawpack46:mwaves                    3
[fclaw]    clawpack46:mthlim                    4 4 4
[fclaw]    clawpack46:mthbc                     1 1 1 1
[fclaw]    clawpack46:ascii-out                 true
[fclaw]    clawpack46:vtk-out                   false
[fclaw]    clawpack5:order                      2 2
[fclaw]    clawpack5:mcapa                      0
[fclaw]    clawpack5:src_term                   false
[fclaw]    clawpack5:use-fwaves                 false
[fclaw]    clawpack5:mwaves                     3
[fclaw]    clawpack5:mthlim                     4 4 4
[fclaw]    clawpack5:mthbc                      1 1 1 1
[fclaw]    clawpack5:ascii-out                  true
[fclaw]    clawpack5:vtk-out                    false
[fclaw]    user:example                         0
[fclaw]    user:g                               1
[fclaw]    user:x0                              0
[fclaw]    user:y0                              0
[fclaw]    user:r0                              0.5
[fclaw]    user:hin                             2
[fclaw]    user:hout                            1
[fclaw]    user:alpha                           0.4
[fclaw]    user:claw-version                    5
[fclaw] Arguments: none
[p4est 1] Local minimum/maximum levels:  1  1
[p4est 1] Patches on level  1:         2
[p4est 0] Local minimum/maximum levels:  1  1
[p4est] Global minimum/maximum levels:  1  1
[p4est 0] Patches on level  1:         2
[fclaw] Max threads set to 0
[fclaw] sum[0] =    2.5785351818898175e+01    0.0000000000000000e+00
[fclaw] sum[1] =    0.0000000000000000e+00    0.0000000000000000e+00
[fclaw] sum[2] =    0.0000000000000000e+00    0.0000000000000000e+00
[fclaw] Time step    0  at time   0.00000000e+00

[fclaw] Level 1 (3-5) step     1 : dt =    6.250e-02; maxcfl (step) =    0.648; Final time =       0.0625
[fclaw] Level 1 (3-5) step     2 : dt =    8.681e-02; maxcfl (step) =    0.905; Final time =       0.1493
[fclaw] Level 1 (3-5) step     3 : dt =    8.628e-02; maxcfl (step) =    0.903; Final time =       0.2356
[fclaw] Level 1 (3-5) step     4 : dt =    8.604e-02; maxcfl (step) =    0.900; Final time =       0.3216
[fclaw] Level 1 (3-5) step     5 : dt =    7.837e-02; maxcfl (step) =    0.814; Final time =       0.4000
[fclaw] sum[0] =    2.5785351818898231e+01    5.6843418860808015e-14
[fclaw] sum[1] =    4.6222643229102135e-17    4.6222643229102135e-17
[fclaw] sum[2] =   -4.9960036108132044e-16    4.9960036108132044e-16
[fclaw] Time step    1  at time   4.00000000e-01

[fclaw] Finalizing run
[libsc] Statistics for   WALLTIME
[libsc]    Global number of values:       2
[libsc]    Mean value (std. dev.):           6.03889 (2.67e-07 = 4.41e-06%)
[libsc]    Minimum attained at rank       0: 6.03889
[libsc]    Maximum attained at rank       1: 6.03889
[libsc] Procs 2 advance 5 3.89875 exchange 83 0.565386 regrid 5 3 0.0257484
[libsc] Max/P 2 advance 5 4.16403 exchange 83 0.592378 regrid 5 3 0.0332874
[libsc] This is libsc 2.8.2.85-1ae8
[libsc] CPP                      mpicc -E
[libsc] CPPFLAGS                 
[libsc] CC                       mpicc
[libsc] CFLAGS                   -O2 -Wall
[libsc] LDFLAGS                  
[libsc] LIBS                     -lz -lm -lgfortran
[p4est] This is p4est 2.8.7-6c4d-dirty
[p4est] CPP                      mpicc -E
[p4est] CPPFLAGS                 
[p4est] CC                       mpicc
[p4est] CFLAGS                   -O2 -Wall
[p4est] LDFLAGS                  
[p4est] LIBS                     -lz -lm -lgfortran
[fclaw] This is ForestClaw 0.1.6358-d618b-dirty
[fclaw] CPP                      mpicc -E
[fclaw] CPPFLAGS                 
[fclaw] F77                      gfortran
[fclaw] FFLAGS                   -O2 -cpp
[fclaw] CC                       mpicc
[fclaw] CFLAGS                   -O2 -Wall
[fclaw] CXX                      mpicxx
[fclaw] CXXFLAGS                 -O2 -Wall
[fclaw] LDFLAGS                  
[fclaw] FLIBS                     -L/usr/lib/gcc/x86_64-linux-gnu/7 -L/usr/lib/gcc/x86_64-linux-gnu/7/../../../x86_64-linux-gnu -L/usr/lib/gcc/x86_64-linux-gnu/7/../../../../lib -L/lib/x86_64-linux-gnu -L/lib/../lib -L/usr/lib/x86_64-linux-gnu -L/usr/lib/../lib -L/usr/lib/gcc/x86_64-linux-gnu/7/../../.. -lgfortran -lm -lquadmath
[fclaw] LIBS                     -lz -lm -lgfortran
[fclaw] Reading file fclaw_options.ini.
[fclaw] Options:
[fclaw]    help                                 false
[fclaw]    version                              false
[fclaw]    print-options                        false
[fclaw]    verbosity                            production
[fclaw]    lib-verbosity                        essential
[fclaw]    initial_dt                           0.015625
[fclaw]    max_cfl                              1
[fclaw]    desired_cfl                          0.9
[fclaw]    reduce-cfl                           true
[fclaw]    use_fixed_dt                         false
[fclaw]    outstyle                             1
[fclaw]    tfinal                               0.1
[fclaw]    nout                                 1
[fclaw]    nstep                                1
[fclaw]    advance-one-step                     false
[fclaw]    outstyle-uses-maxlevel               false
[fclaw]    subcycle                             true
[fclaw]    weighted_partition                   true
[fclaw]    time-sync                            true
[fclaw]    flux-correction                      true
[fclaw]    fluctuation-correction               true
[fclaw]    output                               false
[fclaw]    output-gauges                        false
[fclaw]    gauge-buffer-length                  1
[fclaw]    tikz-out                             false
[fclaw]    tikz-figsize                         8 6
[fclaw]    tikz-plot-prefix                     plot
[fclaw]    tikz-plot-suffix                     png
[fclaw]    tikz-mesh-only                       false
[fclaw]    tikz-plot-fig                        true
[fclaw]    prefix                               fort
[fclaw]    vtkspace                             0
[fclaw]    init_ghostcell                       false
[fclaw]    minlevel                             1
[fclaw]    maxlevel                             5
[fclaw]    regrid_interval                      1
[fclaw]    refratio                             2
[fclaw]    smooth-refine                        true
[fclaw]    smooth-level                         4
[fclaw]    coarsen-delay                        0
[fclaw]    refine_threshold                     1
[fclaw]    coarsen_threshold                    0.5
[fclaw]    run-user-diagnostics                 false
[fclaw]    compute-error                        false
[fclaw]    conservation-check                   false
[fclaw]    report-timing                        true
[fclaw]    report-timing-verbosity              wall
[fclaw]    ghost_patch_pack_area                true
[fclaw]    ghost_patch_pack_extra               false
[fclaw]    ghost_patch_pack_numextrafields      0
[fclaw]    trapfpe                              true
[fclaw]    mpi_debug                            false
[fclaw]    ax                                   -2.5
[fclaw]    bx                                   2.5
[fclaw]    ay                                   -2.5
[fclaw]    by                                   2.5
[fclaw]    az                                   0
[fclaw]    bz                                   1
[fclaw]    manifold                             true
[fclaw]    mi                                   1
[fclaw]    mj                                   1
[fclaw]    periodic_x                           false
[fclaw]    periodic_y                           false
[fclaw]    scale                                2.5 2.5 0
[fclaw]    shift                                0 0 0
[fclaw]    phi                                  0
[fclaw]    theta                                0
[fclaw]    clawpatch:mx                         16
[fclaw]    clawpatch:my                         16
[fclaw]    clawpatch:maux                       16
[fclaw]    clawpatch:mbc                        2
[fclaw]    clawpatch:meqn                       4
[fclaw]    clawpatch:rhs-fields                 1
[fclaw]    clawpatch:interp_stencil_width       3
[fclaw]    clawpatch:ghost_patch_pack_aux       true
[fclaw]    clawpatch:refinement-criteria        gradient
[fclaw]    clawpack46:order                     2 2
[fclaw]    clawpack46:mcapa                     1
[fclaw]    clawpack46:src_term                  false
[fclaw]    clawpack46:use-fwaves                false
[fclaw]    clawpack46:mwaves                    3
[fclaw]    clawpack46:mthlim                    4 4 4
[fclaw]    clawpack46:mthbc                     1 1 1 1
[fclaw]    clawpack46:ascii-out                 true
[fclaw]    clawpack46:vtk-out                   false
[fclaw]    clawpack5:order                      2 2
[fclaw]    clawpack5:mcapa                      1
[fclaw]    clawpack5:src_term                   false
[fclaw]    clawpack5:use-fwaves                 false
[fclaw]    clawpack5:mwaves                     3
[fclaw]    clawpack5:mthlim                     4 4 4
[fclaw]    clawpack5:mthbc                      1 1 1 1
[fclaw]    clawpack5:ascii-out                  true
[fclaw]    clawpack5:vtk-out                    false
[fclaw]    user:example                         1
[fclaw]    user:g                               1
[fclaw]    user:x0                              0
[fclaw]    user:y0                              0
[fclaw]    user:r0                              0.5
[fclaw]    user:hin                             2
[fclaw]    user:hout                            1
[fclaw]    user:alpha                           0.4
[fclaw]    user:claw-version                    4
[fclaw] Arguments: none
[p4est 1] Local minimum/maximum levels:  1  1
[p4est 1] Patches on level  1:        10
[p4est 0] Local minimum/maximum levels:  1  1
[p4est] Global minimum/maximum levels:  1  1
[p4est 0] Patches on level  1:        10
[fclaw] Max threads set to 0
[fclaw] Time step    0  at time   0.00000000e+00

[fclaw] Level 1 (1-5) step     1 : dt =    1.562e-02; maxcfl (step) =    0.405; Final time =       0.0156
[fclaw] Level 1 (1-5) step     2 : dt =    3.473e-02; maxcfl (step) =    0.903; Final time =       0.0504
[fclaw] Level 1 (1-5) step     3 : dt =    3.462e-02; maxcfl (step) =    0.903; Final time =       0.0850
[fclaw] Level 1 (1-5) step     4 : dt =    1.503e-02; maxcfl (step) =    0.392; Final time =       0.1000
[fclaw] Time step    1  at time   1.00000000e-01

[fclaw] Finalizing run
[libsc] Statistics for   WALLTIME
[libsc]    Global number of values:       2
[libsc]    Mean value (std. dev.):           13.2112 (8.26e-07 = 6.25e-06%)
[libsc]    Minimum attained at rank       0: 13.2112
[libsc]    Maximum attained at rank       1: 13.2112
[libsc] Procs 2 advance 4 10.6842 exchange 68 0.859167 regrid 4 4 0.100387
[libsc] Max/P 2 advance 4 10.7643 exchange 68 0.867304 regrid 4 4 0.101058
[libsc] This is libsc 2.8.2.85-1ae8
[libsc] CPP                      mpicc -E
[libsc] CPPFLAGS                 
[libsc] CC                       mpicc
[libsc] CFLAGS                   -O2 -Wall
[libsc] LDFLAGS                  
[libsc] LIBS                     -lz -lm -lgfortran
[p4est] This is p4est 2.8.7-6c4d-dirty
[p4est] CPP                      mpicc -E
[p4est] CPPFLAGS                 
[p4est] CC                       mpicc
[p4est] CFLAGS                   -O2 -Wall
[p4est] LDFLAGS                  
[p4est] LIBS                     -lz -lm -lgfortran
[fclaw] This is ForestClaw 0.1.6358-d618b-dirty
[fclaw] CPP                      mpicc -E
[fclaw] CPPFLAGS                 
[fclaw] F77                      gfortran
[fclaw] FFLAGS                   -O2 -cpp
[fclaw] CC                       mpicc
[fclaw] CFLAGS                   -O2 -Wall
[fclaw] CXX                      mpicxx
[fclaw] CXXFLAGS                 -O2 -Wall
[fclaw] LDFLAGS                  
[fclaw] FLIBS                     -L/usr/lib/gcc/x86_64-linux-gnu/7 -L/usr/lib/gcc/x86_64-linux-gnu/7/../../../x86_64-linux-gnu -L/usr/lib/gcc/x86_64-linux-gnu/7/../../../../lib -L/lib/x86_64-linux-gnu -L/lib/../lib -L/usr/lib/x86_64-linux-gnu -L/usr/lib/../lib -L/usr/lib/gcc/x86_64-linux-gnu/7/../../.. -lgfortran -lm -lquadmath
[fclaw] LIBS                     -lz -lm -lgfortran
[fclaw] Reading file fclaw_options.ini.
[fclaw] Options:
[fclaw]    help                                 false
[fclaw]    version                              false
[fclaw]    print-options                        false
[fclaw]    verbosity                            production
[fclaw]    lib-verbosity                        essential
[fclaw]    initial_dt                           0.015625
[fclaw]    max_cfl                              1
[fclaw]    desired_cfl                          0.9
[fclaw]    reduce-cfl                           true
[fclaw]    use_fixed_dt                         false
[fclaw]    outstyle                             1
[fclaw]    tfinal                               0.1
[fclaw]    nout                                 1
[fclaw]    nstep                                1
[fclaw]    advance-one-step                     false
[fclaw]    outstyle-uses-maxlevel               false
[fclaw]    subcycle                             true
[fclaw]    weighted_partition                   true
[fclaw]    time-sync                            true
[fclaw]    flux-correction                      true
[fclaw]    fluctuation-correction               true
[fclaw]    output                               false
[fclaw]    output-gauges                        false
[fclaw]    gauge-buffer-length                  1
[fclaw]    tikz-out                             false
[fclaw]    tikz-figsize                         8 6
[fclaw]    tikz-plot-prefix                     plot
[fclaw]    tikz-plot-suffix                     png
[fclaw]    tikz-mesh-only                       false
[fclaw]    tikz-plot-fig                        true
[fclaw]    prefix                               fort
[fclaw]    vtkspace                             0
[fclaw]    init_ghostcell                       false
[fclaw]    minlevel                             1
[fclaw]    maxlevel                             5
[fclaw]    regrid_interval                      1
[fclaw]    refratio                             2
[fclaw]    smooth-refine                        true
[fclaw]    smooth-level                         4
[fclaw]    coarsen-delay                        0
[fclaw]    refine_threshold                     1
[fclaw]    coarsen_threshold                    0.5
[fclaw]    run-user-diagnostics                 false
[fclaw]    compute-error                        false
[fclaw]    conservation-check                   false
[fclaw]    report-timing                        true
[fclaw]    report-timing-verbosity              wall
[fclaw]    ghost_patch_pack_area                true
[fclaw]    ghost_patch_pack_extra               false
[fclaw]    ghost_patch_pack_numextrafields      0
[fclaw]    trapfpe                              true
[fclaw]    mpi_debug                            false
[fclaw]    ax                                   -2.5
[fclaw]    bx                                   2.5
[fclaw]    ay                                   -2.5
[fclaw]    by                                   2.5
[fclaw]    az                                   0
[fclaw]    bz                                   1
[fclaw]    manifold                             true
[fclaw]    mi                                   1
[fclaw]    mj                                   1
[fclaw]    periodic_x                           false
[fclaw]    periodic_y                           false
[fclaw]    scale                                2.5 2.5 0
[fclaw]    shift                                0 0 0
[fclaw]    phi                                  0
[fclaw]    theta                                0
[fclaw]    clawpatch:mx                         16
[fclaw]    clawpatch:my                         16
[fclaw]    clawpatch:maux                       16
[fclaw]    clawpatch:mbc                        2
[fclaw]    clawpatch:meqn                       4
[fclaw]    clawpatch:rhs-fields                 1
[fclaw]    clawpatch:interp_stencil_width       3
[fclaw]    clawpatch:ghost_patch_pack_aux       true
[fclaw]    clawpatch:refinement-criteria        gradient
[fclaw]    clawpack46:order                     2 2
[fclaw]    clawpack46:mcapa                     1
[fclaw]    clawpack46:src_term                  false
[fclaw]    clawpack46:use-fwaves                false
[fclaw]    clawpack46:mwaves                    3
[fclaw]    clawpack46:mthlim                    4 4 4
[fclaw]    clawpack46:mthbc                     1 1 1 1
[fclaw]    clawpack46:ascii-out                 true
[fclaw]    clawpack46:vtk-out                   false
[fclaw]    clawpack5:order                      2 2
[fclaw]    clawpack5:mcapa                      1
[fclaw]    clawpack5:src_term                   false
[fclaw]    clawpack5:use-fwaves                 false
[fclaw]    clawpack5:mwaves                     3
[fclaw]    clawpack5:mthlim                     4 4 4
[fclaw]    clawpack5:mthbc                      1 1 1 1
[fclaw]    clawpack5:ascii-out                  true
[fclaw]    clawpack5:vtk-out                    false
[fclaw]    user:example                         1
[fclaw]    user:g                               1
[fclaw]    user:x0                              0
[fclaw]    user:y0                              0
[fclaw]    user:r0                              0.5
[fclaw]    user:hin                             2
[fclaw]    user:hout                            1
[fclaw]    user:alpha                           0.4
[fclaw]    user:claw-version                    5
[fclaw] Arguments: none
[p4est 1] Local minimum/maximum levels:  1  1
[p4est 1] Patches on level  1:        10
[p4est 0] Local minimum/maximum levels:  1  1
[p4est] Global minimum/maximum levels:  1  1
[p4est 0] Patches on level  1:        10
[fclaw] Max threads set to 0
[fclaw] Time step    0  at time   0.00000000e+00

[fclaw] Level 1 (1-5) step     1 : dt =    1.562e-02; maxcfl (step) =    0.405; Final time =       0.0156
[fclaw] Level 1 (1-5) step     2 : dt =    3.473e-02; maxcfl (step) =    0.903; Final time =       0.0504
[nanna:10072] *** Process received signal ***
[nanna:10072] Signal: Floating point exception (8)
[nanna:10072] Signal code: Floating point overflow (4)
[nanna:10072] Failing at address: 0x5644b5045c26
[nanna:10072] [ 0] /lib/x86_64-linux-gnu/libpthread.so.0(+0x128a0)[0x7f47bf0df8a0]
[nanna:10072] [ 1] /home/griesbac/forestclaw_workspace/forestclaw-mpi/applications/clawpack/shallow/2d/radialdam/radialdam(+0x19c26)[0x5644b5045c26]
[nanna:10072] [ 2] /home/griesbac/forestclaw_workspace/forestclaw-mpi/applications/clawpack/shallow/2d/radialdam/radialdam(+0x283aa)[0x5644b50543aa]
[nanna:10072] [ 3] /home/griesbac/forestclaw_workspace/forestclaw-mpi/applications/clawpack/shallow/2d/radialdam/radialdam(+0x2a4ea)[0x5644b50564ea]
[nanna:10072] [ 4] /home/griesbac/forestclaw_workspace/forestclaw-mpi/applications/clawpack/shallow/2d/radialdam/radialdam(+0x28744)[0x5644b5054744]
[nanna:10072] [ 5] /home/griesbac/forestclaw_workspace/forestclaw-mpi/applications/clawpack/shallow/2d/radialdam/radialdam(+0x2325f)[0x5644b504f25f]
[nanna:10072] [ 6] /home/griesbac/forestclaw_workspace/forestclaw-mpi/applications/clawpack/shallow/2d/radialdam/radialdam(+0x551d3)[0x5644b50811d3]
[nanna:10072] [ 7] /home/griesbac/forestclaw_workspace/forestclaw-mpi/applications/clawpack/shallow/2d/radialdam/radialdam(+0x4db6a)[0x5644b5079b6a]
[nanna:10072] [ 8] /home/griesbac/forestclaw_workspace/forestclaw-mpi/applications/clawpack/shallow/2d/radialdam/radialdam(+0x4584c)[0x5644b507184c]
[nanna:10072] [ 9] /home/griesbac/forestclaw_workspace/forestclaw-mpi/applications/clawpack/shallow/2d/radialdam/radialdam(+0x55270)[0x5644b5081270]
[nanna:10072] [10] /home/griesbac/forestclaw_workspace/forestclaw-mpi/applications/clawpack/shallow/2d/radialdam/radialdam(+0x53f41)[0x5644b507ff41]
[nanna:10072] [11] /home/griesbac/forestclaw_workspace/forestclaw-mpi/applications/clawpack/shallow/2d/radialdam/radialdam(+0x54286)[0x5644b5080286]
[nanna:10072] [12] /home/griesbac/forestclaw_workspace/forestclaw-mpi/applications/clawpack/shallow/2d/radialdam/radialdam(+0x49c90)[0x5644b5075c90]
[nanna:10072] [13] /home/griesbac/forestclaw_workspace/forestclaw-mpi/applications/clawpack/shallow/2d/radialdam/radialdam(+0xe201)[0x5644b503a201]
[nanna:10072] [14] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xe7)[0x7f47becfdb97]
[nanna:10072] [15] /home/griesbac/forestclaw_workspace/forestclaw-mpi/applications/clawpack/shallow/2d/radialdam/radialdam(+0xe44a)[0x5644b503a44a]
[nanna:10072] *** End of error message ***
--------------------------------------------------------------------------
mpiexec noticed that process rank 1 with PID 0 on node nanna exited on signal 8 (Floating point exception).
--------------------------------------------------------------------------
FAIL applications/clawpack/shallow/2d/radialdam/regressions.sh (exit status: 1)
donnaaboise commented 2 years ago

Can you try this (just from the command line). 

$ cd $FCLAW/applications/clawpack/shallow/2d/radialdam
$ ln -s $FCLAW_BUILD/applications/clawpack/shallow/2d/radialdam/radialdam   # symlink to executable
$ radialdam -F regression_map.ini --user:claw-version=5 --user:example=1

If you get the same error (which I guess you will), can you compile with "-O0 -g" and run the code in gdb to see where it crashes? 

$ gdb radialdam
(gdb) set args -F regression_map.ini --user:claw-version=5 --user:example=1
(gdb) run
(gdb) bt 5

and report the results?

The code doesn't crash for me - I get : 

[fclaw] Time step    0  at time   0.00000000e+00

[fclaw] Level 1 (1-5) step     1 : dt =    1.562e-02; maxcfl (step) =    0.405; Final time =       0.0156
[fclaw] Level 1 (1-5) step     2 : dt =    3.473e-02; maxcfl (step) =    0.903; Final time =       0.0504
[fclaw] Level 1 (1-5) step     3 : dt =    3.462e-02; maxcfl (step) =    0.903; Final time =       0.0850
[fclaw] Level 1 (1-5) step     4 : dt =    1.503e-02; maxcfl (step) =    0.392; Final time =       0.1000
[fclaw] Time step    1  at time   1.00000000e-01

So there could be a memory issue with the corner cells in the five patch grid.

tim-griesbach commented 2 years ago

Can you try this (just from the command line).

$ cd $FCLAW/applications/clawpack/shallow/2d/radialdam
$ ln -s $FCLAW_BUILD/applications/clawpack/shallow/2d/radialdam/radialdam   # >symlink to executable
$ radialdam -F regression_map.ini --user:claw-version=5 --user:example=1

Yes, I have tried it but it results in the same error.

If you get the same error (which I guess you will), can you compile with "-O0 -g" and run the code in gdb to see where it crashes?

$ gdb radialdam
(gdb) set args -F regression_map.ini --user:claw-version=5 --user:example=1
(gdb) run
(gdb) bt 5

and report the results?

Sure, I obtain

[fclaw] Level 1 (1-5) step     1 : dt =    1.562e-02; maxcfl (step) =    0.405; Final time =       0.0156
[fclaw] Level 1 (1-5) step     2 : dt =    3.473e-02; maxcfl (step) =    0.903; Final time =       0.0504

Thread 1 "radialdam" received signal SIGFPE, Arithmetic exception.
0x000055555556dc26 in clawpack5_rpt2_manifold (ixy=<optimized out>, imp=<optimized out>, maxm=<optimized out>, meqn=<optimized out>, mwaves=<optimized out>, maux=<optimized out>, mbc=2, mx=16, ql=..., qr=..., aux1=..., aux2=..., aux3=..., 
    asdq=<error reading variable: Cannot access memory at address 0x4>, bmasdq=<error reading variable: Cannot access memory at address 0x4>, bpasdq=<error reading variable: Cannot access memory at address 0x4>)
    at ../forestclaw/applications/clawpack/shallow/2d/rp/clawpack5_rpt2_manifold.f90:253
253                 + dmin1(s(mw,i), 0.d0)*wave(m,mw,i)
(gdb) bt 5
#0  0x000055555556dc26 in clawpack5_rpt2_manifold (ixy=<optimized out>, imp=<optimized out>, maxm=<optimized out>, meqn=<optimized out>, mwaves=<optimized out>, maux=<optimized out>, mbc=2, mx=16, ql=..., qr=..., aux1=..., aux2=..., aux3=..., 
    asdq=<error reading variable: Cannot access memory at address 0x4>, bmasdq=<error reading variable: Cannot access memory at address 0x4>, bpasdq=<error reading variable: Cannot access memory at address 0x4>)
    at ../forestclaw/applications/clawpack/shallow/2d/rp/clawpack5_rpt2_manifold.f90:253
#1  0x000055555557c3aa in clawpack5_flux2_ ()
#2  0x000055555557e4ea in clawpack5_step2 (maxm=16, meqn=4, maux=16, mbc=2, mx=<optimized out>, my=16, qold=..., aux=..., dx=0.001953125, dy=0.001953125, dt=0.0021636445380113938, cflgrid=0.81021702105952675, fm=..., fp=..., gm=..., gp=..., 
    rpn2=0x55555556bba0 <clawpack5_rpn2_manifold>, rpt2=0x55555556c950 <clawpack5_rpt2_manifold>, block_corner_count=...) at ../forestclaw/src/solvers/fc2d_clawpack5/fortran_source/clawpack5_step2.f90:158
#3  0x000055555557c744 in clawpack5_step2_wrap_ ()
#4  0x000055555557725f in clawpack5_update(fclaw2d_global*, fclaw2d_patch*, int, int, double, double, void*) ()
(More stack frames follow...)
(gdb) bt
#0  0x000055555556dc26 in clawpack5_rpt2_manifold (ixy=<optimized out>, imp=<optimized out>, maxm=<optimized out>, meqn=<optimized out>, mwaves=<optimized out>, maux=<optimized out>, mbc=2, mx=16, ql=..., qr=..., aux1=..., aux2=..., aux3=..., 
    asdq=<error reading variable: Cannot access memory at address 0x4>, bmasdq=<error reading variable: Cannot access memory at address 0x4>, bpasdq=<error reading variable: Cannot access memory at address 0x4>)
    at ../forestclaw/applications/clawpack/shallow/2d/rp/clawpack5_rpt2_manifold.f90:253
#1  0x000055555557c3aa in clawpack5_flux2_ ()
#2  0x000055555557e4ea in clawpack5_step2 (maxm=16, meqn=4, maux=16, mbc=2, mx=<optimized out>, my=16, qold=..., aux=..., dx=0.001953125, dy=0.001953125, dt=0.0021636445380113938, cflgrid=0.81021702105952675, fm=..., fp=..., gm=..., gp=..., 
    rpn2=0x55555556bba0 <clawpack5_rpn2_manifold>, rpt2=0x55555556c950 <clawpack5_rpt2_manifold>, block_corner_count=...) at ../forestclaw/src/solvers/fc2d_clawpack5/fortran_source/clawpack5_step2.f90:158
#3  0x000055555557c744 in clawpack5_step2_wrap_ ()
#4  0x000055555557725f in clawpack5_update(fclaw2d_global*, fclaw2d_patch*, int, int, double, double, void*) ()
#5  0x00005555555a91d3 in cb_single_step ()
#6  0x00005555555a1b6a in fclaw2d_domain_iterate_level ()
#7  0x000055555559984c in fclaw2d_global_iterate_level ()
#8  0x00005555555a9270 in fclaw2d_update_single_step ()
#9  0x00005555555a7f41 in advance_level ()
#10 0x00005555555a8286 in fclaw2d_advance_all_levels ()
#11 0x000055555559dc90 in fclaw2d_run ()
#12 0x0000555555562201 in main ()
donnaaboise commented 2 years ago

I may have to look at this a bit more, but in the meantime, can you print out the value of i from gdb? It may not let you - since it looks like lots of values have been optimized out.

What OS are you running in?

tim-griesbach commented 2 years ago

I may have to look at this a bit more, but in the meantime, can you print out the value >of i from gdb? It may not let you - since it looks like lots of values have been >optimized out.

Unfortunately, the value of i is not printed by gdb. As you expected I get <optimized out> instead of a value.

What OS are you running in?

I am using Ubuntu 18.04.3 LTS as OS.

donnaaboise commented 2 years ago

One last thing to try is to set --trapfpe=F. Does the code run successfully?

I can't reproduce the error, but am pretty sure I know where it is coming from.

In this particular example, the mapping is done on the "fivepatch" grid - a squared-disk mapping made up of five Cartesian blocks. The interior corner ghost cell values (where three patches meet) may have nan values and although these values may be involved in arithmetic expressions, their results don't normally propagate into the solution. If a calculation touches these values, the compiler may trigger an exception.

Maybe this is compiler dependent? I am compiling with gcc v11. Would you be able to try with a later version of gcc?

tim-griesbach commented 2 years ago

One last thing to try is to set --trapfpe=F. Does the code run successfully?

Yes, this flag fixes the problem and the test runs successfully.

[fclaw] Level 1 (1-5) step     1 : dt =    1.562e-02; maxcfl (step) =    0.405; Final time =       0.0156
[fclaw] Level 1 (1-5) step     2 : dt =    3.473e-02; maxcfl (step) =    0.903; Final time =       0.0504
[fclaw] Level 1 (1-5) step     3 : dt =    3.462e-02; maxcfl (step) =    0.903; Final time =       0.0850
[fclaw] Level 1 (1-5) step     4 : dt =    1.503e-02; maxcfl (step) =    0.392; Final time =       0.1000
[fclaw] Time step    1  at time   1.00000000e-01

[fclaw] Finalizing run
[libsc] Statistics for   WALLTIME
[libsc]    Global number of values:       1
[libsc]    Mean value (std. dev.):           23.9439 (0 = 0%)
[libsc]    Minimum attained at rank       0: 23.9439
[libsc]    Maximum attained at rank       0: 23.9439
[libsc] Procs 1 advance 4 21.2158 exchange 68 1.79465 regrid 4 4 0.328972
[libsc] Max/P 1 advance 4 21.2158 exchange 68 1.79465 regrid 4 4 0.328972

Maybe this is compiler dependent? I am compiling with gcc v11. Would you be able to try with a later version of gcc?

I will also try a newer version of gcc and will then report if --trapfpe=F is still necessary.

donnaaboise commented 2 years ago

Excellent! Go ahead and update the regression test (regressions.sh) by including --trapfpe=F and issuing a PR. I'd also be curious to know if the later compiler also complains about this. These FPEs are also likely hardware dependent, so probably we should just make sure that all the tests run with FPE traps turned off.

tim-griesbach commented 2 years ago

Excellent! Go ahead and update the regression test (regressions.sh) by including --trapfpe=F and issuing a PR.

I created a PR (https://github.com/ForestClaw/forestclaw/pull/164) with the adjusted regression test.

I'd also be curious to know if the later compiler also complains about this.

For gcc 10.2.0 running the tests result in the same problem as before if I run the tests by make check-TESTS. However, using your earlier suggestion

Can you try this (just from the command line).

$ cd $FCLAW/applications/clawpack/shallow/2d/radialdam
$ ln -s $FCLAW_BUILD/applications/clawpack/shallow/2d/radialdam/radialdam   # symlink to executable
$ radialdam -F regression_map.ini --user:claw-version=5 --user:example=1

the radialdam test works even without --trapfpe=F. Do you know why this help?

./radialdam -F regression_map.ini --user:claw-version=5 --user:example=1
[libsc] This is libsc 2.8.2.85-1ae8
[libsc] CPP                      mpicc -E
[libsc] CPPFLAGS                 
[libsc] CC                       mpicc
[libsc] CFLAGS                   -O0 -g -Wall
[libsc] LDFLAGS                  
[libsc] LIBS                     -lz -lm -lgfortran
[p4est] This is p4est 2.8.7-6c4d
[p4est] CPP                      mpicc -E
[p4est] CPPFLAGS                 
[p4est] CC                       mpicc
[p4est] CFLAGS                   -O0 -g -Wall
[p4est] LDFLAGS                  
[p4est] LIBS                     -lz -lm -lgfortran
[fclaw] This is ForestClaw 0.1.6365-a05d-dirty
[fclaw] CPP                      mpicc -E
[fclaw] CPPFLAGS                 
[fclaw] F77                      gfortran
[fclaw] FFLAGS                   -O0 -g -cpp
[fclaw] CC                       mpicc
[fclaw] CFLAGS                   -O0 -g -Wall
[fclaw] CXX                      mpicxx
[fclaw] CXXFLAGS                 -O0 -g -Wall
[fclaw] LDFLAGS                  
[fclaw] FLIBS                     -L/opt/software/easybuild/software/binutils/2.37/lib/../lib64 -L/opt/software/easybuild/software/zlib/1.2.11-GCCcore-10.2.0/lib/../lib64 -L/opt/software/easybuild/software/GCCcore/10.2.0/lib/gcc/x86_64-pc-linux-gnu/10.2.0 -L/opt/software/easybuild/software/GCCcore/10.2.0/lib/gcc/x86_64-pc-linux-gnu/10.2.0/../../../../lib64 -L/lib/x86_64-linux-gnu -L/lib/../lib64 -L/usr/lib/x86_64-linux-gnu -L/opt/software/easybuild/software/binutils/2.37/lib -L/opt/software/easybuild/software/zlib/1.2.11-GCCcore-10.2.0/lib -L/opt/software/easybuild/software/GCCcore/10.2.0/lib/gcc/x86_64-pc-linux-gnu/10.2.0/../../.. -lgfortran -lm -lquadmath
[fclaw] LIBS                     -lz -lm -lgfortran
[fclaw] Reading file fclaw_options.ini.
[fclaw] Options:
[fclaw]    help                                 false
[fclaw]    version                              false
[fclaw]    print-options                        false
[fclaw]    verbosity                            production
[fclaw]    lib-verbosity                        essential
[fclaw]    initial_dt                           0.015625
[fclaw]    max_cfl                              1
[fclaw]    desired_cfl                          0.9
[fclaw]    reduce-cfl                           true
[fclaw]    use_fixed_dt                         false
[fclaw]    outstyle                             1
[fclaw]    tfinal                               0.1
[fclaw]    nout                                 1
[fclaw]    nstep                                1
[fclaw]    advance-one-step                     false
[fclaw]    outstyle-uses-maxlevel               false
[fclaw]    subcycle                             true
[fclaw]    weighted_partition                   true
[fclaw]    time-sync                            true
[fclaw]    flux-correction                      true
[fclaw]    fluctuation-correction               true
[fclaw]    output                               false
[fclaw]    output-gauges                        false
[fclaw]    gauge-buffer-length                  1
[fclaw]    tikz-out                             false
[fclaw]    tikz-figsize                         8 6
[fclaw]    tikz-plot-prefix                     plot
[fclaw]    tikz-plot-suffix                     png
[fclaw]    tikz-mesh-only                       false
[fclaw]    tikz-plot-fig                        true
[fclaw]    prefix                               fort
[fclaw]    vtkspace                             0
[fclaw]    init_ghostcell                       false
[fclaw]    minlevel                             1
[fclaw]    maxlevel                             5
[fclaw]    regrid_interval                      1
[fclaw]    refratio                             2
[fclaw]    smooth-refine                        true
[fclaw]    smooth-level                         4
[fclaw]    coarsen-delay                        0
[fclaw]    refine_threshold                     1
[fclaw]    coarsen_threshold                    0.5
[fclaw]    run-user-diagnostics                 false
[fclaw]    compute-error                        false
[fclaw]    conservation-check                   false
[fclaw]    report-timing                        true
[fclaw]    report-timing-verbosity              wall
[fclaw]    ghost_patch_pack_area                true
[fclaw]    ghost_patch_pack_extra               false
[fclaw]    ghost_patch_pack_numextrafields      0
[fclaw]    trapfpe                              true
[fclaw]    mpi_debug                            false
[fclaw]    ax                                   -2.5
[fclaw]    bx                                   2.5
[fclaw]    ay                                   -2.5
[fclaw]    by                                   2.5
[fclaw]    az                                   0
[fclaw]    bz                                   1
[fclaw]    manifold                             true
[fclaw]    mi                                   1
[fclaw]    mj                                   1
[fclaw]    periodic_x                           false
[fclaw]    periodic_y                           false
[fclaw]    scale                                2.5 2.5 0
[fclaw]    shift                                0 0 0
[fclaw]    phi                                  0
[fclaw]    theta                                0
[fclaw]    clawpatch:mx                         16
[fclaw]    clawpatch:my                         16
[fclaw]    clawpatch:maux                       16
[fclaw]    clawpatch:mbc                        2
[fclaw]    clawpatch:meqn                       4
[fclaw]    clawpatch:rhs-fields                 1
[fclaw]    clawpatch:interp_stencil_width       3
[fclaw]    clawpatch:ghost_patch_pack_aux       true
[fclaw]    clawpatch:refinement-criteria        gradient
[fclaw]    clawpack46:order                     2 2
[fclaw]    clawpack46:mcapa                     1
[fclaw]    clawpack46:src_term                  false
[fclaw]    clawpack46:use-fwaves                false
[fclaw]    clawpack46:mwaves                    3
[fclaw]    clawpack46:mthlim                    4 4 4
[fclaw]    clawpack46:mthbc                     1 1 1 1
[fclaw]    clawpack46:ascii-out                 true
[fclaw]    clawpack46:vtk-out                   false
[fclaw]    clawpack5:order                      2 2
[fclaw]    clawpack5:mcapa                      1
[fclaw]    clawpack5:src_term                   false
[fclaw]    clawpack5:use-fwaves                 false
[fclaw]    clawpack5:mwaves                     3
[fclaw]    clawpack5:mthlim                     4 4 4
[fclaw]    clawpack5:mthbc                      1 1 1 1
[fclaw]    clawpack5:ascii-out                  true
[fclaw]    clawpack5:vtk-out                    false
[fclaw]    user:example                         1
[fclaw]    user:g                               1
[fclaw]    user:x0                              0
[fclaw]    user:y0                              0
[fclaw]    user:r0                              0.5
[fclaw]    user:hin                             2
[fclaw]    user:hout                            1
[fclaw]    user:alpha                           0.4
[fclaw]    user:claw-version                    5
[fclaw] Arguments: none
[p4est 0] Local minimum/maximum levels:  1  1
[p4est] Global minimum/maximum levels:  1  1
[p4est 0] Patches on level  1:        20
[fclaw] Max threads set to 0
[fclaw] Time step    0  at time   0.00000000e+00

[fclaw] Level 1 (1-5) step     1 : dt =    1.562e-02; maxcfl (step) =    0.405; Final time =       0.0156
[fclaw] Level 1 (1-5) step     2 : dt =    3.473e-02; maxcfl (step) =    0.903; Final time =       0.0504
[fclaw] Level 1 (1-5) step     3 : dt =    3.462e-02; maxcfl (step) =    0.903; Final time =       0.0850
[fclaw] Level 1 (1-5) step     4 : dt =    1.503e-02; maxcfl (step) =    0.392; Final time =       0.1000
[fclaw] Time step    1  at time   1.00000000e-01

[fclaw] Finalizing run
[libsc] Statistics for   WALLTIME
[libsc]    Global number of values:       1
[libsc]    Mean value (std. dev.):           24.0652 (0 = 0%)
[libsc]    Minimum attained at rank       0: 24.0652
[libsc]    Maximum attained at rank       0: 24.0652
[libsc] Procs 1 advance 4 21.3319 exchange 68 1.80817 regrid 4 4 0.340671
[libsc] Max/P 1 advance 4 21.3319 exchange 68 1.80817 regrid 4 4 0.340671

These FPEs are also likely hardware dependent, so probably we should just make sure that all the tests run with FPE traps turned off.

Okay, I could add this in the pull request if you like.

donnaaboise commented 2 years ago

The tests are run in parallel, where the order in which ghost cells are filled may be slightly different than in serial.

Does this work? 

$ mpirun -n 2 radialdam -F regression_map.ini --user:claw-version=5 --user:example=1
tim-griesbach commented 2 years ago

Does this work?

$ mpirun -n 2 radialdam -F regression_map.ini --user:claw-version=5 --user:example=1

No, using gcc 10.2.0 I get the same error as before.

donnaaboise commented 2 years ago

I think just turning off the exceptions for this problem is still the best solution (I merged in your PR - thanks). Eventually, though, it would be interesting to track down exactly where the error is happening. I suspect that it only happens with dynamic adaption where new grids are created throughout the run. To test this idea, you could try running with minlevel=maxlevel=6 (might slow down the code considerably, though : 

mpirun -n 2 radialdam -F regression_map.ini --user:claw-version=5 --user:example=1 --minlevel=6 --maxlevel=6

If that takes too long, you can also just try setting minlevel=maxlevel=3 (or 4 or 5).

tim-griesbach commented 2 years ago

I think just turning off the exceptions for this problem is still the best solution (I merged in your PR - thanks). Eventually, though, it would be interesting to track down exactly where the error is happening. I suspect that it only happens with dynamic adaption where new grids are created throughout the run. To test this idea, you could try running with minlevel=maxlevel=6 (might slow down the code considerably, though :

mpirun -n 2 radialdam -F regression_map.ini --user:claw-version=5 --user:example=1 --minlevel=6 --maxlevel=6

If that takes too long, you can also just try setting minlevel=maxlevel=3 (or 4 or 5).

For minlevel=maxlevel=3 the test runs successfully but for minlevel=maxlevel \in {4,5,6} the test results again in a floating point error.

minlevel=maxlevel=3

[libsc] This is libsc 2.8.2.85-1ae8
[libsc] CPP                      mpicc -E
[libsc] CPPFLAGS                 
[libsc] CC                       mpicc
[libsc] CFLAGS                   -O0 -g -Wall
[libsc] LDFLAGS                  
[libsc] LIBS                     -lz -lm -lgfortran
[p4est] This is p4est 2.8.7-6c4d
[p4est] CPP                      mpicc -E
[p4est] CPPFLAGS                 
[p4est] CC                       mpicc
[p4est] CFLAGS                   -O0 -g -Wall
[p4est] LDFLAGS                  
[p4est] LIBS                     -lz -lm -lgfortran
[fclaw] This is ForestClaw 0.1.6366-63b48-dirty
[fclaw] CPP                      mpicc -E
[fclaw] CPPFLAGS                 
[fclaw] F77                      gfortran
[fclaw] FFLAGS                   -O0 -g -cpp
[fclaw] CC                       mpicc
[fclaw] CFLAGS                   -O0 -g -Wall
[fclaw] CXX                      mpicxx
[fclaw] CXXFLAGS                 -O0 -g -Wall
[fclaw] LDFLAGS                  
[fclaw] FLIBS                     -L/opt/software/easybuild/software/binutils/2.35-GCCcore-10.2.0/lib/../lib64 -L/opt/software/easybuild/software/zlib/1.2.11-GCCcore-10.2.0/lib/../lib64 -L/opt/software/easybuild/software/GCCcore/10.2.0/lib/gcc/x86_64-pc-linux-gnu/10.2.0 -L/opt/software/easybuild/software/GCCcore/10.2.0/lib/gcc/x86_64-pc-linux-gnu/10.2.0/../../../../lib64 -L/lib/x86_64-linux-gnu -L/lib/../lib64 -L/usr/lib/x86_64-linux-gnu -L/opt/software/easybuild/software/binutils/2.35-GCCcore-10.2.0/lib -L/opt/software/easybuild/software/zlib/1.2.11-GCCcore-10.2.0/lib -L/opt/software/easybuild/software/GCCcore/10.2.0/lib/gcc/x86_64-pc-linux-gnu/10.2.0/../../.. -lgfortran -lm -lquadmath
[fclaw] LIBS                     -lz -lm -lgfortran
[fclaw] Reading file fclaw_options.ini.
[fclaw] Options:
[fclaw]    help                                 false
[fclaw]    version                              false
[fclaw]    print-options                        false
[fclaw]    verbosity                            production
[fclaw]    lib-verbosity                        essential
[fclaw]    initial_dt                           0.015625
[fclaw]    max_cfl                              1
[fclaw]    desired_cfl                          0.9
[fclaw]    reduce-cfl                           true
[fclaw]    use_fixed_dt                         false
[fclaw]    outstyle                             1
[fclaw]    tfinal                               0.1
[fclaw]    nout                                 1
[fclaw]    nstep                                1
[fclaw]    advance-one-step                     false
[fclaw]    outstyle-uses-maxlevel               false
[fclaw]    subcycle                             true
[fclaw]    weighted_partition                   true
[fclaw]    time-sync                            true
[fclaw]    flux-correction                      true
[fclaw]    fluctuation-correction               true
[fclaw]    output                               false
[fclaw]    output-gauges                        false
[fclaw]    gauge-buffer-length                  1
[fclaw]    tikz-out                             false
[fclaw]    tikz-figsize                         8 6
[fclaw]    tikz-plot-prefix                     plot
[fclaw]    tikz-plot-suffix                     png
[fclaw]    tikz-mesh-only                       false
[fclaw]    tikz-plot-fig                        true
[fclaw]    prefix                               fort
[fclaw]    vtkspace                             0
[fclaw]    init_ghostcell                       false
[fclaw]    minlevel                             3
[fclaw]    maxlevel                             3
[fclaw]    regrid_interval                      1
[fclaw]    refratio                             2
[fclaw]    smooth-refine                        true
[fclaw]    smooth-level                         4
[fclaw]    coarsen-delay                        0
[fclaw]    refine_threshold                     1
[fclaw]    coarsen_threshold                    0.5
[fclaw]    run-user-diagnostics                 false
[fclaw]    compute-error                        false
[fclaw]    conservation-check                   false
[fclaw]    report-timing                        true
[fclaw]    report-timing-verbosity              wall
[fclaw]    ghost_patch_pack_area                true
[fclaw]    ghost_patch_pack_extra               false
[fclaw]    ghost_patch_pack_numextrafields      0
[fclaw]    trapfpe                              true
[fclaw]    mpi_debug                            false
[fclaw]    ax                                   -2.5
[fclaw]    bx                                   2.5
[fclaw]    ay                                   -2.5
[fclaw]    by                                   2.5
[fclaw]    az                                   0
[fclaw]    bz                                   1
[fclaw]    manifold                             true
[fclaw]    mi                                   1
[fclaw]    mj                                   1
[fclaw]    periodic_x                           false
[fclaw]    periodic_y                           false
[fclaw]    scale                                2.5 2.5 0
[fclaw]    shift                                0 0 0
[fclaw]    phi                                  0
[fclaw]    theta                                0
[fclaw]    clawpatch:mx                         16
[fclaw]    clawpatch:my                         16
[fclaw]    clawpatch:maux                       16
[fclaw]    clawpatch:mbc                        2
[fclaw]    clawpatch:meqn                       4
[fclaw]    clawpatch:rhs-fields                 1
[fclaw]    clawpatch:interp_stencil_width       3
[fclaw]    clawpatch:ghost_patch_pack_aux       true
[fclaw]    clawpatch:refinement-criteria        gradient
[fclaw]    clawpack46:order                     2 2
[fclaw]    clawpack46:mcapa                     1
[fclaw]    clawpack46:src_term                  false
[fclaw]    clawpack46:use-fwaves                false
[fclaw]    clawpack46:mwaves                    3
[fclaw]    clawpack46:mthlim                    4 4 4
[fclaw]    clawpack46:mthbc                     1 1 1 1
[fclaw]    clawpack46:ascii-out                 true
[fclaw]    clawpack46:vtk-out                   false
[fclaw]    clawpack5:order                      2 2
[fclaw]    clawpack5:mcapa                      1
[fclaw]    clawpack5:src_term                   false
[fclaw]    clawpack5:use-fwaves                 false
[fclaw]    clawpack5:mwaves                     3
[fclaw]    clawpack5:mthlim                     4 4 4
[fclaw]    clawpack5:mthbc                      1 1 1 1
[fclaw]    clawpack5:ascii-out                  true
[fclaw]    clawpack5:vtk-out                    false
[fclaw]    user:example                         1
[fclaw]    user:g                               1
[fclaw]    user:x0                              0
[fclaw]    user:y0                              0
[fclaw]    user:r0                              0.5
[fclaw]    user:hin                             2
[fclaw]    user:hout                            1
[fclaw]    user:alpha                           0.4
[fclaw]    user:claw-version                    5
[fclaw] Arguments: none
[p4est 0] Local minimum/maximum levels:  3  3
[p4est] Global minimum/maximum levels:  3  3
[p4est 0] Patches on level  3:       160
[p4est 1] Local minimum/maximum levels:  3  3
[p4est 1] Patches on level  3:       160
[fclaw] Max threads set to 0
[fclaw] Time step    0  at time   0.00000000e+00

[fclaw] Level 3 (3-3) step     1 : dt =    1.562e-02; maxcfl (step) =    1.423; Final time =       0.0156
[fclaw]    WARNING : Maximum CFL exceeded; retaking time step
[fclaw] Level 3 (3-3) step     1 : dt =    9.885e-03; maxcfl (step) =    0.900; Final time =       0.0099
[fclaw] Level 3 (3-3) step     2 : dt =    9.885e-03; maxcfl (step) =    0.974; Final time =       0.0198
[fclaw] Level 3 (3-3) step     3 : dt =    9.131e-03; maxcfl (step) =    0.948; Final time =       0.0289
[fclaw] Level 3 (3-3) step     4 : dt =    8.666e-03; maxcfl (step) =    0.903; Final time =       0.0376
[fclaw] Level 3 (3-3) step     5 : dt =    8.634e-03; maxcfl (step) =    0.892; Final time =       0.0462
[fclaw] Level 3 (3-3) step     6 : dt =    8.712e-03; maxcfl (step) =    0.899; Final time =       0.0549
[fclaw] Level 3 (3-3) step     7 : dt =    8.722e-03; maxcfl (step) =    0.901; Final time =       0.0636
[fclaw] Level 3 (3-3) step     8 : dt =    8.714e-03; maxcfl (step) =    0.900; Final time =       0.0724
[fclaw] Level 3 (3-3) step     9 : dt =    8.714e-03; maxcfl (step) =    0.901; Final time =       0.0811
[fclaw] Level 3 (3-3) step    10 : dt =    8.703e-03; maxcfl (step) =    0.901; Final time =       0.0898
[fclaw] Level 3 (3-3) step    11 : dt =    8.698e-03; maxcfl (step) =    0.901; Final time =       0.0985
[fclaw] Level 3 (3-3) step    12 : dt =    1.535e-03; maxcfl (step) =    0.159; Final time =       0.1000
[fclaw] Time step    1  at time   1.00000000e-01

[fclaw] Finalizing run
[libsc] Statistics for   WALLTIME
[libsc]    Global number of values:       2
[libsc]    Mean value (std. dev.):           1.43115 (5.96e-07 = 4.16e-05%)
[libsc]    Minimum attained at rank       0: 1.43115
[libsc]    Maximum attained at rank       1: 1.43115
[libsc] Procs 2 advance 13 1.21158 exchange 13 0.0985478 regrid 12 0 5.56707e-05
[libsc] Max/P 2 advance 13 1.21197 exchange 13 0.0991428 regrid 12 0 5.67436e-05

minlevel=maxlevel=4

[libsc] This is libsc 2.8.2.85-1ae8
[libsc] CPP                      mpicc -E
[libsc] CPPFLAGS                 
[libsc] CC                       mpicc
[libsc] CFLAGS                   -O0 -g -Wall
[libsc] LDFLAGS                  
[libsc] LIBS                     -lz -lm -lgfortran
[p4est] This is p4est 2.8.7-6c4d
[p4est] CPP                      mpicc -E
[p4est] CPPFLAGS                 
[p4est] CC                       mpicc
[p4est] CFLAGS                   -O0 -g -Wall
[p4est] LDFLAGS                  
[p4est] LIBS                     -lz -lm -lgfortran
[fclaw] This is ForestClaw 0.1.6366-63b48-dirty
[fclaw] CPP                      mpicc -E
[fclaw] CPPFLAGS                 
[fclaw] F77                      gfortran
[fclaw] FFLAGS                   -O0 -g -cpp
[fclaw] CC                       mpicc
[fclaw] CFLAGS                   -O0 -g -Wall
[fclaw] CXX                      mpicxx
[fclaw] CXXFLAGS                 -O0 -g -Wall
[fclaw] LDFLAGS                  
[fclaw] FLIBS                     -L/opt/software/easybuild/software/binutils/2.35-GCCcore-10.2.0/lib/../lib64 -L/opt/software/easybuild/software/zlib/1.2.11-GCCcore-10.2.0/lib/../lib64 -L/opt/software/easybuild/software/GCCcore/10.2.0/lib/gcc/x86_64-pc-linux-gnu/10.2.0 -L/opt/software/easybuild/software/GCCcore/10.2.0/lib/gcc/x86_64-pc-linux-gnu/10.2.0/../../../../lib64 -L/lib/x86_64-linux-gnu -L/lib/../lib64 -L/usr/lib/x86_64-linux-gnu -L/opt/software/easybuild/software/binutils/2.35-GCCcore-10.2.0/lib -L/opt/software/easybuild/software/zlib/1.2.11-GCCcore-10.2.0/lib -L/opt/software/easybuild/software/GCCcore/10.2.0/lib/gcc/x86_64-pc-linux-gnu/10.2.0/../../.. -lgfortran -lm -lquadmath
[fclaw] LIBS                     -lz -lm -lgfortran
[fclaw] Reading file fclaw_options.ini.
[fclaw] Options:
[fclaw]    help                                 false
[fclaw]    version                              false
[fclaw]    print-options                        false
[fclaw]    verbosity                            production
[fclaw]    lib-verbosity                        essential
[fclaw]    initial_dt                           0.015625
[fclaw]    max_cfl                              1
[fclaw]    desired_cfl                          0.9
[fclaw]    reduce-cfl                           true
[fclaw]    use_fixed_dt                         false
[fclaw]    outstyle                             1
[fclaw]    tfinal                               0.1
[fclaw]    nout                                 1
[fclaw]    nstep                                1
[fclaw]    advance-one-step                     false
[fclaw]    outstyle-uses-maxlevel               false
[fclaw]    subcycle                             true
[fclaw]    weighted_partition                   true
[fclaw]    time-sync                            true
[fclaw]    flux-correction                      true
[fclaw]    fluctuation-correction               true
[fclaw]    output                               false
[fclaw]    output-gauges                        false
[fclaw]    gauge-buffer-length                  1
[fclaw]    tikz-out                             false
[fclaw]    tikz-figsize                         8 6
[fclaw]    tikz-plot-prefix                     plot
[fclaw]    tikz-plot-suffix                     png
[fclaw]    tikz-mesh-only                       false
[fclaw]    tikz-plot-fig                        true
[fclaw]    prefix                               fort
[fclaw]    vtkspace                             0
[fclaw]    init_ghostcell                       false
[fclaw]    minlevel                             4
[fclaw]    maxlevel                             4
[fclaw]    regrid_interval                      1
[fclaw]    refratio                             2
[fclaw]    smooth-refine                        true
[fclaw]    smooth-level                         4
[fclaw]    coarsen-delay                        0
[fclaw]    refine_threshold                     1
[fclaw]    coarsen_threshold                    0.5
[fclaw]    run-user-diagnostics                 false
[fclaw]    compute-error                        false
[fclaw]    conservation-check                   false
[fclaw]    report-timing                        true
[fclaw]    report-timing-verbosity              wall
[fclaw]    ghost_patch_pack_area                true
[fclaw]    ghost_patch_pack_extra               false
[fclaw]    ghost_patch_pack_numextrafields      0
[fclaw]    trapfpe                              true
[fclaw]    mpi_debug                            false
[fclaw]    ax                                   -2.5
[fclaw]    bx                                   2.5
[fclaw]    ay                                   -2.5
[fclaw]    by                                   2.5
[fclaw]    az                                   0
[fclaw]    bz                                   1
[fclaw]    manifold                             true
[fclaw]    mi                                   1
[fclaw]    mj                                   1
[fclaw]    periodic_x                           false
[fclaw]    periodic_y                           false
[fclaw]    scale                                2.5 2.5 0
[fclaw]    shift                                0 0 0
[fclaw]    phi                                  0
[fclaw]    theta                                0
[fclaw]    clawpatch:mx                         16
[fclaw]    clawpatch:my                         16
[fclaw]    clawpatch:maux                       16
[fclaw]    clawpatch:mbc                        2
[fclaw]    clawpatch:meqn                       4
[fclaw]    clawpatch:rhs-fields                 1
[fclaw]    clawpatch:interp_stencil_width       3
[fclaw]    clawpatch:ghost_patch_pack_aux       true
[fclaw]    clawpatch:refinement-criteria        gradient
[fclaw]    clawpack46:order                     2 2
[fclaw]    clawpack46:mcapa                     1
[fclaw]    clawpack46:src_term                  false
[fclaw]    clawpack46:use-fwaves                false
[fclaw]    clawpack46:mwaves                    3
[fclaw]    clawpack46:mthlim                    4 4 4
[fclaw]    clawpack46:mthbc                     1 1 1 1
[fclaw]    clawpack46:ascii-out                 true
[fclaw]    clawpack46:vtk-out                   false
[fclaw]    clawpack5:order                      2 2
[fclaw]    clawpack5:mcapa                      1
[fclaw]    clawpack5:src_term                   false
[fclaw]    clawpack5:use-fwaves                 false
[fclaw]    clawpack5:mwaves                     3
[fclaw]    clawpack5:mthlim                     4 4 4
[fclaw]    clawpack5:mthbc                      1 1 1 1
[fclaw]    clawpack5:ascii-out                  true
[fclaw]    clawpack5:vtk-out                    false
[fclaw]    user:example                         1
[fclaw]    user:g                               1
[fclaw]    user:x0                              0
[fclaw]    user:y0                              0
[fclaw]    user:r0                              0.5
[fclaw]    user:hin                             2
[fclaw]    user:hout                            1
[fclaw]    user:alpha                           0.4
[fclaw]    user:claw-version                    5
[fclaw] Arguments: none
[p4est 0] Local minimum/maximum levels:  4  4
[p4est] Global minimum/maximum levels:  4  4
[p4est 1] Local minimum/maximum levels:  4  4
[p4est 1] Patches on level  4:       640
[p4est 0] Patches on level  4:       640
[fclaw] Max threads set to 0
[fclaw] Time step    0  at time   0.00000000e+00

[fclaw] Level 4 (4-4) step     1 : dt =    1.562e-02; maxcfl (step) =    2.837; Final time =       0.0156
[fclaw]    WARNING : Maximum CFL exceeded; retaking time step
[fclaw] Level 4 (4-4) step     1 : dt =    4.957e-03; maxcfl (step) =    0.900; Final time =       0.0050
[fclaw] Level 4 (4-4) step     2 : dt =    4.957e-03; maxcfl (step) =    0.979; Final time =       0.0099
[fclaw] Level 4 (4-4) step     3 : dt =    4.556e-03; maxcfl (step) =    0.948; Final time =       0.0145
[fclaw] Level 4 (4-4) step     4 : dt =    4.323e-03; maxcfl (step) =    0.904; Final time =       0.0188
[fclaw] Level 4 (4-4) step     5 : dt =    4.305e-03; maxcfl (step) =    0.894; Final time =       0.0231
[fclaw] Level 4 (4-4) step     6 : dt =    4.334e-03; maxcfl (step) =    0.897; Final time =       0.0274
[fclaw] Level 4 (4-4) step     7 : dt =    4.349e-03; maxcfl (step) =    0.900; Final time =       0.0318
[fclaw] Level 4 (4-4) step     8 : dt =    4.347e-03; maxcfl (step) =    0.900; Final time =       0.0361
[fclaw] Level 4 (4-4) step     9 : dt =    4.345e-03; maxcfl (step) =    0.900; Final time =       0.0405
[fclaw] Level 4 (4-4) step    10 : dt =    4.346e-03; maxcfl (step) =    0.901; Final time =       0.0448
[fclaw] Level 4 (4-4) step    11 : dt =    4.343e-03; maxcfl (step) =    0.900; Final time =       0.0492

===================================================================================
=   BAD TERMINATION OF ONE OF YOUR APPLICATION PROCESSES
=   PID 31987 RUNNING AT nanna
=   EXIT CODE: 8
=   CLEANING UP REMAINING PROCESSES
=   YOU CAN IGNORE THE BELOW CLEANUP MESSAGES
===================================================================================
YOUR APPLICATION TERMINATED WITH THE EXIT STRING: Floating point exception (signal 8)
This typically refers to a problem with your application.
Please see the FAQ page for debugging suggestions

minlevel=maxlevel=5

[libsc] This is libsc 2.8.2.85-1ae8
[libsc] CPP                      mpicc -E
[libsc] CPPFLAGS                 
[libsc] CC                       mpicc
[libsc] CFLAGS                   -O0 -g -Wall
[libsc] LDFLAGS                  
[libsc] LIBS                     -lz -lm -lgfortran
[p4est] This is p4est 2.8.7-6c4d
[p4est] CPP                      mpicc -E
[p4est] CPPFLAGS                 
[p4est] CC                       mpicc
[p4est] CFLAGS                   -O0 -g -Wall
[p4est] LDFLAGS                  
[p4est] LIBS                     -lz -lm -lgfortran
[fclaw] This is ForestClaw 0.1.6366-63b48-dirty
[fclaw] CPP                      mpicc -E
[fclaw] CPPFLAGS                 
[fclaw] F77                      gfortran
[fclaw] FFLAGS                   -O0 -g -cpp
[fclaw] CC                       mpicc
[fclaw] CFLAGS                   -O0 -g -Wall
[fclaw] CXX                      mpicxx
[fclaw] CXXFLAGS                 -O0 -g -Wall
[fclaw] LDFLAGS                  
[fclaw] FLIBS                     -L/opt/software/easybuild/software/binutils/2.35-GCCcore-10.2.0/lib/../lib64 -L/opt/software/easybuild/software/zlib/1.2.11-GCCcore-10.2.0/lib/../lib64 -L/opt/software/easybuild/software/GCCcore/10.2.0/lib/gcc/x86_64-pc-linux-gnu/10.2.0 -L/opt/software/easybuild/software/GCCcore/10.2.0/lib/gcc/x86_64-pc-linux-gnu/10.2.0/../../../../lib64 -L/lib/x86_64-linux-gnu -L/lib/../lib64 -L/usr/lib/x86_64-linux-gnu -L/opt/software/easybuild/software/binutils/2.35-GCCcore-10.2.0/lib -L/opt/software/easybuild/software/zlib/1.2.11-GCCcore-10.2.0/lib -L/opt/software/easybuild/software/GCCcore/10.2.0/lib/gcc/x86_64-pc-linux-gnu/10.2.0/../../.. -lgfortran -lm -lquadmath
[fclaw] LIBS                     -lz -lm -lgfortran
[fclaw] Reading file fclaw_options.ini.
[fclaw] Options:
[fclaw]    help                                 false
[fclaw]    version                              false
[fclaw]    print-options                        false
[fclaw]    verbosity                            production
[fclaw]    lib-verbosity                        essential
[fclaw]    initial_dt                           0.015625
[fclaw]    max_cfl                              1
[fclaw]    desired_cfl                          0.9
[fclaw]    reduce-cfl                           true
[fclaw]    use_fixed_dt                         false
[fclaw]    outstyle                             1
[fclaw]    tfinal                               0.1
[fclaw]    nout                                 1
[fclaw]    nstep                                1
[fclaw]    advance-one-step                     false
[fclaw]    outstyle-uses-maxlevel               false
[fclaw]    subcycle                             true
[fclaw]    weighted_partition                   true
[fclaw]    time-sync                            true
[fclaw]    flux-correction                      true
[fclaw]    fluctuation-correction               true
[fclaw]    output                               false
[fclaw]    output-gauges                        false
[fclaw]    gauge-buffer-length                  1
[fclaw]    tikz-out                             false
[fclaw]    tikz-figsize                         8 6
[fclaw]    tikz-plot-prefix                     plot
[fclaw]    tikz-plot-suffix                     png
[fclaw]    tikz-mesh-only                       false
[fclaw]    tikz-plot-fig                        true
[fclaw]    prefix                               fort
[fclaw]    vtkspace                             0
[fclaw]    init_ghostcell                       false
[fclaw]    minlevel                             5
[fclaw]    maxlevel                             5
[fclaw]    regrid_interval                      1
[fclaw]    refratio                             2
[fclaw]    smooth-refine                        true
[fclaw]    smooth-level                         4
[fclaw]    coarsen-delay                        0
[fclaw]    refine_threshold                     1
[fclaw]    coarsen_threshold                    0.5
[fclaw]    run-user-diagnostics                 false
[fclaw]    compute-error                        false
[fclaw]    conservation-check                   false
[fclaw]    report-timing                        true
[fclaw]    report-timing-verbosity              wall
[fclaw]    ghost_patch_pack_area                true
[fclaw]    ghost_patch_pack_extra               false
[fclaw]    ghost_patch_pack_numextrafields      0
[fclaw]    trapfpe                              true
[fclaw]    mpi_debug                            false
[fclaw]    ax                                   -2.5
[fclaw]    bx                                   2.5
[fclaw]    ay                                   -2.5
[fclaw]    by                                   2.5
[fclaw]    az                                   0
[fclaw]    bz                                   1
[fclaw]    manifold                             true
[fclaw]    mi                                   1
[fclaw]    mj                                   1
[fclaw]    periodic_x                           false
[fclaw]    periodic_y                           false
[fclaw]    scale                                2.5 2.5 0
[fclaw]    shift                                0 0 0
[fclaw]    phi                                  0
[fclaw]    theta                                0
[fclaw]    clawpatch:mx                         16
[fclaw]    clawpatch:my                         16
[fclaw]    clawpatch:maux                       16
[fclaw]    clawpatch:mbc                        2
[fclaw]    clawpatch:meqn                       4
[fclaw]    clawpatch:rhs-fields                 1
[fclaw]    clawpatch:interp_stencil_width       3
[fclaw]    clawpatch:ghost_patch_pack_aux       true
[fclaw]    clawpatch:refinement-criteria        gradient
[fclaw]    clawpack46:order                     2 2
[fclaw]    clawpack46:mcapa                     1
[fclaw]    clawpack46:src_term                  false
[fclaw]    clawpack46:use-fwaves                false
[fclaw]    clawpack46:mwaves                    3
[fclaw]    clawpack46:mthlim                    4 4 4
[fclaw]    clawpack46:mthbc                     1 1 1 1
[fclaw]    clawpack46:ascii-out                 true
[fclaw]    clawpack46:vtk-out                   false
[fclaw]    clawpack5:order                      2 2
[fclaw]    clawpack5:mcapa                      1
[fclaw]    clawpack5:src_term                   false
[fclaw]    clawpack5:use-fwaves                 false
[fclaw]    clawpack5:mwaves                     3
[fclaw]    clawpack5:mthlim                     4 4 4
[fclaw]    clawpack5:mthbc                      1 1 1 1
[fclaw]    clawpack5:ascii-out                  true
[fclaw]    clawpack5:vtk-out                    false
[fclaw]    user:example                         1
[fclaw]    user:g                               1
[fclaw]    user:x0                              0
[fclaw]    user:y0                              0
[fclaw]    user:r0                              0.5
[fclaw]    user:hin                             2
[fclaw]    user:hout                            1
[fclaw]    user:alpha                           0.4
[fclaw]    user:claw-version                    5
[fclaw] Arguments: none
[p4est 0] Local minimum/maximum levels:  5  5
[p4est] Global minimum/maximum levels:  5  5
[p4est 1] Local minimum/maximum levels:  5  5
[p4est[p4est 0] Patches on level  5:      2560
 1] Patches on level  5:      2560
[fclaw] Max threads set to 0
[fclaw] Time step    0  at time   0.00000000e+00

[fclaw] Level 5 (5-5) step     1 : dt =    1.562e-02; maxcfl (step) =    5.665; Final time =       0.0156
[fclaw]    WARNING : Maximum CFL exceeded; retaking time step
[fclaw] Level 5 (5-5) step     1 : dt =    2.482e-03; maxcfl (step) =    0.900; Final time =       0.0025

===================================================================================
=   BAD TERMINATION OF ONE OF YOUR APPLICATION PROCESSES
=   PID 32025 RUNNING AT nanna
=   EXIT CODE: 8
=   CLEANING UP REMAINING PROCESSES
=   YOU CAN IGNORE THE BELOW CLEANUP MESSAGES
===================================================================================
YOUR APPLICATION TERMINATED WITH THE EXIT STRING: Floating point exception (signal 8)
This typically refers to a problem with your application.
Please see the FAQ page for debugging suggestions

minlevel=maxlevel=6

[libsc] This is libsc 2.8.2.85-1ae8
[libsc] CPP                      mpicc -E
[libsc] CPPFLAGS                 
[libsc] CC                       mpicc
[libsc] CFLAGS                   -O0 -g -Wall
[libsc] LDFLAGS                  
[libsc] LIBS                     -lz -lm -lgfortran
[p4est] This is p4est 2.8.7-6c4d
[p4est] CPP                      mpicc -E
[p4est] CPPFLAGS                 
[p4est] CC                       mpicc
[p4est] CFLAGS                   -O0 -g -Wall
[p4est] LDFLAGS                  
[p4est] LIBS                     -lz -lm -lgfortran
[fclaw] This is ForestClaw 0.1.6366-63b48-dirty
[fclaw] CPP                      mpicc -E
[fclaw] CPPFLAGS                 
[fclaw] F77                      gfortran
[fclaw] FFLAGS                   -O0 -g -cpp
[fclaw] CC                       mpicc
[fclaw] CFLAGS                   -O0 -g -Wall
[fclaw] CXX                      mpicxx
[fclaw] CXXFLAGS                 -O0 -g -Wall
[fclaw] LDFLAGS                  
[fclaw] FLIBS                     -L/opt/software/easybuild/software/binutils/2.35-GCCcore-10.2.0/lib/../lib64 -L/opt/software/easybuild/software/zlib/1.2.11-GCCcore-10.2.0/lib/../lib64 -L/opt/software/easybuild/software/GCCcore/10.2.0/lib/gcc/x86_64-pc-linux-gnu/10.2.0 -L/opt/software/easybuild/software/GCCcore/10.2.0/lib/gcc/x86_64-pc-linux-gnu/10.2.0/../../../../lib64 -L/lib/x86_64-linux-gnu -L/lib/../lib64 -L/usr/lib/x86_64-linux-gnu -L/opt/software/easybuild/software/binutils/2.35-GCCcore-10.2.0/lib -L/opt/software/easybuild/software/zlib/1.2.11-GCCcore-10.2.0/lib -L/opt/software/easybuild/software/GCCcore/10.2.0/lib/gcc/x86_64-pc-linux-gnu/10.2.0/../../.. -lgfortran -lm -lquadmath
[fclaw] LIBS                     -lz -lm -lgfortran
[fclaw] Reading file fclaw_options.ini.
[fclaw] Options:
[fclaw]    help                                 false
[fclaw]    version                              false
[fclaw]    print-options                        false
[fclaw]    verbosity                            production
[fclaw]    lib-verbosity                        essential
[fclaw]    initial_dt                           0.015625
[fclaw]    max_cfl                              1
[fclaw]    desired_cfl                          0.9
[fclaw]    reduce-cfl                           true
[fclaw]    use_fixed_dt                         false
[fclaw]    outstyle                             1
[fclaw]    tfinal                               0.1
[fclaw]    nout                                 1
[fclaw]    nstep                                1
[fclaw]    advance-one-step                     false
[fclaw]    outstyle-uses-maxlevel               false
[fclaw]    subcycle                             true
[fclaw]    weighted_partition                   true
[fclaw]    time-sync                            true
[fclaw]    flux-correction                      true
[fclaw]    fluctuation-correction               true
[fclaw]    output                               false
[fclaw]    output-gauges                        false
[fclaw]    gauge-buffer-length                  1
[fclaw]    tikz-out                             false
[fclaw]    tikz-figsize                         8 6
[fclaw]    tikz-plot-prefix                     plot
[fclaw]    tikz-plot-suffix                     png
[fclaw]    tikz-mesh-only                       false
[fclaw]    tikz-plot-fig                        true
[fclaw]    prefix                               fort
[fclaw]    vtkspace                             0
[fclaw]    init_ghostcell                       false
[fclaw]    minlevel                             6
[fclaw]    maxlevel                             6
[fclaw]    regrid_interval                      1
[fclaw]    refratio                             2
[fclaw]    smooth-refine                        true
[fclaw]    smooth-level                         4
[fclaw]    coarsen-delay                        0
[fclaw]    refine_threshold                     1
[fclaw]    coarsen_threshold                    0.5
[fclaw]    run-user-diagnostics                 false
[fclaw]    compute-error                        false
[fclaw]    conservation-check                   false
[fclaw]    report-timing                        true
[fclaw]    report-timing-verbosity              wall
[fclaw]    ghost_patch_pack_area                true
[fclaw]    ghost_patch_pack_extra               false
[fclaw]    ghost_patch_pack_numextrafields      0
[fclaw]    trapfpe                              true
[fclaw]    mpi_debug                            false
[fclaw]    ax                                   -2.5
[fclaw]    bx                                   2.5
[fclaw]    ay                                   -2.5
[fclaw]    by                                   2.5
[fclaw]    az                                   0
[fclaw]    bz                                   1
[fclaw]    manifold                             true
[fclaw]    mi                                   1
[fclaw]    mj                                   1
[fclaw]    periodic_x                           false
[fclaw]    periodic_y                           false
[fclaw]    scale                                2.5 2.5 0
[fclaw]    shift                                0 0 0
[fclaw]    phi                                  0
[fclaw]    theta                                0
[fclaw]    clawpatch:mx                         16
[fclaw]    clawpatch:my                         16
[fclaw]    clawpatch:maux                       16
[fclaw]    clawpatch:mbc                        2
[fclaw]    clawpatch:meqn                       4
[fclaw]    clawpatch:rhs-fields                 1
[fclaw]    clawpatch:interp_stencil_width       3
[fclaw]    clawpatch:ghost_patch_pack_aux       true
[fclaw]    clawpatch:refinement-criteria        gradient
[fclaw]    clawpack46:order                     2 2
[fclaw]    clawpack46:mcapa                     1
[fclaw]    clawpack46:src_term                  false
[fclaw]    clawpack46:use-fwaves                false
[fclaw]    clawpack46:mwaves                    3
[fclaw]    clawpack46:mthlim                    4 4 4
[fclaw]    clawpack46:mthbc                     1 1 1 1
[fclaw]    clawpack46:ascii-out                 true
[fclaw]    clawpack46:vtk-out                   false
[fclaw]    clawpack5:order                      2 2
[fclaw]    clawpack5:mcapa                      1
[fclaw]    clawpack5:src_term                   false
[fclaw]    clawpack5:use-fwaves                 false
[fclaw]    clawpack5:mwaves                     3
[fclaw]    clawpack5:mthlim                     4 4 4
[fclaw]    clawpack5:mthbc                      1 1 1 1
[fclaw]    clawpack5:ascii-out                  true
[fclaw]    clawpack5:vtk-out                    false
[fclaw]    user:example                         1
[fclaw]    user:g                               1
[fclaw]    user:x0                              0
[fclaw]    user:y0                              0
[fclaw]    user:r0                              0.5
[fclaw]    user:hin                             2
[fclaw]    user:hout                            1
[fclaw]    user:alpha                           0.4
[fclaw]    user:claw-version                    5
[fclaw] Arguments: none
[p4est 0] Local minimum/maximum levels:  6  6
[p4est] Global minimum/maximum levels:  6  6
[p4est 1] Local minimum/maximum levels:  6  6
[p4est 0] Patches on level  6:     10240
[p4est 1] Patches on level  6:     10240
[fclaw] Max threads set to 0
[fclaw] Time step    0  at time   0.00000000e+00

===================================================================================
=   BAD TERMINATION OF ONE OF YOUR APPLICATION PROCESSES
=   PID 32044 RUNNING AT nanna
=   EXIT CODE: 8
=   CLEANING UP REMAINING PROCESSES
=   YOU CAN IGNORE THE BELOW CLEANUP MESSAGES
===================================================================================
YOUR APPLICATION TERMINATED WITH THE EXIT STRING: Floating point exception (signal 8)
This typically refers to a problem with your application.
Please see the FAQ page for debugging suggestions
donnaaboise commented 2 years ago

Interesting. One other thing to try is to make sure data is initialized. In file fclaw2d_farraybox.cpp, line 139, you can try setting data to big_number, rather than snan. Then compile with

-DFCLAW_ENABLE_DEBUG=1.

tim-griesbach commented 2 years ago

One other thing to try is to make sure data is initialized. In file fclaw2d_farraybox.cpp, line 139, you can try setting data to big_number, rather than snan. Then compile with

-DFCLAW_ENABLE_DEBUG=1.

I tried it and it results again in a floating point error. The output for minlevel=maxlevel=4 is below.

Becasue of your conjecture regarding parallel runs.

The tests are run in parallel, where the order in which ghost cells are filled may be slightly different than in serial.

I tried to run

mpirun -n 1 ./radialdam -F regression_map.ini --user:claw-version=5 --user:example=1 --minlevel=4 --maxlevel=4

including the big_number adjustestment but it also crashes because of a floating point exception. The serial output is at the end of this comment.

minlevel=maxlevel=4

[libsc] This is libsc 2.8.2.85-1ae8
[libsc] CPP                      mpicc -E
[libsc] CPPFLAGS                 
[libsc] CC                       mpicc
[libsc] CFLAGS                   -O0 -g -Wall
[libsc] LDFLAGS                  
[libsc] LIBS                     -lz -lm -lgfortran
[p4est] This is p4est 2.8.7-6c4d
[p4est] CPP                      mpicc -E
[p4est] CPPFLAGS                 
[p4est] CC                       mpicc
[p4est] CFLAGS                   -O0 -g -Wall
[p4est] LDFLAGS                  
[p4est] LIBS                     -lz -lm -lgfortran
[fclaw] This is ForestClaw 0.1.6366-63b48-dirty
[fclaw] CPP                      mpicc -E
[fclaw] CPPFLAGS                 
[fclaw] F77                      gfortran
[fclaw] FFLAGS                   -O0 -g -DFCLAW_ENABLE_DEBUG=1 -cpp
[fclaw] CC                       mpicc
[fclaw] CFLAGS                   -O0 -g -Wall
[fclaw] CXX                      mpicxx
[fclaw] CXXFLAGS                 -O0 -g -Wall
[fclaw] LDFLAGS                  
[fclaw] FLIBS                     -L/opt/software/easybuild/software/binutils/2.35-GCCcore-10.2.0/lib/../lib64 -L/opt/software/easybuild/software/zlib/1.2.11-GCCcore-10.2.0/lib/../lib64 -L/opt/software/easybuild/software/GCCcore/10.2.0/lib/gcc/x86_64-pc-linux-gnu/10.2.0 -L/opt/software/easybuild/software/GCCcore/10.2.0/lib/gcc/x86_64-pc-linux-gnu/10.2.0/../../../../lib64 -L/lib/x86_64-linux-gnu -L/lib/../lib64 -L/usr/lib/x86_64-linux-gnu -L/opt/software/easybuild/software/binutils/2.35-GCCcore-10.2.0/lib -L/opt/software/easybuild/software/zlib/1.2.11-GCCcore-10.2.0/lib -L/opt/software/easybuild/software/GCCcore/10.2.0/lib/gcc/x86_64-pc-linux-gnu/10.2.0/../../.. -lgfortran -lm -lquadmath
[fclaw] LIBS                     -lz -lm -lgfortran
[fclaw] Reading file fclaw_options.ini.
[libsc] This is libsc 2.8.2.85-1ae8
[libsc] CPP                      mpicc -E
[libsc] CPPFLAGS                 
[libsc] CC                       mpicc
[libsc] CFLAGS                   -O0 -g -Wall
[libsc] LDFLAGS                  
[libsc] LIBS                     -lz -lm -lgfortran
[p4est] This is p4est 2.8.7-6c4d
[p4est] CPP                      mpicc -E
[p4est] CPPFLAGS                 
[p4est] CC                       mpicc
[p4est] CFLAGS                   -O0 -g -Wall
[p4est] LDFLAGS                  
[p4est] LIBS                     -lz -lm -lgfortran
[fclaw] This is ForestClaw 0.1.6366-63b48-dirty
[fclaw] CPP                      mpicc -E
[fclaw] CPPFLAGS                 
[fclaw] F77                      gfortran
[fclaw] FFLAGS                   -O0 -g -DFCLAW_ENABLE_DEBUG=1 -cpp
[fclaw] CC                       mpicc
[fclaw] CFLAGS                   -O0 -g -Wall
[fclaw] CXX                      mpicxx
[fclaw] CXXFLAGS                 -O0 -g -Wall
[fclaw] LDFLAGS                  
[fclaw] FLIBS                     -L/opt/software/easybuild/software/binutils/2.35-GCCcore-10.2.0/lib/../lib64 -L/opt/software/easybuild/software/zlib/1.2.11-GCCcore-10.2.0/lib/../lib64 -L/opt/software/easybuild/software/GCCcore/10.2.0/lib/gcc/x86_64-pc-linux-gnu/10.2.0 -L/opt/software/easybuild/software/GCCcore/10.2.0/lib/gcc/x86_64-pc-linux-gnu/10.2.0/../../../../lib64 -L/lib/x86_64-linux-gnu -L/lib/../lib64 -L/usr/lib/x86_64-linux-gnu -L/opt/software/easybuild/software/binutils/2.35-GCCcore-10.2.0/lib -L/opt/software/easybuild/software/zlib/1.2.11-GCCcore-10.2.0/lib -L/opt/software/easybuild/software/GCCcore/10.2.0/lib/gcc/x86_64-pc-linux-gnu/10.2.0/../../.. -lgfortran -lm -lquadmath
[fclaw] LIBS                     -lz -lm -lgfortran
[fclaw] Options:
[fclaw]    help                                 false
[fclaw]    version                              false
[fclaw]    print-options                        false
[fclaw]    verbosity                            production
[fclaw]    lib-verbosity                        essential
[fclaw]    initial_dt                           0.015625
[fclaw]    max_cfl                              1
[fclaw]    desired_cfl                          0.9
[fclaw]    reduce-cfl                           true
[fclaw]    use_fixed_dt                         false
[fclaw]    outstyle                             1
[fclaw]    tfinal                               0.1
[fclaw]    nout                                 1
[fclaw]    nstep                                1
[fclaw]    advance-one-step                     false
[fclaw]    outstyle-uses-maxlevel               false
[fclaw]    subcycle                             true
[fclaw]    weighted_partition                   true
[fclaw]    time-sync                            true
[fclaw]    flux-correction                      true
[fclaw]    fluctuation-correction               true
[fclaw]    output                               false
[fclaw]    output-gauges                        false
[fclaw]    gauge-buffer-length                  1
[fclaw]    tikz-out                             false
[fclaw]    tikz-figsize                         8 6
[fclaw]    tikz-plot-prefix                     plot
[fclaw]    tikz-plot-suffix                     png
[fclaw]    tikz-mesh-only                       false
[fclaw]    tikz-plot-fig                        true
[fclaw]    prefix                               fort
[fclaw]    vtkspace                             0
[fclaw]    init_ghostcell                       false
[fclaw]    minlevel                             4
[fclaw]    maxlevel                             4
[fclaw]    regrid_interval                      1
[fclaw]    refratio                             2
[fclaw]    smooth-refine                        true
[fclaw]    smooth-level                         4
[fclaw]    coarsen-delay                        0
[fclaw]    refine_threshold                     1
[fclaw]    coarsen_threshold                    0.5
[fclaw]    run-user-diagnostics                 false
[fclaw]    compute-error                        false
[fclaw]    conservation-check                   false
[fclaw]    report-timing                        true
[fclaw]    report-timing-verbosity              wall
[fclaw]    ghost_patch_pack_area                true
[fclaw]    ghost_patch_pack_extra               false
[fclaw]    ghost_patch_pack_numextrafields      0
[fclaw]    trapfpe                              true
[fclaw]    mpi_debug                            false
[fclaw]    ax                                   -2.5
[fclaw]    bx                                   2.5
[fclaw]    ay                                   -2.5
[fclaw]    by                                   2.5
[fclaw]    az                                   0
[fclaw]    bz                                   1
[fclaw]    manifold                             true
[fclaw]    mi                                   1
[fclaw]    mj                                   1
[fclaw]    periodic_x                           false
[fclaw]    periodic_y                           false
[fclaw]    scale                                2.5 2.5 0
[fclaw]    shift                                0 0 0
[fclaw]    phi                                  0
[fclaw]    theta                                0
[fclaw]    clawpatch:mx                         16
[fclaw]    clawpatch:my                         16
[fclaw]    clawpatch:maux                       16
[fclaw]    clawpatch:mbc                        2
[fclaw]    clawpatch:meqn                       4
[fclaw]    clawpatch:rhs-fields                 1
[fclaw]    clawpatch:interp_stencil_width       3
[fclaw]    clawpatch:ghost_patch_pack_aux       true
[fclaw]    clawpatch:refinement-criteria        gradient
[fclaw]    clawpack46:order                     2 2
[fclaw]    clawpack46:mcapa                     1
[fclaw]    clawpack46:src_term                  false
[fclaw]    clawpack46:use-fwaves                false
[fclaw]    clawpack46:mwaves                    3
[fclaw]    clawpack46:mthlim                    4 4 4
[fclaw]    clawpack46:mthbc                     1 1 1 1
[fclaw]    clawpack46:ascii-out                 true
[fclaw]    clawpack46:vtk-out                   false
[fclaw]    clawpack5:order                      2 2
[fclaw]    clawpack5:mcapa                      1
[fclaw]    clawpack5:src_term                   false
[fclaw]    clawpack5:use-fwaves                 false
[fclaw]    clawpack5:mwaves                     3
[fclaw]    clawpack5:mthlim                     4 4 4
[fclaw]    clawpack5:mthbc                      1 1 1 1
[fclaw] Reading file fclaw_options.ini.
[fclaw]    clawpack5:ascii-out                  true
[fclaw]    clawpack5:vtk-out                    false
[fclaw]    user:example                         1
[fclaw]    user:g                               1
[fclaw]    user:x0                              0
[fclaw]    user:y0                              0
[fclaw]    user:r0                              0.5
[fclaw]    user:hin                             2
[fclaw]    user:hout                            1
[fclaw]    user:alpha                           0.4
[fclaw]    user:claw-version                    5
[fclaw] Arguments: none
[fclaw] Options:
[fclaw]    help                                 false
[fclaw]    version                              false
[fclaw]    print-options                        false
[fclaw]    verbosity                            production
[fclaw]    lib-verbosity                        essential
[fclaw]    initial_dt                           0.015625
[fclaw]    max_cfl                              1
[fclaw]    desired_cfl                          0.9
[fclaw]    reduce-cfl                           true
[fclaw]    use_fixed_dt                         false
[fclaw]    outstyle                             1
[fclaw]    tfinal                               0.1
[fclaw]    nout                                 1
[fclaw]    nstep                                1
[fclaw]    advance-one-step                     false
[fclaw]    outstyle-uses-maxlevel               false
[fclaw]    subcycle                             true
[fclaw]    weighted_partition                   true
[fclaw]    time-sync                            true
[fclaw]    flux-correction                      true
[fclaw]    fluctuation-correction               true
[fclaw]    output                               false
[fclaw]    output-gauges                        false
[fclaw]    gauge-buffer-length                  1
[fclaw]    tikz-out                             false
[fclaw]    tikz-figsize                         8 6
[fclaw]    tikz-plot-prefix                     plot
[fclaw]    tikz-plot-suffix                     png
[fclaw]    tikz-mesh-only                       false
[fclaw]    tikz-plot-fig                        true
[fclaw]    prefix                               fort
[fclaw]    vtkspace                             0
[fclaw]    init_ghostcell                       false
[fclaw]    minlevel                             4
[fclaw]    maxlevel                             4
[fclaw]    regrid_interval                      1
[fclaw]    refratio                             2
[fclaw]    smooth-refine                        true
[fclaw]    smooth-level                         4
[fclaw]    coarsen-delay                        0
[fclaw]    refine_threshold                     1
[fclaw]    coarsen_threshold                    0.5
[fclaw]    run-user-diagnostics                 false
[fclaw]    compute-error                        false
[fclaw]    conservation-check                   false
[fclaw]    report-timing                        true
[fclaw]    report-timing-verbosity              wall
[fclaw]    ghost_patch_pack_area                true
[fclaw]    ghost_patch_pack_extra               false
[fclaw]    ghost_patch_pack_numextrafields      0
[fclaw]    trapfpe                              true
[fclaw]    mpi_debug                            false
[fclaw]    ax                                   -2.5
[fclaw]    bx                                   2.5
[fclaw]    ay                                   -2.5
[fclaw]    by                                   2.5
[fclaw]    az                                   0
[fclaw]    bz                                   1
[fclaw]    manifold                             true
[fclaw]    mi                                   1
[fclaw]    mj                                   1
[fclaw]    periodic_x                           false
[fclaw]    periodic_y                           false
[fclaw]    scale                                2.5 2.5 0
[fclaw]    shift                                0 0 0
[fclaw]    phi                                  0
[fclaw]    theta                                0
[fclaw]    clawpatch:mx                         16
[fclaw]    clawpatch:my                         16
[fclaw]    clawpatch:maux                       16
[fclaw]    clawpatch:mbc                        2
[fclaw]    clawpatch:meqn                       4
[fclaw]    clawpatch:rhs-fields                 1
[fclaw]    clawpatch:interp_stencil_width       3
[fclaw]    clawpatch:ghost_patch_pack_aux       true
[fclaw]    clawpatch:refinement-criteria        gradient
[fclaw]    clawpack46:order                     2 2
[fclaw]    clawpack46:mcapa                     1
[fclaw]    clawpack46:src_term                  false
[fclaw]    clawpack46:use-fwaves                false
[fclaw]    clawpack46:mwaves                    3
[fclaw]    clawpack46:mthlim                    4 4 4
[fclaw]    clawpack46:mthbc                     1 1 1 1
[fclaw]    clawpack46:ascii-out                 true
[fclaw]    clawpack46:vtk-out                   false
[fclaw]    clawpack5:order                      2 2
[fclaw]    clawpack5:mcapa                      1
[fclaw]    clawpack5:src_term                   false
[fclaw]    clawpack5:use-fwaves                 false
[fclaw]    clawpack5:mwaves                     3
[fclaw]    clawpack5:mthlim                     4 4 4
[fclaw]    clawpack5:mthbc                      1 1 1 1
[fclaw]    clawpack5:ascii-out                  true
[fclaw]    clawpack5:vtk-out                    false
[fclaw]    user:example                         1
[fclaw]    user:g                               1
[fclaw]    user:x0                              0
[fclaw]    user:y0                              0
[fclaw]    user:r0                              0.5
[fclaw]    user:hin                             2
[fclaw]    user:hout                            1
[fclaw]    user:alpha                           0.4
[fclaw]    user:claw-version                    5
[fclaw] Arguments: none
[p4est 0] Local minimum/maximum levels:  4  4
[p4est] Global minimum/maximum levels:  4  4
[p4est 0] Patches on level  4:      1280
[fclaw] Max threads set to 0
[p4est 0] Local minimum/maximum levels:  4  4
[p4est] Global minimum/maximum levels:  4  4
[p4est 0] Patches on level  4:      1280
[fclaw] Max threads set to 0
[fclaw] Time step    0  at time   0.00000000e+00

[fclaw] Time step    0  at time   0.00000000e+00

[fclaw] Level 4 (4-4) step     1 : dt =    1.562e-02; maxcfl (step) =    2.837; Final time =       0.0156
[fclaw]    WARNING : Maximum CFL exceeded; retaking time step
[fclaw] Level 4 (4-4) step     1 : dt =    1.562e-02; maxcfl (step) =    2.837; Final time =       0.0156
[fclaw]    WARNING : Maximum CFL exceeded; retaking time step
[fclaw] Level 4 (4-4) step     1 : dt =    4.957e-03; maxcfl (step) =    0.900; Final time =       0.0050
[fclaw] Level 4 (4-4) step     1 : dt =    4.957e-03; maxcfl (step) =    0.900; Final time =       0.0050
[fclaw] Level 4 (4-4) step     2 : dt =    4.957e-03; maxcfl (step) =    0.979; Final time =       0.0099
[fclaw] Level 4 (4-4) step     2 : dt =    4.957e-03; maxcfl (step) =    0.979; Final time =       0.0099
[fclaw] Level 4 (4-4) step     3 : dt =    4.556e-03; maxcfl (step) =    0.948; Final time =       0.0145
[fclaw] Level 4 (4-4) step     3 : dt =    4.556e-03; maxcfl (step) =    0.948; Final time =       0.0145
[fclaw] Level 4 (4-4) step     4 : dt =    4.323e-03; maxcfl (step) =    0.904; Final time =       0.0188
[fclaw] Level 4 (4-4) step     4 : dt =    4.323e-03; maxcfl (step) =    0.904; Final time =       0.0188
[fclaw] Level 4 (4-4) step     5 : dt =    4.305e-03; maxcfl (step) =    0.894; Final time =       0.0231
[fclaw] Level 4 (4-4) step     5 : dt =    4.305e-03; maxcfl (step) =    0.894; Final time =       0.0231
--------------------------------------------------------------------------
mpirun noticed that process rank 0 with PID 0 on node nanna exited on signal 8 (Floating point exception).
--------------------------------------------------------------------------

minlevel=maxlevel=4 serial

[libsc] This is libsc 2.8.2.85-1ae8
[libsc] CPP                      mpicc -E
[libsc] CPPFLAGS                 
[libsc] CC                       mpicc
[libsc] CFLAGS                   -O0 -g -Wall
[libsc] LDFLAGS                  
[libsc] LIBS                     -lz -lm -lgfortran
[p4est] This is p4est 2.8.7-6c4d
[p4est] CPP                      mpicc -E
[p4est] CPPFLAGS                 
[p4est] CC                       mpicc
[p4est] CFLAGS                   -O0 -g -Wall
[p4est] LDFLAGS                  
[p4est] LIBS                     -lz -lm -lgfortran
[fclaw] This is ForestClaw 0.1.6366-63b48-dirty
[fclaw] CPP                      mpicc -E
[fclaw] CPPFLAGS                 
[fclaw] F77                      gfortran
[fclaw] FFLAGS                   -O0 -g -DFCLAW_ENABLE_DEBUG=1 -cpp
[fclaw] CC                       mpicc
[fclaw] CFLAGS                   -O0 -g -Wall
[fclaw] CXX                      mpicxx
[fclaw] CXXFLAGS                 -O0 -g -Wall
[fclaw] LDFLAGS                  
[fclaw] FLIBS                     -L/opt/software/easybuild/software/binutils/2.35-GCCcore-10.2.0/lib/../lib64 -L/opt/software/easybuild/software/zlib/1.2.11-GCCcore-10.2.0/lib/../lib64 -L/opt/software/easybuild/software/GCCcore/10.2.0/lib/gcc/x86_64-pc-linux-gnu/10.2.0 -L/opt/software/easybuild/software/GCCcore/10.2.0/lib/gcc/x86_64-pc-linux-gnu/10.2.0/../../../../lib64 -L/lib/x86_64-linux-gnu -L/lib/../lib64 -L/usr/lib/x86_64-linux-gnu -L/opt/software/easybuild/software/binutils/2.35-GCCcore-10.2.0/lib -L/opt/software/easybuild/software/zlib/1.2.11-GCCcore-10.2.0/lib -L/opt/software/easybuild/software/GCCcore/10.2.0/lib/gcc/x86_64-pc-linux-gnu/10.2.0/../../.. -lgfortran -lm -lquadmath
[fclaw] LIBS                     -lz -lm -lgfortran
[fclaw] Reading file fclaw_options.ini.
[fclaw] Options:
[fclaw]    help                                 false
[fclaw]    version                              false
[fclaw]    print-options                        false
[fclaw]    verbosity                            production
[fclaw]    lib-verbosity                        essential
[fclaw]    initial_dt                           0.015625
[fclaw]    max_cfl                              1
[fclaw]    desired_cfl                          0.9
[fclaw]    reduce-cfl                           true
[fclaw]    use_fixed_dt                         false
[fclaw]    outstyle                             1
[fclaw]    tfinal                               0.1
[fclaw]    nout                                 1
[fclaw]    nstep                                1
[fclaw]    advance-one-step                     false
[fclaw]    outstyle-uses-maxlevel               false
[fclaw]    subcycle                             true
[fclaw]    weighted_partition                   true
[fclaw]    time-sync                            true
[fclaw]    flux-correction                      true
[fclaw]    fluctuation-correction               true
[fclaw]    output                               false
[fclaw]    output-gauges                        false
[fclaw]    gauge-buffer-length                  1
[fclaw]    tikz-out                             false
[fclaw]    tikz-figsize                         8 6
[fclaw]    tikz-plot-prefix                     plot
[fclaw]    tikz-plot-suffix                     png
[fclaw]    tikz-mesh-only                       false
[fclaw]    tikz-plot-fig                        true
[fclaw]    prefix                               fort
[fclaw]    vtkspace                             0
[fclaw]    init_ghostcell                       false
[fclaw]    minlevel                             4
[fclaw]    maxlevel                             4
[fclaw]    regrid_interval                      1
[fclaw]    refratio                             2
[fclaw]    smooth-refine                        true
[fclaw]    smooth-level                         4
[fclaw]    coarsen-delay                        0
[fclaw]    refine_threshold                     1
[fclaw]    coarsen_threshold                    0.5
[fclaw]    run-user-diagnostics                 false
[fclaw]    compute-error                        false
[fclaw]    conservation-check                   false
[fclaw]    report-timing                        true
[fclaw]    report-timing-verbosity              wall
[fclaw]    ghost_patch_pack_area                true
[fclaw]    ghost_patch_pack_extra               false
[fclaw]    ghost_patch_pack_numextrafields      0
[fclaw]    trapfpe                              true
[fclaw]    mpi_debug                            false
[fclaw]    ax                                   -2.5
[fclaw]    bx                                   2.5
[fclaw]    ay                                   -2.5
[fclaw]    by                                   2.5
[fclaw]    az                                   0
[fclaw]    bz                                   1
[fclaw]    manifold                             true
[fclaw]    mi                                   1
[fclaw]    mj                                   1
[fclaw]    periodic_x                           false
[fclaw]    periodic_y                           false
[fclaw]    scale                                2.5 2.5 0
[fclaw]    shift                                0 0 0
[fclaw]    phi                                  0
[fclaw]    theta                                0
[fclaw]    clawpatch:mx                         16
[fclaw]    clawpatch:my                         16
[fclaw]    clawpatch:maux                       16
[fclaw]    clawpatch:mbc                        2
[fclaw]    clawpatch:meqn                       4
[fclaw]    clawpatch:rhs-fields                 1
[fclaw]    clawpatch:interp_stencil_width       3
[fclaw]    clawpatch:ghost_patch_pack_aux       true
[fclaw]    clawpatch:refinement-criteria        gradient
[fclaw]    clawpack46:order                     2 2
[fclaw]    clawpack46:mcapa                     1
[fclaw]    clawpack46:src_term                  false
[fclaw]    clawpack46:use-fwaves                false
[fclaw]    clawpack46:mwaves                    3
[fclaw]    clawpack46:mthlim                    4 4 4
[fclaw]    clawpack46:mthbc                     1 1 1 1
[fclaw]    clawpack46:ascii-out                 true
[fclaw]    clawpack46:vtk-out                   false
[fclaw]    clawpack5:order                      2 2
[fclaw]    clawpack5:mcapa                      1
[fclaw]    clawpack5:src_term                   false
[fclaw]    clawpack5:use-fwaves                 false
[fclaw]    clawpack5:mwaves                     3
[fclaw]    clawpack5:mthlim                     4 4 4
[fclaw]    clawpack5:mthbc                      1 1 1 1
[fclaw]    clawpack5:ascii-out                  true
[fclaw]    clawpack5:vtk-out                    false
[fclaw]    user:example                         1
[fclaw]    user:g                               1
[fclaw]    user:x0                              0
[fclaw]    user:y0                              0
[fclaw]    user:r0                              0.5
[fclaw]    user:hin                             2
[fclaw]    user:hout                            1
[fclaw]    user:alpha                           0.4
[fclaw]    user:claw-version                    5
[fclaw] Arguments: none
[p4est 0] Local minimum/maximum levels:  4  4
[p4est] Global minimum/maximum levels:  4  4
[p4est 0] Patches on level  4:      1280
[fclaw] Max threads set to 0
[fclaw] Time step    0  at time   0.00000000e+00

[fclaw] Level 4 (4-4) step     1 : dt =    1.562e-02; maxcfl (step) =    2.837; Final time =       0.0156
[fclaw]    WARNING : Maximum CFL exceeded; retaking time step
[fclaw] Level 4 (4-4) step     1 : dt =    4.957e-03; maxcfl (step) =    0.900; Final time =       0.0050
[fclaw] Level 4 (4-4) step     2 : dt =    4.957e-03; maxcfl (step) =    0.979; Final time =       0.0099
[fclaw] Level 4 (4-4) step     3 : dt =    4.556e-03; maxcfl (step) =    0.948; Final time =       0.0145
[fclaw] Level 4 (4-4) step     4 : dt =    4.323e-03; maxcfl (step) =    0.904; Final time =       0.0188
[fclaw] Level 4 (4-4) step     5 : dt =    4.305e-03; maxcfl (step) =    0.894; Final time =       0.0231
--------------------------------------------------------------------------
mpirun noticed that process rank 0 with PID 0 on node nanna exited on signal 8 (Floating point exception).
--------------------------------------------------------------------------
donnaaboise commented 1 year ago

This doesn't seem to be causing problems with CI anymore.