Skip to content
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

Refactoring pH and matrix solvers #2184

Merged
merged 14 commits into from Nov 28, 2023
Merged

Refactoring pH and matrix solvers #2184

merged 14 commits into from Nov 28, 2023

Conversation

steff-mueller
Copy link
Contributor

@steff-mueller steff-mueller commented Sep 21, 2023

  • [riccati] Add solve_pos_ricc_dense and bindings for scipy and slycot
  • [iosys] add pr_c_dense and pr_o_dense Gramians
  • [spectral_factor] use pr_o_dense Gramian instead of pr_o_lrcf because we need the full dense solution anyway
  • [iosys] use dense solvers in from_passive_LTIModel instead of lrcf ones
  • [iosys] add to_abcde_matrices to LTIModel such that the standard ABCDE matrices can still be exported if a child class overrides to_matrices (like for instance PHLTIModel)
  • [iosys] add format parameter to to_matrices methods
  • [spectral_factor] use LTIModel.to_abcde_matrices instead of converting the matrices individually

@steff-mueller steff-mueller added the pr:change Change in existing functionality label Sep 22, 2023
@steff-mueller
Copy link
Contributor Author

pre-commit.ci autofix

@steff-mueller steff-mueller marked this pull request as ready for review September 22, 2023 09:37
src/pymor/bindings/scipy.py Outdated Show resolved Hide resolved
src/pymordemos/phlti.py Show resolved Hide resolved
src/pymortests/algorithms/riccati.py Outdated Show resolved Hide resolved
src/pymortests/algorithms/riccati.py Show resolved Hide resolved
@pmli pmli added this to the 2023.2 milestone Nov 2, 2023
@pmli pmli self-assigned this Nov 2, 2023
@pmli
Copy link
Member

pmli commented Nov 10, 2023

Hi @steff-mueller, do you have time to work on this before the hard freeze (Nov 30)?

steff-mueller and others added 8 commits November 22, 2023 15:13
Then we can assume that the solution X to the Riccati equation is nonsingular.
…atrices can still be exported if a child class overrides `to_matrices` (like for instance `PHLTIModel`)

- add `format` parameter to to_matrices methods
Co-authored-by: Petar Mlinarić <mlinaric@vt.edu>
@steff-mueller
Copy link
Contributor Author

pre-commit.ci autofix

Copy link
Member

@pmli pmli left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Looks good to me, just found a typo.

src/pymor/models/iosys.py Outdated Show resolved Hide resolved
Co-authored-by: Petar Mlinarić <mlinaric@vt.edu>
@pmli pmli added this pull request to the merge queue Nov 27, 2023
Merged via the queue into main with commit f13dc29 Nov 28, 2023
18 checks passed
@pmli pmli deleted the refactor-ph branch November 28, 2023 00:18
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
pr:change Change in existing functionality
Projects
None yet
Development

Successfully merging this pull request may close these issues.

None yet

3 participants