PCL; is there a string formatting which can produce 8 digits Nastran exponentials?

Question

Arthur Hickson · Accepted Answer

There is no inbuilt function that does exactly what you want.

The translator codes for writing the nastran input files would be written in Fortran or C etc and will be implementing this in their native code. It is a seemingly simple but in practice quite involved exercise to get the most accuracy into a Nastran 8 digit real string from a real number of arbitrary magnitude. The challenge is to change the format conversion depending on the magnitude and also to account for rounding.

The following function will give you an idea how you might do this in PCL though performance may be poor.

The example does not extract the maximum accuracy for small numbers like -0.098765 which will come out as -0.09877 when the most accurate for NASTRAN is -.098765 , to get this accuracy requires doing the real to string conversion with another digit of accuracy and then trimming the leading 0 from the resultant string, but leaving any negative sign in place.

Also remember that for small numbers with negative exponential data it is possible to change format like "-0.1234E-5" to "-.1234-5" which gains 2 more digits compared with truncating the number as "-0.12E-5", again this example code does not do this.

Hopefully the example will help but working with the wide field 16 digits is of course much easier and more accurate though people do like the 8 character wide fields despite their lack of significant digits.

FUNCTION r2nas(invalue)

REAL value

REAL absval

REAL invalue

STRING sval[8]

INTEGER sdig

INTEGER digbefore

INTEGER digaft

value = invalue

absval=abs(value)

digaft=1

sdig = 0

IF (value < 0.) then

sdig =-1

END IF

IF ( absval < .99999995e-9 ) then

sdig += 3

digbefore = 1

ELSE IF ( absval < .0099999995) then

sdig += 4

digaft=2

digbefore = 1

ELSE IF ( absval < .099999995) then

sdig += 5

digaft = 2

digbefore = 1

ELSE IF ( absval < 0.99999995) then

sdig += 6

digbefore = 1

ELSE IF ( absval < 999999.95) then

sdig += 7

digbefore = sdig - 1

ELSE IF ( absval < 9999999500.) then

sdig += 5

digbefore = sdig - 1

ELSE

sdig += 4

digbefore = sdig - 1

END IF

str_from_real_signif_digit(value, sdig, digbefore, digaft, 0, sval)

RETURN sval

END FUNCTION