Trev's profile - activity

This is the code:

var('a_1')  
var('z')   
var('b_1')   
var('c_1')   
a = vector([a_1,sqrt(1-a_1^2-z^2),z])   
b = vector([b_1,sqrt(1-b_1^2-z^2),z])   
c = vector([c_1,sqrt(1-c_1^2-z^2),z])   
Circum = vector([0,0,z])   
H_1 = a+b+c-3*Circum   
def f(x,y):   
     return x.dot_product(y)   
H_d = f(H_1,H_1)   
d = H_1 - sqrt (1-H_d) * vector([0,0,1])   
def p(x,y):   
     return (f(x,x)*f(y,y)-f(x,y)*f(y,y))*x/(2*f(x,x)*f(y,y)-2*f(x,y)*f(x,y))+(f(x,x)*f(y,y)-f(x,y)*f(x,x))*y/(2*f(x,x)*f(y,y)-2*f(x,y)*f(x,y))   
def Circum2d(r,s,t):   
     return p((r-t),(s-t))+t   
G = 1/4*(a+b+c+d)   
H = Circum+Circum2d(a,b,d)+Circum(a,c,d)+Circum(b,c,d)   
f(H,G)/sqrt(f(H,H)*f(G,G)).simplify_rational()

For a test, I have something around 500 digits long when not simplified, and that takes less than 1 minute to simplify

I have a big chunk of expression that needs simplification, but it seems to take ages. I am not sure if Sage is not working while my computer rests for a while, or just that the simplification really takes that long time.

A range would be enough. I have roughly 30000 - 40000 digits long when not simplified.