Examples

We'll use a few of the functions used for testing the package to demonstrate its usage. These are

• simplices_sharing_vertices. Generate a set of simplices which intersect in some arbitrary way, but sharing at least one vertex.
• nontrivially_intersecting_simplices. Generate a set of non-trivially intersecting simplices (i.e. intersections are not only along boundaries or vertices).
• childsimplex. Generate a simplex completely contained within a parent simplex.

Note that these functions take as inputs simplices of shape (dim + 1, dim). This will be fixed in a future release.

In the examples, some of the functions used to generate simplices return the simplex arrays with rows as columns. Therefore, we transpose the simplices before calling simplexintersection.

Simplices sharing vertices

Let's compute the intersection between a set of simplices that share at least one vertex.

using Simplices
s₁, s₂ = transpose.(simplices_sharing_vertices(3))

simplexintersection(s₁, s₂)

0.014958296370133663

Nontrivially intersecting simplices

Simplices can also intersect in nontrivial ways, meaning that they have an intersection beyond a common boundary or vertex.

using Simplices
s₁, s₂ = transpose.(nontrivially_intersecting_simplices(3))
simplexintersection(s₁, s₂)

0.0017481765842190034

One simplex fully contained within the other

We'll generate a random simplex s₁, then generate a simplex s₂ fully contained within that simplex. If s₂ is fully contained, the intersection volume should be the volume of s₂.

using Simplices

Ds =  2:7;
intersection_vols = zeros(Float64, length(Ds));
analytical_vols   = zeros(Float64, length(Ds));
for i = 1:length(Ds)
    s₁ = rand(Ds[i] + 1, Ds[i])
    s₂ = childsimplex(s₁)
    intersection_vols[i] = simplexintersection(transpose(s₁), transpose(s₂))
    analytical_vols[i] = volume(s₂)
end

# Within numerical error, the results should be the same.
all([isapprox(intersection_vols[i], analytical_vols[i]; atol = 1e-9)
    for i = 1:length(Ds)])

true

Simplices are identical

If simplices are identical, the intersection volume should equal the volume of either simplex:

using Simplices

s₁ = rand(4, 3); s₂ = s₁;

simplexintersection(transpose(s₁), transpose(s₂)) .≈ volume(s₁) .≈ volume(s₂)

true