This function checks the isomorphism of two balanced covering stars of
PG(n-1,2). If they are isomorphic, it returns the list of isomorphism establishing collineations. The option is provided to enumerate all isomorphism establishing collineations or to terminate after the first one is found.
A star of
A star of
An indicator of whether all isomorphism establishing collineations (IECs) should be returned (default), or terminate only after the first one is found.
This code considers all possible collineations of
PG(n-1,2) to search for isomorphism establishing collineations from star1 to star2. The search is conducted by first projecting onto a lower dimensional space described as Algorithm 2 in Spencer et al. (2019). Equivalence is assessed using the bitstring comparison method described in Spencer et al. (2019). Both input stars should be formatted as 3-dimensional arrays with
star1[i,j,k] indicating whether or not the
ith basic factor is present in the
jth effect of the
kth flat of star1.
A list containing two objects. The first object is a Boolean indicating whether or not star1 is isomorphic to star2. If isomorphic, the second object is a list of isomorphism establishing collineation matrices. If not isomorphic, the second object is NA.
Neil Spencer, Pritam Ranjan, Franklin Mendivil
Spencer, N.A., Ranjan, P., and Mendivil, F., (2019), "Isomorphism Check for 2^n Factorial Designs with Randomization Restrictions", Journal of Statistical Theory and Practice, 13(60),1-24 [https://doi.org/10.1007/s42519-019-0064-5]
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
## Example 1: Two stars of PG(4,2) consisting of 4-flats data(starn5t3a) data(starn5t3b) test1 <- checkStarIsomorphism(starn5t3a, starn5t3b, returnfirstIEC = TRUE) test1$result # the test indicates that they are isomorphic (IECstar <- test1$IECs[]) # the first IEC ## Example 2: Two stars of PG(7,2) consisting of 6-flats data(starn8t5a) data(starn8t5b) test2 <- checkStarIsomorphism(starn8t5a, starn8t5b, returnfirstIEC = TRUE) test2$result
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.