Skip to content

Commit

Permalink
Change image layout
Browse files Browse the repository at this point in the history
  • Loading branch information
annatartaglia authored Oct 25, 2024
1 parent 5a5dd35 commit bb02992
Showing 1 changed file with 4 additions and 2 deletions.
6 changes: 4 additions & 2 deletions docs/src/benchmarks.md
Original file line number Diff line number Diff line change
Expand Up @@ -8,13 +8,15 @@ The Julia implementation of the EHT scattering framework takes advantage of Juli

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.

![](assets/speed_dipole.png) ![](assets/speed_boxcar.png) ![](assets/speed_vonmises.png)
| ![](assets/speed_dipole.png) | ![](assets/speed_boxcar.png) | ![](assets/speed_vonmises.png) |
|:----------------------------:|:----------------------------:|:------------------------------:|

With this significant speed up, the Julia implementation paves the pathway for the joint modeling of scattering parameters and reconstructed images in self-consistent way.

## 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}$, which is well under scattering kernel uncertainty limits.

![](assets/kernel87.png) ![](assets/kernel36.png)
| ![](assets/kernel87.png) | ![](assets/kernel36.png) |
|:------------------------:|:------------------------:|

0 comments on commit bb02992

Please sign in to comment.