#19872: regular symmetric Hadamard matrices for n=324
-------------------------+-------------------------------------------------
       Reporter:         |        Owner:
  dimpase                |       Status:  needs_review
           Type:         |    Milestone:  sage-7.0
  enhancement            |   Resolution:
       Priority:  major  |    Merged in:
      Component:  graph  |    Reviewers:  Nathann Cohen
  theory                 |  Work issues:
       Keywords:         |       Commit:
        Authors:  Dima   |  e123de45de77715fb84f03697e52d3d2b3c6e8a0
  Pasechnik              |     Stopgaps:
Report Upstream:  N/A    |
         Branch:         |
  public/JKandJKT        |
   Dependencies:         |
  #19861                 |
-------------------------+-------------------------------------------------

Comment (by dimpase):

 I'm working on speeding it up. Already using libgap gives a 3-fold
 speedup...
 {{{
 diff --git a/src/sage/graphs/generators/smallgraphs.py
 b/src/sage/graphs/generators/smallgraphs
 index 8ea232d..840c346 100644
 --- a/src/sage/graphs/generators/smallgraphs.py
 +++ b/src/sage/graphs/generators/smallgraphs.py
 @@ -4941,7 +4941,7 @@ def JankoKharaghaniTonchevGraph():

      EXAMPLES::

 -        sage: Gamma=graphs.JankoKharaghaniTonchevGraph()  # long time
 +        sage: Gamma=graphs.JankoKharaghaniTonchevGraph()  # long time (14
 sec)
          sage: Gamma.is_strongly_regular(parameters=True)  # long time
          (324, 153, 72, 72)

 @@ -4956,7 +4956,7 @@ def JankoKharaghaniTonchevGraph():
      from itertools import product
      from sage.misc.misc_c import prod
      from sage.combinat.permutation import Permutation as P
 -    from sage.groups.perm_gps.permgroup import PermutationGroup
 +    from sage.libs.gap.libgap import libgap

      m1=prod([P((9*x+k,9*x+k+3,9*x+k+6)) for k,x in
 product(xrange(1,4),xrange(36))])
      m2=prod([P((3*x+1,3*x+2,3*x+3)) for x in xrange(108)])
 @@ -4975,7 +4975,7 @@ def JankoKharaghaniTonchevGraph():
 (18*x+4,18*x+13),(18*x+5,18*x+14),(18*x+6,18*x+15),(18*x+7,18*x+16),
                  (18*x+8,18*x+17),(18*x+9,18*x+18)]))
                   for x in xrange(18))
 -    G=PermutationGroup([m1,m2,t,n1,n2,s,k])
 +    G=libgap.Group(map(lambda p: libgap.PermList(p),
 [m1,m2,t,n1,n2,s,k]))
 B1=(19,22,25,29,30,31,33,34,35,37,40,43,47,48,49,51,52,53,55,56,57,65,
 66,67,68,70,72,76,77,78,79,80,81,82,86,90,92,93,95,96,98,99,100,105,107,
 109,110,111,119,120,121,122,124,126,128,129,131,132,134,135,136,141,143,
 @@ -4995,6 +4995,6 @@ def JankoKharaghaniTonchevGraph():
      Gamma=Graph(multiedges=False,name='Janko-Kharaghani-Tonchev')
      for i,b in ((1,B1),(163,B163)):
          for j in b:
 -            Gamma.add_edges(map(tuple,G.orbit((i,j),action="OnSets")))
 +            Gamma.add_edges(map(tuple,G.Orbit(libgap.Set([i,j]),
 libgap.OnSets)))
      Gamma.relabel()
      return Gamma
 }}}

--
Ticket URL: <http://trac.sagemath.org/ticket/19872#comment:11>
Sage <http://www.sagemath.org>
Sage: Creating a Viable Open Source Alternative to Magma, Maple, Mathematica, 
and MATLAB

-- 
You received this message because you are subscribed to the Google Groups 
"sage-trac" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to [email protected].
To post to this group, send email to [email protected].
Visit this group at https://groups.google.com/group/sage-trac.
For more options, visit https://groups.google.com/d/optout.

Reply via email to