SciML / SciML/OrdinaryDiffEq.jl

iterative linear solvers in the context of conservation laws

Open
#2,139 0 comments 0 reactions 0 assignees View on GitHub

Nobody has claimed this yet.

Dominant language
Julia
Stars
684
Forks
276
Avg merge
1d 5h
Merged PRs (30d)
198

Description

https://link.springer.com/article/10.1007/s10543-023-00992-w

Q: Do iterative methods (in this case Newton) preserve nonlinear invariants? A: No. Q: What should one do about it? A: Relax the time step but leave the iterations alone.

https://link.springer.com/article/10.1007/s10915-022-01923-7

Q: Do iterative methods preserve basic properties of PDE discretizations such as conservation? A: Some do, some don’t. Q: So if we pick a conservative one, then all is fine? A: No, because consistency is lost (time is slowed down).

https://epubs.siam.org/doi/abs/10.1137/22M1503348

An extension of the previous paper. A bit harder to read but contains some new interesting results, not least on the connection between Krylov methods and pseudotime stepping.

https://onlinelibrary.wiley.com/doi/abs/10.1002/pamm.202200157

A deeper dive into the inconsistency (time dilation) caused by Krylov methods. I wrote this paper in about 48 hours so it misses a crucial simplification in the main result that I can tell you about if you want to know more.

Contributor guide

Open the contributing guide

First steps

  1. Read the whole issue, then the project's contributing guide.
  2. Comment on the issue to say you are picking it up — it saves two people doing the same work.
  3. Fork the repository and make your change on a branch.
  4. Open a pull request that references the issue number.

Research direction

Start by reading the four linked papers, especially the discussion of Newton, conservative iterative methods, Krylov methods, and pseudotime stepping. The issue does not name files, tests, entry points, or a concrete acceptance condition, so the first step is to clarify which solver behavior should change and how it should be validated.

Written by the indexing model from the issue text.

Assessment

Tech stack
julia
Domain
backend
Issue type
Feature
Difficulty
5/5
Estimated time
Over a week
Activity status
Stale
Clarity
Needs clarification
Newbie friendliness
15/100

Get new issues in your inbox

A short digest of beginner-friendly GitHub issues.