From: SourceForge.net <no...@so...> - 2004-02-25 15:37:41
|
Bugs item #904295, was opened at 2004-02-25 10:26 Message generated for change (Tracker Item Submitted) made by Item Submitter You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=904295&group_id=4933 Category: None Group: None Status: Open Resolution: None Priority: 5 Submitted By: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: ratweight expt coredump/FIX Initial Comment: ratweights(x,1,y,1); rat((1+x+y)^10000),ratwtlvl:3; coredumps (too deep recursion) The fix is simple, and improves the algorithm to boot: In rat3b.lisp: (defun wtpexpt (x n) (cond ((= n 0) 1) ((= n 1) x) ((evenp n) (let ((xn2 (wtpexpt x (/ n 2)))) (wtptimes xn2 xn2 0))) (t (wtptimes x (wtpexpt x (1- n)) 0)))) Now even rat((1+x+y)^(10^30)),ratwtlvl:5; works. Maxima 5.9.0 GCL 2.5.0 mingw32 W2k Athlon ---------------------------------------------------------------------- You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=904295&group_id=4933 |
From: SourceForge.net <no...@so...> - 2004-11-23 16:56:30
|
Bugs item #904295, was opened at 2004-02-25 10:26 Message generated for change (Comment added) made by rtoy You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=904295&group_id=4933 Category: None Group: None >Status: Closed >Resolution: Fixed Priority: 5 Submitted By: Stavros Macrakis (macrakis) Assigned to: Nobody/Anonymous (nobody) Summary: ratweight expt coredump/FIX Initial Comment: ratweights(x,1,y,1); rat((1+x+y)^10000),ratwtlvl:3; coredumps (too deep recursion) The fix is simple, and improves the algorithm to boot: In rat3b.lisp: (defun wtpexpt (x n) (cond ((= n 0) 1) ((= n 1) x) ((evenp n) (let ((xn2 (wtpexpt x (/ n 2)))) (wtptimes xn2 xn2 0))) (t (wtptimes x (wtpexpt x (1- n)) 0)))) Now even rat((1+x+y)^(10^30)),ratwtlvl:5; works. Maxima 5.9.0 GCL 2.5.0 mingw32 W2k Athlon ---------------------------------------------------------------------- >Comment By: Raymond Toy (rtoy) Date: 2004-11-23 11:56 Message: Logged In: YES user_id=28849 Assuming you meant wtpexpt in ratout.lisp, this all works nicely. Suggested fix applied. ---------------------------------------------------------------------- You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=904295&group_id=4933 |