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

Sun Jan 3 09:48:03 CET 2016

Thank you for your comments.

It is indeed the RcppFaddeeva package. I had also tried to put 
RcppFaddeeva on the LinkingTo line with no success.
But I will try to ask on the Rcpp mailing list.

Best regards,
Pierre L.

Le 02/01/2016 21:09, Duncan Murdoch a écrit :
> On 02/01/2016 2:39 PM, Pierre Lafaye de Micheaux wrote:
>> 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.
>
> I don't know of an R package called Faddeeva.  There's an RcppFaddeeva 
> 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.
>
> You shouldn't need to copy the code if you're using RcppFaddeeva, but 
> I can't tell you exactly what you need to do:  the Rcpp mailing list 
> would be a better place to ask.
>
>>
>> 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 ...
>
> Depends is not the way to link at the C++ level.  If this is possible, 
> you'll need to use LinkingTo as well as Imports.  I'm not sure if it's 
> possible, but people on the Rcpp list would know...
>
> Duncan Murdoch
>
>>
>> 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