Re: [sage-combinat-devel] Re: the unfortunate case of affine type A twisted
Hi Anne, On Fri, Mar 09, 2012 at 04:09:13PM -0800, Anne Schilling wrote: Impressive patch! Thanks for working on this. :-) Thanks so much for your very useful and quick review. I really appreciate that! Followup on http://trac.sagemath.org/sage_trac/ticket/6588 Cheers, Nicolas -- Nicolas M. Thiéry Isil nthi...@users.sf.net http://Nicolas.Thiery.name/ -- You received this message because you are subscribed to the Google Groups sage-combinat-devel group. To post to this group, send email to sage-combinat-devel@googlegroups.com. To unsubscribe from this group, send email to sage-combinat-devel+unsubscr...@googlegroups.com. For more options, visit this group at http://groups.google.com/group/sage-combinat-devel?hl=en.
Re: [sage-combinat-devel] Re: the unfortunate case of affine type A twisted
Christian: I had to slightly rebase trac_11187-finite_reflection_groups-cs.patch for the doctests I just added to apply_simple_reflection and friends. Cheers, Nicolas -- Nicolas M. Thiéry Isil nthi...@users.sf.net http://Nicolas.Thiery.name/ -- You received this message because you are subscribed to the Google Groups sage-combinat-devel group. To post to this group, send email to sage-combinat-devel@googlegroups.com. To unsubscribe from this group, send email to sage-combinat-devel+unsubscr...@googlegroups.com. For more options, visit this group at http://groups.google.com/group/sage-combinat-devel?hl=en.
Re: [sage-combinat-devel] Re: the unfortunate case of affine type A twisted
Hi -- as you are currently discussing root systems (including 6588, thanks a lot for working on that!): here are some issues that made it very hard for me to work with root systems for finite reflection groups: # the first and the second seem to behave very differently: sage: x = CartanType([['A',4]]) sage: type(x) class 'sage.combinat.root_system.type_reducible.CartanType_with_superclass' sage: x.is_irreducible() False sage: x = CartanType(['A',4]) sage: type(x) class 'sage.combinat.root_system.type_A.CartanType' sage: x.is_irreducible() True # there is a typo here (a double ll in the middle): sage: x.is_crystalographic() True # the standard index set is not python/sage convention, is that on purpose? sage: x.index_set() [1, 2, 3, 4] # the convention for the exceptional node in types B/C is the other choice that in chevie (compare http://www.math.jussieu.fr/~jmichel/gap3/htm/chap071.htm): sage: x = CartanType(['B',4]) sage: x.dynkin_diagram() O---O---O==O 1 2 3 4 B4 # the roots are accessible only in the root space, my impression is that a root system has roots, why not having a method there calling them? # the order in which the roots are given comes from the backtracker: sage: x = CartanType(['A',4]) sage: r = x.root_system().root_space() sage:. r.positive_roots() sage.combinat.backtrack.TransitiveIdeal sage: [ beta for beta in r.positive_roots() ] [alpha[2], alpha[2] + alpha[3], alpha[1], alpha[3], alpha[1] + alpha[2], alpha[1] + alpha[2] + alpha[3], alpha[2] + alpha[3] + alpha[4], alpha[4], alpha[3] + alpha[4], alpha[1] + alpha[2] + alpha[3] + alpha[4]] # it would be great if we were able to get the same order as obtained from chevie, as this gives directly the permutation representation of the corresponding Coxeter group (acting on positions of roots). finally, the non-crystallographic types are missing. using the universal cyclotomic field we should also be able to get these root systems into sage. I would be happy to work on (parts of) the things above, but first wanted to hear what might be good to get into 6588, and your opinion on others. Best, Christian -- You received this message because you are subscribed to the Google Groups sage-combinat-devel group. To post to this group, send email to sage-combinat-devel@googlegroups.com. To unsubscribe from this group, send email to sage-combinat-devel+unsubscr...@googlegroups.com. For more options, visit this group at http://groups.google.com/group/sage-combinat-devel?hl=en.
[sage-combinat-devel] B-matrix class of finite type cluster algebras
Dear all, I am running Sage Version 4.8 with combinat compiled 3 days ago. I was trying to compute the class of extended B-matices for the principal coefficients cluster algebra of type $E_8$ and I got the following error. sage: S=ClusterSeed(['E',8]) sage: T=S.principal_extension() sage: T.b_matrix_class() --- IndexErrorTraceback (most recent call last) /home/plutone/staff/stella/sage-4.8/ipython console in module() /home/plutone/staff/stella/sage-4.8/local/lib/python2.6/site-packages/sage/combinat/cluster_algebra_quiver/cluster_seed.pyc in b_matrix_class(self, depth, up_to_equivalence) 1474 assert self.is_mutation_finite(), 'The B-matrix class can - for infinite mutation types - only be computed up to a given depth' 1475 - 1476 return [ M for M in self.b_matrix_class_iter( depth=depth, up_to_equivalence=up_to_equivalence ) ] 1477 1478 def variable_class_iter(self, depth=infinity, ignore_bipartite_belt=False): /home/plutone/staff/stella/sage-4.8/local/lib/python2.6/site-packages/sage/combinat/cluster_algebra_quiver/cluster_seed.pyc in b_matrix_class_iter(self, depth, up_to_equivalence) 1450 1451 Q = self.quiver() - 1452 for M in Q.mutation_class_iter( depth=depth, up_to_equivalence=up_to_equivalence, data_type='matrix' ): 1453 yield M 1454 /home/plutone/staff/stella/sage-4.8/local/lib/python2.6/site-packages/sage/combinat/cluster_algebra_quiver/quiver.pyc in mutation_class_iter(self, depth, show_depth, return_paths, data_type, up_to_equivalence, sink_source) 917 dg = DiGraph( self._digraph ) 918 MC_iter = _mutation_class_iter( dg, self._n, self._m, depth=depth, return_dig6=return_dig6, show_depth=show_depth, up_to_equivalence=up_to_equivalence, sink_source=sink_source ) -- 919 for data in MC_iter: 920 if data_type == quiver: 921 Q = Quiver( data[0] ) /home/plutone/staff/stella/sage-4.8/local/lib/python2.6/site-packages/sage/combinat/cluster_algebra_quiver/mutation_class.pyc in _mutation_class_iter(dg, n, m, depth, return_dig6, show_depth, up_to_equivalence, sink_source) 167 nr_mutations = 1 168 if up_to_equivalence: -- 169 iso, orbits = _dg_canonical_form( dg, n, m ) 170 iso_inv = dict( (iso[a],a) for a in iso ) 171 /home/plutone/staff/stella/sage-4.8/local/lib/python2.6/site-packages/sage/combinat/cluster_algebra_quiver/mutation_class.pyc in _dg_canonical_form(dg, n, m) 143 if v = n: 144 del iso[v] -- 145 v1,v2,label1 = [ edge for edge in dg._backend.iterator_in_edges([v],True) ][0] 146 w1,w2,label2 = [ edge for edge in dg._backend.iterator_out_edges([v],True) ][0] 147 dg._backend.del_edge(v1,v2,label1,True) IndexError: list index out of range The same error reproduces in may other cases of principal extensions. Is there some issue with dg._backend.iterator_in_edges for rectangular matrices? Adding the flag up_to_equivalence=False to b_matrix_class I get the expected answer (I assume this flag prevent _dg_canonical_form from running) but my hardware is too small to handle the case $E_8$ without the reduction up to equivalece. Is there any workaround for this problem? Thanks S. -- You received this message because you are subscribed to the Google Groups sage-combinat-devel group. To post to this group, send email to sage-combinat-devel@googlegroups.com. To unsubscribe from this group, send email to sage-combinat-devel+unsubscr...@googlegroups.com. For more options, visit this group at http://groups.google.com/group/sage-combinat-devel?hl=en.
Re: [sage-combinat-devel] B-matrix class of finite type cluster algebras
Dear Salvatore, sage: S=ClusterSeed(['E',8]) sage: T=S.principal_extension() sage: T.b_matrix_class() that was a bug in our code (that was there actually from the beginning on). I hope I fixed it; I pushed the changes, could you please recheck that you get the right numbers (if that is actually possible!) in some small examples like S=ClusterSeed(X) T=S.principal_extension() T.b_matrix_class() for X in types A and B in rank 2 and 3? Thanks for reporting, Christian -- You received this message because you are subscribed to the Google Groups sage-combinat-devel group. To post to this group, send email to sage-combinat-devel@googlegroups.com. To unsubscribe from this group, send email to sage-combinat-devel+unsubscr...@googlegroups.com. For more options, visit this group at http://groups.google.com/group/sage-combinat-devel?hl=en.
Re: [sage-combinat-devel] B-matrix class of finite type cluster algebras
Dear Christian, at the moment I am unable to pull: the update fails when applying trac_6588-categories-root_systems-review-nt.patch which was pushed to the queue some hours ago. I will try again shortly and report back. Thank you for the quick fix S. * Christian Stump christian.st...@gmail.com [2012-03-10 18:43:30]: Dear Salvatore, sage: S=ClusterSeed(['E',8]) sage: T=S.principal_extension() sage: T.b_matrix_class() that was a bug in our code (that was there actually from the beginning on). I hope I fixed it; I pushed the changes, could you please recheck that you get the right numbers (if that is actually possible!) in some small examples like S=ClusterSeed(X) T=S.principal_extension() T.b_matrix_class() for X in types A and B in rank 2 and 3? Thanks for reporting, Christian -- You received this message because you are subscribed to the Google Groups sage-combinat-devel group. To post to this group, send email to sage-combinat-devel@googlegroups.com. To unsubscribe from this group, send email to sage-combinat-devel+unsubscr...@googlegroups.com. For more options, visit this group at http://groups.google.com/group/sage-combinat-devel?hl=en. -- You received this message because you are subscribed to the Google Groups sage-combinat-devel group. To post to this group, send email to sage-combinat-devel@googlegroups.com. To unsubscribe from this group, send email to sage-combinat-devel+unsubscr...@googlegroups.com. For more options, visit this group at http://groups.google.com/group/sage-combinat-devel?hl=en.
Re: [sage-combinat-devel] B-matrix class of finite type cluster algebras
Sorry I meant there are problem applying trac_6588-categories-root_systems-review-nt.patch S. * Christian Stump christian.st...@gmail.com [2012-03-10 18:43:30]: Dear Salvatore, sage: S=ClusterSeed(['E',8]) sage: T=S.principal_extension() sage: T.b_matrix_class() that was a bug in our code (that was there actually from the beginning on). I hope I fixed it; I pushed the changes, could you please recheck that you get the right numbers (if that is actually possible!) in some small examples like S=ClusterSeed(X) T=S.principal_extension() T.b_matrix_class() for X in types A and B in rank 2 and 3? Thanks for reporting, Christian -- You received this message because you are subscribed to the Google Groups sage-combinat-devel group. To post to this group, send email to sage-combinat-devel@googlegroups.com. To unsubscribe from this group, send email to sage-combinat-devel+unsubscr...@googlegroups.com. For more options, visit this group at http://groups.google.com/group/sage-combinat-devel?hl=en. -- You received this message because you are subscribed to the Google Groups sage-combinat-devel group. To post to this group, send email to sage-combinat-devel@googlegroups.com. To unsubscribe from this group, send email to sage-combinat-devel+unsubscr...@googlegroups.com. For more options, visit this group at http://groups.google.com/group/sage-combinat-devel?hl=en.
Re: [sage-combinat-devel] Re: the unfortunate case of affine type A twisted
Hi Christian, as you are currently discussing root systems (including 6588, thanks a lot for working on that!): here are some issues that made it very hard for me to work with root systems for finite reflection groups: # the first and the second seem to behave very differently: sage: x = CartanType([['A',4]]) sage: type(x) class 'sage.combinat.root_system.type_reducible.CartanType_with_superclass' sage: x.is_irreducible() False sage: x = CartanType(['A',4]) sage: type(x) class 'sage.combinat.root_system.type_A.CartanType' sage: x.is_irreducible() True This is definitely strange. # there is a typo here (a double ll in the middle): sage: x.is_crystalographic() True Please fix this! # the standard index set is not python/sage convention, is that on purpose? sage: x.index_set() [1, 2, 3, 4] I think this is intentional. The nodes in the Dynkin diagram are labeled 1,2,...,rank of g. # the convention for the exceptional node in types B/C is the other choice that in chevie (compare http://www.math.jussieu.fr/~jmichel/gap3/htm/chap071.htm): sage: x = CartanType(['B',4]) sage: x.dynkin_diagram() O---O---O==O 1 2 3 4 B4 There are many different conventions. I think in sage the default is the one used in Kac. There is also a relabeling option if you want to use a different convention. # the roots are accessible only in the root space, my impression is that a root system has roots, why not having a method there calling them? I am confused. This also exists for other spaces: sage: R = RootSystem(['A',4]) sage: R.weight_lattice().roots() [-Lambda[2] + 2*Lambda[3] - Lambda[4], -Lambda[1] + Lambda[2] + Lambda[3] - Lambda[4], 2*Lambda[1] - Lambda[2], -Lambda[1] + 2*Lambda[2] - Lambda[3], Lambda[1] + Lambda[2] - Lambda[3], -Lambda[1] + Lambda[2] + Lambda[4], Lambda[1] + Lambda[3] - Lambda[4], -Lambda[2] + Lambda[3] + Lambda[4], -Lambda[3] + 2*Lambda[4], Lambda[1] + Lambda[4], Lambda[2] - 2*Lambda[3] + Lambda[4], Lambda[1] - Lambda[2] - Lambda[3] + Lambda[4], -2*Lambda[1] + Lambda[2], Lambda[1] - 2*Lambda[2] + Lambda[3], -Lambda[1] - Lambda[2] + Lambda[3], Lambda[1] - Lambda[2] - Lambda[4], -Lambda[1] - Lambda[3] + Lambda[4], Lambda[2] - Lambda[3] - Lambda[4], Lambda[3] - 2*Lambda[4], -Lambda[1] - Lambda[4]] sage: R.root_lattice().roots() [alpha[2], alpha[2] + alpha[3], alpha[1], alpha[3], alpha[1] + alpha[2], alpha[1] + alpha[2] + alpha[3], alpha[2] + alpha[3] + alpha[4], alpha[4], alpha[3] + alpha[4], alpha[1] + alpha[2] + alpha[3] + alpha[4], -alpha[2], -alpha[2] - alpha[3], -alpha[1], -alpha[3], -alpha[1] - alpha[2], -alpha[1] - alpha[2] - alpha[3], -alpha[2] - alpha[3] - alpha[4], -alpha[4], -alpha[3] - alpha[4], -alpha[1] - alpha[2] - alpha[3] - alpha[4]] I would be happy to work on (parts of) the things above, but first wanted to hear what might be good to get into 6588, and your opinion on others. As far as I understand Nicolas would like to finish this patch soon. If there are some things that can be implemented quickly, why don't you write a quick patch on top of Nicolas? Otherwise it might be better to open a new ticket. Best, Anne -- You received this message because you are subscribed to the Google Groups sage-combinat-devel group. To post to this group, send email to sage-combinat-devel@googlegroups.com. To unsubscribe from this group, send email to sage-combinat-devel+unsubscr...@googlegroups.com. For more options, visit this group at http://groups.google.com/group/sage-combinat-devel?hl=en.
Re: [sage-combinat-devel] Re: sage.combinat.sf.sf into the reference manual?
On Fri, Mar 09, 2012 at 04:47:50PM +0100, Nicolas M. Thiery wrote: On Fri, Mar 09, 2012 at 06:51:15AM -0800, Simon King wrote: Still: Is there a list of available markups? Florent? Standard ReST markup: http://docutils.sourceforge.net/docs/ref/rst/directives.html Sphinx specific markup: http://sphinx.pocoo.org/markup/index.html Florent -- You received this message because you are subscribed to the Google Groups sage-combinat-devel group. To post to this group, send email to sage-combinat-devel@googlegroups.com. To unsubscribe from this group, send email to sage-combinat-devel+unsubscr...@googlegroups.com. For more options, visit this group at http://groups.google.com/group/sage-combinat-devel?hl=en.
Re: [sage-combinat-devel] B-matrix class of finite type cluster algebras
Dear Christian, I can confirm that your patch does the job. It took me some time to veryfy because I had to install sage 5 (the patch trac_6588-categories-root_systems-review-nt.patch does not apply to 4.8 (even to a clean install) Thank you S. * Christian Stump christian.st...@gmail.com [2012-03-10 18:43:30]: Dear Salvatore, sage: S=ClusterSeed(['E',8]) sage: T=S.principal_extension() sage: T.b_matrix_class() that was a bug in our code (that was there actually from the beginning on). I hope I fixed it; I pushed the changes, could you please recheck that you get the right numbers (if that is actually possible!) in some small examples like S=ClusterSeed(X) T=S.principal_extension() T.b_matrix_class() for X in types A and B in rank 2 and 3? Thanks for reporting, Christian -- You received this message because you are subscribed to the Google Groups sage-combinat-devel group. To post to this group, send email to sage-combinat-devel@googlegroups.com. To unsubscribe from this group, send email to sage-combinat-devel+unsubscr...@googlegroups.com. For more options, visit this group at http://groups.google.com/group/sage-combinat-devel?hl=en. -- You received this message because you are subscribed to the Google Groups sage-combinat-devel group. To post to this group, send email to sage-combinat-devel@googlegroups.com. To unsubscribe from this group, send email to sage-combinat-devel+unsubscr...@googlegroups.com. For more options, visit this group at http://groups.google.com/group/sage-combinat-devel?hl=en.
[sage-combinat-devel] Re: sage.combinat.sf.sf into the reference manual?
Hi Florent, On 10 Mrz., 10:07, Florent Hivert florent.hiv...@lri.fr wrote: Standard ReST markup: http://docutils.sourceforge.net/docs/ref/rst/directives.html Sphinx specific markup: http://sphinx.pocoo.org/markup/index.html Thank you! Simon -- You received this message because you are subscribed to the Google Groups sage-combinat-devel group. To post to this group, send email to sage-combinat-devel@googlegroups.com. To unsubscribe from this group, send email to sage-combinat-devel+unsubscr...@googlegroups.com. For more options, visit this group at http://groups.google.com/group/sage-combinat-devel?hl=en.