Fix Type in Aowdt (#137)

* Fix Type in Aowdt

* To v0.19

* To v0.18.2

Author: JeffFessler

Project.toml

@@ -1,7 +1,7 @@
 name = "MIRT"
 uuid = "7035ae7a-3787-11e9-139a-5545ed3dc201"
 authors = ["fessler <fessler@umich.edu>"]
-version = "0.18.1"
+version = "0.18.2"
 
 [deps]
 AVSfldIO = "b6189060-daf9-4c28-845a-cc0984b81781"

src/regularize/Aodwt.jl

@@ -33,7 +33,7 @@ which is useful when imposing scale-dependent regularization
 """
 function Aodwt(
     dims::Dims ;
-    T::Type{<:Complex{<:AbstractFloat}} = ComplexF32,
+    T::Type{<:Union{AbstractFloat,Complex{<:AbstractFloat}}} = Float32,
     level::Int = 3,
     wt = wavelet(WT.haar),
     operator::Bool = true, # !
@@ -52,7 +52,7 @@ function Aodwt(
             mfun = (A, x) -> A * x
             return mfun, LinearMapAA(forw!, back!,
                 (prod(dims), prod(dims)) ;
-                prop = (wt=wt, level=lev), T=T,
+                prop = (wt, level=lev), T,
                 operator = true, idim=dims, odim=dims,
             )
         else
@@ -63,7 +63,7 @@ function Aodwt(
             #    x -> vec(dwt(reshape(x, dims), wt, level)),
             #    y -> vec(idwt(reshape(y, dims), wt, level)),
                 (prod(dims), prod(dims)) ;
-                prop = (wt=wt, level=lev), T=T,
+                prop = (wt, level=lev), T,
             )
         end
     end
@@ -72,7 +72,7 @@ function Aodwt(
 
     scales = zeros(dims)
     for il in 1:level
-        _,Al = mfunA(il)
+        _, Al = mfunA(il)
         tmp = mfun(Al, ones(dims)) .== 0
         scales += il * (tmp .& (scales .== 0))
     end