add benchmark

Author: mileslucas

bench/Project.toml

@@ -0,0 +1,12 @@
+[deps]
+BenchmarkTools = "6e4b80f9-dd63-53aa-95a3-0cdb28fa8baf"
+CSV = "336ed68f-0bac-5ca0-87d4-7b16caf5d00b"
+DataFrames = "a93c6f00-e57d-5684-b7b6-d8193f3e46c0"
+PSFModels = "9ba017d1-7760-46cd-84a3-1e79e9ae9ddc"
+PyCall = "438e738f-606a-5dbb-bf0a-cddfbfd45ab0"
+
+[compat]
+BenchmarkTools = "0.5"
+CSV = "0.8"
+PSFModels = "0.1.1"
+PyCall = "1"

bench/bench.jl

@@ -0,0 +1,42 @@
+using BenchmarkTools
+using CSV
+using DataFrames
+using PSFModels: Gaussian, AiryDisk, Moffat
+using PyCall
+
+apm = pyimport("astropy.modeling.functional_models")
+
+@info "Starting PSF evaluation benchmark"
+
+julia_models = [
+    Gaussian((5.3, 4.7)), 
+    AiryDisk(12.7), 
+    Moffat(8.1)
+]
+astropy_models = [
+    apm.Gaussian2D(x_stddev=5.3, y_stddev=4.7),
+    apm.AiryDisk2D(radius=12.7),
+    apm.Moffat2D(gamma=8.1, alpha=-1)
+]
+
+names = ("Gaussian", "AiryDisk", "Moffat")
+jl_times = []
+py_times = []
+
+for (jl_mod, py_mod, name) in zip(julia_models, astropy_models, names)
+    jl_time = @belapsed $jl_mod(0.5, 0.7)
+    push!(jl_times, jl_time)
+
+    py_time = @belapsed $py_mod(0.5, 0.7)
+    push!(py_times, py_time)
+    
+    @info "$name" PSFModels=jl_time astropy=py_time
+end
+
+filename = joinpath(@__DIR__, "results.csv")
+DataFrame(name=names, psfmodels=jl_times, astropy=py_times) |> CSV.write(filename)
+
+
+@info "Finished PSF evaluation benchmark. Results saved to $filename"
+
+nothing

bench/results.csv

@@ -0,0 +1,4 @@
+name,psfmodels,astropy
+"(""Gaussian"", ""AiryDisk"", ""Moffat"")",1.5927855711422846e-8,0.000165969
+"(""Gaussian"", ""AiryDisk"", ""Moffat"")",3.7001009081735617e-8,0.000135154
+"(""Gaussian"", ""AiryDisk"", ""Moffat"")",1.192685370741483e-8,0.000134448