The Meson Build System
http://mesonbuild.com/
You can not select more than 25 topics
Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
34 lines
895 B
34 lines
895 B
#include <stdio.h> |
|
|
|
// #include <mkl.h> |
|
// #include <mkl_scalapack.h> |
|
// #include <mkl_blacs.h> |
|
|
|
extern float pslamch_(const int *, const char *); |
|
extern void blacs_pinfo_(int *, int *); |
|
extern void blacs_get_(const int *, const int *, int *); |
|
extern void blacs_gridinit_(int *, const char *, const int *, const int *); |
|
extern void blacs_gridinfo_(const int *, int *, int *, int *, int *); |
|
extern void blacs_gridexit_(const int *); |
|
extern void blacs_exit_(const int *); |
|
|
|
int main(void){ |
|
|
|
int myid, nprocs, ictxt, mycol, myrow, npcol=2, nprow=2; |
|
const int i0=0, i1=1, in1=-1; |
|
|
|
blacs_pinfo_(&myid, &nprocs); |
|
blacs_get_(&in1, &i0, &ictxt); |
|
blacs_gridinit_(&ictxt, "C", &nprocs, &i1); |
|
|
|
blacs_gridinfo_(&ictxt, &nprow, &npcol, &myrow, &mycol); |
|
|
|
float eps = pslamch_(&ictxt, "E"); |
|
|
|
if (myrow == mycol) printf("OK: Scalapack C: eps= %f\n", eps); |
|
|
|
blacs_gridexit_(&ictxt); |
|
blacs_exit_(&i0); |
|
|
|
return 0; |
|
} |