function rinv(
R::AbstractVector,
δ::DifferentialMobilityAnalyzer;
λ₁ = 1e-2,
λ₂ = 1e1,
n = 1,
)
setupRegularization(δ.𝐀, δ.𝐈, R, inv(δ.𝐒) * R, n) # setup the system
λopt = lcorner(λ₁, λ₂; n = 10, r = 3) # compute the optimal λ
N = clean(Ninv(λopt)) # find the inverted size
return SizeDistribution([], δ.De, δ.Dp, δ.ΔlnD, N ./ δ.ΔlnD, N, :regularized)
end
It seems that there is no S value returned. But in [Inversion of Size Distribution From Data.ipynb], df2 = DataFrame(Dp = 𝕟ⁱⁿᵛ².Dp, S = 𝕟ⁱⁿᵛ².S, Dist = ["𝕟ⁱⁿᵛ²" for i = 1:length(𝕟ⁱⁿᵛ².Dp)]), includes 𝕟ⁱⁿᵛ².S