From: SourceForge.net <no...@so...> - 2009-07-02 12:57:24
|
Bugs item #2815766, was opened at 2009-07-02 12:57 Message generated for change (Tracker Item Submitted) made by nobody You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=2815766&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: Lisp Core - Solving equations Group: None Status: Open Resolution: None Priority: 5 Private: No Submitted By: Nobody/Anonymous (nobody) Assigned to: Nobody/Anonymous (nobody) Summary: strange results from eigenvalues Initial Comment: large matrices like matrix([72126.5,-25007.8,7670.53],[-25007.8,8670.74,-2659.53],[7670.53,-2659.53,818.809]) produce strange outputs: (%i32) eigenvalues(T); rat: replaced 25007.8 by 125039/5 = 25007.8 rat: replaced 2.0400004650900003E+7 by 20400005/1 = 2.0400005E+7 rat: replaced -25007.8 by -125039/5 = -25007.8 rat: replaced 818.809 by 72874/89 = 818.8089887640449 rat: replaced 7670.53 by 767053/100 = 7670.53 rat: replaced 6.6508994334000006E+7 by 66508994/1 = 6.6508994E+7 rat: replaced -7670.53 by -767053/100 = -7670.53 rat: replaced 8670.74 by 433537/50 = 8670.74 rat: replaced -7073099.8209 by -35365499/5 = -7073099.8 rat: replaced 818.809 by 72874/89 = 818.8089887640449 rat: replaced 8670.74 by 433537/50 = 8670.74 rat: replaced 72126.5 by 144253/2 = 72126.5 (%o32) [[(-(sqrt(3)*%i)/2-1/2)*((3*sqrt(12608488167242884134811801138679920231)*%i)/(3960500000*sqrt(2))+2838509476210882261900497/140993800000)^(1/3)+(58619247149883207*((sqrt(3)*%i)/2-1/2))/(79210000*((3*sqrt(12608488167242884134811801138679920231)*%i)/(3960500000*sqrt(2))+2838509476210882261900497/140993800000)^(1/3))+60531903/2225,((sqrt(3)*%i)/2-1/2)*((3*sqrt(12608488167242884134811801138679920231)*%i)/(3960500000*sqrt(2))+2838509476210882261900497/140993800000)^(1/3)+(58619247149883207*(-(sqrt(3)*%i)/2-1/2))/(79210000*((3*sqrt(12608488167242884134811801138679920231)*%i)/(3960500000*sqrt(2))+2838509476210882261900497/140993800000)^(1/3))+60531903/2225,((3*sqrt(12608488167242884134811801138679920231)*%i)/(3960500000*sqrt(2))+2838509476210882261900497/140993800000)^(1/3)+58619247149883207/(79210000*((3*sqrt(12608488167242884134811801138679920231)*%i)/(3960500000*sqrt(2))+2838509476210882261900497/140993800000)^(1/3))+60531903/2225],[1,1,1]] ovtave gives followin: octave.exe:3> t = >[72126.5, -25007.8, 7670.53; > -25007.8, 8670.74, -2659.53; > 7670.53, -2659.53, 818.809]; octave.exe:4> eig(t) ans = 8.2521e-004 3.0311e+000 8.1613e+004 ---------------------------------------------------------------------- You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=2815766&group_id=4933 |