Skip to content

Matrix-free evaluator for monodomain problem #150

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.

Sign up for GitHub

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

Jump to bottom
Merged
fdrmrc merged 3 commits into from
Jan 28, 2026
Merged

Matrix-free evaluator for monodomain problem #150

fdrmrc merged 3 commits into from
Jan 28, 2026

Conversation

@fdrmrc
Copy link
Owner

@fdrmrc fdrmrc commented Dec 30, 2025 •
edited
Loading

This PR allows to exploit, on the finest level of the multigrid hierarchy (and within smoothers application), matrix-free evaluators for the monodomain problem.

  • Add model parameters in Utils::Physics namespace
  • Implement matrix-free evaluator for interior penalty DG form associated with monodomain problem
  • Use matrix-free action when applying smoothers and system matrix
  • Use matrix-free action in rhs evaluations

TODOs:

  • Adapt 3D test with matrix-free
  • Test realistic example on cluster(s)
  • Remove code duplication.

We should parametrize all solvers/models through ParameterHandlers. Let's keep everything hardcoded for the time being.

@fdrmrc fdrmrc mentioned this pull request Dec 30, 2025
Closed
@fdrmrc fdrmrc force-pushed the matrix_free_monodomain branch from 028957d to 7ad71f9 Compare December 31, 2025 12:48
@fdrmrc fdrmrc added enhancement New feature or request design Related to code design parallel Related to distributed computing Applications Related to applications/examples snippets multigrid Related to multigrid labels Jan 7, 2026
@fdrmrc fdrmrc marked this pull request as ready for review January 26, 2026 15:46
@fdrmrc
Use matrix-free action on finer level (and within Chebyshev smoother)…
d08ac21 
…. Apply to mondomain problems.

Use exact same penalty for matrix-free and matrix-based variant. Assemble finest multigrid matrix properly. Use <filesystem> to handle output results.

Improve rhs evaluation using matrix-free

Improve timers for setup phase

Improve hierarchy setup and relative timers
@fdrmrc fdrmrc force-pushed the matrix_free_monodomain branch from 4b63cf0 to d08ac21 Compare January 28, 2026 17:46
@fdrmrc
Copy link
Owner Author

fdrmrc commented Jan 28, 2026

I’ll merge this as-is and handle solver parametrization in a follow-up

@fdrmrc fdrmrc merged commit 5bd676e into main Jan 28, 2026
4 checks passed
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Reviewers

No reviews

Assignees

No one assigned

Labels

Applications Related to applications/examples snippets design Related to code design enhancement New feature or request multigrid Related to multigrid parallel Related to distributed computing

Projects

None yet

Milestone

No milestone

Development

Successfully merging this pull request may close these issues.

1 participant

@fdrmrc