X-Git-Url: http://gitweb.michael.orlitzky.com/?a=blobdiff_plain;f=mjo%2Feja%2Feja_utils.py;h=8422fbff3c3a3f1523a84708ee659bd605da7ffe;hb=71ed5e9dc86ef368e81e13122aad6046bf056a28;hp=79d8ecfce61c555375deffd40501b2fb100c0379;hpb=40850626cb85d115363995ad6beaee8cb17d83da;p=sage.d.git diff --git a/mjo/eja/eja_utils.py b/mjo/eja/eja_utils.py index 79d8ecf..8422fbf 100644 --- a/mjo/eja/eja_utils.py +++ b/mjo/eja/eja_utils.py @@ -1,41 +1,5 @@ from sage.structure.element import is_Matrix -def _charpoly_sage_input(s): - r""" - Helper function that you can use on the string output from sage - to convert a charpoly coefficient into the corresponding input - to be cached. - - SETUP:: - - sage: from mjo.eja.eja_algebra import JordanSpinEJA - sage: from mjo.eja.eja_utils import _charpoly_sage_input - - EXAMPLES:: - - sage: J = JordanSpinEJA(4,QQ) - sage: a = J._charpoly_coefficients() - sage: a[0] - X1^2 - X2^2 - X3^2 - X4^2 - sage: _charpoly_sage_input(str(a[0])) - 'X[0]**2 - X[1]**2 - X[2]**2 - X[3]**2' - - """ - import re - - exponent_out = r"\^" - exponent_in = r"**" - - digit_out = r"X([0-9]+)" - - def replace_digit(m): - # m is a match object - return "X[" + str(int(m.group(1)) - 1) + "]" - - s = re.sub(exponent_out, exponent_in, s) - return re.sub(digit_out, replace_digit, s) - - def _scale(x, alpha): r""" Scale the vector, matrix, or cartesian-product-of-those-things @@ -301,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()