Update finitedifference_2d.jl

Author: dlfivefifty

examples/finitedifference_2d.jl

@@ -3,7 +3,7 @@
 #
 ###########
 
-using BlockBandedMatrices, BandedMatrices, ArrayLayouts, FillArrays
+using BlockBandedMatrices, BandedMatrices, ArrayLayouts, FillArrays, LazyBandedMatrices, LazyArrays
 
 
 function finitedifference_2d(n)
@@ -58,14 +58,14 @@ n = 1000;
     @time _gaussseidel(L,U, b, x, y, 20) # 0.4s
 
     @time L̃ = sparse(BandedBlockBandedMatrix(L))
-    @time Ũ = sparse(U)
+    @time Ũ = sparse(U)
     x = copy(b)
     y = copy(x)
-    @time _gaussseidel2(L̃,Ũ, b, x, y, 20)
+    @time _gaussseidel2(L̃,Ũ, b, x, y, 20)
 
-Ã = sparse(A)
-@time qr(Ã);
-@time à \b ;
+Ã = sparse(A)
+@time qr(Ã);
+@time à \b ;
 
 A = randn(n,n)
 
@@ -82,98 +82,3 @@ n = 1000;
         lu(A)
     end
 
-0.02*n
-
-A*x - b
-
-
-
-x̃ = copy(x)
-ỹ = copy(y)
-
-@view(y[1:end-1]) .= Mul(U , @view(x[2:end]))
-
-mul!(@view(ỹ[1:end-1]) , Ũ , @view(x̃[2:end]))
-y[end] = 0
-ỹ[end] = 0
-x = copy(b)
-    y = copy(x)
-    _gaussseidel(L,U, b, x, y, 1) - x_bb
-
-x .= b .- y
-ỹ .= b .- ỹ
-x .= Ldiv(L, x)
-
-x̃ .= L̃ \ ỹ
-
-
-x - x_bb
-
-x̃-x
-
-
-
-(L̃ \ x) - (L \ x)
-
-ỹ - y |> norm
-x - x_bb
-
-A*x - b
-
-
-L̃ == L
-Ũ == U
-
-norm(A*x - b)
-@time
-
-
-@time copy(x)
-L̃ = sparse(L)
-Ũ = sparse(U)
-
-
-
-
-y = x[1:end-1]
-z = similar(y)
-@time z .= Mul(U,y)
-y = similar(x)
-@time (@view(y[1:end-1]) .= Mul(U , @view(x[2:end])))
-y[end] = 0
-@time x .= b .- y
-@time x .= Ldiv(L, x)
-A = randn(9,40_000);
-    @time A*x;
-
-A
-S = sparse(A);
-
-400^2
-
-@time L\x
-
-
-h = 1/n
-    @time D² = BandedMatrix(0 => Fill(-2,n), 1 => Fill(1,n-1), -1 => Fill(1,n-1))/h^2
-    @time D_xx = BandedBlockBandedMatrix(Kron(D², Eye(n)))
-    @time D_yy = BandedBlockBandedMatrix(Kron(Eye(n), D²))
-    @time D_xx .+ D_yy
-    4
-
-
-
-using Profile
-dest = Δ;
-    @time BlockBandedMatrices.blockbanded_copyto!(dest, D_yy)
-    view(dest, Block(1,1)) .= view(dest, Block(1,1)) .+ view(D_yy, Block(1,1))
-D_yy
-@time copyto!(dest, Broadcast.broadcasted(+, D_xx , D_yy))
-
-
-
-
-using FFTW
-A = randn(400, 400);
-    @time fft(A);
-dest[Block(1,1)]