coarser level matrix update during resetup process
Nobody has claimed this yet.
- Dominant language
- Cuda
- Stars
- 692
- Forks
- 197
- PR merge metrics
- No merged PRs in 30d
Description
Hi,
I notice that if I update the matrix entries value and keep the matrix sparsity pattern the same, I may able to use the resetup function to reuse the previous step coarsening structure. The structure_reuse_levels option needs to be set to "-1".
My question is although the structure is reused, computeAOperator function will still re-compute the sparsity before fill the Ac matrix, is it possible to skip this sparsity recomputing?
Thanks
Frank
Contributor guide
No contributing guide indexed for this repository
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
Start by tracing the resetup process with structure_reuse_levels set to "-1", then inspect the computeAOperator step that recomputes sparsity before filling the Ac matrix. Determine whether unchanged matrix sparsity can safely bypass that recomputation, and verify the behavior for updated matrix values with the same sparsity pattern.
Written by the indexing model from the issue text.
Assessment
- Domain
- hpc, performance
- Issue type
- Feature
- Difficulty
- 4/5
- Estimated time
- 3-5 days
- Activity status
- Stale
- Clarity
- Mostly clear
- Newbie friendliness
- 30/100