Dear developers, I am conducting a simulation of chemical reactions in the DEM. I encountered a problem. The calculation could proceed normally at the beginning, but when the calculation time reached around 8.5 seconds, an error occurred, indicating “Program received signal SIGFPE: Floating-point exception - erroneous arithmetic operation.
Backtrace for this error:
#0 leq_bicgs_mod_MOD_leq_bicgs0
at leq_bicgs.f:404
#1 leq_bicgs_mod_MOD_leq_bicgs
at leq_bicgs.f:91
#2 solve_lin_eq_mod_MOD_solve_lin_eq
at solve_lin_eq.f:165
#3 solve_species_eq_mod_MOD_solve_species_eq
at solve_species_eq.f:192
#4 iterate_MOD_do_iteration
at iterate.f:367
#5 run_fluid
at mfix.f:187
#6 run_mfix
at mfix.f:142
#7 main_MOD_run_mfix0
at main.f:83
#8 start_thread
At /build/glibc-SzIz7B/glibc-2.31/nptl/pthread_create.c:477 ”, The error persists even after I tried recalculating, and it still occurs around this simulation time. The attachment contains the materials . Thank you for your assistance.
dem_2026-05-08T103610.780643.zip (45.5 MB)