testing Jacobi FEM

Author: dlfivefifty

appveyor.yml

@@ -0,0 +1,43 @@
+environment:
+  matrix:
+  - julia_version: 0.7
+  - julia_version: 1
+  - julia_version: nightly
+
+platform:
+  - x86 # 32-bit
+  - x64 # 64-bit
+
+# # Uncomment the following lines to allow failures on nightly julia
+# # (tests will run but not make your overall status red)
+# matrix:
+#   allow_failures:
+#   - julia_version: nightly
+
+branches:
+  only:
+    - master
+    - /release-.*/
+
+notifications:
+  - provider: Email
+    on_build_success: false
+    on_build_failure: false
+    on_build_status_changed: false
+
+install:
+  - ps: iex ((new-object net.webclient).DownloadString("https://raw.githubusercontent.com/JuliaCI/Appveyor.jl/version-1/bin/install.ps1"))
+
+build_script:
+  - echo "%JL_BUILD_SCRIPT%"
+  - C:\julia\bin\julia -e "%JL_BUILD_SCRIPT%"
+
+test_script:
+  - echo "%JL_TEST_SCRIPT%"
+  - C:\julia\bin\julia -e "%JL_TEST_SCRIPT%"
+
+# # Uncomment to support code coverage upload. Should only be enabled for packages
+# # which would have coverage gaps without running on Windows
+# on_success:
+#   - echo "%JL_CODECOV_SCRIPT%"
+#   - C:\julia\bin\julia -e "%JL_CODECOV_SCRIPT%"

test/runtests.jl

@@ -157,14 +157,6 @@ end
 S = Jacobi(true,true)
 W = Diagonal(JacobiWeight(true,true))
 D = Derivative(axes(W,1))
-D*S
-
-D*S[:,1:10]
-
-materialize(D*W*S)
-
-J = Jacobi(false,false)
-J'J
 
 
 S'W'W*S
@@ -176,21 +168,59 @@ L = D*W*S[:,1:N]
 
 Δ = L'L # weak second derivative
 
-f = Legendre() * Vcat(1:10, Zeros(∞))
+f = Legendre() * Vcat(1:2, Zeros(∞))
+L'f
 
+A = (L').factors[2]*f.factors[1]
 
-L'f
+A*f.factors[2]
 
 axes(L')
 
 Legendre()'f
 
 A = Diagonal(1:∞)
-x = Vcat(1:10, Zeros(∞))
+
 
 @which A*x
+1
+using InfiniteArrays, BandedMatrices, InteractiveUtils, Test
+
+A = BandedMatrices._BandedMatrix((1:∞)', ∞, -1,1)
+x = Vcat(1:3, Zeros(∞))
+
+ Zeros(∞) ==  Zeros(∞)
+
+
+@test A*x == Vcat([4.0,9.0], Zeros(∞))
+
+y = similar(Mul
+
+y = Vcat(randn(9), Zeros(∞))
+copyto!(y , Mul(A,x))
+
+MemoryLayout(x)
+
+similar(Mul(A,x), Float64)
+
+
+typeof( Mul(B,x))
+
+Mul2{<:BandedColumnMajor,<:Any,<:AbstractMatrix,
+        <:Vcat{<:AbstractVector,<:Zeros}}
+
+
+
+
+
+typeof(V)
+V = view(bandeddata(B),1:1,3:∞)
+    _BandedMatrix(V, ∞, 2,-2)
+
+using InteractiveUtils
+
+B
 
-B = BandedMatrices._BandedMatrix((1:∞)', ∞, -1,1)
 
 B*x