Continued fraction/Arithmetic/G(matrix ng, continued fraction n1, continued fraction n2): Difference between revisions

Content added Content deleted
Line 3,768: Line 3,768:
else:
else:
term = x[env[ix]]
term = x[env[ix]]
env[ix] += 1
if term is not None:
if term is not None:
new_ng = (a2 + (a12 * term), a + (a1 * term), a12, a1,
new_ng = (a2 + (a12 * term), a + (a1 * term), a12, a1,
Line 3,773: Line 3,774:
if not too_big (new_ng):
if not too_big (new_ng):
env[ng] = new_ng
env[ng] = new_ng
env[ix] += 1
else:
else:
env[ng] = (a12, a1, a12, a1, b12, b1, b12, b1)
env[ng] = (a12, a1, a12, a1, b12, b1, b12, b1)
Line 3,779: Line 3,779:
else:
else:
env[ng] = (a12, a1, a12, a1, b12, b1, b12, b1)
env[ng] = (a12, a1, a12, a1, b12, b1, b12, b1)
env[ix] += 1
return
return


Line 3,788: Line 3,787:
else:
else:
term = y[env[iy]]
term = y[env[iy]]
env[iy] += 1
if term is not None:
if term is not None:
new_ng = (a1 + (a12 * term), a12, a + (a2 * term), a2,
new_ng = (a1 + (a12 * term), a12, a + (a2 * term), a2,
Line 3,793: Line 3,793:
if not too_big (new_ng):
if not too_big (new_ng):
env[ng] = new_ng
env[ng] = new_ng
env[iy] += 1
else:
else:
env[ng] = (a12, a12, a2, a2, b12, b12, b2, b2)
env[ng] = (a12, a12, a2, a2, b12, b12, b2, b2)
Line 3,799: Line 3,798:
else:
else:
env[ng] = (a12, a12, a2, a2, b12, b12, b2, b2)
env[ng] = (a12, a12, a2, a2, b12, b12, b2, b2)
env[iy] += 1
return
return