[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[Help-gsl] Embedded GSL in Fortran 90
From: |
Eric Germaneau |
Subject: |
[Help-gsl] Embedded GSL in Fortran 90 |
Date: |
Mon, 06 Oct 2008 06:11:08 -0400 |
User-agent: |
Thunderbird 2.0.0.17 (X11/20080914) |
Dear gsl users,
I wish to use gsl in my Fortran 90 program.
So far, I've written 2 C routines as the follwing:
*#include <gsl/gsl_rng.h>
#include <gsl/gsl_randist.h>
#include "gslwr.h"
static gsl_rng* r;
void rng_initialise__(int* s)
{
r = gsl_rng_alloc(gsl_rng_taus);
gsl_rng_set(r, (unsigned long int)(*s));
void rng__(double* zetha)
{
*zetha = (double)gsl_rng_uniform_pos(r) ;
}*
And my program in Fortran 90 looks like that:
*PROGRAM random
IMPLICIT NONE
integer :: seed = 0
double precision :: zetha
call system_clock(seed)
call rng_initialise(seed)
call rng(zetha)
PRINT*, "zetha=",zetha
END PROGRAM random
*
I'm able to creat *.o file doing *icc -c gslwr.c* and *ifc -c gsl.f90*.
But when I do *ifc gsl.o gslw.o* I get
/*gsl.o: In function `MAIN__':
gsl.f90:(.text+0x2c): undefined reference to `rng_initialise_'
gsl.f90:(.text+0x36): undefined reference to `rng_'*/
Does someone understand what's wrong with this?
On the other hand I've heard about fgsl
<http://www.lrz-muenchen.de/services/software/mathematik/gsl/fortran/index.html>.
Dose someone use it?
What do you think about it?
I thank you in advance,
Best,
Eric.
--
/Be the change you wish to see in the world
/ — GANDHI —
Dr. Éric Germaneau <mailto:address@hidden>
Chemical Engineering Department
University of Massachusetts Amherst
159 Gœssmann Laboratory
686 North Pleasant Street
Amherst, MA 01003-9303
USA
Tel: 413-545-6209 (Office)
/msn/: address@hidden
/skype/: aihaike
/ Please consider the environment before printing this email.
Considérez svp l'environnement avant d'imprimer cet email. /