diff options
Diffstat (limited to 'packages/base/src/Internal/Algorithms.hs')
| -rw-r--r-- | packages/base/src/Internal/Algorithms.hs | 2 |
1 files changed, 1 insertions, 1 deletions
diff --git a/packages/base/src/Internal/Algorithms.hs b/packages/base/src/Internal/Algorithms.hs index 1ca373f..f5bddc6 100644 --- a/packages/base/src/Internal/Algorithms.hs +++ b/packages/base/src/Internal/Algorithms.hs | |||
| @@ -522,7 +522,7 @@ eig = {-# SCC "eig" #-} eig' | |||
| 522 | -- Eigenvalues are represented as pairs of alpha, beta, where eigenvalue = alpha / beta. Alpha is always | 522 | -- Eigenvalues are represented as pairs of alpha, beta, where eigenvalue = alpha / beta. Alpha is always |
| 523 | -- complex, but betas has the same type as the input matrix. | 523 | -- complex, but betas has the same type as the input matrix. |
| 524 | -- | 524 | -- |
| 525 | -- If @(alphas, betas, v) = geig a b@, then @a \<> v == diag (alphas / betas) \<> b \<> v@ | 525 | -- If @(alphas, betas, v) = geig a b@, then @a \<> v == b \<> v \<> diag (alphas / betas)@ |
| 526 | -- | 526 | -- |
| 527 | -- Note that beta can be 0 and that has reasonable interpretation. | 527 | -- Note that beta can be 0 and that has reasonable interpretation. |
| 528 | geig :: Field t => Matrix t -> Matrix t -> (Vector (Complex Double), Vector t, Matrix (Complex Double)) | 528 | geig :: Field t => Matrix t -> Matrix t -> (Vector (Complex Double), Vector t, Matrix (Complex Double)) |