Changeset f4a60f


Ignore:
Timestamp:
03/13/11 04:19:11 (4 years ago)
Author:
Raymond Toy <toy.raymond@…>
Branches:
master
Children:
8ade17
Parents:
5bd5df
Message:

Fix typo; update accuracy requirements.

qd-elliptic.lisp:
o Fix missing paren.

rt-tests.lisp:
o Reduce accuracy requirements so the tests pass for the random

arguments.

Files:
2 edited

Legend:

Unmodified
Added
Removed
  • qd-elliptic.lisp

    r5bd5df rf4a60f  
    715715         (sin-phi (sin phi))
    716716         (cos-phi (cos phi))
    717          (m-sin2 (- 1 (* m sin-phi sin-phi)))
     717         (m-sin2 (- 1 (* m sin-phi sin-phi))))
    718718    (- (* sin-phi (carlson-rf (expt cos-phi 2) m-sin2 1))
    719719       (* (/ nn 3) (expt sin-phi 3)
  • rt-tests.lisp

    r5bd5df rf4a60f  
    2626(in-package #:oct)
    2727
     28(eval-when (:compile-toplevel :load-toplevel :execute)
     29  (setf *readtable* *oct-readtable*))
     30
    2831;; For the tests, we need to turn off underflow for clisp.
    2932#+clisp
     
    941944       for epi = (elliptic-pi 0 phi m)
    942945       for ef = (elliptic-f phi m)
    943        for result = (check-accuracy 53 epi ef)
     946       for result = (check-accuracy 51 epi ef)
    944947       unless (eq nil result)
    945948       append (list (list phi m) result))
     
    9991002       for true = (/ (atan (* (tan phi) (sqrt (- 1 n))))
    10001003                     (sqrt (- 1 n)))
    1001        for result = (check-accuracy 50 epi true)
     1004       for result = (check-accuracy 48 epi true)
    10021005       unless (eq nil result)
    10031006       append (list (list (list k n phi) result)))
     
    10091012       for epi = (elliptic-pi 1 phi 0)
    10101013       for true = (tan phi)
    1011        for result = (check-accuracy 43 epi true)
     1014       for result = (check-accuracy 37 epi true)
    10121015       unless (eq nil result)
    10131016       append (list (list (list k phi) result)))
     
    10461049       for epi = (elliptic-pi 1 phi 0)
    10471050       for true = (tan phi)
    1048        for result = (check-accuracy 205 epi true)
     1051       for result = (check-accuracy 200 epi true)
    10491052       unless (eq nil result)
    10501053       append (list (list (list k phi) result)))
     
    10581061       for true = (/ (atanh (* (tan phi) (sqrt (- n 1))))
    10591062                     (sqrt (- n 1)))
    1060        for result = (check-accuracy 208 epi true)
     1063       for result = (check-accuracy 207 epi true)
    10611064       ;; Not sure if this formula holds when atanh gives a complex
    10621065       ;; result.  Wolfram doesn't say
Note: See TracChangeset for help on using the changeset viewer.