All Questions
5 questions
1
vote
0
answers
129
views
Kahan Compensated Summation
Is there an example/library that implements Kahan summation in Mathematica for matrix-matrix/vector-matrix multiplication?
Having some trouble adapting something like this Compilation of Total with ...
6
votes
2
answers
275
views
Catastrophic loss of accuracy in Orthogonalize
Context
In connection to this question
I am interested in orthogonalizing known matrices.
As a test case, let us consider the definite positive $15 \times 15$ matrix ...
16
votes
2
answers
5k
views
more numerically accurate inverse matrix
I encountered the following matrix
mat = {{2,
2.161209223472559` + 1.682941969615793` I}, {2.161209223472559` -
1.682941969615793` I, 2}}
and ...
0
votes
1
answer
778
views
Precision of Eigensystem? [closed]
I was using Eigensystem to obtain the rotation matrix. However, I find out Mathematica does not fully diagonalize my matrix (or say not precise enough). My matrix ...
4
votes
1
answer
488
views
Wrong generalized eigenvalues of large numeric matrices
I've been working on a toy model of a variational method of a Hamiltonian H with an overlap matrix S. The point is that when H and S are small enough, Mathematica gives the correct eigenvalues (≈ -1,-...