make it work for all combinations of precisions

Author: ranjanan

src/classical.jl

@@ -82,7 +82,7 @@ function rs_direct_interpolation_pass1(T, A, splitting)
                                                 splitting::Vector{Ti},
                                                 Bp::Vector{Ti}) where {Tv,Ti}
 
-
+                                                
     Bx = zeros(Tv, Bp[end] - 1)
     Bj = zeros(Ti, Bp[end] - 1)
 

src/multilevel.jl

@@ -66,8 +66,10 @@ function solve{T}(ml::MultiLevel, b::Vector{T},
                                     tol = 1e-5;
                                     verbose = false,
                                     log = false)
-    x = zeros(T, size(b))
-    residuals = Vector{T}()
+    V = promote_type(eltype(ml.levels[1].A), eltype(b))
+    x = zeros(V, size(b))
+    tol = eltype(b)(tol)
+    residuals = Vector{V}()
     A = length(ml) == 1 ? ml.final_A : ml.levels[1].A
     normb = norm(b)
     if normb != 0
@@ -91,14 +93,14 @@ function solve{T}(ml::MultiLevel, b::Vector{T},
         return x
     end
 end
-function __solve{T}(v::V, ml, x::Vector{T}, b::Vector{T}, lvl)
+function __solve(v::V, ml, x, b, lvl)
 
     A = ml.levels[lvl].A
     presmoother!(ml.presmoother, A, x, b)
 
     res = b - A * x
     coarse_b = ml.levels[lvl].R * res
-    coarse_x = zeros(T, size(coarse_b))
+    coarse_x = zeros(eltype(coarse_b), size(coarse_b))
 
     if lvl == length(ml.levels)
         coarse_x = coarse_solver(ml.coarse_solver, ml.final_A, coarse_b)
@@ -113,5 +115,4 @@ function __solve{T}(v::V, ml, x::Vector{T}, b::Vector{T}, lvl)
     x
 end
 
-coarse_solver{Tv,Ti}(::Pinv, A::SparseMatrixCSC{Tv,Ti}, b::Vector{Tv}) =
-                                                        pinv(full(A)) * b
+coarse_solver(::Pinv, A, b) = pinv(full(A)) * b

src/smoother.jl

@@ -32,7 +32,7 @@ smoother!(s::GaussSeidel, ::BackwardSweep, A, x, b) =
     gs!(A, b, x, size(A,1), -1, 1)
 
 
-function gs!{T,Ti}(A::SparseMatrixCSC{T,Ti}, b::Vector{T}, x::Vector{T}, start, step, stop)
+function gs!(A, b, x, start, step, stop)
     n = size(A, 1)
     z = zero(eltype(A))
     for i = start:step:stop

test/runtests.jl

@@ -51,6 +51,7 @@ end
 # Direct Interpolation
 using AMG
 A = poisson(5)
+A = Float64.(A)
 splitting = [1,0,1,0,1]
 P, R = AMG.direct_interpolation(A, copy(A), splitting)
 @test P ==  [ 1.0  0.0  0.0