[R-pkg-devel] Problem (with complex.h?) under Microsoft Windows

[Pierre Lafaye de Micheaux]

Dear all,

Thank you very much for your answers. And sorry for the last reply.
Here is how I solved my problem.

I incorporated some files from the src/ directory of the R package 
Faddeeva to
the subdirectory /libraries of the src/ directory of my own package.

This is because I needed to use some functions (in C)  from this package.
Then I used something as:
#include <R.h>
#include "Rmath.h"
#include "libraries/callFaddeeva.cpp"
#include "libraries/Faddeeva.cpp"
#include "libraries/RcppExports.cpp"

It is then possible, to initialise and manipulate complex numbers in C 
in my own code. For example:

cmplx I = C(0.0, 1.0); # For a + ib where a = 0 and b = 1.
cmplx res;
res = Faddeeva::erfi(I, 0.0);

And this compiles well under Windows now.

Note that I first tried to add the Faddeeva package in the Depends field 
of my own package thinking that this might prevent me from having to put 
the Faddeeva source files in my own package but it did not work ...

Best regards,

Pierre L.




Le 13/11/2015 21:36, Duncan Murdoch a écrit :
> On 13/11/2015 11:24 AM, Pierre Lafaye de Micheaux wrote:
>> Dear all,
>>
>> I created a new version of the package IndependenceTests where I call
>> some C and Fortran functions (using the .C() or .Fortran() interface).
>> It works perfectly under Linux, with gcc. No problem with R CMD check. I
>> would like to upload this package on the CRAN.
>>
>> But when I first tried to use winbuilder, I got an error with the
>> following message:
>>
>> * installing *source* package 'IndependenceTests' ...
>> ** libs
>>
>> *** arch - i386
>> g++  -I"D:/RCompile/recent/R/include" -DNDEBUG
>> -I"d:/RCompile/r-compiling/local/local320/include"     -O2 -Wall
>> -mtune=core2 -c Cnhat.cpp -o Cnhat.o
>> Cnhat.cpp:9:89: error: 'complex' has not been declared
>> Cnhat.cpp:9:105: error: two or more data types in declaration of 'res'
>> ....
>> .....
>> etc.
>>
>>
>> I think I understand that there is a problem (only under Windows though)
>> with the complex type.
>> The beginning of my file Cnhat.cpp is as follows:
>>
>> 1. #include <R.h>
>> 2. #include "Rmath.h"
>> 3. #include <complex.h>
>> 4. #include <iostream>
>> 5. using namespace std;
>> 6.
>> 7. extern"C" {
>> 8.
>> 9.   void CnhatC(double *vecs, double *vect, double *X, int *n, int *q,
>> int *p, int *vecd, complex double *res) {
>> 10.
>> 11.    void phinhatReturn(double *vect1, double *vect2, double *vect3,
>> double *X, int *q, int *n, double _Complex *res1, double _Complex *res2,
>> double _Complex *res3);
>> 12.    double _Complex tmp2, prod1 = 1.0 + 0.0*_Complex_I, prod2 = 1.0 +
>> 0.0*_Complex_I, somme = 0.0 + 0.0*_Complex_I;
>> 13.    int indbeginblocl, l, i, *vecdl, k, j;
>> 14.    double *vecsl, *mvectl, *diffvecsvectl, *Xl;
>> 15.    double _Complex *res1, *res2, *res3;
>> 16.    vecdl = new int[1];
>> 17.    res1 =  new _Complex double[1];
>> ....
>> .....
>> etc.
>>
>> As you can see, there seem to be some problem (on line 9.) with: complex
>> double *res
>>
>> Could you please point me to something (e.g., another package) that
>> might help me solve the problem? If possible, I would like to avoid
>> invest too much time using the new .Call() interface because I am quite
>> familiar
>> with the old .C(). And also I am not even sure if it would help solving
>> the problem.
>
> You should be using Rcomplex in the type declaration. (See Writing R 
> Extensions for a discussion of this.)  It is defined as
>
> typedef struct {
>     double r;
>     double i;
> } Rcomplex;
>
> which is likely compatible with your "complex double" on those other 
> systems (which are using newer gcc versions), but it's not guaranteed 
> to be, even there.
>
> Duncan Murdoch

-- 
Pierre Lafaye de Micheaux

Adresse courrier:
Université de Montréal
Pavillon André-Aisenstadt
Département de Mathématiques et Statistique
CP 6128 Succursale Centre-ville
Montréal Qc H3C 3J7
CANADA

Adresse physique:
Département de Mathématiques et Statistique
Bureau 4249, Pavillon André-Aisenstadt
2920, chemin de la Tour
Montréal, Québec H3T 1J4
CANADA

Tél.: (00-1) 514-343-6607 / Fax: (00-1) 514-343-5700
lafaye at dms.umontreal.ca
http://www.biostatisticien.eu