From dmbates at gmail.com  Tue Jul  5 15:58:25 2005
From: dmbates at gmail.com (Douglas Bates)
Date: Tue, 5 Jul 2005 08:58:25 -0500
Subject: [Rd] calling fortran functions CHOL and DPOTRF form Fortran
In-Reply-To: <200507051312.48997.gilles.guillot@inapg.inra.fr>
References: <200507051312.48997.gilles.guillot@inapg.inra.fr>
Message-ID: <40e66e0b05070506582d4702ec@mail.gmail.com>

On 7/5/05, Gilles GUILLOT <gilles.guillot at inapg.inra.fr> wrote:
> Hi all,
> 
> I'm working out some Fortran code for which
> I want to compute the Choleski decomposition of a covariance matrix
> in Fortran.
> 
> I tried to do it by two methods :
> 
> 1) Calling the lapack function DPOTRF.
> I can see the source code and check that my call is correct,
> but it does not compile with:
> system("R CMD SHLIB ~/main.f")
> dyn.load("~/main.so")
> 
> I get:
> Error in dyn.load(x, as.logical(local), as.logical(now)) :
>         unable to load shared library
> [...]
>  undefined symbol: dpotrf_
> 
> Could anyone point out how to compile with such a call to lapack?

Make a small package from your source code then follow the instructions in

http://tolstoy.newcastle.edu.au/R/devel/05/06/1382.html

> 
> 2) Calling the Fortran function CHOL
> The previous combination of R CMD SHLIB  and dyn.load
> compiles fine but  I would like to see the source code of this function
> to check if my call is correct. I've not been able to find it in  the R source
> distribution.
> Could anyone point the place of that file?

${RSRC}/src/appl/chol.f

> 
> Thanks in advance,
> 
> Gilles
> 
> 
> --
> _________________________________________________________________________
> Gilles GUILLOT
> 
> INRA -D?partement Math?matiques et Informatique Appliqu?es
> Unit? Mixte de Recherche INRA - INAPG - ENGREF
> 
> Institut National Agronomique de Paris-Grignon
> 16, rue Claude Bernard
> 75231 Paris cedex 5 France
> 
> phone 33 1 44 08 18 42
> fax   33 1 44 08 16 66
> http://www.inapg.fr/ens_rech/mathinfo/personnel/guillot/welcome.html
> 
> ______________________________________________
> R-devel at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-devel
>



[Rd] calling fortran functions CHOL and DPOTRF form Fortran
Douglas Bates