I get this:
Code: Select all
rumba-jj) ty test.f90
program test
complex*16 x
complex*8 y
complex z
x = (1.0,0.0)
y = (1.0,0.0)
z = (1.0,0.0)
write(*,*) sqrt( x ) , sqrt( y ) , sqrt( z )
end
13-JUN-2023 12:19:10
rumba-jj) f90/noop test.f90
rumba-jj) link test
rumba-jj) run test
(0.000000000000000E+000,2.00000000000000) (1.000000,2.000000)
(1.000000,2.000000)
rumba-jj) f90/vers
VSI Fortran X8.5-0002 (GEM 50X2N) for X86 systems
rumba-jj)
Added in 40 minutes 25 seconds:
I tried to rewrite the sqrt function using the exp function, but that one seems to crash. (try to run the above sample where sqrt is replaced by exp)
Added in 1 hour 5 minutes 49 seconds:
Looks like the sin and cos functions on complex arguments are also giving wrong results or crash.