From c16c447979811e8e6449a3329abb9d294d31cb78 Mon Sep 17 00:00:00 2001
From: Jon Drobny <drobny2@illinois.edu>
Date: Wed, 27 Oct 2021 20:18:41 -0500
Subject: [PATCH] Added context to error when particle energy goes NaN
 following collision.

---
 src/bca.rs | 1 +
 1 file changed, 1 insertion(+)

diff --git a/src/bca.rs b/src/bca.rs
index 9c70511..4f911bd 100644
--- a/src/bca.rs
+++ b/src/bca.rs
@@ -447,6 +447,7 @@ pub fn update_particle_energy<T: Geometry>(particle_1: &mut particle::Particle,
     recoil_energy: f64, x0: f64, strong_collision_Z: f64, strong_collision_index: usize, options: &Options) {
 
     //If particle energy  drops below zero before electronic stopping calcualtion, it produces NaNs
+    assert!(!recoil_energy.is_nan(), "Numerical error: recoil Energy is NaN. E0 = {} Za = {} Ma = {} x0 = {} Zb = {} delta-x = {}", particle_1.E, particle_1.Z, particle_1.m, x0, strong_collision_Z, distance_traveled);
     particle_1.E -= recoil_energy;
     assert!(!particle_1.E.is_nan(), "Numerical error: particle energy is NaN following collision.");
     if particle_1.E < 0. {