@@ -9,7 +9,7 @@ import BlockBandedMatrices: isblockbanded, _BlockBandedMatrix
9
9
import MatrixFactorizations: QLPackedQ
10
10
import BandedMatrices: bandeddata, _BandedMatrix, BandedStyle
11
11
import LazyArrays: colsupport, ApplyStyle, MemoryLayout, ApplyLayout, LazyArrayStyle, arguments
12
- import InfiniteArrays: OneToInf
12
+ import InfiniteArrays: OneToInf, oneto
13
13
import LazyBandedMatrices: BroadcastBandedBlockBandedLayout, BroadcastBandedLayout
14
14
15
15
40
40
41
41
@testset " ∞-block arrays" begin
42
42
@testset " fixed block size" begin
43
- k = Base. OneTo .(Base . oneto (∞))
44
- n = Fill .(Base . oneto (∞),Base . oneto (∞))
43
+ k = Base. OneTo .(oneto (∞))
44
+ n = Fill .(oneto (∞),oneto (∞))
45
45
@test broadcast (length,k) ≡ map (length,k) ≡ OneToInf ()
46
46
@test broadcast (length,n) ≡ map (length,n) ≡ OneToInf ()
47
47
b = mortar (Fill ([1 ,2 ],∞))
52
52
end
53
53
54
54
@testset " 1:∞ blocks" begin
55
- a = blockedrange (Base . oneto (∞))
55
+ a = blockedrange (oneto (∞))
56
56
@test axes (a,1 ) == a
57
57
o = Ones ((a,))
58
58
@test Base. BroadcastStyle (typeof (a)) isa LazyArrayStyle{1 }
102
102
103
103
@testset " triangle recurrences" begin
104
104
@testset " n and k" begin
105
- n = mortar (Fill .(Base . oneto (∞),Base . oneto (∞)))
106
- k = mortar (Base. OneTo .(Base . oneto (∞)))
105
+ n = mortar (Fill .(oneto (∞),oneto (∞)))
106
+ k = mortar (Base. OneTo .(oneto (∞)))
107
107
108
108
@test n[Block (5 )] ≡ layout_getindex (n, Block (5 )) ≡ view (n,Block (5 )) ≡ Fill (5 ,5 )
109
109
@test k[Block (5 )] ≡ layout_getindex (k, Block (5 )) ≡ view (k,Block (5 )) ≡ Base. OneTo (5 )
135
135
end
136
136
137
137
@testset " BlockHcat copyto!" begin
138
- n = mortar (Fill .(Base . oneto (∞),Base . oneto (∞)))
139
- k = mortar (Base. OneTo .(Base . oneto (∞)))
138
+ n = mortar (Fill .(oneto (∞),oneto (∞)))
139
+ k = mortar (Base. OneTo .(oneto (∞)))
140
140
141
141
a = b = c = 0.0
142
142
dat = BlockHcat (
160
160
161
161
@testset " BlockBanded" begin
162
162
a = b = c = 0.0
163
- n = mortar (Fill .(Base . oneto (∞),Base . oneto (∞)))
164
- k = mortar (Base. OneTo .(Base . oneto (∞)))
163
+ n = mortar (Fill .(oneto (∞),oneto (∞)))
164
+ k = mortar (Base. OneTo .(oneto (∞)))
165
165
Dy = BlockBandedMatrices. _BandedBlockBandedMatrix ((k .+ (b+ c))' , axes (k,1 ), (- 1 ,1 ), (- 1 ,1 ))
166
166
N = 100 ;
167
167
@test Dy[Block .(1 : N), Block .(1 : N)] == BlockBandedMatrices. _BandedBlockBandedMatrix ((k .+ (b+ c))[Block .(1 : N)]' , axes (k,1 )[Block .(1 : N)], (- 1 ,1 ), (- 1 ,1 ))
0 commit comments