forums.silverfrost.com

The NSWC Fortran subroutine library contains the routine CGAMMA, which appears to be quite suitable as a replacement for your GAMCAR:

https://wp.csiro.au/alanmiller/NSWC/cgamma.f90

Thanks, Paul. The use of '*' as a repeat count for format items was introduced in F2008. There is a description of the feature at https://stevelionel.com/drfortran/2009/07/01/doctor-fortran-in-revert- ...

Given the following program:

program testg0
implicit none
character(*), parameter :: cfmt = '(*(g0,1x))'
integer :: i,j,k
i=2; j=3; k=i*j
print cfmt,i,j ...

The following short program came out of a discussion at Fortran.Discourse.Group. When run, it reproduces the result from the 1966 paper by Lander and Parkin in 1966, which disproved the long-standing ...

Thanks, Paul. Updating to the latest DLLs fixed the problem.

Thanks, Kenneth, and I have no excuse for not being aware of those additions to F2008 in 2025, other than not running into an instance of one of them being beneficial in a Fortran program.

Paul, I took Bill's source code and removed the lines for the module, leaving me a small source file, with invocations of SHIFTR and SHIFTL, which should be available in the current version of FTN95. ...

Paul, the new intrinsics SHIFTL and SHIFTR in Version 9.10 work with /x64, but fail in 32-bit EXEs. The problem is not present in Version 9.14

program tst
implicit none
print *,shiftr& ...

Gfortran, Intel Fortran, Absoft Fortran and NAG Fortran provide TRAILZ as an intrinsic function.

Another application for TRAILZ is in rapid calculation of Collatz sequences. I am curious to know i ...

Thanks, Paul.
The trailz function is also useful in deciding whether a four digit year is a leap year.

integer:: yyyy, ntz
logical :: cent, leap

ntz = trailz(yyy ...

Paul, please consider adding new intrinsic functions that make the X64 instructions LEADZ, TRAILZ and POPCNT available to the Fortran programmer. There are several situations where such a function cal ...

Applying the DOW function to dates from decades ago can be beset by many problems. For instance, consider the Alaska problem:

the purchase date of Friday, October 6 (Julian) was followed once agai ...

See https://www.nottingham.ac.uk/manuscriptsandspecialcollections/researchguidance/datingdocuments/juliangregorian.aspx

where you can even see an image of part of the orders and text with a partial ...

We could also use the MODULO() intrinsic function instead of MOD():
dow = modulo(v,7)+1

This may not suit users of Fortran 77.
I, too, have found myself making errors with code involving signed ...

Here is a little puzzle that can be solved using a DOW() function.

Find years in the range 1901 to 2099 for which:

(i) the weekdays of January 1 and December 31 are the same, and
(ii) Decembe ...