Skip to content

Commit

Permalink
clean up subcell ops
Browse files Browse the repository at this point in the history
  • Loading branch information
jlchan committed Apr 19, 2024
1 parent afa5d52 commit 0eca724
Showing 1 changed file with 1 addition and 4 deletions.
5 changes: 1 addition & 4 deletions src/RefElemData_SBP.jl
Original file line number Diff line number Diff line change
Expand Up @@ -419,14 +419,11 @@ function subcell_limiting_operators(Qr::AbstractMatrix; tol = 100 * eps())
end
end

d = (ones(size(D, 1))' * D)'
d = vec(ones(size(D, 1))' * D)
ids = findall(@. abs(d) > 1e2 * eps())

Z = nullspace(D)[ids, :]
A = pinv(D)[ids,:]
X = randn(size(Z, 2), size(D, 1))
y = randn(length(ids))
e = ones(size(D, 1))

# compute R via linear algebra
m, n = size(A)
Expand Down

0 comments on commit 0eca724

Please sign in to comment.