Update for QuasiArrays.jl (#13)

* add ultraspherical bases

* fixes

* add Conversion

* Evaluation and collocation work now

* Update runtests.jl

* Update Project.toml

* Move out QuasiArrays

* Updating for QuasiArrays v0.0.1

* Tests pass again!

* Update runtests.jl

* Update Project.toml

* more tests pass

* quasimulapplystyle

* Remove AdjointStyle, BasisStyle

* work on materialize

* Ldiv improvements

Author: dlfivefifty

.gitignore

@@ -2,3 +2,4 @@
 *.jl.*.cov
 *.jl.mem
 deps/deps.jl
+.DS_Store

Project.toml

@@ -0,0 +1,27 @@
+name = "ContinuumArrays"
+uuid = "7ae1f121-cc2c-504b-ac30-9b923412ae5c"
+
+[deps]
+BandedMatrices = "aae01518-5342-5314-be14-df237901396f"
+ContinuumArrays = "7ae1f121-cc2c-504b-ac30-9b923412ae5c"
+DomainSets = "5b8099bc-c8ec-5219-889f-1d9e522a28bf"
+FillArrays = "1a297f60-69ca-5386-bcde-b61e274b549b"
+InfiniteArrays = "4858937d-0d70-526a-a4dd-2d5cb5dd786c"
+InfiniteLinearAlgebra = "cde9dba0-b1de-11e9-2c62-0bab9446c55c"
+IntervalSets = "8197267c-284f-5f27-9208-e0e47529a953"
+LazyArrays = "5078a376-72f3-5289-bfd5-ec5146d43c02"
+LinearAlgebra = "37e2e46d-f89d-539d-b4ee-838fcccc9c8e"
+QuasiArrays = "c4ea9172-b204-11e9-377d-29865faadc5c"
+
+[extras]
+Test = "8dfed614-e22c-5e08-85e1-65c5234f0b40"
+
+[compat]
+BandedMatrices = "0.10"
+DomainSets = "0.1"
+FillArrays = "0.6"
+InfiniteArrays = "0.1.1"
+IntervalSets = "0.3.1"
+LazyArrays = "0.10"
+julia = "1.0"
+QuasiArrays = "0.0.2"

REQUIRE

@@ -1,5 +1,13 @@
-using ContinuumArrays, FillArrays, InfiniteArrays
-import ContinuumArrays.QuasiArrays: Inclusion, QuasiDiagonal
+using ContinuumArrays, QuasiArrays, FillArrays, InfiniteArrays
+import QuasiArrays: Inclusion, QuasiDiagonal
+
+using Plots; pyplot();
+
+S = LinearSpline(range(0,1; length=10))
+xx = range(0,1; length=1000)
+
+S = Jacobi(true,true)
+
 
 P = Legendre()
 X = QuasiDiagonal(Inclusion(-1..1))

examples/collocation.jl

@@ -1,19 +1,19 @@
 module ContinuumArrays
-using IntervalSets, LinearAlgebra, LazyArrays, BandedMatrices, InfiniteArrays, DomainSets
+using IntervalSets, LinearAlgebra, LazyArrays, FillArrays, BandedMatrices, InfiniteArrays, DomainSets, InfiniteLinearAlgebra, QuasiArrays
 import Base: @_inline_meta, axes, getindex, convert, prod, *, /, \, +, -,
-                IndexStyle, IndexLinear, ==, OneTo, tail
+                IndexStyle, IndexLinear, ==, OneTo, tail, similar, copyto!
 import Base.Broadcast: materialize
-import LazyArrays: Mul2, MemoryLayout, Applied, ApplyStyle, flatten, _flatten
+import LazyArrays: MemoryLayout, Applied, ApplyStyle, flatten, _flatten, colsupport, adjointlayout, LdivApplyStyle
 import LinearAlgebra: pinv
 import BandedMatrices: AbstractBandedLayout, _BandedMatrix
+import FillArrays: AbstractFill, getindex_value
 
-include("QuasiArrays/QuasiArrays.jl")
-using .QuasiArrays
-import .QuasiArrays: cardinality, checkindex, QuasiAdjoint, QuasiTranspose, slice, Inclusion, SubQuasiArray,
+import QuasiArrays: cardinality, checkindex, QuasiAdjoint, QuasiTranspose, Inclusion, SubQuasiArray,
                     QuasiDiagonal, MulQuasiArray, MulQuasiMatrix, MulQuasiVector, QuasiMatMulMat,
-                    ApplyQuasiArray, ApplyQuasiMatrix, LazyQuasiArrayApplyStyle
+                    ApplyQuasiArray, ApplyQuasiMatrix, LazyQuasiArrayApplyStyle, AbstractQuasiArrayApplyStyle,
+                    LazyQuasiArray, LazyQuasiVector, LazyQuasiMatrix, quasimulapplystyle, LazyLayout
 
-export Spline, LinearSpline, HeavisideSpline, DiracDelta, Derivative, JacobiWeight, Jacobi, Legendre,
+export Spline, LinearSpline, HeavisideSpline, DiracDelta, Derivative, JacobiWeight, Jacobi, Legendre, Chebyshev, Ultraspherical,
             fullmaterialize
 
 ####
@@ -24,8 +24,8 @@ struct AlephInfinity{N} <: Integer end
 const ℵ₁ = AlephInfinity{1}()
 
 
-const QMul2{A,B} = Mul{<:Any, <:Tuple{A,B}}
-const QMul3{A,B,C} = Mul{<:Any, <:Tuple{A,B,C}}
+const QMul2{A,B} = Mul{<:AbstractQuasiArrayApplyStyle, <:Tuple{A,B}}
+const QMul3{A,B,C} = Mul{<:AbstractQuasiArrayApplyStyle, <:Tuple{A,B,C}}
 
 cardinality(::AbstractInterval) = ℵ₁
 *(ℵ::AlephInfinity) = ℵ