Skip to main content

All 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 ...
chunckz7's user avatar
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 ...
chris's user avatar
  • 23.4k
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 ...
matheorem's user avatar
  • 17.7k
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 ...
Dhwister's user avatar
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,-...
FMulder's user avatar
  • 43