paul=Cauchy, remove last printlns

Author: dsweber2

src/ContinuousWavelets.jl

@@ -42,7 +42,7 @@ are:
 
     ``\\psi\\hat(\\omega) \\propto \\textrm{e}^{-\\frac{\\mu^2}{2}}\\big(\\textrm{e}^{-(\\mu - \\omega)^2} -\\textrm{e}^{\\frac{\\omega^2-\\mu^2}{2}}\\big)``
 
-- `Paul{N}`: A complex analytic wavelet. `pauln` for n in `1:20` e.g. `paul5`
+- `Paul{N}`: A complex analytic wavelet, also known as Cauchy wavelets. `pauln` for n in `1:20` e.g. `paul5`
 
     ``\\psi\\hat(\\omega) \\propto \\chi_{\\omega \\geq 0} \\omega^N\\textrm{e}^{-\\omega}``
 - `Dog{N}`: Derivative of a Gaussian, where N is the number of

src/apply.jl

@@ -298,7 +298,6 @@ function icwt(res::AbstractArray, cWav::CWT, inverseStyle::PenroseDelta)
     Ŵ = computeWavelets(size(res,1), cWav)[1]
     β = computeDualWeights(Ŵ, cWav)
     testDualCoverage(β, Ŵ)
-    println(size(res),size(β))
     compXRecon = sum(res .* β, dims=2)
     imagXRecon = irfft(im*rfft(imag.(compXRecon),1), size(compXRecon,1)) # turns out the dual frame for the imaginary part is rather gross in the time domain
     return imagXRecon + real.(compXRecon)

src/utils.jl

@@ -276,7 +276,6 @@ compute the weight given to each wavelet so that in the Fourier domain, the sum
 """
 function computeDualWeights(Ŵ, wav)
     @views lastReasonableFreq = computeLastFreq(Ŵ[:,end], wav)
-    println(size(Ŵ))
     Wdag = pinv(Ŵ[1:lastReasonableFreq,:])
     β = conj.(Wdag * ones(size(Wdag,2)))'
     return β