Infinite loop formatting a very (very) small number

Came across an infinite loop in Format.exp_format() when the double value is smaller than 1.0e-309.
This arose in the PCA calculation after duplicating uniref50.fa 5 times to make 480 sequences.
Matrix.printD() calls Format.print(ps, "%15.4e", 1.0E-323);