Changeset bba9f8940c9f904bf14adc405d795a38ac333c24 for rttests.lisp
 Timestamp:
 04/11/12 08:47:28 (2 years ago)
 Children:
 7c5a3186070096ee93e16f2ddf51b2c84e7c5895
 Parents:
 104efdeae6ef52c5d370e0b3a048a5087cdb1ea2
 gitcommitter:
 Raymond Toy <rtoy@…> (04/11/12 08:47:28)
 Files:

 1 modified
Legend:
 Unmodified
 Added
 Removed

rttests.lisp
r104efd rbba9f8 1115 1115 for true = (/ (atanh (* (tan phi) (sqrt ( n 1)))) 1116 1116 (sqrt ( n 1))) 1117 for result = (checkaccuracy 4 7epi true)1117 for result = (checkaccuracy 45.85 epi true) 1118 1118 ;; Not sure if this formula holds when atanh gives a complex 1119 1119 ;; result. Wolfram doesn't say … … 1121 1121 append (list (list (list k n phi) result))) 1122 1122 nil) 1123 1124 ;; Failed test case: 1125 ;; ((89 66.68551748022054d0 0.12266024127708153d0) 1126 ;; (45.868614757480834d0 47 0.47787458521306514d0 1127 ;; 0.4778745852130726d0)) 1128 ;; New threshold is 45.85 bits. 1129 (rt:deftest oct.ellipticpi.n2.d1 1130 (let* ((n 66.68551748022054d0) 1131 (phi 0.12266024127708153d0) 1132 (epi (ellipticpi n phi 0)) 1133 (true (/ (atanh (* (tan phi) (sqrt ( n 1)))) 1134 (sqrt ( n 1))))) 1135 (checkaccuracy 45.8686d0 epi true)) 1136 nil) 1137 1123 1138 1124 1139 (rt:deftest oct.ellipticpi.n0.q … … 1242 1257 (let ((g (cfincompletegamma 3d0 5d0)) 1243 1258 (true ( 2 (* 37 (exp 5d0))))) 1244 (checkaccuracy 5 3g true))1259 (checkaccuracy 51.2 g true)) 1245 1260 nil) 1246 1261