diff --git a/.gitattributes b/.gitattributes index e9f7119..46d8444 100644 --- a/.gitattributes +++ b/.gitattributes @@ -1,3 +1,5 @@ input/AB->ABBBBBBBBB.txt filter=lfs diff=lfs merge=lfs -text input/AB->ABBBBBBB.txt filter=lfs diff=lfs merge=lfs -text -*.zip filter=lfs diff=lfs merge=lfs -text +*.zip filter=lfs diff=lfs merge=lfs +*.gif filter=lfs diff=lfs merge=lfs +*.jld2 filter=lfs diff=lfs merge=lfs diff --git a/images/contour_plot_congruent_in_photons.gif b/images/contour_plot_congruent_in_photons.gif new file mode 100644 index 0000000..caa9700 --- /dev/null +++ b/images/contour_plot_congruent_in_photons.gif @@ -0,0 +1,3 @@ +version https://git-lfs.github.com/spec/v1 +oid sha256:ad271f33b35e3e637f0614171201d3b31f81e651a0f699ffa7b529f2e83058ca +size 5781188 diff --git a/notebooks/congruent_ph_eval.ipynb b/notebooks/congruent_ph_eval.ipynb new file mode 100644 index 0000000..9baef6d --- /dev/null +++ b/notebooks/congruent_ph_eval.ipynb @@ -0,0 +1,143 @@ +{ + "cells": [ + { + "cell_type": "code", + "execution_count": 1, + "metadata": {}, + "outputs": [], + "source": [ + "using JLD2\n", + "using NamedDims\n", + "using QEDbase\n", + "using Plots" + ] + }, + { + "cell_type": "code", + "execution_count": 2, + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "3-element Vector{Symbol}:\n", + " :omegas\n", + " :phis\n", + " :thetas" + ] + }, + "metadata": {}, + "output_type": "display_data" + } + ], + "source": [ + "data = Dict()\n", + "for i in 1:4\n", + " @load \"../results/$(i)_congruent_photons_grid.jld2\" results\n", + " data[i] = results\n", + "end\n", + "\n", + "@load \"../results/1_congruent_photons_grid.jld2\" omegas phis thetas" + ] + }, + { + "cell_type": "code", + "execution_count": 3, + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "5-element Vector{Float64}:\n", + " 0.0\n", + " 1.5707963267948966\n", + " 3.141592653589793\n", + " 4.71238898038469\n", + " 6.283185307179586" + ] + }, + "metadata": {}, + "output_type": "display_data" + } + ], + "source": [ + "# https://stackoverflow.com/questions/58962826/plotting-with-a-custom-scale-0-%CF%80-2-%CF%80-3%CF%80-2-2%CF%80-in-julia\n", + "isintegralmultipleof(f, x) = f != 0 && x / f ≈ round(x / f)\n", + "\n", + "\n", + "function pifrac(x)\n", + " if isintegralmultipleof(pi / 2, x)\n", + " n = Int(round(2 * x / π))\n", + " return iseven(n) ? \"$(n ÷ 2) \\\\pi\" : \"$(n) \\\\pi / 2\"\n", + " end\n", + " return string(x)\n", + "end\n", + "\n", + "pirange = [0:π/2:2π;]" + ] + }, + { + "cell_type": "code", + "execution_count": 5, + "metadata": {}, + "outputs": [ + { + "name": "stderr", + "output_type": "stream", + "text": [ + "┌ Info: Saved animation to /tmp/jl_Jyze3K1gdA.gif\n", + "└ @ Plots /home/antonr/.julia/packages/Plots/ju9dp/src/animation.jl:156\n" + ] + }, + { + "data": { + "text/html": [ + "" + ], + "text/plain": [ + "Plots.AnimatedGif(\"/tmp/jl_Jyze3K1gdA.gif\")" + ] + }, + "metadata": {}, + "output_type": "display_data" + } + ], + "source": [ + "@gif for i in 1:length(omegas)\n", + " contours = Vector()\n", + " for j in 1:4\n", + " push!(\n", + " contours, \n", + " contour(\n", + " thetas, \n", + " phis,\n", + " data[j][omegas=i], \n", + " xticks = (pirange, pifrac.(pirange)),\n", + " yticks = (pirange, pifrac.(pirange)),\n", + " xlabel=\"θ\",\n", + " ylabel=\"ϕ\",\n", + " title=\"$j incoming photons\\nω = e$(round(log10(omegas[i]), digits=2))\"\n", + " )\n", + " )\n", + " end\n", + " plot(contours..., layout = (2, 2), size = (1280, 720))\n", + "end fps=16" + ] + } + ], + "metadata": { + "kernelspec": { + "display_name": "Julia 1.10.4", + "language": "julia", + "name": "julia-1.10" + }, + "language_info": { + "file_extension": ".jl", + "mimetype": "application/julia", + "name": "julia", + "version": "1.10.4" + } + }, + "nbformat": 4, + "nbformat_minor": 2 +} diff --git a/results/1_congruent_photons_grid.jld2 b/results/1_congruent_photons_grid.jld2 new file mode 100644 index 0000000..cafd995 --- /dev/null +++ b/results/1_congruent_photons_grid.jld2 @@ -0,0 +1,3 @@ +version https://git-lfs.github.com/spec/v1 +oid sha256:24392cae9b598470b053349e421731d7bf5e7e6be0b160e26d8fae947df1cd01 +size 536894196 diff --git a/results/2_congruent_photons_grid.jld2 b/results/2_congruent_photons_grid.jld2 new file mode 100644 index 0000000..d8c793c --- /dev/null +++ b/results/2_congruent_photons_grid.jld2 @@ -0,0 +1,3 @@ +version https://git-lfs.github.com/spec/v1 +oid sha256:e8eb849bfcacdeac9630db9aabfcc7523607c78f98567959f86d75b9a05c2a69 +size 536894196 diff --git a/results/3_congruent_photons_grid.jld2 b/results/3_congruent_photons_grid.jld2 new file mode 100644 index 0000000..a7ccde7 --- /dev/null +++ b/results/3_congruent_photons_grid.jld2 @@ -0,0 +1,3 @@ +version https://git-lfs.github.com/spec/v1 +oid sha256:6855909a3c4435e20063151fab49599f19d79f4b8e55ace15d207d38963eda2a +size 536894196 diff --git a/results/4_congruent_photons_grid.jld2 b/results/4_congruent_photons_grid.jld2 new file mode 100644 index 0000000..91e731a --- /dev/null +++ b/results/4_congruent_photons_grid.jld2 @@ -0,0 +1,3 @@ +version https://git-lfs.github.com/spec/v1 +oid sha256:bfb4461b365c05c596c62a120b083ce5947c59749c03ef2a6bf71f1f42128e88 +size 536894196 diff --git a/results/5_congruent_photons_grid.jld2 b/results/5_congruent_photons_grid.jld2 new file mode 100644 index 0000000..d134e4c --- /dev/null +++ b/results/5_congruent_photons_grid.jld2 @@ -0,0 +1,3 @@ +version https://git-lfs.github.com/spec/v1 +oid sha256:33e553e072962a2657348a7d69d848a70469aa39e340f9f685f9ae61a32f9cde +size 536894196