forums.silverfrost.com

Kenneth_Smith · Posted: Fri May 16, 2025 8:58 am Post subject: Fortran 2023 trig functions

The fortran 2023 standard includes new trig related functions, namely cosd, sind, tand, acosd, asind, atand, atan2d, acospi, asinpi, atanpi, atan2pi, cospi, sinpi, and tanpi.

**D functions which use degrees rather than radians were implemented in gFortran and Ifort as extensions long before the 2023 standard, so they do appear fairly regularly in code I see written by others.

**PI functions operate on multiples of half revolutions. For example:: COSPI(x) = COS(PI*X), so COSPI(1.0) = -1.0 These functions have obvious applications in signal processing etc.

Perhaps these new intrinsics should be added to the FTN95 wish list, and added to the library when time permits? This would keep the core mathematical intrinsics in FTN95 in synch with the latest fortran standard.

The functions should be easy to implement - no gotchas like asin(Z) Very Happy

Further details can be found here:

https://wg5-fortran.org/N2201-N2250/N2212.pdf

PaulLaidler · Posted: Sat May 17, 2025 6:46 am Post subject:

Many thanks for the feedback. I have added this to the wish list.

PaulLaidler · Posted: Thu Jun 05, 2025 8:24 am Post subject:

These functions have now been implemented for the next release of FTN95 and the associated DLLs.

wahorger · Joined: 13 Oct 2014 Posts: 1257 Location: Morrison, CO, USA

Paul, this will be a nice addition. I have some complex equations that will be simplified (to the eye, at least) using these new features.

Thanks for all you do!

PaulLaidler · Posted: Thu Jun 05, 2025 5:04 pm Post subject:

Bill

If "complex" means passing and returning complex numbers then this addition will not help.

As I understand it, the Standard does not include complex numbers in this context. Please let me know if I am wrong.

wahorger · Joined: 13 Oct 2014 Posts: 1257 Location: Morrison, CO, USA

Oh, sorry! I "misspoke". I meant that there are a number of calculations that would be better done using *D or *PI functions, making the reading of them simpler, versus complex. Mostly involving spherical trig and datum conversions...

I can certainly see using that term would cause some concern!

Kenneth_Smith · Posted: Fri Jun 06, 2025 8:23 am Post subject:

Good work Paul, I appreciate this addition to the library.

Just to clarify, the aforementioned functions are indeed intended to operate on real and not complex variables.

Kenneth_Smith · Posted: Wed Jun 25, 2025 12:22 pm Post subject:

Paul, there's something not correct with the *d functions. The following code should return two columns of data with equal entries in each row.

PaulLaidler · Posted: Wed Jun 25, 2025 2:02 pm Post subject:

Thanks Ken. I will take a look at this.

PaulLaidler · Posted: Thu Jun 26, 2025 9:15 am Post subject:

This has now been fixed for the next release of FTN95.

Kenneth_Smith · Posted: Thu Jun 26, 2025 2:03 pm Post subject:

Thanks Paul,

Did you also look at the *pi functions? I see today that they have a similar issue.

PaulLaidler · Posted: Thu Jun 26, 2025 3:07 pm Post subject:

Yes and also a failure when used in an initialisation (PARAMETER) statement.

LitusSaxonicum · Posted: Tue Jul 01, 2025 7:44 pm Post subject:

I am fascinated by the introduction of new functions like this. What did people do before? What did they name their self-developed functions? Is there any conflict - and do the system functions perform better (faster? more accurate?)?

PaulLaidler · Posted: Tue Jul 01, 2025 8:59 pm Post subject:

As I understand it, in general there is no reason to expect faster or more accurate results. Any differences will normally be in the round-off error.

But there are some interesting possibilities. For example cosd(60.0) could be exactly 0.5 with no round-off error, whilst tand(90.0) could print as INFINITY.