X-Git-Url: http://gitweb.michael.orlitzky.com/?a=blobdiff_plain;f=mjo%2Feja%2Feja_utils.py;h=8422fbff3c3a3f1523a84708ee659bd605da7ffe;hb=71ed5e9dc86ef368e81e13122aad6046bf056a28;hp=7c6c581329dd64aaa8f3d958e00b2f9eb3221fcb;hpb=f38e111cd026f1ecdab6dd2d2ed194a8745252a8;p=sage.d.git diff --git a/mjo/eja/eja_utils.py b/mjo/eja/eja_utils.py index 7c6c581..8422fbf 100644 --- a/mjo/eja/eja_utils.py +++ b/mjo/eja/eja_utils.py @@ -1,6 +1,4 @@ -from sage.functions.other import sqrt -from sage.matrix.constructor import matrix -from sage.modules.free_module_element import vector +from sage.structure.element import is_Matrix def _scale(x, alpha): r""" @@ -111,9 +109,16 @@ def _all2list(x): # first needing to convert them to a list of octonions and # then recursing down into the list. It also avoids the wonky # list(x) when x is an element of a CFM. I don't know what it - # returns but it aint the coordinates. This will fall through - # to the iterable case the next time around. - return _all2list(x.to_vector()) + # returns but it aint the coordinates. We don't recurse + # because vectors can only contain ring elements as entries. + return x.to_vector().list() + + if is_Matrix(x): + # This sucks, but for performance reasons we don't want to + # call _all2list recursively on the contents of a matrix + # when we don't have to (they only contain ring elements + # as entries) + return x.list() try: xl = list(x) @@ -124,15 +129,8 @@ def _all2list(x): # Avoid the retardation of list(QQ(1)) == [1]. return [x] - return sum(list( map(_all2list, xl) ), []) - - + return sum( map(_all2list, xl) , []) -def _mat2vec(m): - return vector(m.base_ring(), m.list()) - -def _vec2mat(v): - return matrix(v.base_ring(), sqrt(v.degree()), v.list()) def gram_schmidt(v, inner_product=None): """ @@ -267,8 +265,6 @@ def gram_schmidt(v, inner_product=None): # cool return v - R = v[0].base_ring() - # Our "zero" needs to belong to the right space for sum() to work. zero = v[0].parent().zero()