you might try to find a combinatorial invariant distinguishing graphs here, 
e.g. the 4-vertex condition
(introduced in http://www.ams.org/mathscinet-getitem?mr=340088
more references 
here: http://link.springer.com/article/10.1007/s10801-014-0554-1)


On Thursday, November 3, 2016 at 12:13:52 PM UTC, Paul Leopardi wrote:
>
>
>
> On Thursday, 3 November 2016 04:58:11 UTC+11, Dima Pasechnik wrote:
>>
>> On Wednesday, November 2, 2016 at 11:57:20 AM UTC, Paul Leopardi wrote:
>>>
>>> On Wednesday, 2 November 2016 20:08:20 UTC+11, Dima Pasechnik wrote:
>>>>
>>>> I must say I don't know why these libgap.set_global and get_global are 
>>>> needed.
>>>>
>>>>>
>>> I used them to get my large matrices into libgap as globals. Trying to 
>>> pass them through a text based interface does not work.
>>>
>>
>> ld you show an example of this? 
>> IMHO you might indeed have a problem if you use libgap.eval(). But you 
>> should not...
>> The other thought is that you appear to work with Cayley graphs, and for 
>> them you only need to pass one row of the adjacency matrix,
>> and the group; passing the whole matrix is an overkill.
>>
>>
> Hi Dima,
> Thanks for your patience and for all your help. I finally got the 
> following to work:
>
>
> gap.load_package('grape')
>
>
> def check_graphs_using_gap(g_0, g_1):
>     m_0 = g_0.adjacency_matrix()
>     m_1 = g_1.adjacency_matrix()
>     dim = str(m_0.dimensions()[0])
>     gap.eval('M_0 := '+gap(m_0).str())
>     gap.eval('M_1 := '+gap(m_1).str())
>     gap.eval('GR := Group( () )')
>     gap.eval('G_0 := Graph( GR, [1..'+dim+'], OnPoints, function(x,y) 
> return M_0[x][y] = 1; end,true );')
>     gap.eval('G_1 := Graph( GR, [1..'+dim+'], OnPoints, function(x,y) 
> return M_1[x][y] = 1; end,true );')
>     return gap.eval('IsIsomorphicGraph(G_0, G_1);')
>
> Assuming that the algorithm used is Nauty rather than Bliss, running this 
> code against my graphs confirms the cases of non-isomorphism that I 
> previously identified.
>
>

-- 
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