Add profiling script

Author: langestefan

examples/Project.toml

@@ -7,6 +7,7 @@ ModelingToolkit = "961ee093-0014-501f-94e3-6117800e7a78"
 OhMyThreads = "67456a42-1dca-4109-a031-0a68de7e3ad5"
 OrdinaryDiffEq = "1dea7af3-3e70-54e6-95c3-0bf5283fa5ed"
 Plots = "91a5bcdd-55d7-5caf-9e0b-520d859cae80"
+ProfileView = "c46f51b8-102a-5cf2-8d2c-8597cb0e0da7"
 SolarPosition = "5b9d1343-a731-5a90-8730-7bf8d89bf3eb"
 StructArrays = "09ab397b-f2b6-538f-b94a-2f83cf4a842a"
 

examples/profiling.jl

@@ -0,0 +1,22 @@
+"""Profiling example for SolarPosition.jl using ProfileView."""
+
+using Pkg
+Pkg.activate(@__DIR__)
+Pkg.instantiate()
+
+using Dates
+using SolarPosition
+
+# define observer location (latitude, longitude, altitude in meters)
+obs = Observer(52.358134610343214, 4.881269505489815, 0.0)  # Van Gogh Museum
+
+# a whole year of hourly timestamps for profiling
+times = collect(DateTime(2023):Minute(1):DateTime(2024))
+
+# warm up the function to compile
+solar_position(obs, times, PSA())
+
+# profile the PSA algorithm with multiple iterations for better sampling
+@profview for _ = 1:1000
+    solar_position(obs, times, PSA())
+end