HOW TO CALL THE IMSL FORTRAN LIBRARIES FROM C by ...

More documents

Recommendations

Info

$Intel(R) Math Kernel Library for Linux* OS User's Guide$

3.4 Example /* This example solves a linear system of 100 by 100 using the IMSL routine LSLRG. */ #include #include #include #include #define MAXN 100 main() { int long int long int double double i, nerr; n=MAXN, nn=MAXN*MAXN, lda=MAXN; ipath=2, seed=123457; ermax=l.e-7; a[MAXN][MAXN], b[MAXN], x[MAXN], ans[MAXN]; /* generate matrix a with uniform random numbers */ RNSET (&seed); RNUN (&nn, a); /* generate a pre-determined answer, ans */ for(i=0;i
* print the result */ if (nerr) printf ("IMSL: %2 differences occurred.\n", nerr); else printf ("IMSL: LSLRG is correct with ermax = %g.\n", ermax); } The output from CRAY is: IMSL: LSLRG is correct with ermax = 1e-07. 18
Page 1 and 2: HOW TO CALL THE IMSL FORTRAN LIBRAR
Page 3 and 4: How to Call the IMSL Fortran Librar
Page 5 and 6: 1. Syntax Differences Between C and
Page 7 and 8: double a[2] [3] = {11., 12., 13., 2
Page 9 and 10: 2. Unix Platforms The topics discus
Page 11 and 12: * * myfunc calculates and returns *
Page 13 and 14: #include #define MAXN 100 main ()
Page 15 and 16: 3. Cray The Cray UNICOS environment
Page 17: main () { double exp(), fabs(), myf
Page 21 and 22: 2x3 matrix 1 2 3 1 11 12 13 2 21 22
Page 23 and 24: drnun (&nn, a); /* generate a pre-d
Page 25 and 26: 5.1 Passing a Character String Argu
Page 27 and 28: Exact answer = 8.38906 IMSL dqdng =
Page 29 and 30: if (nerr) printf ("IMSL: %2 differe
Page 31 and 32: The user should compile C and Fortr

HOW TO CALL THE IMSL FORTRAN LIBRARIES FROM C by ...

Create successful ePaper yourself

Delete template?

Save as template?