Right, this looks like a recent regression. AFAIR, it used to work.

If you use gap.* rather than libgap.*, things seem to work for me, e.g.

sage: gap.load_package('grape')
sage: gap.eval('GG:=JohnsonGraph(5,2);')
'rec( adjacencies := [ [ 2, 3, 4, 5, 6, 7 ] ], group := Group([ 
(1,5,8,10,4)\n  (2,6,9,3,7), (2,5)(3,6)(4,7) ]), isGraph := true, isSimple 
:= true, \n  names := [ [ 1, 2 ], [ 1, 3 ], [ 1, 4 ], [ 1, 5 ], [ 2, 3 ], [ 
2, 4 ], \n      [ 2, 5 ], [ 3, 4 ], [ 3, 5 ], [ 4, 5 ] ], order := 10, \n 
 representatives := [ 1 ], schreierVector := [ -1, 2, 2, 1, 1, 1, 2, 1, 1, 
1 \n     ] )'
sage: gap.eval('AutGroupGraph(GG);')
'Group([ (3,4)(6,7)(8,9), (2,3)(5,6)(9,10), (2,5)(3,6)(4,7), 
(1,2)(6,8)(7,9) ])'
sage: 

(whereas if I use libgap.* I end up with the same error as you)

This is now
https://trac.sagemath.org/ticket/21794


On Tuesday, November 1, 2016 at 12:48:05 PM UTC, Paul Leopardi wrote:
>
>
>
> On Tuesday, 25 October 2016 06:02:06 UTC+11, Dima Pasechnik wrote:
>>
>> I would do a check using GAP's Grape package, which allows for checking 
>> isomorphisms
>> (it uses nauty as the backend)
>> While there is no ready function to call Grape from Sage, this should be 
>> easy to write using e.g.
>> libgap.function_factory()
>>
>
> I finally managed to figure out how to get graphs from Sage to Grape, but 
> now Grape can't find nauty. What do I need to do to configure Grape 
> properly?
>
> sage: GG=libgap.eval('GG := Graph( G, [1..16], OnPoints, function(x,y) 
> return MCG[x][y] = 1; end,true );')
> sage: GG
> rec( adjacencies := [ [ 2, 3, 5, 7, 10, 11 ], [ 1, 4, 6, 7, 10, 12 ], [ 1, 
> 4, 5, 8, 11, 13 ], [ 2, 3, 6, 8, 12, 13 ], [ 1, 3, 6, 9, 11, 14 ], [ 2, 4, 
> 5, 9, 12, 14 ], [ 1, 2, 8, 9, 10, 15 ], [ 3, 4, 7, 9, 13, 15 ], [ 5, 6, 7, 
> 8, 14, 15 ], [ 1, 2, 7, 13, 14, 16 ], [ 1, 3, 5, 12, 15, 16 ], [ 2, 4, 6, 
> 11, 15, 16 ], [ 3, 4, 8, 10, 14, 16 ], [ 5, 6, 9, 10, 13, 16 ], [ 7, 8, 9, 
> 11, 12, 16 ], [ 10, 11, 12, 13, 14, 15 ] ], group := Group(()), isGraph := 
> true, names := [ 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16 ], 
> order := 16, representatives := [ 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 
> 13, 14, 15, 16 ], schreierVector := [ -1, -2, -3, -4, -5, -6, -7, -8, -9, 
> -10, -11, -12, -13, -14, -15, -16 ] )
> sage: HH=libgap.eval('HH := Graph( G, [1..16], OnPoints, function(x,y) 
> return MLG[x][y] = 1; end,true );')
> sage: HH
> rec( adjacencies := [ [ 2, 3, 5, 7, 10, 11 ], [ 1, 4, 6, 7, 10, 12 ], [ 1, 
> 4, 5, 8, 11, 13 ], [ 2, 3, 6, 8, 12, 13 ], [ 1, 3, 6, 9, 11, 14 ], [ 2, 4, 
> 5, 9, 12, 14 ], [ 1, 2, 8, 9, 10, 15 ], [ 3, 4, 7, 9, 13, 15 ], [ 5, 6, 7, 
> 8, 14, 15 ], [ 1, 2, 7, 13, 14, 16 ], [ 1, 3, 5, 12, 15, 16 ], [ 2, 4, 6, 
> 11, 15, 16 ], [ 3, 4, 8, 10, 14, 16 ], [ 5, 6, 9, 10, 13, 16 ], [ 7, 8, 9, 
> 11, 12, 16 ], [ 10, 11, 12, 13, 14, 15 ] ], group := Group(()), isGraph := 
> true, names := [ 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16 ], 
> order := 16, representatives := [ 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 
> 13, 14, 15, 16 ], schreierVector := [ -1, -2, -3, -4, -5, -6, -7, -8, -9, 
> -10, -11, -12, -13, -14, -15, -16 ] )
> sage: 
> libgap.eval('GraphIsomorphism(GG,HH)')                                        
>                      
>  
> ---------------------------------------------------------------------------
> ValueError                                Traceback (most recent call last)
> <ipython-input-80-325cd83b9cf4> in <module>()
> ----> 1 libgap.eval('GraphIsomorphism(GG,HH)')
>
> /home/leopardi/opt/sage/sage-7.2/src/sage/libs/gap/libgap.pyx in 
> sage.libs.gap.libgap.Gap.eval 
> (/home/leopardi/opt/sage/sage-7.2/src/build/cythonized/sage/libs/gap/libgap.c:4183)()
>     429         if not isinstance(gap_command, basestring):
>     430             gap_command = str(gap_command._gap_init_())
> --> 431         return make_any_gap_element(self, gap_eval(gap_command))
>     432 
>     433     @cached_method
>
> /home/leopardi/opt/sage/sage-7.2/src/sage/libs/gap/util.pyx in 
> sage.libs.gap.util.gap_eval 
> (/home/leopardi/opt/sage/sage-7.2/src/build/cythonized/sage/libs/gap/util.c:4947)()
>     286             sig_off()
>     287         except RuntimeError as msg:
> --> 288             raise ValueError('libGAP: '+str(msg).strip())
>     289 
>     290         if libGAP_Symbol != libGAP_S_SEMICOLON:
>
> ValueError: libGAP: d
> $1n16g
> 2
> 3
> 5
> 7
> 10
> 11;
> 1
> 4
> 6
> 7
> 10
> 12;
> 1
> 4
> 5
> 8
> 11
> 13;
> 2
> 3
> 6
> 8
> 12
> 13;
> 1
> 3
> 6
> 9
> 11
> 14;
> 2
> 4
> 5
> 9
> 12
> 14;
> 1
> 2
> 8
> 9
> 10
> 15;
> 3
> 4
> 7
> 9
> 13
> 15;
> 5
> 6
> 7
> 8
> 14
> 15;
> 1
> 2
> 7
> 13
> 14
> 16;
> 1
> 3
> 5
> 12
> 15
> 16;
> 2
> 4
> 6
> 11
> 15
> 16;
> 3
> 4
> 8
> 10
> 14
> 16;
> 5
> 6
> 9
> 10
> 13
> 16;
> 7
> 8
> 9
> 11
> 12
> 16;
> 10
> 11
> 12
> 13
> 14
> 15.
> > /tmp/tmd2qUpq/ftmp1 p,cx
> >> /tmp/tmd2qUpq/ftmp2 bq
> 16
> Error, cannot find output produced by `dreadnaut'
>
>

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

Reply via email to