Hi, I am trying to seed some SQP particle in some given geometry to generate particle output dat file (to be used in another simulation). However I am continuously getting floating point error after a certain time period.

It says:

Program received signal SIGFPE: Floating-point exception - erroneous arithmetic operation.

Backtrace for this error:

#0 sq_contact_newton_dpmethod_mod_MOD_func_dp_a

at des/sq_contact_detection_newton.f:964

#1 sq_contact_newton_dpmethod_mod_MOD_sq_contact_newton_dp_a

at des/sq_contact_detection_newton.f:267

#2 sq_contact_evolution_mod_MOD_sq_evolution_newton_method_a

at des/sq_contact_detection_evolution.f:115

#3 sq_calc_force_mod_MOD_calc_force_superdem

at des/sq_calc_force_superdem.f:354

#4 des_time_march_MOD_des_time_step

#5 run_dem

at mfix.f:214

#6 run_mfix

at mfix.f:149

#7 main_MOD_run_mfix0

at main.f:85

Floating point exception (core dumped).

I have also attached the mfx file for reference. How can i possibly resolve this issue

MI_SQP.mfx (12.2 KB)