Continued fraction/Arithmetic/G(matrix ng, continued fraction n1, continued fraction n2): Difference between revisions
Content added Content deleted
m (→{{header|D}}) |
|||
Line 3,768: | Line 3,768: | ||
else: |
else: |
||
term = x[env[ix]] |
term = x[env[ix]] |
||
⚫ | |||
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) |
||
⚫ | |||
return |
return |
||
Line 3,788: | Line 3,787: | ||
else: |
else: |
||
term = y[env[iy]] |
term = y[env[iy]] |
||
⚫ | |||
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) |
||
⚫ | |||
return |
return |
||