Skip to content

Conversation

Shreyas-Ekanathan
Copy link

builds on the previous PR to implement a full butterfly factorization method

@ChrisRackauckas
Copy link
Member

@Shreyas-Ekanathan can you push your version?

Comment on lines +73 to +92
#=
if (N % 2 == 1) # N odd
n = N2
for m in 1:M
A[m, n] = u[m] * A[m, n] * v[n]
end
end

if (M % 2 == 1) # M odd
m = M2
for n in 1:N
A[m, n] = u[m] * A[m, n] * v[n]
end
end

if (M % 2 == 1) && (N % 2 == 1)
m = M2
n = N2
A[m, n] /= (u[m] * v[n])
end =#

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

delete this

Comment on lines +42 to +45
Mh = M >>> 1
Nh = N >>> 1
M2 = M - Mh
N2 = N - Nh

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

now that we're back to padding, we know that all these are equal

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

3 participants