  Computation of the extended gcd of two quadratic integers. The ring of
integers considered is principal but could be euclidean or not euclidean ring.
This method rely on principal ideal ring and reduction of binary quadratic
forms.
