Therefore

Therefore is a deterministic testbed that couples ~higher~ order discretization schemes and two separate iteration methods for transient radiation transport to various hardware acceleration backends in a sparse linear algebra framework. We started with a full Python based implementation using Numba, SciPy and CuPy (1D, transient mono-energetic) then moved into a C++ (energy-dependent) based implementation using cuSOLVER, hipSOLVER, and LAPACK.

We have published various posters, papers, and presentations about this work at:

• 2024 18th Conference on Iterative Methods (SIAM) (presentation)
• 2023 UCSD HEDS (poster)
• 2023 International ANS M&C (paper & presentation)
• 2022 CEMeNT AST Review meeting (poster)

This transient reeds problem which was solved on a GPU is an example of the type of transient transport we can do!

All of our code is open source and available on github