version 1.0.0

Author: Hinamoooon

.gitignore

@@ -1,5 +1,4 @@
 /Manifest.toml
 /docs/build/
-/developper_notes
 
-/example_dev
+/result

Project.toml

@@ -1,7 +1,7 @@
 name = "jlmie"
 uuid = "83f16f31-eeab-45e4-be0e-7fc4231b6e5b"
 authors = ["Tatsuki Hinamoto"]
-version = "0.1.0"
+version = "1.0.0"
 
 [deps]
 SpecialFunctions = "276daf66-3868-5448-9aa4-cd146d93841b"

developper_notes.txt

@@ -0,0 +1,69 @@
+# This code separately visualize the 3D radiation pattern
+# of electric and magnetic resonances.
+# Because the amplitude of scattering intensity is normalized, 
+# the setting of wavelength is not very important...
+
+# import needed packages
+# for jlmie.jl
+include("..\\src\\jlmie.jl")
+using SpecialFunctions
+# for this script
+using Plots
+
+# settings
+# structure
+nmat = 3.5  # refractive index of a sphere
+radius = 100*1e-9
+
+# environment
+nenv = 1.00
+
+# wavelenth (not very important)
+lbd0ED = 618*1e-9  # ED
+lbd0MD = 831*1e-9  # MD
+lbd0 = [lbd0ED lbd0MD]
+
+# angular range
+theta = range(0,π,length=37) # per 5 degrees with length 37
+phi = range(0,2π,length=73) # per 5 degrees with length 73
+
+# other settings
+n = 1
+EorH = 2  # 1: Eelectric, 2: Magnemtic
+Isff = zeros(length(theta),length(phi))
+for i = 1:length(theta)
+    for j = 1:length(phi)
+        Isff[i,j] = jlmie_Isff_EM_n(nmat,radius,lbd0[EorH],nenv,theta[i],phi[j],n)[EorH]
+    end
+end
+Isff = Isff/maximum(Isff)
+mat_x = Isff.* (sin.(theta) .* cos.(phi)')
+mat_y = Isff.* (sin.(theta) .* sin.(phi)')
+mat_z = Isff.* (cos.(theta) .* cos.(phi*0)')
+
+# Just a sphere
+# mat_x = 1 .* (sin.(theta) * cos.(phi)')
+# mat_y = 1 .* (sin.(theta) * sin.(phi)')
+# mat_z = 1 .* (cos.(theta) * cos.(phi*0)')
+
+# show results
+pyplot()
+plt = plot(mat_x,mat_y,mat_z,
+           st=:surface,
+           camera=(30,30), # azimuth, elevate
+           xlabel="x",
+           ylabel="y",
+           zlabel="z",
+           xlims=(-1,1),
+           ylims=(-1,1),
+           zlims=(-1,1),
+)
+display(plt)
+
+# save
+if false
+    runningfilename = splitext(splitpath(@__FILE__)[end])[1]
+    outfilename = ".\\result\\" * runningfilename
+    savefig(plt,outfilename)
+    println("Graph has been saved in result directory")
+end

example_dev/test.jl

@@ -0,0 +1 @@
+using jlmie

example_dev/test2.jl

@@ -0,0 +1,35 @@
+using jlmie
+using Plots
+
+n = 3
+# angular range
+theta = range(0,2π,length=73) # per 5 degrees with length 73
+mu = cos.(theta)
+pin,taun = jlmie_pt(mu,n)
+# show results
+# validation can be done by comparing to ITMO Mie Calculator
+# https://physics.itmo.ru/en/mie
+plt = plot(theta, pin, proj=:polar,
+    # xlabel       = "Wavelength (nm)",
+    # ylabel       = "Scattering efficiency",
+    # legend       = false,
+    label        = "pi",
+    # xlims        =(-3,3),
+    # ylims        =(-3,3),
+    # aspect_ratio =0.5,
+    # title        ="TITLE",
+    # linecolor    =:blue,
+    # linewidth    =5,
+    # linestyle    =:dot,
+    # size         =(400,300),
+) # https://qiita.com/I_ppp/items/dca3552affa6a672e4bd
+plot!(theta, taun, proj=:polar,label="tau")
+display(plt)
+
+# save
+if true
+    runningfilename = splitext(splitpath(@__FILE__)[end])[1]
+    outfilename = ".\\result\\" * runningfilename
+    savefig(plt,outfilename)
+    println("Graph has been saved in result directory")
+end

example_dev/xx01_AngularFunctions.jl

@@ -0,0 +1,65 @@
+# This code separately visualize the 3D radiation pattern
+# of electric and magnetic resonances.
+# Because the amplitude of scattering intensity is normalized, 
+# the setting of wavelength is not very important...
+
+using jlmie
+using Plots
+
+# settings
+# structure
+nmat = 3.5  # refractive index of a sphere
+radius = 100*1e-9
+
+# environment
+nenv = 1.00
+
+# wavelenth (not very important)
+lbd0ED = 618*1e-9  # ED
+lbd0MD = 831*1e-9  # MD
+lbd0 = [lbd0ED lbd0MD]
+
+# angular range
+theta = range(0,π,length=37) # per 5 degrees with length 37
+phi = range(0,2π,length=73) # per 5 degrees with length 73
+
+# other settings
+n = 1
+EorH = 2  # 1: Eelectric, 2: Magnemtic
+Isff = zeros(length(theta),length(phi))
+for i = 1:length(theta)
+    for j = 1:length(phi)
+        Isff[i,j] = jlmie_Isff_EM_n(nmat,radius,lbd0[EorH],nenv,theta[i],phi[j],n)[EorH]
+    end
+end
+Isff = Isff/maximum(Isff)
+mat_x = Isff.* (sin.(theta) * cos.(phi)')
+mat_y = Isff.* (sin.(theta) * sin.(phi)')
+mat_z = Isff.* (cos.(theta) * cos.(phi*0)')
+
+# Just a sphere
+# mat_x = 1 .* (sin.(theta) * cos.(phi)')
+# mat_y = 1 .* (sin.(theta) * sin.(phi)')
+# mat_z = 1 .* (cos.(theta) * cos.(phi*0)')
+
+# show results
+pyplot()
+plt = plot(mat_x,mat_y,mat_z,
+           st=:surface,
+           camera=(30,30), # azimuth, elevate
+           xlabel="x",
+           ylabel="y",
+           zlabel="z",
+           xlims=(-1,1),
+           ylims=(-1,1),
+           zlims=(-1,1),
+)
+display(plt)
+
+# save
+if true
+    runningfilename = splitext(splitpath(@__FILE__)[end])[1]
+    outfilename = ".\\result\\" * runningfilename
+    savefig(plt,outfilename)
+    println("Graph has been saved in result directory")
+end

example_dev/xx02_3DRadiationPatternEM_n4r100.jl

@@ -0,0 +1,54 @@
+using jlmie
+using Plots
+
+# # xz plot range
+# x = (-150:5:150)
+# y = 0.0
+# z = (-150:5:150)
+# # calculation using broadcast
+# r, theta, phi = cart2sph(x',y,z)
+# # heatmap(x,z,r)
+# # heatmap(x,z,theta)
+# # heatmap(x,z,phi)
+
+# meshgrid type calculation
+# mat_x = x' .* ones(length(z))
+# mat_z = ones(length(x))' .* z
+# mat_r,mat_theta,mat_phi = cart2sph(mat_x,y,mat_z)
+# calculation using broadcast
+# r, theta, phi = cart2sph(x',y,z)
+# heatmap(x,z,r)
+# heatmap(x,z,theta)
+# heatmap(x,z,phi)
+
+x = (-150:1:150)
+y = (-150:1:150)
+z = (-150:1:150)
+
+xyz = 2  # 1:x, 2:y, 3:z
+if xyz == 1
+    # xz plot range
+    mat_x = x' .* ones(length(z))
+    mat_z = ones(length(x))' .* z
+    mat_r,mat_theta,mat_phi = cart2sph(mat_x,0,mat_z)
+    p1 = heatmap(x,z,mat_r)
+    p2 = heatmap(x,z,mat_theta)
+    p3 = heatmap(x,z,mat_phi)
+elseif xyz == 2
+    # yz plot range
+    mat_y = y' .* ones(length(y))
+    mat_z = ones(length(y))' .* z
+    mat_r,mat_theta,mat_phi = cart2sph(0,mat_y,mat_z)
+    p1 = heatmap(y,z,mat_r)
+    p2 = heatmap(y,z,mat_theta)
+    p3 = heatmap(y,z,mat_phi)
+elseif xyz == 3
+    # xy plot range
+    mat_x = x' .* ones(length(y))
+    mat_y = ones(length(x))' .* y
+    mat_r,mat_theta,mat_phi = cart2sph(mat_x,mat_y,0)
+    p1 = heatmap(x,y,mat_r)
+    p2 = heatmap(x,y,mat_theta)
+    p3 = heatmap(x,y,mat_phi)
+end
+plot(p1,p2,p3)