diff --git a/src/image.jl b/src/image.jl index 20ca3732a7620b2e53984b739e391d275ca3814d..0011d6ea5b3f97ff220032ac1b9f74ad19ea303c 100644 --- a/src/image.jl +++ b/src/image.jl @@ -268,33 +268,6 @@ function project_qi_lagrange!(u::FeFunction, img) return u 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} "cell vertices" A::A