Ask Your Question
2

memory error when computing restricted_automorphism_group

asked 2020-10-09 13:28:02 +0100

selia gravatar image

I am trying to compute the restricted_automorphism_groups of permutahedra. Starting with the 5 dimensional permutahredon in R^6, I get the memory error shown below. Is there anything I can do to work around this? For what I am doing, I am really just interested in a small subgroup of the restricted_automorphism_group, but I don't know how to generate it without first generating the entire group. Thank you for any feedback!

sage: P6 = polytopes.permutahedron(6,backend = 'normaliz')                                                          
sage: G6 = P6.restricted_automorphism_group(output = 'permutation')                                                 
Python(1290,0x1086055c0) malloc: can't allocate region
*** mach_vm_map(size=18446744068325527552) failed (error code=3)
Python(1290,0x1086055c0) malloc: *** set a breakpoint in malloc_error_break to debug
---------------------------------------------------------------------------
KeyError                                  Traceback (most recent call last)
~/sage/local/lib/python3.7/site-packages/sage/misc/cachefunc.pyx in sage.misc.cachefunc.CachedMethodCaller.__call__ (build/cythonized/sage/misc/cachefunc.c:10304)()
   1942             try:
-> 1943                 return cache[k]
   1944             except TypeError:  # k is not hashable

KeyError: (('permutation',), ())

During handling of the above exception, another exception occurred:

MemoryError                               Traceback (most recent call last)
<ipython-input-2-3ba71fcc8be0> in <module>
----> 1 G6 = P6.restricted_automorphism_group(output = 'permutation')

~/sage/local/lib/python3.7/site-packages/sage/misc/cachefunc.pyx in sage.misc.cachefunc.CachedMethodCaller.__call__ (build/cythonized/sage/misc/cachefunc.c:10438)()
   1946                 return cache[k]
   1947         except KeyError:
-> 1948             w = self._instance_call(*args, **kwds)
   1949             cache[k] = w
   1950             return w

~/sage/local/lib/python3.7/site-packages/sage/misc/cachefunc.pyx in sage.misc.cachefunc.CachedMethodCaller._instance_call (build/cythonized/sage/misc/cachefunc.c:9917)()
   1822             True
   1823         """
-> 1824         return self.f(self._instance, *args, **kwds)
   1825 
   1826     cdef fix_args_kwds(self, tuple args, dict kwds):

~/sage/local/lib/python3.7/site-packages/sage/geometry/polyhedron/base.py in restricted_automorphism_group(self, output)
   9616                 G.add_edge(index0+i, index0+j, edge_label(i, j, c_ij))
   9617 
-> 9618         permgroup = G.automorphism_group(edge_labels=True)
   9619         if output == "permutation":
   9620             return permgroup

~/sage/local/lib/python3.7/site-packages/sage/graphs/generic_graph.py in automorphism_group(self, partition, verbosity, edge_labels, order, return_group, orbits, algorithm)
  21994             GC = HB.c_graph()[0]
  21995             partition = [[G_to[vv] for vv in cell] for cell in partition]
> 21996             A = search_tree(GC, partition, lab=False, dict_rep=True, dig=dig, verbosity=verbosity, order=order)
  21997             if order:
  21998                 a,b,c = A

~/sage/local/lib/python3.7/site-packages/sage/groups/perm_gps/partn_ref/refinement_graphs.pyx in sage.groups.perm_gps.partn_ref.refinement_graphs.search_tree (build/cythonized/sage/groups/perm_gps/partn_ref/refinement_graphs.c:13503)()
    458 
    459     lab_new = lab or certificate
--> 460     output = get_aut_gp_and_can_lab(<void *>GS, part, G.num_verts, &all_children_are_equivalent, &refine_by_degree, &compare_graphs, lab, NULL, NULL, NULL)
    461     sig_free( GS.scratch )
    462     # prepare output

~/sage/local/lib/python3.7/site-packages/sage/groups/perm_gps/partn_ref/automorphism_group_canonical_label.pyx in sage.groups.perm_gps.partn_ref.automorphism_group_canonical_label.get_aut_gp_and_can_lab (build/cythonized/sage/groups/perm_gps/partn_ref/automorphism_group_canonical_label.c:11657)()
    503             if output_prealloc is NULL:
    504                 deallocate_agcl_output(output)
--> 505             raise MemoryError
    506 
    507     # Allocate:

MemoryError:
edit retag flag offensive close merge delete

1 Answer

Sort by ยป oldest newest most voted
1

answered 2020-10-10 14:35:53 +0100

vdelecroix gravatar image

On Sage 9.2.beta14 with both bliss and nauty installed I obtain the group without an error.

edit flag offensive delete link more

Comments

Turned comment into an answer so it can be accepted to mark the question as solved.

slelievre gravatar imageslelievre ( 2021-01-02 15:49:39 +0100 )edit

Your Answer

Please start posting anonymously - your entry will be published after you log in or create a new account.

Add Answer

Question Tools

1 follower

Stats

Asked: 2020-10-09 13:28:02 +0100

Seen: 263 times

Last updated: Oct 09 '20