diff options
Diffstat (limited to 'src/Equations.c')
| -rw-r--r-- | src/Equations.c | 4 |
1 files changed, 4 insertions, 0 deletions
diff --git a/src/Equations.c b/src/Equations.c index 1c9c5f4..4bbcbef 100644 --- a/src/Equations.c +++ b/src/Equations.c @@ -88,6 +88,10 @@ BY_Aijofxyz (double x, double y, double z, double Aij[3][3]) r2_plus = (x - par_b) * (x - par_b) + y * y + z * z; r2_minus = (x + par_b) * (x + par_b) + y * y + z * z; + if (r2_plus < TP_Tiny) + r2_plus = TP_Tiny; + if (r2_minus < TP_Tiny) + r2_minus = TP_Tiny; r_plus = sqrt (r2_plus); r_minus = sqrt (r2_minus); r3_plus = r_plus * r2_plus; |