From 5ae4b63595fbf2ebe2349a6c86a66702a1e37274 Mon Sep 17 00:00:00 2001
From: Stephan Hilb <stephan@ecshi.net>
Date: Sun, 12 Dec 2021 14:49:35 +0100
Subject: [PATCH] remove unused method
---
src/image.jl | 27 ---------------------------
1 file changed, 27 deletions(-)
diff --git a/src/image.jl b/src/image.jl
index 20ca373..0011d6e 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
--
GitLab