Changeset bba9f8 for rttests.lisp
 Timestamp:
 04/11/12 15:47:28 (3 years ago)
 Branches:
 master
 Children:
 7c5a31
 Parents:
 104efd
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

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