Get a human-readable description of the current benchmark mode.

Apply dark theme styling to a matplotlib axis.

Apply dark theme styling to a matplotlib axis.

Compute gravitational acceleration for all particles in-place (full version). Non-symmetric loop: compute all N² forces explicitly.

Compute gravitational acceleration with blocked loops for better cache utilization. For N ≳ a few thousand, blocking improves L1/L2 cach

Compute gravitational acceleration for all particles in-place (symmetric version). This is the computationally intensive O(N²) part of t

Apply light theme styling to a matplotlib axis.

Apply light theme styling to a matplotlib axis.

Compute total kinetic energy of the system. KE = Σ(½ × m × v²) for all particles Parameters: ----------- m, u, v, w : n

Compute total gravitational potential energy of the system. PE = -G × Σ(mᵢ × mⱼ / rᵢⱼ) for all particle pairs Parameters: -

Generate initial conditions for a rotating stellar disk. Creates a 2D disk with small thickness in the z-direction. Particles are di

Find a subset of 'target_size' runs that are consistent with each other. This is the key fix: we look for ANY subset of 5 runs that a

Save a two-panel snapshot matching the live visualization. Creates a standalone figure with the same layout as the live visualization:

Save a two-panel snapshot matching the live visualization. Creates a standalone figure with the same layout as the live visualization:

Run the Julia script and measure execution time.

Run the MATLAB script and measure execution time.

Run the Python script and measure execution time.

Generate initial conditions for a spiral galaxy disk. Creates an exponential disk with logarithmic spiral arm overdensities. This pr

Generate initial conditions for a spiral galaxy disk. Creates an exponential disk with logarithmic spiral arm overdensities. This pr

Compute gravitational acceleration for all particles. This is the computationally intensive O(N²) part of the simulation. Uses Plumm

Generate initial conditions for a simple galaxy disk. Creates a disk with exponential-like density profile and circular orbits. This

Generate initial conditions for a rotating stellar disk.

Generate initial conditions for a spiral galaxy disk.

Generate initial conditions for the N-body simulation. Creates different particle distributions based on the selected scenario. All

Generate initial conditions for the N-body simulation. Creates different particle distributions based on the selected scenario. All

Compute total kinetic energy of the system. KE = Σ(½ × m × v²) for all particles Parameters: ----------- m, u, v, w : n

Compute total gravitational potential energy of the system. PE = -G × Σ(mᵢ × mⱼ / rᵢⱼ) for all particle pairs Parameters: -

Build script paths based on user's language selection. Returns: -------- dict : Dictionary mapping language

Check if the latest run is consistent with the existing pattern. This is used to provide real-time feedback during data collection.

Generate a bar plot comparing execution times.

Generate a detailed benchmark report.

Get version information for each programming language.

Collect system information for the benchmark report.

Conditional Numba JIT decorator. If Numba is available, applies JIT compilation for performance. If not available, returns the origi

Conditional Numba JIT decorator. If Numba is available, applies JIT compilation for performance. If not available, returns the origi

Run benchmark for a specific language implementation.

Selects the fastest 'target_size' runs from a larger set. This is used for the 'best_of_10' mode. Parameters:

Selects the most consistent subset of runs from a larger set. This is useful when outlier detection fails to achieve consistency.

Create and configure the main visualization figure. Sets up a 3-panel layout: positions (top 2/3) and energy (bottom 1/3). Retu

Create and configure the main visualization figure. Sets up a 3-panel layout: positions (top 2/3) and energy (bottom 1/3). Retu

Apply velocity boundary conditions at inlet and outlet. Maintains constant velocity at domain boundaries.

Apply velocity boundary conditions at inlet and outlet with OPTIMIZED vectorized operations. Computes all 9 directions at once using broadcas

Apply periodic boundary conditions to the distribution functions. This ensures particles leaving one edge of the domain re-enter from the opp

Perform the collision step using the BGK (Bhatnagar-Gross-Krook) approximation. Distribution functions relax toward their equilibrium values.

Perform collision step in-place using 2D scratch buffer to avoid 3D Feq allocation. OPTIMIZED VERSION with precomputed constants for better p

Compute macroscopic fluid variables (density and velocity) from the distribution functions using moment integrals.

Compute macroscopic fluid variables using OPTIMIZED single-pass computation. Fuses density and momentum calculations into one loop over direc

Create a boolean mask representing the cylinder in the computational domain.

Create a boolean mask representing the cylinder in the computational domain.

Export initial conditions to a portable JSON file.

Generate initial conditions for the specified scenario.

Generate initial conditions for a Plummer sphere.

Handle the no-slip boundary condition at the cylinder surface. Implements the bounce-back method where particles hitting the cylinder rev

Import initial conditions from a JSON file generated by nbody_ic.py. Parameters: ----------- filename : str Path to the

Import initial conditions from a JSON file generated by nbody_ic.py. Parameters: ----------- filename : str Path to the

Initialize the flow field with uniform density and inlet velocity. Sets up the initial distribution functions F based on equilibrium conditio

Initialize the flow field with uniform density and inlet velocity. Sets up the initial distribution functions F based on equilibrium conditio

Initialize all simulation parameters including domain size, lattice properties, cylinder geometry, and fluid properties.

Initialize all simulation parameters including domain size, lattice properties, cylinder geometry, and fluid properties.

Main simulation loop implementing the Lattice Boltzmann Method. This function orchestrates the complete LBM simulation by: 1. Settin

Main simulation loop implementing the Lattice Boltzmann Method. This function orchestrates the complete LBM simulation by: 1. Setting up

Main function to run the benchmarking.

Main simulation function. Implements the leapfrog integration scheme: 1. Half-kick: update velocities by half a time step 2. Dri

Main simulation function (optimized version). Implements the leapfrog integration scheme using optimized kernels: 1. Preallocates ac

Main function to generate and export initial conditions.

Precompute incoming particle masks for bounce-back boundary conditions. Uses the CORRECTED algorithm that matches the baseline exactly.

Print a summary of the benchmark results.

Linear blend in RdBu-like diverging scheme (blue->white->red) */

read whole file to string (NUL-terminated) */