Overloading the REAL intrinsic

I was reading about Fortran’s intrinsic function REAL recently and discovered that it behaved inconsistently, so I fixed it by overloading the function in a module that I generally use everywhere.

I use REAL to extract the real part of a complex number. The KIND of the return value is equal to that of it’s argument, but only if the argument is COMPLEX. If the argument itself is of type REAL, then it returns a value using the compiler’s default REAL KIND. For me, since I program in double precision, this can cause me to accidentally lose precision, if I’m not careful.

I see three ways to fix this:

  1. Always include the KIND when calling REAL, for example REAL(foo,8). I’m likely to forget, and it clutters my code without adding information to the reader.
  2. Change the compiler default KIND. In gfortran, I just need a command-line argument: gfortran -fdefault-real-8. However, this option can change from compiler to compiler.
  3. Overload the intrinsic function to return a REAL value with the same KIND as the argument.

So, I went with the last option. To get this to work, I put the function in my module prec_def, which is where I set the precision for the entire code. I used an interface for the overloading, so that I only overrode the intrinsic for the argument names and types that I specify.

Here’s the result:

module prec_def
use iso_fortran_env, only : REAL64 !to get the KIND parameter for 64-bit precision
implicit none

integer, parameter :: long=REAL64

interface real
module procedure real_real
end interface

contains

real(Long) function real_real(num)
real(Long), intent(in) :: num

!call the intrinsic function, specifying KIND
real_real = real(num,Long)

end function real_real

end module prec_def

Let me know what you think!

2 thoughts on “Overloading the REAL intrinsic

  1. Pingback: [redacted]

    • Not sure if you are a person, but you were asking an offtopic question and linking to your site. I call that spam. Since you sound like a person, I’ll edit your comment to remove advertising, but answer your question. In the future, if the comment is offtopic, please direct it to my email address.

      So, one offtopic answer: You can do either one. Modern blog software is easily installable to your own webserver, or you can get free or paid hosting online. It’s mostly WYSIWYG at that point, unless you want to do some custom tweaks to the site.

Leave a Reply to Brent W. Barker Cancel reply

Your email address will not be published. Required fields are marked *

     

You may use these HTML tags and attributes: <a href="" title=""> <abbr title=""> <acronym title=""> <b> <blockquote cite=""> <cite> <code> <del datetime=""> <em> <i> <q cite=""> <strike> <strong>