[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Permutations in expm
From: |
Marco Caliari |
Subject: |
Permutations in expm |
Date: |
Fri, 25 Jan 2008 13:47:47 +0100 (CET) |
Hi all.
I investigated how the reverse permutation works in expm. The
enclosed patches (against 3.0.0) are much simpler than my previous ones. In
practice, there was an error just in the order of the computation of
iperm. I tested the patches on random permutations of balanced
matrices, both real and complex, up to order 100.
The ChangeLog could be the same
2008-01-25 Marco Caliari <marco.caliari at univr.it>
* dMatrix.cc (Matrix::expm): Correctly perform reverse permutation.
* CMatrix.cc (ComplexMatrix::expm): Likewise.
Best regards,
Marco
dMatrix.cc.patch
Description: Text document
CMatrix.cc.patch
Description: Text document
- Permutations in expm,
Marco Caliari <=