From: SourceForge.net <no...@so...> - 2009-05-01 17:28:59
|
Bugs item #2784167, was opened at 2009-04-30 06:41 Message generated for change (Comment added) made by nobody You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=2784167&group_id=4933 Please note that this message will contain a full copy of the comment thread, including the initial issue submission, for this request, not just the latest update. Category: Tests Group: None Status: Open Resolution: None Priority: 5 Private: No Submitted By: Kenneth Sloan (kennethrsloan) Assigned to: Nobody/Anonymous (nobody) Summary: 101 & 139 fail MacOSX Initial Comment: ------------------------------------------------------------- Maxima version: 5.18.1 Maxima build date: 14:34 4/18/2009 host type: i686-apple-darwin8.11.1 lisp-implementation-type: CMU Common Lisp lisp-implementation-version: 19f (19F) ------------------------------------------------------------- ********************** Problem 101 *************** Input: ev(e5, au = 0, omega = 2) Result: [0.4, 2.104312740630307e-12, 225, 0] This differed from the expected result: [.4000000000000001, 2.216570948815925e-11, 175, 0] ********************** Problem 139 *************** Input: expintegral_chi(z) - conjugate(expintegral_chi(conjugate(z))) Result: 0.0339413106060743 - .07647309974404348 %i This differed from the expected result: 0.0 ---------------------------------------------------------------------- Comment By: Nobody/Anonymous (nobody) Date: 2009-05-01 17:28 Message: I'm using the MacOSX build - not building from source. I ran the entire test suite - these were the only two errors. Working on the numerical results in response to the previous comment, now. ---------------------------------------------------------------------- Comment By: Raymond Toy (rtoy) Date: 2009-05-01 17:01 Message: I ran the entire testsuite using CMUCL 2009-04 on Mac OS X 10.5 with the current CVS code. I do not get any test failures either. Did you run the entire testsuite or just parts of it? ---------------------------------------------------------------------- Comment By: Dieter Kaiser (crategus) Date: 2009-05-01 14:36 Message: Hello Kenneth, I cannot verify the failures with GCL 2.6.8 and CLISP 2.44. There might be problem with the numericial calculations of expintegral_chi which is special to your Lisp implementation. For the problem 139 in rtest_expintegral.mac I get the following numerical results with CLISP 2.44: (%i12) z:0.5+%i; (%o12) %i + 0.5 (%i13) expintegral_chi(z); (%o13) 1.341686846905066 %i + 0.49714695079962 (%i14) conjugate(expintegral_chi(conjugate(z))); (%o14) 1.341686846905066 %i + 0.49714695079962 Internally, the following formula for the calculation is used: (%i20) -0.5*(expintegral_e(1,z)+expintegral_e(1,-z)+log(-z)-log(z)); (%o20) - 0.5 (- 2.683373693810132 %i - 0.99429390159924) (%i21) expand(%); (%o21) 1.341686846905066 %i + 0.49714695079962 Can you post your numerical results. Dieter Kaiser ---------------------------------------------------------------------- You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=2784167&group_id=4933 |