******* Branson ******* This is the documentation for the FCR Benchmark Branson - 3D hohlraum multi-node using domain decomposition. https://github.com/lanl/branson Purpose ======= From [Branson]_: Branson is not an acronym. Branson is a proxy application for parallel Monte Carlo transport. It contains a particle passing method for domain decomposition. Characteristics =============== Problems -------- The benchmark performance problem is a multi-node 3D hohlraum problem that is meant to be run with a 30 group build of Branson. It is in domain decomposition mode which means there is both MPI communication for particle movement and end of cycle reductions. Four problem configurations are provided: #. CPU, decomposed, history, SoA #. GPU, decomposed, history, SoA #. GPU, decomposed, event, SoA #. GPU, decomposed, event, AoS Figure of Merit --------------- The Figure of Merit is defined as particles/second and is obtained by dividing the number of particles in the problem divided by the `Total transport` value. This value is labeled "Photons Per Second (FOM):" in Branson's output. Source code modifications ========================= Please see :ref:`GlobalRunRules` for general guidance on allowed modifications. Building ======== Accessing the sources * Clone the FCR branch/tag? from the branson github https://github.com/lanl/branson.git .. code-block:: bash git clone https://github.com/lanl/branson.git cd branson git checkout FCR .. Build requirements: * C/C++ compiler(s) with support for C11 and C++14. * `CMake 3.9X `_ * MPI 3.0+ * `OpenMPI 1.10+ `_ * `mpich `_ * There is only one CMake user option right now: ``CMAKE_BUILD_TYPE`` which can be set on the command line with ``-DCMAKE_BUILD_TYPE=`` and the default is Release. * If cmake has trouble finding your installed TPLs, you can try * appending their locations to ``CMAKE_PREFIX_PATH``, * try running ``ccmake .`` from the build directory and changing the values of build system variables related to TPL locations. * If building a CUDA enabled version of Branson use the ``CUDADIR`` environment variable to specify your CUDA directory. * If building for multi-node runs Metis should be used for mesh partitioning. See README.md from Branson for more details. Single node CPU and single node GPU runs for SSNI should not use Metis. To build metis: .. code-block:: bash cd make config cc= prefix= shared=1 make install .. To build branson: .. code-block:: bash export CXX=`which g++` cd mkdir build cd build cmake -DCMAKE_BUILD_TYPE=Release -DCMAKE_INSTALL_PREFIX= make -j .. Testing the build: .. code-block:: bash cd $build_dir ctest -j 32 .. Running ======= The ``inputs`` folder contains the 3D hohlraum input file. 3D hohlraums and should be run with a 30 group build of Branson (see Special builds section above). The ``3D_hohlraum_multi_node.xml`` problem is meant to be run on multiple nodes. It is run with: .. code-block:: bash mpirun -n .. Memory footprint is the sum of all Branson processes resident set size (or equivalent) on the node. This can be obtained on a CPU system using the following (while the application is in step 2): .. code-block:: bash ps -C BRANSON -o euser,c,pid,ppid,cmd,%cpu,%mem,rss --sort=-rss ps -C BRANSON -o rss | awk '{sum+=$1;} END{print sum/1024/1024;}' .. Validation ========== Branson has two main checks on correctness. The first is a looser check that's meant as a "smoke test" to see if a code change has introduced an error. After every timestep, a summary block is printed: .. code-block:: bash ******************************************************************************** Step: 5 Start Time: 0.04 End Time: 0.05 dt: 0.01 source time: 0.166658 WARNING: use_gpu_transporter set to true but GPU kernel not available, running transport on CPU Total Photons transported: 10632225 Emission E: 4.43314e-05, Source E: 0, Absorption E: 4.1747e-05, Exit E: 2.59802e-06 Pre census E: 3.5321e-07 Post census E: 3.396e-07 Post census Size: 219902 Pre mat E: 0.0130731 Post mat E: 0.0130705 Radiation conservation: -5.83707e-17 Material conservation: -5.8599e-15 Sends posted: 0, sends completed: 0 Receives posted: 0, receives completed: 0 Transport time max/min: 7.31594/7.20329 .. Two lines in the block specifically relate to conservation: .. code-block:: bash Radiation conservation: -5.83707e-17 Material conservation: -5.8599e-15 .. The radiation conservation should capture roughly half of the range of the floating point type compared to the amount of radiation energy in the problem. The standard version of Branson uses double precision for all floating point values in both CPU and GPU versions. For the timestep shown above, there's 4.43314e-5 jerks of energy being emitted and the conservation quantity is -5.837e-17, so the relative accuracy is about 1.0e-12, which is well above half the range of a double. The same check can be done for the material energy conservation: here the total energy in the material at the end of the timestep is 0.0130705 jerks, and the conservation value is -5.8599e-15, representing relative precision of 1.0e-13. As mentioned above, conservation is a relatively loose check as more particles and more cells represent more summmations and more opportunities for loss of precision. This is further complicated by MPI reductions. Still, this check is accurate enough to clearly detect particles that may havbe been lost in a modified MPI scheme (for example). The second check on correctness is much simpler. For any changes to Branson, the code should produce the same temperature in a standard marshak wave problem after 100 cycles. For the `marshak wave input `_ file, the following temperature profile should be reproduced to 3% after 100 cycles, as shown below: .. code-block:: bash Step: 100 Start Time: 0.99 End Time: 1 dt: 0.01 source time: 0.094371 -------- VERBOSE PRINT BLOCK: CELL TEMPERATURE -------- cell T_e T_r abs_E 0 0.9864821 0.98624394 2.3231089e-05 1 0.97376231 0.97335755 2.2986719e-05 2 0.95987812 0.95921396 2.2604072e-05 3 0.94448294 0.94359619 2.223203e-05 4 0.92838247 0.92729361 2.1860113e-05 5 0.91059797 0.90933099 2.1487142e-05 6 0.89041831 0.88903414 2.1098101e-05 7 0.86713097 0.86559489 2.0554045e-05 8 0.83972062 0.83807018 1.9926467e-05 9 0.80754477 0.80583439 1.9216495e-05 10 0.76586319 0.76409724 1.8223846e-05 11 0.71065544 0.70892379 1.6994308e-05 12 0.6190012 0.61733211 1.5009059e-05 13 0.36540211 0.35970671 1.1687053e-05 14 0.016821133 0.016162407 6.3406719e-07 15 0.01 0.0099763705 2.356755e-07 16 0.010000399 0.0099766379 2.3568489e-07 17 0.0099989172 0.0099752306 2.3564998e-07 18 0.010000684 0.0099769858 2.3569162e-07 19 0.009999951 0.0099762996 2.3567434e-07 20 0.0099997415 0.0099761208 2.356694e-07 21 0.010000476 0.0099768182 2.3568672e-07 22 0.0099993136 0.0099756288 2.3565932e-07 23 0.010000237 0.0099765577 2.3568109e-07 24 0.010000281 0.0099765314 2.3568212e-07 ------------------------------------------------------- .. This output is expected as long as the spatial, boundary and region blocks are kept the same in the input file. The IMC method that Branson uses is stochastic so changing the random number seed or the number of particles will produce a slightly different answer, but the difference should not be more than 3% if one million or more particles are used. This test is sensitive to precision changes in Branson as propagating the energy correctly involves many small summations as particle's slowly lose their energy into the material. Example Scalability Results =========================== Memory Usage ============ Strong Scaling on El Capitan ============================ Weak Scaling on El Capitan ========================== References ========== .. [Branson] Alex R. Long, 'Branson', 2026. [Online]. Available: https://github.com/lanl/branson. [Accessed: 06- Jan- 2026]