From kristian at buchardt.net  Wed Jun  4 11:53:41 2014
From: kristian at buchardt.net (Kristian Buchardt)
Date: Wed, 04 Jun 2014 11:53:41 +0200
Subject: [Rd] Lapack / dgesvx function declaration
Message-ID: <538EECA5.9070602@buchardt.net>

Hello there

I'm writing some c-code to solve a numerically tough problem for me in 
R. Looking in Lapack.h, i find the following line

F77_NAME(dgesvx)(const int* fact, const char* trans, const int* n,

and I believe that "fact" should've been char instead of int, i.e.

F77_NAME(dgesvx)(const char* fact, const char* trans, const int* n,


My reasoning: In the R-source:
src/modules/lapack/dlapack.f

I find the following lines:

*       SUBROUTINE DGESVX( FACT, TRANS, N, NRHS, A, LDA, AF, LDAF, IPIV,
*                          EQUED, R, C, B, LDB, X, LDX, RCOND, FERR, BERR,
*                          WORK, IWORK, INFO )
*
*       .. Scalar Arguments ..
*       CHARACTER          EQUED, FACT, TRANS
*       INTEGER            INFO, LDA, LDAF, LDB, LDX, N, NRHS
*       DOUBLE PRECISION   RCOND

where FACT is specified to a char.


Might this simply be an error in Lapack.h? And can I just edit my 
Lapack.h and compile?

-- 
Best Regards
Kristian Buchardt

Ph.d.-student
http://math.ku.dk/~buchardt

Department of Mathematical Sciences, University of Copenhagen
Universitetsparken 5, DK-2100 Copenhagen O

PFA Pension
Sundkrogsgade 4, DK-2100 Copenhagen O



[Rd] Lapack / dgesvx function declaration
Kristian Buchardt