-
Notifications
You must be signed in to change notification settings - Fork 16
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Mumps solver updates #39
Commits on Aug 15, 2023
-
Configuration menu - View commit details
-
Copy full SHA for 6126a45 - Browse repository at this point
Copy the full SHA 6126a45View commit details -
Configuration menu - View commit details
-
Copy full SHA for 12aeda2 - Browse repository at this point
Copy the full SHA 12aeda2View commit details -
Configuration menu - View commit details
-
Copy full SHA for afcef62 - Browse repository at this point
Copy the full SHA afcef62View commit details -
Configuration menu - View commit details
-
Copy full SHA for 0c4617f - Browse repository at this point
Copy the full SHA 0c4617fView commit details -
Configuration menu - View commit details
-
Copy full SHA for a7a7670 - Browse repository at this point
Copy the full SHA a7a7670View commit details -
Configuration menu - View commit details
-
Copy full SHA for 645bfde - Browse repository at this point
Copy the full SHA 645bfdeView commit details -
Refactor unit tests - now pytest format
The Mumps test_singular test is available now. The Mumps test_1to5_T test is still failing.
Configuration menu - View commit details
-
Copy full SHA for e0e410e - Browse repository at this point
Copy the full SHA e0e410eView commit details -
Configuration menu - View commit details
-
Copy full SHA for 8f94db0 - Browse repository at this point
Copy the full SHA 8f94db0View commit details -
Configuration menu - View commit details
-
Copy full SHA for b3d44b2 - Browse repository at this point
Copy the full SHA b3d44b2View commit details
Commits on Aug 16, 2023
-
Handle unknown MUMPS errors and warnings
Also correctly handle warnings as flags.
Configuration menu - View commit details
-
Copy full SHA for 6fcd2a7 - Browse repository at this point
Copy the full SHA 6fcd2a7View commit details -
Configuration menu - View commit details
-
Copy full SHA for 24c9935 - Browse repository at this point
Copy the full SHA 24c9935View commit details -
Fix misplaced call to Mumps destroy function
There was an issue where if you copied a matrix and the copy passed out of scope, the original matrix could no longer be used. This issue was the reason that the test_1to5_T test was failing.
Configuration menu - View commit details
-
Copy full SHA for 17e5126 - Browse repository at this point
Copy the full SHA 17e5126View commit details -
Configuration menu - View commit details
-
Copy full SHA for 539c330 - Browse repository at this point
Copy the full SHA 539c330View commit details -
Configuration menu - View commit details
-
Copy full SHA for 7cec0a1 - Browse repository at this point
Copy the full SHA 7cec0a1View commit details -
Configuration menu - View commit details
-
Copy full SHA for b22147a - Browse repository at this point
Copy the full SHA b22147aView commit details
Commits on Aug 17, 2023
-
Configuration menu - View commit details
-
Copy full SHA for 26f3cb8 - Browse repository at this point
Copy the full SHA 26f3cb8View commit details -
Update shape, then check accuracy
Otherwise you get ValueError: dimension mismatch when the solver doesn't do the reshaping
Configuration menu - View commit details
-
Copy full SHA for 67351f8 - Browse repository at this point
Copy the full SHA 67351f8View commit details -
Full attribute parity and result parity between Pardiso and Mumps
It's now possible to use the same solver_opts values in SimPEG for both the Pardiso and MUMPS solvers. We now run all the Pardiso test cases for Mumps as well. Running these test cases uncovered a Mumps issue where matrix transposes didn't commute when operating on complex matrices. If you asked the Mumps solver: Ainv = Mumps(A) AinvT_1 = Mumps(A).T AinvT_2 = Mumps(A.T) then the two inverse transposes weren't equal when A is a complex matrix. For small matrices, you can look at the inverse transposes by solving this trivial system: AinvT_1 * np.identity(A.size) != AinvT_2 * np.identity(A.size) This issue occurred because pymatsolver treated the complex case as the conjugate transpose -- but we weren't asking for the conjugate transpose; we were explicitly asking for just the transpose. There was no issue with Mumps itself, just a few extra lines of Fortran code in the Mumps interface. Now that those lines are removed, the Mumps solver and the Pardiso solver operate identically with these transposes.
Configuration menu - View commit details
-
Copy full SHA for 47ef682 - Browse repository at this point
Copy the full SHA 47ef682View commit details -
Configuration menu - View commit details
-
Copy full SHA for 8f5cb6e - Browse repository at this point
Copy the full SHA 8f5cb6eView commit details
Commits on Aug 18, 2023
-
Configuration menu - View commit details
-
Copy full SHA for 3d6be06 - Browse repository at this point
Copy the full SHA 3d6be06View commit details -
Configuration menu - View commit details
-
Copy full SHA for 4836604 - Browse repository at this point
Copy the full SHA 4836604View commit details -
Configuration menu - View commit details
-
Copy full SHA for c848dd8 - Browse repository at this point
Copy the full SHA c848dd8View commit details
Commits on Aug 21, 2023
-
Configuration menu - View commit details
-
Copy full SHA for cde8692 - Browse repository at this point
Copy the full SHA cde8692View commit details