From b2a572b7904a7dd620f1057aeed3bd9d296782e2 Mon Sep 17 00:00:00 2001
From: nemu <nemu>
Date: Fri, 2 Oct 2009 15:18:59 +0000
Subject: [PATCH] fixed a bug in the spin glass function for negative times

---
 src/classes/PTheory.cpp | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/src/classes/PTheory.cpp b/src/classes/PTheory.cpp
index 1deba070..980d5007 100644
--- a/src/classes/PTheory.cpp
+++ b/src/classes/PTheory.cpp
@@ -1497,8 +1497,8 @@ double PTheory::SpinGlass(register double t, const PDoubleVector& paramValues, c
   double lambda_t_2_q = tt*tt*lambda_2*val[2];
   double rate_2       = 4.0*lambda_2*(1.0-val[2])*tt/val[1];
 
-  double rateL = TMath::Sqrt(rate_2);
-  double rateT = TMath::Sqrt(rate_2+lambda_t_2_q);
+  double rateL = TMath::Sqrt(fabs(rate_2));
+  double rateT = TMath::Sqrt(fabs(rate_2)+lambda_t_2_q);
 
   return 0.333333333333333*(TMath::Exp(-rateL) + 2.0*(1.0-lambda_t_2_q/rateT)*TMath::Exp(-rateT));
 }