Changeset f4a60f
 Timestamp:
 03/13/11 04:19:11 (4 years ago)
 Branches:
 master
 Children:
 8ade17
 Parents:
 5bd5df
 Files:

 2 edited
Legend:
 Unmodified
 Added
 Removed

qdelliptic.lisp
r5bd5df rf4a60f 715 715 (sinphi (sin phi)) 716 716 (cosphi (cos phi)) 717 (msin2 ( 1 (* m sinphi sinphi))) 717 (msin2 ( 1 (* m sinphi sinphi)))) 718 718 ( (* sinphi (carlsonrf (expt cosphi 2) msin2 1)) 719 719 (* (/ nn 3) (expt sinphi 3) 
rttests.lisp
r5bd5df rf4a60f 26 26 (inpackage #:oct) 27 27 28 (evalwhen (:compiletoplevel :loadtoplevel :execute) 29 (setf *readtable* *octreadtable*)) 30 28 31 ;; For the tests, we need to turn off underflow for clisp. 29 32 #+clisp … … 941 944 for epi = (ellipticpi 0 phi m) 942 945 for ef = (ellipticf phi m) 943 for result = (checkaccuracy 5 3epi ef)946 for result = (checkaccuracy 51 epi ef) 944 947 unless (eq nil result) 945 948 append (list (list phi m) result)) … … 999 1002 for true = (/ (atan (* (tan phi) (sqrt ( 1 n)))) 1000 1003 (sqrt ( 1 n))) 1001 for result = (checkaccuracy 50epi true)1004 for result = (checkaccuracy 48 epi true) 1002 1005 unless (eq nil result) 1003 1006 append (list (list (list k n phi) result))) … … 1009 1012 for epi = (ellipticpi 1 phi 0) 1010 1013 for true = (tan phi) 1011 for result = (checkaccuracy 43epi true)1014 for result = (checkaccuracy 37 epi true) 1012 1015 unless (eq nil result) 1013 1016 append (list (list (list k phi) result))) … … 1046 1049 for epi = (ellipticpi 1 phi 0) 1047 1050 for true = (tan phi) 1048 for result = (checkaccuracy 20 5epi true)1051 for result = (checkaccuracy 200 epi true) 1049 1052 unless (eq nil result) 1050 1053 append (list (list (list k phi) result))) … … 1058 1061 for true = (/ (atanh (* (tan phi) (sqrt ( n 1)))) 1059 1062 (sqrt ( n 1))) 1060 for result = (checkaccuracy 20 8epi true)1063 for result = (checkaccuracy 207 epi true) 1061 1064 ;; Not sure if this formula holds when atanh gives a complex 1062 1065 ;; result. Wolfram doesn't say
Note: See TracChangeset
for help on using the changeset viewer.