default(realprecision,38);
default(echo,1);
test(D) =
{
  for (d = D, D+10^3,
    if (!isfundamental(d),next);
    print(d, " ", quadclassunit(d).cyc)
  );
}
setrand(1); test(10^15); test(-10^15)

\\ #1195
setrand(11);quadclassunit(-8419588).cyc
setrand(2);quadclassunit(-1459008).cyc
setrand(7);quadclassunit(-3799812).cyc
setrand(1); quadclassunit(-13163208).cyc

\\ #1195 with non-fundamental discriminants [oo loop]
setrand(38);quadclassunit(-29920).cyc

quadclassunit(-13163208,,[0.1]).cyc
setrand(1); quadclassunit((2^70+25)).cyc
setrand(1); quadclassunit(8*3*5*7).cyc

\\ #1700
setrand(1);quadclassunit(-612556842419).cyc

\\ #2015
setrand(2); quadclassunit(-699,,[6,6]).cyc
