PGI User Forum

[stegan]

I'm using a fortran code (f90) for MonteCarlo calculations. During the execution the program uses a lot of times the intrinsic function "anint(x)" to get the integer value (returned as real) nearer to x.
I tried to use the same program, written in C, and the performance is better than fortran (about 30%). After some tests, I realized that the function "anint" in fortran is slower than the correspond "rint" in C; in fact if I don't use these functions in two programs, the fortran code is faster.
So, there is a solution to use the "rint" C function in the fortran code?
Thank's a lot.

[mkcolg]

Chapter 10 of the PGI user's guide gives a good introduction to inter-language calling conventions. Basically, you need to create a C wrapper function which calls rint. Note the extra underscore after the function name and the use of pointer arguements. Also note that rint is a macro that complains if its compiled with pgcc, so you'll need to compile this file with gcc.

Example:

myrint.c

[stegan]

Thank's for help!
I tried to do this, but the needed time for execution is longer.
In fact by using ANINT (in fortran), the required time is about 33% less than by using RINT (in C).
Probably the problem is another!
I ask you another question: why by declaring DOUBLE the myrint function is very and very slower than with FLOAT?
Thank's a lot.
Stefano

[mkcolg]

Off hand I don't know why the double version is so much slower. When you compile the with "gcc -S" and view the generated assembly file the float version calls "rintf" and the double version calls "rint". My guess is that is simply takes more instructions to convert a double to and int than it does a float.

Have you tried compiling with a profiling flag and using the PGI profiler pgprof? This tool can really help in determining performance bottlenecks. A complete description of pgrof can be found in the PGI tools guide. http://www.pgroup.com/resources/docs.htm

- Mat