Update rescale option in InelasticNeutronScatteringSpectra.

Quantum-Many-Body · Dec 4, 2024 · f375167 · f375167
1 parent d1ca9f6
commit f375167
Show file tree

Hide file tree

Showing 3 changed files with 3 additions and 3 deletions.
diff --git a/docs/src/examples/Phonon.md b/docs/src/examples/Phonon.md
@@ -53,7 +53,7 @@ The inelastic neutron scattering spectra of phonons can also be computed:
 # scale: the scale of the intensity
 spectra = phonon(
     :INSS,
-    InelasticNeutronScatteringSpectra(path, range(0.0, 2.5, length=501); fwhm=0.05, rescale=log)
+    InelasticNeutronScatteringSpectra(path, range(0.0, 2.5, length=501); fwhm=0.05, rescale=x->log(1+x))
 )
 plt = plot()
 plot!(plt, spectra)

diff --git a/src/Core.jl b/src/Core.jl
@@ -946,7 +946,7 @@ function run!(tba::Algorithm{<:CompositeTBA{Phononic, <:AbstractLattice}}, inss:
             end
         end
     end
-    inss.data[3][:, :] = get(inss.action.options, :rescale, identity).(inss.data[3].+1)
+    inss.data[3][:, :] = get(inss.action.options, :rescale, identity).(inss.data[3])
 end
 @inline function checkpolarizations(qs₁::AbstractMatrix, qs₂::AbstractMatrix, momentum)
     inner = mapreduce((e₁, e₂)->norm(conj(e₁)*e₂), +, qs₁, qs₂)/norm(qs₁)/norm(qs₂)

diff --git a/test/TightBindingApproximation.jl b/test/TightBindingApproximation.jl
@@ -274,7 +274,7 @@ end
     energybands = phonon(:EB, EnergyBands(path))
     savefig(plot(energybands), "phonon.png")
 
-    inelastic = phonon(:INSS, InelasticNeutronScatteringSpectra(path, range(0.0, 2.5, length=501); fwhm=0.05, rescale=log))
+    inelastic = phonon(:INSS, InelasticNeutronScatteringSpectra(path, range(0.0, 2.5, length=501); fwhm=0.05, rescale=x->log(1+x)))
     plt = plot()
     plot!(plt, inelastic)
     plot!(plt, energybands, color=:white, linestyle=:dash)