From 6f38d9e346f2e5a613ee063f8dab47ef88e14a97 Mon Sep 17 00:00:00 2001 From: Anna Tartaglia <111704070+annatartaglia@users.noreply.github.com> Date: Wed, 23 Oct 2024 10:50:24 -0400 Subject: [PATCH] Update benchmarks.md --- docs/src/benchmarks.md | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/docs/src/benchmarks.md b/docs/src/benchmarks.md index 973502f..e31cede 100644 --- a/docs/src/benchmarks.md +++ b/docs/src/benchmarks.md @@ -1,10 +1,10 @@ ```@meta CurrentModule = ScatteringOptics ``` -#Benchmarks +# Benchmarks The Julia implementation of the EHT scattering framework takes advantage of Julia's fast Auto-Differentiation (AD) capabilities, which accelerate the simulation of interstellar scattering. As the overall imaging process can consist of thousands to tens of thousands of iterations in order to optimize scattering parameters and the final output image, even a small imprdovements in simulation speed can result in significant cuts in runtime. It is then crucial that our implementation performs faster than the existing python framework, `eht-imaging`, and produces consistent results. Here, we present benchmarks between the two implementations. -##Speed +## Speed With Julia's speed improvements, the `ScatteringOptics.jl` scattering kernel loads 100 times faster than the `eht-imaging` kernel. The scattering kernel may be called upon to compute visibilities for a given sample of Fourier space points. In these computations, `ScatteringOptics.jl` also exhibits speed improvements up to 100 times over the three included models (Dipole, Boxcar, and Von Mises; see [Use Non-default Models](@ref) for more info) and varying observing wavelengths. These computation runtimes are compared over a range of sample in the plots below. @@ -12,7 +12,7 @@ With Julia's speed improvements, the `ScatteringOptics.jl` scattering kernel loa With this significant speed up, the Julia implementation allows for the joint modeling of scattering parameters and reconstructed images in self-consistent way. -##Consistency +## Consistency Our scattering kernel produces visibilities that are consistent with those computed by `eht-imaging` across all 3 kernel models as well as the range of standard observing wavelengths in radio astronomy. Below, we demonstrate this by plotting the Dipole kernel over the Fourier space major and minor axis for both the Julia and Python implementations in one plot per observing wavelength. The lines completely overlap, indicating complete consistency in results. All fractional errors between the two implementations range on the order of $10^{-6}$ and $10^{-7}$ (Figure \ref{kernelerrors}), which is well under scattering kernel uncertainty limits.