forums.silverfrost.com

Bill's post of March 4 reminds me of the "Cray Pointers" feature, which is an extension provided by several vendors of Fortran 77. See, for example, http://www.icl.utk.edu/~mucci/MPPopt/addi ...

The expected precision of 32-bit reals is seven significant decimal digits. You are taking the difference between two floating point numbers that differ only in their eighth digits; that difference in ...

Is there a standard fortran statement that will declare all reals as kind=2 (double precision) ?
First of all, note that kind numbers are NON-PORTABLE.

For real variables, the default kind number ...

IDNINT is the specific name for the generic NINT when the argument is double precision. Either name produces the same inline code when FTN95 is used to compile for X87.

The argument, if an express ...

The program
program repcnt
implicit none
logical :: b1,b2,b3,b4,b5,b6
character(len=9) :: line = ' 4*T 2*F '
!
read(line,*)b1,b2,b3,b4,b5,b6
print '(1x,6L3& ...

Optimisation should only touch the efficiency of execution not accuracy of the results surely.
Only if exact arithmetic is used (i.e., integer, character and boolean variables) and no integer overfl ...

I noticed the name "HST2D" in your post. It so happens that I worked with the successor USGS code HST3D 2.2.16 in 2016. I found several bugs (at least 15 of them) in the source code, using F ...

Thanks for your comments, John.

I agree that we are working with a corner case. However, these findings struck me as strange:

(i) X87 arithmetic (as used by programs compiled with FTN77, FTN95- ...

FTN95 (32-bit, many versions up to 8.3) was giving me a result off by 1.0 when evaluating NINT(double precision expression). The numbers are such that each term in the expression, as well as the sum o ...

And if you can't find an old Windows 95 machine or VM, you can do the installation of a 32-bit package that uses a 16-bit installer on a Linux PC with Wine installed. Once the installation is done, co ...

In many cases, it is not necessary to hook a profiler to a library function entry point to count the number of function calls, especially if all the calls occur in source code that you have and can ch ...

Eddie, before working on an EINS catalogue, you will probably need to have an Enigma machine (from eBay?) or an Enigma simulator, just to test out some ideas.

If you relax the requirement from &qu ...

Paul's work-around succeeds if the Fortran source file is saved as UTF-8 without BOM (byte order mark).

I tried with UTF-8 + BOM and with UTF-16 LE, and the compiler rejected these versions.

Si ...

well that worked. I didn't try that (from earlier posts) because I saw I didn't have the c:\windows\sysnative folder and so I didn't bother trying it.

Out of curiosity, what does it stand for?
Th ...

Fortran compilers lag behind C compilers in regard to their support for multi-byte and Unicode characters. Without such support (i.e., character types with more than just one KIND), what you can do in ...