Computing the Discrete Compactness of Orthogonal Pseudo-Polytopes via Their 𝑛D-EVM Representation

Algorithm 5

Computing L_C(p) for an $n$D-OPP expressed in the EVM.

Input:   An $n$D-EVM p and the number n of dimensions.

Output: The number of internal contacts, perpendicular to X_i-axis,

         i=1,2,…,n, between the hypervoxels that originally composed p.

Procedure TotalInternalContacts(EVM p, int n)

      int Lc = 0

       for each sorting in $\{$X₁X₂$\dots$X_n-1X_n,X₂X₃$\dots$X_nX₁,X₃X₄$\dots$X₁X₂,$\dots$,X_nX₁$\dots$X_n-2X_n-1$\}$  do

         SortEVM(p, n, sorting)

         Lc = Lc + InternalContacts(p, n) // Call to Algorithm 4.

       end-of-for

       return Lc

end-of-procedure