New answers tagged

1

This is just a long comment really. Out of curiosity, I compared the result of three packages: Mathematica, MATLAB and eig_sym from Armadillo (compiled on OS X). You said that your C++ code uses Armadillo. I get very close but not identical results: Mathematica MachinePrecision: {1.08568*10^15, -1.08568*10^15, 2.04979*10^13, 1.98037*10^13, ...


13

Numerics in Mathematica can be as precise as you like. However, precision comes at price; you pay for it in computation time and in additional coding effort. In Mathematica there are several computational classes of non-complex numbers, which form a tree like this. The computation you made was made with machine reals because you included 0.5 as a term. ...



Top 50 recent answers are included