Added BLAS/LAPACK support for matrix inverse inv() #2
Conversation
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
only implemented if a BLAS library isused since a solver for linear equations is needed. The Pade method is used if you use expm() instead of the standard exp() for matrices or the super operator. It is not obvious when the exp() function which uses blocking of the matrix and diagonalization or the expm() function using the Pade method is faster. This seems to vary strongly. prop() for operators is still implemented using diagonalization.
Merge branch 'main' of https://github.com/pygamma-mrs/gamma
…floq2_op. No standard functionality has been changed. These are just additional functions that can be used as drop-in replacements for the exp function if desired. They are only different from exp if a BLAS library is used.
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
Add this suggestion to a batch that can be applied as a single commit.
This suggestion is invalid because no changes were made to the code.
Suggestions cannot be applied while the pull request is closed.
Suggestions cannot be applied while viewing a subset of changes.
Only one suggestion per line can be applied in a batch.
Add this suggestion to a batch that can be applied as a single commit.
Applying suggestions on deleted lines is not supported.
You must change the existing code in this line in order to create a valid suggestion.
Outdated suggestions cannot be applied.
This suggestion has been applied or marked resolved.
Suggestions cannot be applied from pending reviews.
Suggestions cannot be applied on multi-line comments.
Suggestions cannot be applied while the pull request is queued to merge.
Suggestion cannot be applied right now. Please check back later.
These are the changes to support BLAS/LAPACK for matrix inversion. This only changes n_matrix.cc and h_matrix.cc and is fully compatible with the BLAS/LAPACK support as it was done for diagonalization.
Also fixed a parsing error in MultiSys where Exch() produced in error for more than 10 components.