application: topaz

The TOPology Application Zoo deals with abstract simplicial complexes. A complex is given as a list of facets. You can ask for its global properties (manifold recognition, homology groups, etc.), explore the local vertex environment (stars, links, etc.), and make a lot of constructions.

The visualization means are constrained, as they are mostly based on the GRAPH (1-skeleton) of a complex.

imports from: common, graph
uses: ideal

Objects

User Functions

  •  

    These functions compare two SimplicialComplex

    •  
      find_facet_vertex_permutations (complex1, complex2) → Pair<Array<Int>, Array<int>>

      Find the permutations of facets and vertices which maps the first complex to the second one. The facet permutation is the first component of the return value. If the complexes are not isomorphic, an exception is thrown.

    •  
      isomorphic (complex1, complex2) → Bool

      Determine whether two given complexes are combinatorially isomorphic. The problem is reduced to graph isomorphism of the vertex-facet incidence graphs.

      Parameters
      SimplicialComplexcomplex1
      SimplicialComplexcomplex2
      Returns
      Bool
    •  
      pl_homeomorphic (complex1, complex2) → Bool

      Tries to determine whether two complexes are pl-homeomorphic by using bistellar flips and a simulated annealing strategy.

      You may specify the maximal number of rounds, how often the system may relax before heating up and how much heat should be applied. The function stops computing, once the size of the triangulation has not decreased for rounds iterations. If the abs flag is set, the function stops after rounds iterations regardless of when the last improvement took place. Additionally, you may set the threshold min_n_facets for the number of facets when the simplification ought to stop. Default is d+2 in the CLOSED_PSEUDO_MANIFOLD case and 1 otherwise.

      If you want to influence the distribution of the dimension of the moves when warming up you may do so by specifying a distribution. The number of values in distribution determines the dimensions used for heating up. The heating and relaxing parameters decrease dynamically unless the constant flag is set. The function prohibits to execute the reversed move of a move directly after the move itself unless the allow_rev_move flag is set. Setting the allow_rev_move flag might help solve a particular resilient problem.

      If you are interested in how the process is coming along, try the verbose option. It specifies after how many rounds the current best result is displayed.

      The obj determines the objective function used for the optimization. If obj is set to 0, the function searches for the triangulation with the lexicographically smallest f-vector, if obj is set to 1, the function searches for the triangulation with the reversed-lexicographically smallest f-vector and if obj is set to 2 the sum of the f-vector entries is used. The default is 1.

      Parameters
      SimplicialComplexcomplex1
      SimplicialComplexcomplex2
      Options
      Introunds
      Boolabs
      Intobj
      Intrelax
      Intheat
      Boolconstant
      Boolallow_rev_move
      Intmin_n_facets
      Intverbose
      Intseed
      Boolquiet
      Array<Int>distribution
      Returns
      Bool
  •  

    Special purpose functions.

    •  
      is_generalized_shelling (FaceList) → Bool

      Check if a given sequence of faces of a simplicial complex is a generalized shelling.

      Parameters
      Array<Set>FaceList
      Options
      Boolverbose
      Returns
      Bool
    •  
      is_vertex_decomposition (complex, vertices) → Bool

      Check whether a given ordered subset of the vertex set is a vertex decomposition. Works for 1-, 2- and 3-manifolds only!

      Parameters
      SimplicialComplexcomplex
      Array<Int>vertices
      shedding vertices
      Options
      Boolverbose
      Returns
      Bool
    •  
      mixed_graph (complex)

      Produces the mixed graph of a complex.

      Parameters
      SimplicialComplexcomplex
      Options
      Floatedge_weight
    •  
      morse_matching (complex) → EdgeMap

      Compute a Morse matching. Two heuristics are implemented:

      (1) A simple greedy algorithm:
      The arcs are visited in lexicographical order, i.e.:
      we proceed by levels from top to bottom,
      visit the faces in each dimension in lexicographical order,
      and visited the faces covered by these faces in lexicographical order.
      This heuristic is used by default and with heuristic => 1.
      (2) A Morse matching can be improved by canceling critical cells
      along unique alternating paths, see function
      processAlternatingPaths() in file morse_matching_tools.h .
      This idea is due to Robin Forman:
      Morse Theory for Cell-Complexes,
      Advances in Math., 134 (1998), pp. 90-145.
      This heuristic is used by default and with heuristic => 2.

      The default setting is to use both, i.e., to run the greedy algorithm and then improve the result by the canceling algorithm.

      Morse matchings for the bottom level can be found optimally by spanning tree techniques. This can be enabled by the option levels => 1. If the complex is a pseudo-manifold the same can be done for the top level (option levels => 2). By specifying option levels => 0, both levels can be computed by spanning trees. For 2-dim pseudo-manifolds this computes an optimal Morse matching.

      Parameters
      SimplicialComplexcomplex
      given by its Hasse diagram
      Options
      Intheuristic
      (1=greedy, 2=cancel, 0=both (default))
      Intlevels
      (1=bottom, 2=top, 0=both (default))
      Returns
      EdgeMap
      matching a labelling of the edges of the Hasse diagram with integer values, where 1 means that the edge is in the matching
    •  
      morse_matching_size (complex) → Int

      Compute the number of edges in a Morse matching.

      Parameters
      SimplicialComplexcomplex
      a complex with a Morse matching
      Returns
      Int
      the number of edges in the matching.
    •  
      stanley_reisner (complex) → ideal::Ideal

      Creates the Stanley-Reisner ideal of a simplicial complex. Optional ring parameter is required to have precisely as many variables as C has vertices.

      Parameters
      SimplicialComplexcomplex
      Options
      Ring<Rational,int>ring
      Returns
      ideal::Ideal
    •  
      stiefel_whitney (facets) → Array<PowerSet<Int>>

      Computes Stiefel-Whitney classes of mod 2 Euler space (in particular, closed manifold). Use option verbose to show regular pairs and cycles. A narrower dimension range of interest can be specified. Negative values are treated as co-dimension - 1

      Parameters
      Array<Set<Int>>facets
      the facets of the simplicial complex
      Options
      Inthigh_dim
      Intlow_dim
      Boolverbose
      Returns
      Array<PowerSet<Int>>
  •  

    These functions construct a new SimplicialComplex from other objects of the same type.

    •  
      alexander_dual (complex) → SimplicialComplex

      Computes the Alexander dual complex, that is, the complements of all non-faces. The vertex labels are preserved unless the nol flag is specified.

      Parameters
      SimplicialComplexcomplex
      Options
      Boolnol
      Returns
      SimplicialComplex
    •  
      barycentric_subdivision (complex) → SimplicialComplex

      Computes the barycentric subdivision of complex.

      Parameters
      SimplicialComplexcomplex
      Options
      Stringpin_hasse_section
      default: HASSE_DIAGRAM
      Stringlabel_section
      default: VERTEX_LABELS
      Stringcoord_section
      default: VERTICES
      Boolgeometric_realization
      set to 1 to obtain a GeometricSimplicialComplex, default: 0
      Returns
      SimplicialComplex
    •  
      bistellar_simplification (complex) → SimplicialComplex

      Heuristic for simplifying the triangulation of the given manifold without changing its PL-type. The function uses bistellar flips and a simulated annealing strategy.

      You may specify the maximal number of rounds, how often the system may relax before heating up and how much heat should be applied. The function stops computing, once the size of the triangulation has not decreased for rounds iterations. If the abs flag is set, the function stops after rounds iterations regardless of when the last improvement took place. Additionally, you may set the threshold min_n_facets for the number of facets when the simplification ought to stop. Default is d+2 in the CLOSED_PSEUDO_MANIFOLD case and 1 otherwise.

      If you want to influence the distribution of the dimension of the moves when warming up you may do so by specifying a distribution. The number of values in distribution determines the dimensions used for heating up. The heating and relaxing parameters decrease dynamically unless the constant flag is set. The function prohibits to execute the reversed move of a move directly after the move itself unless the allow_rev_move flag is set. Setting the allow_rev_move flag might help solve a particular resilient problem.

      If you are interested in how the process is coming along, try the verbose option. It specifies after how many rounds the current best result is displayed.

      The obj determines the objective function used for the optimization. If obj is set to 0, the function searches for the triangulation with the lexicographically smallest f-vector, if obj is set to any other value the sum of the f-vector entries is used. The default is 1.

      Parameters
      SimplicialComplexcomplex
      Options
      Introunds
      Boolabs
      Intobj
      Intrelax
      Intheat
      Boolconstant
      Boolallow_rev_move
      Intmin_n_facets
      Intverbose
      Intseed
      Boolquiet
      Array<Int>distribution
      Returns
      SimplicialComplex
    •  
      bs2quotient (P, complex) → SimplicialComplex

      Create a simplicial complex from a simplicial subdivision of a given complex by identifying vertices on the boundary of the original complex according to a group that acts on vertices.

      Contained in extension bundled:group.
      Parameters
      PolytopeP
      the underlying polytope
      SimplicialComplexcomplex
      a sufficiently fine subdivision of P, for example the second barycentric subdivision
      Returns
      SimplicialComplex
    •  
      cone (complex, k) → SimplicialComplex

      Produce the k-cone over a given simplicial complex.

      Parameters
      SimplicialComplexcomplex
      intk
      default is 1
      Options
      Array<String>apex_labels
      labels of the apex vertices. Default labels have the form apex_0, apex_1, .... In the case the input complex has already vertex labels of this kind, the duplicates are avoided.
      Boolnol
      don't generate any vertex labels.
      Returns
      SimplicialComplex
    •  
      connected_sum (complex1, complex2, f1, f2) → SimplicialComplex

      Compute the connected sum of two complexes.

      Parameters f_1 and f_2// specify which facet of the first and second complex correspondingly are glued together. Default is the 0-th facet of both.

      The vertices in the selected facets are identified with each other according to their order in the facet (that is, in icreasing index order). The option permutation allows to get an alternative identification. It should specify a permutation of the vertices of the second facet.

      The vertices of the new complex get the original labels with _1 or _2 appended, according to the input complex they came from. If you set the no_labels flag, the label generation will be omitted.

      Parameters
      SimplicialComplexcomplex1
      SimplicialComplexcomplex2
      intf1
      default: 0
      intf2
      default: 0
      Options
      Array<int>permutation
      Boolno_lables
      Returns
      SimplicialComplex
    •  
      deletion (complex, face) → SimplicialComplex

      Remove the given face and all the faces containing it.

      Parameters
      SimplicialComplexcomplex
      Set<Int>face
      specified by vertex indices. Please use labeled_vertices if you want to specify the face by vertex labels.
      Options
      Boolno_labels
      do not write vertex labels.
      Returns
      SimplicialComplex
    •  
      disjoint_union (complex1, complex2) → SimplicialComplex

      Produce the disjoint union of the two given complexes.

      Parameters
      SimplicialComplexcomplex1
      SimplicialComplexcomplex2
      Options
      labelscreates
      VERTEX_LABELS. The vertex labels are built from the original labels with a suffix _1 or _2 appended.
      Returns
      SimplicialComplex
    •  
      edge_contraction (complex) → SimplicialComplex

      Heuristic for simplifying the triangulation of the given manifold without changing its PL-type. Choosing a random order of the vertices, the function tries to contract all incident edges.

      Parameters
      SimplicialComplexcomplex
      Options
      Intseed
      Returns
      SimplicialComplex
    •  
      foldable_prism (complex) → GeometricSimplicialComplex

      Produce a prism over a given SimplicialComplex.

      Parameters
      GeometricSimplicialComplexcomplex
      Options
      Boolgeometric_realization
      Returns
      GeometricSimplicialComplex
    •  
      h_induced_quotient (C, vertices) → SimplicialComplex

      Let C be the given simplicial and A the subcomplex induced by the given vertices. Then this function produces a simplicial complex homotopy equivalent to C mod A by adding the cone over A with apex a to C. The label of the apex my be specified via the option apex.

      Parameters
      SimplicialComplexC
      Set<Int>vertices
      Options
      Boolno_labels
      tells the client not to write any labels.
      Stringapex
      Returns
      SimplicialComplex
    •  
      induced_subcomplex (complex, vertices) → SimplicialComplex

      Produce the subcomplex consisting of all faces which are contained in the given set of vertices.

      Parameters
      SimplicialComplexcomplex
      Set<Int>vertices
      Options
      Boolno_label
      tells the client not to create any labels.
      Boolgeom_real
      tells the client to inherit the COORDINATES.
      Returns
      SimplicialComplex
    •  
      iterated_barycentric_subdivision (complex, k) → SimplicialComplex

      Computes the k-th barycentric subdivision of complex by iteratively calling barycentric_subdivision.

      Parameters
      SimplicialComplexcomplex
      Intk
      Options
      Stringpin_hasse_section
      default: HASSE_DIAGRAM
      Stringlabel_section
      default: VERTEX_LABELS
      Stringcoord_section
      default: VERTICES
      Boolgeometric_realization
      set to 1 to obtain a GeometricSimplicialComplex, default: 0
      Returns
      SimplicialComplex
    •  
      join_complexes (complex1, complex2) → SimplicialComplex

      Creates the join of complex1 and complex2.

      Parameters
      SimplicialComplexcomplex1
      SimplicialComplexcomplex2
      Options
      Boollabels
      creates VERTEX_LABELS. The vertex labels are built from the original labels with a suffix _1 or _2 appended.
      Returns
      SimplicialComplex
    •  
      k_skeleton (complex, k) → SimplicialComplex

      Produce the k-skeleton.

      Parameters
      SimplicialComplexcomplex
      intk
      Options
      Boolvertex_labels
      whether to create VERTEX_LABELS
      Returns
      SimplicialComplex
    •  
      k_skeleton (complex, k) → GeometricSimplicialComplex

      Produce the k-skeleton.

      Parameters
      GeometricSimplicialComplexcomplex
      intk
      Options
      Boolvertex_labels
      whether to create VERTEX_LABELS
      Returns
      GeometricSimplicialComplex
    • link_complex (complex, face) → SimplicialComplex
    •  
      p_sum_triangulation (P, Q, WebOfStars) → GeometricSimplicialComplex

      Produce a specific P-sum-triangulation of two given triangulations and a WebOfStars.

      Parameters
      GeometricSimplicialComplexP
      First complex which will be favoured
      GeometricSimplicialComplexQ
      second complex
      IncidenceMatrixWebOfStars
      Every row corresponds to a full dimensional simplex in P and every column to a full dimensional simplex in Q.
      Returns
      GeometricSimplicialComplex
    •  
      simplicial_product (complex1, complex2) → SimplicialComplex

      Computes the simplicial product of two complexes. Vertex orderings may be given as options.

      Parameters
      SimplicialComplexcomplex1
      SimplicialComplexcomplex2
      Options
      Array<Int>vertex_order1
      Array<Int>vertex_order2
      Boolgeometric_realization
      default 0
      Boolcolor_cons
      Boolno_labels
      Returns
      SimplicialComplex
    •  
      simplicial_product <Scalar> (complex1, complex2) → GeometricSimplicialComplex<Scalar>

      Computes the simplicial product of two complexes. Vertex orderings may be given as options.

      Type Parameters
      Scalar
      Parameters
      GeometricSimplicialComplexcomplex1
      GeometricSimplicialComplexcomplex2
      Options
      Array<Int>vertex_order1
      Array<Int>vertex_order2
      Boolgeometric_realization
      default 1
      Boolcolor_cons
      Boolno_labels
      Returns
      GeometricSimplicialComplex<Scalar>
    •  
      star (complex, face) → SimplicialComplex

      Produce the star of the face of the complex.

      Parameters
      SimplicialComplexcomplex
      Set<int>face
      Options
      Boollabels
      creates VERTEX_LABELS.
      Returns
      SimplicialComplex
    •  
      star_deletion (complex, face) → SimplicialComplex

      Remove the star of a given face.

      Parameters
      SimplicialComplexcomplex
      Set<Int>face
      specified by vertex indices. Please use labeled_vertices if you want to specify the face by vertex labels.
      Options
      Boolno_labels
      do not write vertex labels.
      Returns
      SimplicialComplex
    •  
      stellar_subdivision (complex, faces) → SimplicialComplex

      Computes the complex obtained by stellar subdivision of the given faces of the complex.

      Parameters
      SimplicialComplexcomplex
      Array<Set<Int>>faces
      Options
      Boolno_labels
      Boolgeometric_realization
      default 0
      Returns
      SimplicialComplex
    •  
      stellar_subdivision (complex, face) → SimplicialComplex

      Computes the complex obtained by stellar subdivision of the given face of the complex.

      Parameters
      SimplicialComplexcomplex
      Set<Int>face
      Options
      Boolno_labels
      Boolgeometric_realization
      default 0
      Returns
      SimplicialComplex
    •  
      suspension (complex, k) → SimplicialComplex

      Produce the k-suspension over a given simplicial complex.

      Parameters
      SimplicialComplexcomplex
      Intk
      default value is 1
      Options
      Array<String>labels
      for the apices. By default apices are labeled with apex_0+, apex_0-, apex_1+, etc. If one of the specified labels already exists, a unique one is made by appending _i where i is some small number.
      Boolnol
      do not produce any labels.
      Returns
      SimplicialComplex
    •  
      union (complex1, complex2) → SimplicialComplex

      Produce the union of the two given complexes, identifying vertices with equal labels.

      Parameters
      SimplicialComplexcomplex1
      SimplicialComplexcomplex2
      Options
      Boollabels
      creates VERTEX_LABELS.
      Returns
      SimplicialComplex
  •  

    These functions construct a new SimplicialComplex from other objects.

  •  

    With these clients you can create special examples of simplicial complexes and complexes belonging to parameterized families.

    •  
      ball (d) → GeometricSimplicialComplex

      A d-dimensional ball, realized as the d-simplex.

      Parameters
      Intd
      dimension
      Returns
      GeometricSimplicialComplex
    •  
      complex_projective_plane () → SimplicialComplex

      The complex projective plane with the vertex-minimal triangulation by Kühnel and Brehm.

    •  
      cube_complex ()

      Produces a triangulated pile of hypercubes: Each cube is split into d! tetrahedra, and the tetrahedra are all grouped around one of the diagonal axes of the cube. DOC_FIXME args: x_1, ... , x_d

    •  
      klein_bottle () → SimplicialComplex

      The Klein bottle.

    •  
      rand_knot (n_edges) → SimplicialComplex

      Produce a random knot (or link) as a polygonal closed curve in 3-space. The knot (or each connected component of the link) has n_edges edges.

      The vertices are uniformly distributed in [-1,1]3, unless the on_sphere option is set. In the latter case the vertices are uniformly distributed on the 3-sphere. Alternatively the brownian option produces a knot by connecting the ends of a simulated brownian motion.

      Parameters
      Intn_edges
      Options
      Intn_comp
      number of components, default is 1.
      Boolon_sphere
      Boolbrownian
      Intseed
      Returns
      SimplicialComplex
    •  
      real_projective_plane () → SimplicialComplex

      The real projective plane with its unique minimal triangulation on six vertices.

    •  
      simplex (d) → SimplicialComplex

      A simplex of dimension d.

      Parameters
      Intd
      dimension
      Returns
      SimplicialComplex
    •  
      sphere (d) → GeometricSimplicialComplex

      The d-dimensional sphere, realized as the boundary of the (d+1)-simplex.

      Parameters
      Intd
      dimension
      Returns
      GeometricSimplicialComplex
    •  
      surface (g) → SimplicialComplex

      Produce a surface of genus g. For g >= 0 the client produces an orientable surface, otherwise it produces a non-orientable one.

      Parameters
      Intg
      genus
      Returns
      SimplicialComplex
    •  
      torus () → SimplicialComplex

      The Császár Torus. Geometric realization by Frank Lutz, Electronic Geometry Model No. 2001.02.069

    •  
      unknot (m, n) → GeometricSimplicialComplex

      Produces a triangulated 3-sphere with the particularly NASTY embedding of the unknot in its 1-skeleton. The parameters m >= 2 and n >= 1 determine how entangled the unknot is. eps determines the COORDINATES.

      Parameters
      Intm
      Intn
      Options
      Rationaleps
      Returns
      GeometricSimplicialComplex
  •  

    The following functions compute topological invariants.

    •  
      cap_product (complex)

      Compute and print all cap products of cohomology and homology cycles.

      Parameters
      SimplicialComplexcomplex
    •  
      homology (complex, co)

      Calculate the (co-)homology groups of a simplicial complex.

      Parameters
      Array<Set<int>>complex
      Boolco
      set to true for cohomology
      Options
      Intdim_low
      narrows the dimension range of interest, with negative values being treated as co-dimensions
      Intdim_high
      see dim_low
    •  
      homology_and_cycles (complex, co)

      Calculate the (co-)homology groups and cycle representatives of a simplicial complex.

      Parameters
      Array<Set<int>>complex
      Boolco
      set to true for cohomology
      Options
      Intdim_low
      narrows the dimension range of interest, with negative values being treated as co-dimensions
      Intdim_high
      see dim_low

Property Types

  •  

    The following property_types are topological invariants.

    •  
      CycleGroup <Scalar>

      A group is encoded as a pair of an integer matrix and a vector of faces. The elements of the group can be obtained by symbolic multiplication of both.

      Access methods: coeff delivers the integer matrix, faces the vector of faces.

      Type Parameters
      Scalar
      integer type of matrix elements
    •  
      HomologyGroup

      A group is encoded as a sequence ( { (t1 m1) ... (tn mn) } f) of non-negative integers, with t1 > t2 > ... > tn > 1, plus an extra non-negative integer f.

      The group is isomorphic to (Z/t1)m1 × ... × (Z/tn)mn × Zf, where Z0 is the trivial group.

      Access methods: torsion delivers the list of Z-groups, betti_number the number f.

    •  
      IntersectionForm

      TODO: explain the meaning of the fields: parity, positive, negative