Of course I could now try to analyze the Pari sources, but maybe you know a link to a high level description of the algorithm used in quadunit?
(Similar questions could be asked about the algorithms used in other Pari functions: I didn't see any such documentation on the Pari website. Maybe I didn't dig deep enough, or just overlooked it?)
Regards, Michael