Incorrect display of engineering format

[dwannetts]

Hi,

I'm trying to use Julia, and NumericIO in particular, to refactor some Fortran.

Julia's output is

      10000.00     1   -500.00000E0    0.0000000E0    1.0000000E0   -510.00000E0
      10000.00     2   -400.00000E0    0.0000000E0    1.0000000E0   -410.00000E0
      10000.00     3   -300.00000E0    0.0000000E0    1.0000000E0   -310.00000E0
      10000.00     4   -200.00000E0    0.0000000E0    1.0000000E0   -210.00000E0
      10000.00     5   -100.00000E0    0.0000000E0    1.0000000E0   -110.00000E0
      10000.00     6    0.0000000E0    0.0000000E0    1.0000000E0   -10.000000E0
      10000.00     7    100.00000E0    0.0000000E0    1.0000000E0    90.000000E0
      10000.00     8    200.00000E0    0.0000000E0    1.0000000E0    190.00000E0
      10000.00     9    300.00000E0    0.0000000E0    1.0000000E0    290.00000E0
      10000.00    10    400.00000E0    0.0000000E0    1.0000000E0    390.00000E0
      10000.00    11    500.00000E0    0.0000000E0    1.0000000E0    490.00000E0

from the code fragment

oline =  @sprintf("  %12.2f  %4d", line, tndx)
	    s = formatted(txdp.l.e, :ENG, ndigits=8, charset=:ASCII); oline *= lpad(s, 15, ' ')
	    s = formatted(txdp.l.n, :ENG, ndigits=8, charset=:ASCII); oline *= lpad(s, 15, ' ')
	    s = formatted(txdp.l.z, :ENG, ndigits=8, charset=:ASCII); oline *= lpad(s, 15, ' ')
	    s = formatted(rxdp[jr].l.e, :ENG, ndigits=8, charset=:ASCII); oline *= lpad(s, 15, ' ')
	    s = formatted(rxdp[jr].l.n, :ENG, ndigits=8, charset=:ASCII); oline *= lpad(s, 15, ' ')
	    s = formatted(rxdp[jr].l.z, :ENG, ndigits=8, charset=:ASCII); oline *= lpad(s, 15, ' ')

I believe that if engineering format is working correctly, the output should be

      10000.00     1   -500.00000E+00    0.0000000E+00    1.0000000E+00   -510.00000E+00
      10000.00     2   -400.00000E+00    0.0000000E+00    1.0000000E+00   -410.00000E+00
      10000.00     3   -300.00000E+00    0.0000000E+00    1.0000000E+00   -310.00000E+00
      10000.00     4   -200.00000E+00    0.0000000E+00    1.0000000E+00   -210.00000E+00
      10000.00     5   -100.00000E+00    0.0000000E+00    1.0000000E+00   -110.00000E+00
      10000.00     6      0.00000E+00    0.0000000E+00    1.0000000E+00    -10.00000E+00
      10000.00     7    100.00000E+00    0.0000000E+00    1.0000000E+00     90.00000E+00
      10000.00     8    200.00000E+00    0.0000000E+00    1.0000000E+00    190.00000E+00
      10000.00     9    300.00000E+00    0.0000000E+00    1.0000000E+00    290.00000E+00
      10000.00    10    400.00000E+00    0.0000000E+00    1.0000000E+00    390.00000E+00
      10000.00    11    500.00000E+00    0.0000000E+00    1.0000000E+00    490.00000E+00

ie. decimal points should align down the columns as they do in Fortran.

[ma-laforge]

Hi @dwannetts,

For context:

• This main feature of NumericIO was purposely written to generate outputs with a given number of significant digits.
• It was not designed to generate outputs with a given number of post-decimal digits (thus not expected to match the Fortran output).
• Specifically, the output in :SI or :ENG modes is designed to make numbers easy for humans to read non-column aligned formats. For example, 100n is preferred over 0.1u (decimal point is positioned to ensure there are 1, 2, or 3 digits to the left).

Here is some test code:

using NumericIO
using Printf

#Convenience functions:
SI(v) = formatted(v, :SI, ndigits=3, charset=:ASCII)
ENG(v) = formatted(v, :ENG, ndigits=3, charset=:ASCII)

for i in -15:1:15
	sistr = SI(3.1*(10.0^i))
	engstr = ENG(3.1*(10.0^i))
	@printf("%8s => %s\n", sistr, engstr)
end

And the resultant output:

   3.10f => 3.10E-15
   31.0f => 31.0E-15
    310f => 310E-15
   3.10p => 3.10E-12
   31.0p => 31.0E-12
    310p => 310E-12
   3.10n => 3.10E-9
   31.0n => 31.0E-9
    310n => 310E-9
   3.10u => 3.10E-6
   31.0u => 31.0E-6
    310u => 310E-6
   3.10m => 3.10E-3
   31.0m => 31.0E-3
    310m => 310E-3
    3.10 => 3.10E0
    31.0 => 31.0E0
     310 => 310E0
   3.10k => 3.10E3
   31.0k => 31.0E3
    310k => 310E3
   3.10M => 3.10E6
   31.0M => 31.0E6
    310M => 310E6
   3.10G => 3.10E9
   31.0G => 31.0E9
    310G => 310E9
   3.10T => 3.10E12
   31.0T => 31.0E12
    310T => 310E12
   3.10P => 3.10E15

(You can also use UTF characters instead of specifying :ASCII):

   3.10μ => 3.10×10⁻⁶
   31.0μ => 31.0×10⁻⁶
    310μ => 310×10⁻⁶

That said:
Although this particular feature is working as expected, It would be possible to augment this package with your more functionality.

At this moment however, I don't personally have much time to work on issues at this time. Sorry.

[dwannetts]

Hi @ma-laforge,

Many thanks for your reply -- much appreciated!

It is a pity that the design wasn't expected to match Fortran output, and I do appreciate time constraints, especially for volunteers.

Perhaps the package can be augmented in the future ....

D.