1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21
|
#include <stdio.h> /* I/O lib ISOC */
#include <stdlib.h> /* Standard Lib ISOC */
#include <mkl.h>
#include "mkl_scalapack.h"
int main()
{
double A[9]={0, 1,2,3,4,5,6,7,8};
double B[9];
int m=3, n=2, d=1;
char *uplo = "All";
pdlacpy(uplo, &m, &n, A, 0, 0, &n, B, 0, &d, &n);
for(i=0; i<3; i++){
for(j=0; j<3; j++){
printf("%lf \t", B[i*3+j]);
}
printf("\n");
}
|