== for Hcat (#64)

* == for Hcat

* D^2 * b

Author: dlfivefifty

Project.toml

@@ -1,7 +1,7 @@
 name = "QuasiArrays"
 uuid = "c4ea9172-b204-11e9-377d-29865faadc5c"
 authors = ["Sheehan Olver <solver@mac.com>"]
-version = "0.5.0"
+version = "0.5.1"
 
 [deps]
 ArrayLayouts = "4c555306-a7a7-4459-81d9-ec55ddd5c99a"

src/lazyquasiarrays.jl

@@ -224,7 +224,7 @@ axes(A::Applied{<:Any,typeof(^),<:Tuple{<:AbstractQuasiMatrix,<:Number}}) = axes
 eltype(::Applied{<:Any,typeof(^),<:Tuple{<:AbstractQuasiMatrix{T},<:Integer}}) where T = T
 eltype(::Applied{<:Any,typeof(^),<:Tuple{<:AbstractQuasiMatrix{T},<:Number}}) where T = complex(T)
 
-function *(App::ApplyQuasiMatrix{<:Any,typeof(^),<:Tuple{<:AbstractQuasiMatrix{T},<:Integer}}, b::AbstractQuasiMatrix) where T
+function *(App::ApplyQuasiMatrix{<:Any,typeof(^),<:Tuple{<:AbstractQuasiMatrix{T},<:Integer}}, b::AbstractQuasiArray) where T
     A,p = arguments(App)
     p < 0 && return ApplyQuasiMatrix(^,inv(A),-p)*b
     p == 0 && return copy(b)

src/quasiconcat.jl

@@ -25,6 +25,11 @@ end
 
 \(A::AbstractQuasiArray, H::ApplyQuasiMatrix{<:Any,typeof(hcat)}) = hcat((Ref(A) .\ H.args)...)
 
+function ==(A::ApplyQuasiMatrix{<:Any,typeof(hcat)}, B::ApplyQuasiMatrix{<:Any,typeof(hcat)})
+    axes.(A.args,2) == axes.(B.args,2) && return all(A.args .== B.args)
+    Base.invoke(==, Tuple{AbstractQuasiArray,AbstractQuasiArray}, A, B)
+end
+
 """
     UnionVcat
 

test/test_quasiconcat.jl

@@ -8,6 +8,9 @@ import QuasiArrays: ApplyQuasiMatrix, UnionVcat
         @test B isa ApplyQuasiMatrix{Float64,typeof(hcat)}
         @test axes(B) == (axes(A,1), Base.OneTo(4))
         @test B[0.0,1] == B[0.0,3] == A[0.0,1]
+        @test B == B
+        @test B ≠ A
+        @test B ≠ [A A A]
         @test_throws BoundsError B[0.1,1]
         @test_throws BoundsError B[0.0,5]
         @test QuasiArray(B) == B