Regina Calculation Engine
|
A utility class for searching through all possible gluing permutation sets that correspond to a given triangle edge pairing. More...
#include <census/dim2gluingpermsearcher.h>
Public Types | |
typedef DimTraits< dim > ::FacetPairing | FacetPairing |
typedef DimTraits< dim >::Perm | Perm |
typedef DimTraits< dim >::Simplex | Simplex |
typedef DimTraits< dim > ::Triangulation | Triangulation |
Public Member Functions | |
Dim2GluingPermSearcher (const Dim2EdgePairing *pairing, const Dim2EdgePairing::IsoList *autos, bool orientableOnly, UseDim2GluingPerms use, void *useArgs=0) | |
Initialises a new search for gluing permutation sets. More... | |
Dim2GluingPermSearcher (std::istream &in, UseDim2GluingPerms use, void *useArgs=0) | |
Initialises a new search manager based on data read from the given input stream. More... | |
virtual | ~Dim2GluingPermSearcher () |
Destroys this search manager and all supporting data structures. More... | |
virtual void | runSearch (long maxDepth=-1) |
Generates all possible gluing permutation sets that satisfy the current search criteria. More... | |
bool | completePermSet () const |
Determines whether this search manager holds a complete gluing permutation set or just a partially completed search state. More... | |
void | dumpTaggedData (std::ostream &out) const |
Dumps all internal data in a plain text format, along with a marker to signify which precise class the data belongs to. More... | |
virtual void | dumpData (std::ostream &out) const |
Dumps all internal data in a plain text format to the given output stream. More... | |
bool | inputError () const |
Was an error found during construction from an input stream? More... | |
unsigned | size () const |
Returns the total number of simplices under consideration. More... | |
const FacetPairing * | getFacetPairing () const |
Returns the specific pairing of simplex facets that this set of gluing permutations complements. More... | |
Perm | gluingPerm (const NFacetSpec< dim > &source) const |
Returns the gluing permutation associated with the given simplex facet. More... | |
Perm | gluingPerm (unsigned simp, unsigned facet) const |
Returns the gluing permutation associated with the given simplex facet. More... | |
Triangulation * | triangulate () const |
Returns a newly created triangulation as modelled by this set of gluing permutations and the associated simplex facet pairing. More... | |
Static Public Member Functions | |
static void | findAllPerms (const Dim2EdgePairing *pairing, const Dim2EdgePairing::IsoList *autos, bool orientableOnly, UseDim2GluingPerms use, void *useArgs=0) |
The main entry routine for running a search for all gluing permutation sets that complement a given edge pairing. More... | |
static Dim2GluingPermSearcher * | bestSearcher (const Dim2EdgePairing *pairing, const Dim2EdgePairing::IsoList *autos, bool orientableOnly, UseDim2GluingPerms use, void *useArgs=0) |
Constructs a search manager of the best possible class for the given search parameters. More... | |
static Dim2GluingPermSearcher * | readTaggedData (std::istream &in, UseDim2GluingPerms use, void *useArgs=0) |
Creates a new search manager based on tagged data read from the given input stream. More... | |
Static Public Attributes | |
static const char | dataTag_ |
A character used to identify this class when reading and writing tagged data in text format. More... | |
Protected Member Functions | |
bool | isCanonical () const |
Compares the current set of gluing permutations with its preimage under each automorphism of the underlying edge pairing, in order to see whether the current set is in canonical form (i.e., is lexicographically smallest). More... | |
virtual char | dataTag () const |
Returns the character used to identify this class when storing tagged data in text format. More... | |
int & | permIndex (const NFacetSpec< dim > &source) |
Returns the index into array Perm::Sn_1 describing how the the given facet is joined to its partner. More... | |
int & | permIndex (unsigned simp, unsigned facet) |
Returns the index into array Perm::Sn_1 describing how the the given facet is joined to its partner. More... | |
const int & | permIndex (const NFacetSpec< dim > &source) const |
Returns the index into array Perm::Sn_1 describing how the the given facet is joined to its partner. More... | |
const int & | permIndex (unsigned simp, unsigned facet) const |
Returns the index into array Perm::Sn_1 describing how the the given facet is joined to its partner. More... | |
int | gluingToIndex (const NFacetSpec< dim > &source, const Perm &gluing) const |
Returns the index into array Perm::Sn_1 corresponding to the given gluing permutation from the given facet to its partner. More... | |
int | gluingToIndex (unsigned simp, unsigned facet, const Perm &gluing) const |
Returns the index into array Perm::Sn_1 corresponding to the given gluing permutation from the given facet to its partner. More... | |
Perm | indexToGluing (const NFacetSpec< dim > &source, int index) const |
Returns the gluing permutation from the given facet to its partner that corresponds to the given index into array Perm::Sn_1. More... | |
Perm | indexToGluing (unsigned simp, unsigned facet, int index) const |
Returns the gluing permutation from the given facet to its partner that corresponds to the given index into array Perm::Sn_1. More... | |
Protected Attributes | |
const Dim2EdgePairing::IsoList * | autos_ |
The set of isomorphisms that define equivalence of gluing permutation sets. More... | |
bool | autosNew |
Did we create the isomorphism list autos_ ourselves (in which case we must destroy it also)? More... | |
bool | orientableOnly_ |
Are we only searching for gluing permutations that correspond to orientable triangulations? More... | |
UseDim2GluingPerms | use_ |
A routine to call each time a gluing permutation set is found during the search. More... | |
void * | useArgs_ |
Additional user-supplied data to be passed as the second argument to the use_ routine. More... | |
bool | started |
Has the search started yet? This helps distinguish between a new search and the resumption of a partially completed search. More... | |
int * | orientation |
Keeps track of the orientation of each triangle in the underlying triangulation. More... | |
Dim2TriangleEdge * | order |
Describes the order in which gluing permutations are assigned to edges. More... | |
int | orderSize |
The total number of edges in the edge pairing graph, i.e., the number of elements of interest in the order[] array. More... | |
int | orderElt |
Marks which element of order[] we are currently examining at this stage of the search. More... | |
const FacetPairing * | pairing_ |
The facet pairing that this permutation set complements. More... | |
int * | permIndices_ |
The index into array Perm::Sn_1 describing how each simplex facet is glued to its partner. More... | |
bool | inputError_ |
Has an error occurred during construction from an input stream? More... | |
A utility class for searching through all possible gluing permutation sets that correspond to a given triangle edge pairing.
In the future, there may be subclasses of Dim2GluingPermSearcher that correspond to specialised search algorithms for use in certain scenarios. The main class Dim2GluingPermSearcher offers a default search algorithm that may be used in a general context.
The simplest way of performing a search through all possible gluing permutations is by calling the static method findAllPerms(). This will examine the search parameters and ensure that the best possible algorithm is used. For finer control over the program flow, the static method bestSearcher() can be used to create a search manager of the most suitable class and then runSearch() can be called on this object directly. For absolute control, a specific algorithm can be forced by explicitly constructing an object of the corresponding class (and again calling runSearch() on that object directly).
Note that this class derives from Dim2GluingPerms. The search will involve building and repeatedly modifying the inherited Dim2GluingPerms data in-place.
regina::Dim2GluingPermSearcher::Dim2GluingPermSearcher | ( | const Dim2EdgePairing * | pairing, |
const Dim2EdgePairing::IsoList * | autos, | ||
bool | orientableOnly, | ||
UseDim2GluingPerms | use, | ||
void * | useArgs = 0 |
||
) |
Initialises a new search for gluing permutation sets.
The search is started by calling runSearch(). Note that the static method findAllPerms() handles both construction and searching, and is the preferred entry point for end users.
The arguments to this constructor describe the search parameters in detail, as well as what should be done with each gluing permutation set that is found.
pairing | the specific pairing of triangle edges that the generated permutation sets will complement. |
autos | the collection of isomorphisms that define equivalence of permutation sets. These are used by runSearch(), which produces each permutation set precisely once up to equivalence. These isomorphisms must all be automorphisms of the given edge pairing, and will generally be the set of all such automorphisms. This parameter may be 0, in which case the set of all automorphisms of the given edge pairing will be generated and used. |
orientableOnly | true if only gluing permutations corresponding to orientable triangulations should be generated, or false if no such restriction should be imposed. |
use | the function to call upon each permutation set that is found. The first parameter passed to this function will be a gluing permutation set. The second parameter will be parameter useArgs as was passed to this routine. |
useArgs | the pointer to pass as the final parameter for the function use which will be called upon each permutation set found. |
regina::Dim2GluingPermSearcher::Dim2GluingPermSearcher | ( | std::istream & | in, |
UseDim2GluingPerms | use, | ||
void * | useArgs = 0 |
||
) |
Initialises a new search manager based on data read from the given input stream.
This may be a new search or a partially completed search.
This routine reads data in the format written by dumpData(). If you wish to read data whose precise class is unknown, consider using dumpTaggedData() and readTaggedData() instead.
If the data found in the input stream is invalid or incorrectly formatted, the routine inputError() will return true
but the contents of this object will be otherwise undefined.
in | the input stream from which to read. |
use | as for the main Dim2GluingPermSearcher constructor. |
useArgs | as for the main Dim2GluingPermSearcher constructor. |
|
virtual |
Destroys this search manager and all supporting data structures.
|
static |
Constructs a search manager of the best possible class for the given search parameters.
Different subclasses of Dim2GluingPermSearcher provide optimised search algorithms for different types of search.
Calling this routine and then calling runSearch() on the result has the same effect as the all-in-one routine findAllPerms(). Unless you have specialised requirements (such as partial searching), you are probably better calling findAllPerms() instead.
The resulting object is newly created, and must be destroyed by the caller of this routine.
See the Dim2GluingPermSearcher constructor for documentation on the arguments to this routine.
|
inline |
Determines whether this search manager holds a complete gluing permutation set or just a partially completed search state.
This may assist the use_ routine when running partial depth-based searches. See runSearch() for further details.
true
if a complete gluing permutation set is held, or false
otherwise.
|
inlineprotectedvirtual |
Returns the character used to identify this class when storing tagged data in text format.
|
virtual |
Dumps all internal data in a plain text format to the given output stream.
This object can be recreated from this text data by calling the input stream constructor for this class.
This routine may be useful for transferring objects from one processor to another.
Note that subclass data is written after superclass data, so it is safe to dump data from a subclass and then recreate a new superclass object from that data (though subclass-specific information will of course be lost).
out | the output stream to which the data should be written. |
Reimplemented from regina::NGenericGluingPerms< 2 >.
void regina::Dim2GluingPermSearcher::dumpTaggedData | ( | std::ostream & | out | ) | const |
Dumps all internal data in a plain text format, along with a marker to signify which precise class the data belongs to.
This routine can be used with readTaggedData() to transport objects from place to place whose precise class is unknown.
out | the output stream to which the data should be written. |
|
static |
The main entry routine for running a search for all gluing permutation sets that complement a given edge pairing.
This routine examines the search parameters, chooses the best possible search algorithm, constructs an object of the corresponding subclass of Dim2GluingPermSearcher and then calls runSearch().
See the Dim2GluingPermSearcher constructor for documentation on the arguments to this routine. See the runSearch() method for documentation on how the search runs and returns its results.
|
inherited |
Returns the specific pairing of simplex facets that this set of gluing permutations complements.
|
inherited |
Returns the gluing permutation associated with the given simplex facet.
source | the simplex facet under investigation. |
|
inherited |
Returns the gluing permutation associated with the given simplex facet.
simp | the simplex under investigation (this must be strictly less than the total number of simplices under consideration). |
facet | the facet of the given simplex under investigation (between 0 and dim inclusive). |
|
protectedinherited |
Returns the index into array Perm::Sn_1 corresponding to the given gluing permutation from the given facet to its partner.
This need not be the index into Perm::Sn_1 that is currently stored for the given facet.
Indices into array Perm::Sn_1 are stored internally in the array permIndices_. Full gluing permutations on the other hand are used in constructing triangulations.
source | the simplex facet under investigation. |
gluing | a possible gluing permutation from the given simplex facet to its partner according to the underlying facet pairing. |
|
protectedinherited |
Returns the index into array Perm::Sn_1 corresponding to the given gluing permutation from the given facet to its partner.
This need not be the index into Perm::Sn_1 that is currently stored for the given facet.
Indices into array Perm::Sn_1 are stored internally in the array permIndices_. Full gluing permutations on the other hand are used in constructing triangulations.
simp | the simplex under investigation; this must be strictly less than the total number of simplices under consideration. |
facet | the facet of the given simplex under investigation; this must be between 0 and dim inclusive. |
gluing | a possible gluing permutation from the given simplex facet to its partner according to the underlying facet pairing. |
|
protectedinherited |
Returns the gluing permutation from the given facet to its partner that corresponds to the given index into array Perm::Sn_1.
This index into Perm::Sn_1 need not be the index that is currently stored for the given facet.
Indices into array Perm::Sn_1 are stored internally in the array permIndices_. Full gluing permutations on the other hand are used in constructing triangulations.
If the given simplex facet and its partner according to the underlying facet pairing are facets x and y of their respective simplices, then the resulting gluing permutation will map x to y.
source | the simplex facet under investigation. |
index | an index into Perm::Sn_1; this must be between 0 and dim!-1 inclusive. |
|
protectedinherited |
Returns the gluing permutation from the given facet to its partner that corresponds to the given index into array Perm::Sn_1.
This index into Perm::Sn_1 need not be the index that is currently stored for the given facet.
Indices into array Perm::Sn_1 are stored internally in the array permIndices_. Full gluing permutations on the other hand are used in constructing triangulations.
If the given simplex facet and its partner according to the underlying facet pairing are facets x and y of their respective simplices, then the resulting gluing permutation will map x to y.
simp | the simplex under investigation; this must be strictly less than the total number of simplices under consideration. |
facet | the facet of the given simplex under investigation; this must be between 0 and dim inclusive. |
index | an index into Perm::Sn_1; this must be between 0 and dim!-1 inclusive. |
|
inherited |
Was an error found during construction from an input stream?
This routine returns true
if an input stream constructor was used to create this object but the data in the input stream was invalid or incorrectly formatted.
If a different constructor was called (i.e., no input stream was used), then this routine will always return false
.
true
if an error occurred during construction from an input stream, or false
otherwise.
|
protected |
Compares the current set of gluing permutations with its preimage under each automorphism of the underlying edge pairing, in order to see whether the current set is in canonical form (i.e., is lexicographically smallest).
true
if the current set is in canonical form, or false
otherwise.
|
protectedinherited |
Returns the index into array Perm::Sn_1 describing how the the given facet is joined to its partner.
Note that this permutation is not a gluing permutation as such, but rather a permutation of 0,...,dim-1 only. For a real facet gluing permutation, see routine gluingPerm().
source | the simplex facet under investigation. |
|
protectedinherited |
Returns the index into array Perm::Sn_1 describing how the the given facet is joined to its partner.
Note that this permutation is not a gluing permutation as such, but rather a permutation of 0,...,dim-1 only. For a real facet gluing permutation, see routine gluingPerm().
simp | the simplex under investigation (this must be strictly less than the total number of simplices under consideration). |
facet | the facet of the given simplex under investigation (between 0 and dim inclusive). |
|
protectedinherited |
Returns the index into array Perm::Sn_1 describing how the the given facet is joined to its partner.
Note that this permutation is not a gluing permutation as such, but rather a permutation of 0,...,dim-1 only. For a real facet gluing permutation, see routine gluingPerm().
source | the simplex facet under investigation. |
|
protectedinherited |
Returns the index into array Perm::Sn_1 describing how the the given facet is joined to its partner.
Note that this permutation is not a gluing permutation as such, but rather a permutation of 0,...,dim-1 only. For a real facet gluing permutation, see routine gluingPerm().
simp | the simplex under investigation (this must be strictly less than the total number of simplices under consideration). |
facet | the facet of the given simplex under investigation (between 0 and dim inclusive). |
|
static |
Creates a new search manager based on tagged data read from the given input stream.
This may be a new search or a partially completed search.
The tagged data should be in the format written by dumpTaggedData(). The precise class of the search manager will be determined from the tagged data, and does not need to be known in advance. This is in contrast to dumpData() and the input stream constructors, where the class of the data being read must be known at compile time.
If the data found in the input stream is invalid or incorrectly formatted, a null pointer will be returned. Otherwise a newly constructed search manager will be returned, and it is the responsibility of the caller of this routine to destroy it after use.
The arguments use and useArgs are the same as for the Dim2GluingPermSearcher constructor.
in | the input stream from which to read. |
|
virtual |
Generates all possible gluing permutation sets that satisfy the current search criteria.
The search criteria are specified in the class constructor, or through the static method findAllPerms().
Each set of gluing permutations will be produced precisely once up to equivalence, where equivalence is defined by the given set of automorphisms of the given edge pairing.
For each permutation set that is generated, routine use_ (as passed to the class constructor) will be called with that permutation set as an argument.
Once the generation of permutation sets has finished, routine use_ will be called once more, this time with null
as its first (permutation set) argument.
Subclasses corresponding to more specialised search criteria should override this routine to use a better optimised algorithm where possible.
It is possible to run only a partial search, branching to a given depth but no further. In this case, rather than producing complete gluing permutation sets, the search will produce a series of partially-complete Dim2GluingPermSearcher objects. These partial searches may then be restarted by calling runSearch() once more (usually after being frozen or passed on to a different processor). If necessary, the use_ routine may call completePermSet() to distinguish between a complete set of gluing permutations and a partial search state.
Note that a restarted search will never drop below its initial depth. That is, calling runSearch() with a fixed depth can be used to subdivide the overall search space into many branches, and then calling runSearch() on each resulting partial search will complete each of these branches without overlap.
maxDepth | the depth of the partial search to run, or a negative number if a full search should be run (the default). |
|
inherited |
Returns the total number of simplices under consideration.
|
inherited |
Returns a newly created triangulation as modelled by this set of gluing permutations and the associated simplex facet pairing.
Each matched pair of facets and their associated permutations will be realised as two simplex facets in the triangulation glued together with the corresponding gluing permutation. Each unmatched facet will be realised as a boundary facet in the triangulation.
It is the responsibility of the caller of this routine to delete this triangulation once it is no longer required.
|
protected |
The set of isomorphisms that define equivalence of gluing permutation sets.
Generally this is the set of all automorphisms of the underlying edge pairing.
|
protected |
Did we create the isomorphism list autos_ ourselves (in which case we must destroy it also)?
|
static |
A character used to identify this class when reading and writing tagged data in text format.
|
protectedinherited |
Has an error occurred during construction from an input stream?
|
protected |
Describes the order in which gluing permutations are assigned to edges.
Specifically, this order is order[0], order[1], ..., order[orderSize-1].
Note that each element of this array corresponds to a single edge of the underlying edge pairing graph, which in turn represents a triangle edge and its image under the given edge pairing.
The specific triangle edge stored in this array for each edge of the underlying edge pairing graph will be the smaller of the two identified triangle edges (unless otherwise specified by a subclass that uses a specialised search algorithm.
|
protected |
Marks which element of order[] we are currently examining at this stage of the search.
|
protected |
The total number of edges in the edge pairing graph, i.e., the number of elements of interest in the order[] array.
|
protected |
Are we only searching for gluing permutations that correspond to orientable triangulations?
|
protected |
Keeps track of the orientation of each triangle in the underlying triangulation.
Orientation is positive/negative, or 0 if unknown. Note that in some algorithms the orientation is simply +/-1, and in some algorithms the orientation counts forwards or backwards from 0 according to how many times the orientation has been set or verified.
|
protectedinherited |
The facet pairing that this permutation set complements.
This is guaranteed to be the minimal representative of its facet pairing isomorphism class.
|
protectedinherited |
The index into array Perm::Sn_1 describing how each simplex facet is glued to its partner.
Note that this is not a gluing permutation as such but rather a permutation of 0,...,dim-1 only (see the routines gluingToIndex() and indexToGluing() for conversions). If a permutation has not yet been selected (e.g., if this permutation set is still under construction) then this index is -1.
|
protected |
Has the search started yet? This helps distinguish between a new search and the resumption of a partially completed search.
|
protected |
A routine to call each time a gluing permutation set is found during the search.
|
protected |
Additional user-supplied data to be passed as the second argument to the use_ routine.