[petsc-users] Auto sparsity detection?
Zou, Ling
lzou at anl.gov
Thu Jan 16 20:50:34 CST 2025
Hi all,
Does PETSc has some automatic matrix sparsity detection algorithm available?
Something like: https://urldefense.us/v3/__https://docs.sciml.ai/NonlinearSolve/stable/basics/sparsity_detection/__;!!G_uCfscf7eWS!ccEx6zmuNrVADqtN50hO2N0k4Qs-A70nztAjMLu-JElnjhK5w84BpYC8CAINd6KihSxaS2rx_LgpqUVM49U$
The background is that I use finite differencing plus matrix coloring to (efficiently) get the Jacobian.
For the matrix coloring part, I color the matrix based on mesh connectivity and variable dependencies, which is not bad, but just try to be lazy to even eliminating this part.
A related but different question, how much does PETSc support automatic differentiation?
I see some old paper:
https://ftp.mcs.anl.gov/pub/tech_reports/reports/P922.pdf
and discussion in the roadmap:
https://urldefense.us/v3/__https://petsc.org/release/community/roadmap/__;!!G_uCfscf7eWS!ccEx6zmuNrVADqtN50hO2N0k4Qs-A70nztAjMLu-JElnjhK5w84BpYC8CAINd6KihSxaS2rx_Lgpw6v6hKE$
I am thinking that if AD works so I don’t even need to do finite differencing Jacobian, or have it as another option.
Best,
-Ling
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20250117/ba6891e5/attachment-0001.html>
More information about the petsc-users
mailing list