Alex Karenin's profile - activity

Thank you! It would work just fine for me.

Overflow in list_plot? The code: l = [1.4312706585e12, -1.4575266189e12, 4.9459419278e11, -4.8146621292e12] list_plot(l

Does sympy break CyclotomicField? One can define a cyclotomic field as CyclotomicField(16) for example. But the code bel

Are rings of integers incompatible with multiprocessing? The code below crashen on sage 9.5: from multiprocessing impor

Sage fails to check if NumberField if UFD The field Q(sqrt(-5)) if not a unique factorization domain. Yet the code below

Sage fails to check if NumberField if UFD The field Q(sqrt(-5)) if not a unique factorization domain. Yet the code below

Note that if you define OK = K.ring_of_integers() L, U = matrix(OK,L), matrix(OK,U) then it magically works. Also

Note that if you define OK = K.ring_of_integers() L, U = matrix(OK,L), matrix(OK,U) then it magically works. Also

UPD: the code works as intended if I define matrices L and U over the ring of integers of the field K.

UPD: the code works as intended if I definine matrices L and U over the ring of integers of the field K.

Determinants over cyclotomic fields are broken? I launch the code below which generates two unimodular matrices L and U

NameError: name '...' is not defined I got the preparsed with sage --preparse command keflll.py file with the function i

MAC M1 Unhandled SIGILL: An illegal instruction occurred. I got MAC OS 12.4 installed and after attempting to do matrix

That's my suspition as well. Is there I can do anything about this issue? I'd like to be able to handle 63 bit numbers.

FractionFieldElement works incorrectly over large fields The code below returns a*b/a instead of b if characteristic of

This code return proof that f^8(x)=f(f(f(f(f(f(f(f(x))))))))=x if f(x)=(x(sqrt(2)-1)-1)/(x-1):

def yx(x): 
    cst=sqrt(2)-1
    if(x==1):
        return -Infinity
    return (cst*x-1)/(x-1)
tmp=x 
for t in range(1,9): 
    tmp=(yx(tmp)).canonicalize_radical() 
    print(t, tmp.canonicalize_radical())

But this code looks directly at orbit of a point and shows that the statement above is false:

tmp=1.5
lst=[]
for t in range(1,9): 
    tmp=(yx(tmp)).canonicalize_radical() 
    print(t, tmp.canonicalize_radical().n(80))
    lst.append((t,tmp.n()))

By giving the next output:

1 -0.75735931288071485359493 2 0.74754689570642838213503 3 2.7345908033901357840028 4 0.076504843704676803069591 5 1.0485281374238889523706 6 0.13254046199725214558851 7 -0.090964954340477885594991 8 0.35391568118499587774957

Why the output is so different? Bit precision don't seem to be a problem.

This code return proof that f^8(x)=f(f(f(f(f(f(f(f(x))))))))=x if f(x)=(x(sqrt(2)-1)-1)/(x-1):

def yx(x): 
    cst=sqrt(2)-1
    if(x==1):
        return -Infinity
    return (cst*x-1)/(x-1)
tmp=x 
for t in range(1,9): 
    tmp=(yx(tmp)).canonicalize_radical() 
    print(t, tmp.canonicalize_radical())

But this code looks directly at orbit of a point and shows that the statement above is false:

tmp=1.5
lst=[]
for t in range(1,9): 
    tmp=(yx(tmp)).canonicalize_radical() 
    print(t, tmp.canonicalize_radical().n(80))
    lst.append((t,tmp.n()))

By giving the next output:

1 -0.75735931288071485359493
2 0.74754689570642838213503
3 2.7345908033901357840028
4 0.076504843704676803069591
5 1.0485281374238889523706
6 0.13254046199725214558851
7 -0.090964954340477885594991
8 0.35391568118499587774957

Why the output is so different? Bit precision don't seem to be a problem.

I am new to SAGE. While attempting to substitute finite field element into polynomial over finite field instead of variable i get positive characteristic not allowed in symbolic computations error. How am i supposed to substitute in this case? Sage notebook 8.5. Not working code:

x,y=var('x'),var('y')
fx=PolynomialRing(GF(32),x)(x^2+x)
fx.substitute(x==1)

Thanks in advanse fot answer.

I am new to Sage. While attempting to substitute finite field element into polynomial over finite field, instead of variable i get positive characteristic not allowed in symbolic computations error. How am i supposed to substitute in this case? Sage notebook 8.5. Not working code:

x, y = var('x'), var('y')
fx = PolynomialRing(GF(32), x)(x^2 + x)
fx.substitute(x == 1)

Thanks in advance for answer.