From 7e2364208b59c1ab04fd2a7622169d21eef4e776 Mon Sep 17 00:00:00 2001 From: Sheehan Olver Date: Sun, 21 Jan 2024 09:52:13 +0000 Subject: [PATCH] Use transform plans instead of factorizations --- docs/src/index.md | 3 --- src/bases/bases.jl | 24 ------------------------ 2 files changed, 27 deletions(-) diff --git a/docs/src/index.md b/docs/src/index.md index b2c1f47..ea66139 100644 --- a/docs/src/index.md +++ b/docs/src/index.md @@ -52,9 +52,6 @@ ContinuumArrays.MulPlan ContinuumArrays.PiecewiseBasis ``` ```@docs -ContinuumArrays.MappedFactorization -``` -```@docs ContinuumArrays.basis ``` ```@docs diff --git a/src/bases/bases.jl b/src/bases/bases.jl index 6b1c04d..f0514df 100644 --- a/src/bases/bases.jl +++ b/src/bases/bases.jl @@ -248,30 +248,6 @@ _sub_factorize(::Tuple{Any,Any}, (kr,jr), L, dims...; kws...) = _factorize(::SubBasisLayout, L, dims...; kws...) = _sub_factorize(size(parent(L)), parentindices(L), L, dims...; kws...) -""" - MappedFactorization(F, map) - -remaps a factorization to a different domain. That is, if `M` is a `MappedFactorization` -then `M \\ f` is equivalent to `F \\ f[map]` -""" -struct MappedFactorization{T, FAC<:Factorization{T}, MAP} <: Factorization{T} - F::FAC - map::MAP -end - -\(a::MappedFactorization, b::AbstractQuasiVector) = a.F \ view(b, a.map) -\(a::MappedFactorization, b::AbstractVector) = a.F \ b -\(a::MappedFactorization, b::AbstractQuasiMatrix) = a.F \ view(b, a.map, :) - - -function invmap end - -function _factorize(::MappedBasisLayout, L, dims...; kws...) - kr, jr = parentindices(L) - P = parent(L) - MappedFactorization(factorize(view(P,:,jr), dims...; kws...), invmap(parentindices(L)[1])) -end - plan_ldiv(A, B::AbstractQuasiVector) = factorize(A) plan_ldiv(A, B::AbstractQuasiMatrix) = factorize(A, size(B,2))