Re: [Rd] Calling FORTRAN function from R issue?

From: Berend Hasselman <bhh_at_xs4all.nl>
Date: Tue, 06 Mar 2012 15:07:39 +0100

On 06-03-2012, at 14:37, Berwin A Turlach wrote:

> G'day Berend,
>
> On Tue, 6 Mar 2012 13:06:34 +0100
> Berend Hasselman <bhh_at_xs4all.nl> wrote:
>
> [... big snip ...]
>

>> But I would really like to hear from an Rexpert why you
>> shouldn't/can't use external here in the Fortran.

>
> Probably less a question for an Rexpert but for a Fortran expert....
>
> If you insert "implicit none" (one of my favourite extensions that I
> always use) as the first statement after
> subroutine callzdotc(retval,n, zx, incx, zy, incy)
> you will see what is going on. The compiler should refuse compilation
> and complain that the type of zdotc was not declared (at least my
> compiler does). For FORTRAN to know that zdotc returns a double
> complex you need the
> double complex zdotc
> declaration in callzdotc.
>

I usually use -fimplicit-none or a similar option as argument for a fortran compiler. I didn't use it in a syntax checking pre compiling run.

> An
> external double complex zdotc
> would be necessary if you were to call another subroutine/function, say
> foo, that accepts functions as formal arguments and you want to pass
> zdotc via such an argument. Something like
>
> subroutine callzdotc(....)
> ....
> external double complex zdotc
> ....
> call foo(a, b, zdotc)
> ....
> return
> end
>
> subroutine(a, b, h)
> double complex h, t
> ....
> t = h(..,..,..,..)
> ....
> return
> end
>
> In C, the qualifier (or whatever the technical term is) "external" is
> used to indicate that the function/variable/symbol is defined in
> another compilation unit. In FORTRAN77, "external" is used to tell the
> compiler that you are passing a function to another
> function/subroutine. At least that is my understanding from what I
> re-read in my FORTRAN documentation.

Not quite true.
It is required to use external if you are passing the name of a function/subroutine to another routine. Otherwise it is just a statement telling the compile that the <external-name> is external to the function/subroutine/.. being compiled.

See http://www.fortran.com/F77_std/rjcnf0001-sh-8.html#sh-8.7

>
> Thus, perhaps strangely, if there is only a
> external double complex zdotc
> declaration in your subroutine, the compiler doesn't know that a call
> ....

If I am reading the above reference correctly, the syntax

        external double complex zdotc

is simply incorrect (correct syntax uses commas to separate names e.g. A,B,C). The compiler (gfortran) simply seems to ignore the "double" and the "complex" if you don't use implicit none.

Berend



R-devel_at_r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel Received on Tue 06 Mar 2012 - 14:09:10 GMT

This quarter's messages: by month, or sorted: [ by date ] [ by thread ] [ by subject ] [ by author ]

All messages

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Tue 06 Mar 2012 - 14:50:24 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-devel. Please read the posting guide before posting to the list.

list of date sections of archive