(in-package :division)


#+nil
(let ((f (string->poly "x^2" '(x y)))
      ;;(y-sq (string->poly "y^2" '(x y)))
      (fl (cdr (string->poly "[x+y,x-2*y]" '(x y))))
      ;;(quotients (cdr (string->poly "[x-y,0]" '(x y))))
      )
  (multiple-value-bind (f c division-count)
      (step (normal-form f fl))
    (format t "Quotient:        ~S~%" (->sexp f))
    (format t "Coefficient:     ~S~%" c)
    (format t "Division count:  ~S~%" division-count)))

(let ((f (string->poly "x^2" '(x y)))
      ;;(y-sq (string->poly "y^2" '(x y)))
      (fl (cdr (string->poly "[x+y,x-2*y]" '(x y))))
      ;;(quotients (cdr (string->poly "[x-y,0]" '(x y))))
      )
  (multiple-value-bind (a r c division-count)
      (step (poly-pseudo-divide f fl))
    (format t "Quotients:        ~S~%" a)
    (format t "Remainder:        ~S~%" (->sexp r))
    (format t "Coefficient:     ~S~%" c)
    (format t "Division count:  ~S~%" division-count)))