We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
There was an error while loading. Please reload this page.
1 parent a657582 commit 3686d11Copy full SHA for 3686d11
src/butterflylu.jl
@@ -26,7 +26,7 @@ struct 🦋workspace{T}
26
out::Vector{T}
27
n::Int
28
function 🦋workspace(A, b, ::Val{SEED} = Val(888)) where {SEED}
29
- N = size(A, 1)
+ len = size(A, 1)
30
out = similar(b, N)
31
if (M % 4 != 0)
32
A = pad!(A)
@@ -36,20 +36,20 @@ struct 🦋workspace{T}
36
U, V = (similar(A), similar(A))
37
ws = 🦋generate_random!(A)
38
materializeUV(U, V, ws)
39
- new{eltype(A)}(A, b, ws, U, V, out, N)
+ new{eltype(A)}(A, b, ws, U, V, out, n)
40
end
41
42
43
function 🦋solve!(workspace::🦋workspace, thread)
44
- (;A, b, ws, U, V, out, N) = workspace
+ (;A, b, ws, U, V, out, n) = workspace
45
🦋mul!(A, ws)
46
F = RecursiveFactorization.lu!(A, Val(false), thread)
47
48
mul!(b, U', b)
49
ldiv!(b, UnitLowerTriangular(F.factors), b, thread)
50
ldiv!(b, UpperTriangular(F.factors), b, thread)
51
mul!(b, V, b)
52
- out .= @view b[1:N]
+ out .= @view b[1:n]
53
out
54
55
0 commit comments