Skip to content
Snippets Groups Projects
Commit 5ae4b635 authored by Stephan Hilb's avatar Stephan Hilb
Browse files

remove unused method

parent df151355
No related branches found
No related tags found
No related merge requests found
...@@ -268,33 +268,6 @@ function project_qi_lagrange!(u::FeFunction, img) ...@@ -268,33 +268,6 @@ function project_qi_lagrange!(u::FeFunction, img)
return u return u
end end
# FIXME: unfinished / unused -> remove?
function quadrature_cell_pixels(mesh, cell; m0)
d = ndims_domain(mesh)
A = SArray{Tuple{ndims_space(mesh), nvertices_cell(mesh)}}(
view(mesh.vertices, :, view(mesh.cells, :, cell)))
# TODO: cleanup when minimum/maximum with dims on StaticArrays becomes
# type-stable
I0 = round.(Int, SVector(minimum(A[1, :]), minimum(A[2, :])) .- m0) .+ 1
I1 = round.(Int, SVector(maximum(A[1, :]), maximum(A[2, :])) .- m0) .+ 1
#quadrature() = SA[1/6, 1/6, 1/6], SA[1/6 4/6 1/6; 1/6 1/6 4/6]
vertices = Float64[]
for I in CartesianIndex(I0[1], I0[2]):CartesianIndex(I1[1], I1[2])
x = SVector(Tuple(I) .- 1 .+ m0)
# TODO: move to global-to-local function or inside-triangle test
xloc = (A[:, SUnitRange(2, 3)] .- A[:, 1])::SArray \
(x .- A[:, 1])
ε = eps()
if all(xloc .>= -ε) && sum(xloc) <= 1 + ε
append!(vertices, xloc)
end
end
qw = ones(length(vertices) ÷ d) ./ length(vertices)
qx = reshape!(vertices, (d, :))
end
struct PixelIterator{A, G} struct PixelIterator{A, G}
"cell vertices" "cell vertices"
A::A A::A
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment