Commit 4d02b3ec authored by rtoy's avatar rtoy
Browse files

f2cl can't handle equivalences of different data types. Replace the

implementation with a new Fortran version that doesn't use
equivalences.
parent de7fc70e
double precision function dcabs1(z)
double complex z,zz
double precision t(2)
equivalence (zz,t(1))
zz = z
dcabs1 = dabs(t(1)) + dabs(t(2))
C ORIGINAL:
c double complex z,zz
c double precision t(2)
c equivalence (zz,t(1))
c zz = z
c dcabs1 = dabs(t(1)) + dabs(t(2))
c NEW
double complex z
dcabs1 = dabs(dble(z)) + dabs(dimag(z))
return
end
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment