*** empty log message ***

quick.cc

@@ -30,7 +30,7 @@ int main(int argc, char *argv[])
     OUT(tmp2);
 
     OUT(T2z::getInstance().converttologp(tmp,tmp2));
-
+    OUT(std::exp(T2z::getInstance().converttologp(tmp,tmp2)));
     OUT(T2z::getInstance().convert(tmp,tmp2));
 
   }
@@ -41,14 +41,3 @@ int main(int argc, char *argv[])
   return 0;
 }
 
-
-
-
-
-
-
-
-
-
-
-

t2z.cc

@@ -114,8 +114,8 @@ namespace MISCMATHS {
 	return (logp < -14.5);
       }
 
-  float T2z::convert(float t, int dof) 
-    {
+  float T2z::convert(float t, int dof) {
+
       float z = 0.0, logp=0.0;
       
       if(!islarget(t,dof,logp)) {
@@ -125,11 +125,16 @@ namespace MISCMATHS {
 	z = MISCMATHS::ndtri(p);
       }
       else {
+	
 	z = logp2largez(logp);
+
+	//	cerr<<endl<<"logp="<<logp<<endl;
+
 	if (t<0) z=-z;
       }
 
       return z;
+
     }
 
 
@@ -140,6 +145,12 @@ namespace MISCMATHS {
       if(!islarget(t,dof,logp)) {
 	logp = log(1-MISCMATHS::stdtr(dof, t));
       }
+      else if(t<0) {
+	// t < 0 and abs(t) is large enough to require asymptotic approx.
+	// but t to logp is not available for negative t 
+	// so just hardcode it to be -1e-12
+	logp=-1e-12;
+      }
 
 //       cerr << "logp = " << logp << endl;
 //       cerr << "exp(logp) = " << std::exp(logp) << endl;