Propag source code excerpts
memory.pdf contains the code for memory allocation, which is quite an issue in a large, experimental, parallelized program.
status.pdf contains the code for saving and loading checkpoint files ("status dumps"). These files are written at regular intervals by the program and are used to continue a simulation if it is stopped unexpectedly. This happens relatively often in this kind of software, because it typically runs for hours or days, and on bleeding-edge hardware that breaks more often than mainstream computers. Checkpoint files can also be useful to start several simulations from a common starting point, or to interrupt a simulation and continue with different parameter settings. A useful example is switching between monodomain and bidomain mode. Checkpointing safely and efficiently on a parallel computer requires some extra attention, as you may see in this document.
bernus.pdf is our implementation of the membrane model by Bernus et al, AJP-H 282:H2296-2308, 2002. It uses Rush-Larsen integration for gating variables, and tabulates some of the computationally expensive functions that are involved.
bicgstab.pdf implements parallel BiCGStab, BiCGStab-p, QMRCGStab, QMRCGStab-2, and CG solvers. The code is only slightly dependent on other parts of the program. The algorithm contains some novelties, such as restarting and constraining of the residuals to avoid divergence when solving very large models without grounding. (Grounding very large bidomain models is somewhat like immobilizing a grizzly by holding its ear - it hurts a lot but it won't immobilize very well).