SciML / SciML/ExponentialUtilities.jl
Classical Gram-Schmidt with DGKS re-orthogonalization criterion
Nobody has claimed this yet.
- Dominant language
- Julia
- Stars
- 107
- Forks
- 35
- Avg merge
- 22h 52m
- Merged PRs (30d)
- 11
Description
In Gram-Schmidt iterations, re-orthogonalization is often required to enforce the obtained spanning space is correct up to some precision in many practical applications.
https://www.cs.cornell.edu/~bindel/class/cs6210-f16/lec/2016-11-16.pdf
I didn't find such procedure in this repo yet. Is it because it is not needed in implementing expmv, or is it simply not implemented yet or have I missed it? I am happy to submit a PR if re-orthogonalization is really needed.
Contributor guide
First steps
- Read the whole issue, then the project's contributing guide.
- Comment on the issue to say you are picking it up — it saves two people doing the same work.
- Fork the repository and make your change on a branch.
- Open a pull request that references the issue number.
Research direction
Review the repository's expmv entry point and its Gram-Schmidt iterations, then compare them with the linked DGKS reference. Determine whether re-orthogonalization is required for the intended precision and spanning-space guarantees. Done means the need is resolved and, if required, the behavior is implemented with appropriate validation.
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
- 25/100