(pgcd) ( b a y x -- a d )

Vocabulary

math.polynomials.private

Inputs

b	an object
a	an object
y	an object
x	an object

Outputs

a	an object
d	an object

Definition

USING: kernel math.polynomials ;

IN: math.polynomials.private

: (pgcd) ( b a y x -- a d )
    dup V{ 0 } p=
    [ drop nip ] [
        [ nip ] [ p/mod ] 2bi
        [ pick p* swap [ swapd p- ] dip ] dip (pgcd)
    ] if ;