From 0eca7246392982246c9e8d959d0d4a5c82159001 Mon Sep 17 00:00:00 2001 From: Jesse Chan Date: Fri, 19 Apr 2024 08:06:48 -0500 Subject: [PATCH] clean up subcell ops --- src/RefElemData_SBP.jl | 5 +---- 1 file changed, 1 insertion(+), 4 deletions(-) diff --git a/src/RefElemData_SBP.jl b/src/RefElemData_SBP.jl index f8ef4620..8fe47bc1 100644 --- a/src/RefElemData_SBP.jl +++ b/src/RefElemData_SBP.jl @@ -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)