Continued fraction/Arithmetic/G(matrix ng, continued fraction n1, continued fraction n2): Difference between revisions
Content added Content deleted
Line 4,223: | Line 4,223: | ||
block |
block |
||
! |
|||
! Let us start with a test of the long division routine, with some |
! Let us start with a test of the long division routine, with some |
||
! numbers known to trigger a bug in the first version I |
! numbers known to trigger a bug in the first version I |
||
! posted. That version had a buggy "add_back" routine. |
! posted. That version had a buggy "add_back" routine. |
||
! |
|||
! (How I found such numbers is easy: I used random search.) |
|||
! |
|||
type(big_integer), allocatable :: a, b, q, r |
type(big_integer), allocatable :: a, b, q, r |