From pgilbert@bank-banque-canada.ca  Thu Aug 28 17:01:05 1997
From: pgilbert@bank-banque-canada.ca (Paul Gilbert)
Date: Thu, 28 Aug 1997 12:01:05 -0400
Subject: R-alpha: pow_ii
Message-ID: <97Aug28.120241edt.26501@mailgate.bank-banque-canada.ca>


>> I'm a bit confused about this discussion, but since I don't have f2c on
...

>A lot of Fortran code doesn't require any libf2c routines.  When routines
>are needed the most common ones are the exponentiation routines pow_dd,
>pow_ii and pow_di (pow_ri is for single precision reals and so is less
>needed in R, which uses double precision).  However, *some* fortran code,
...

Thomas

>From your previous response to my confusion I guessed I would be ok
with pow_ii, but that does not seem to be the case:

R
...
>  dsecurve.function.tests(verbose=T) 
Random number generator tests ...ok
DSE curvature test 1 ... ld.so.1: /R/bin/R.binary: fatal: relocation error: symbol not found: pow_ii: referenced in /dse/dsefor.so
Killed

[4] /R : nm bin/R.binary | grep pow
[146]	|    244556|     312|FUNC |LOCL |0    |7      |complex_pow
[2557]	|   1462800|       0|FUNC |GLOB |0    |UNDEF  |pow
[2795]	|   1053360|      68|FUNC |GLOB |0    |7      |pow_dd
[2800]	|   1053124|     236|FUNC |GLOB |0    |7      |pow_di
[377]	|   1474036|       4|OBJT |LOCL |0    |15     |powerSymbol
[5] /R : 
[5] /R : R


This leads me to wonder how R calculates powers and exp for complex
numbers. I presume that if I can arrange to call the same routines then
I can avoid the extra burden of having to include parts of libf2c with
my distributed code. Do you know what R uses and can I call these from
fortran?

Paul Gilbert
=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-
r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-devel-request@stat.math.ethz.ch
=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-


R-alpha: pow_ii
Paul Gilbert