Actual source code: example3.edp
 1 /*******************************************************
 2  *  This file is part of the FreeFEM tutorial made by  *
 3  *      Pierre Jolivet <pierre@joliv.et>               *
 4  *                                                     *
 5  *  See https://joliv.et/FreeFem-tutorial for more     *
 6  *                                                     *
 7  *   Description: multithreaded BLAS                   *
 8  *******************************************************/
  
10 load "lapack"
 
12 complex[int, int] A = [[-4*1i,-1+2*1i,3+1i], [2,11*1i,3], [0,-1+2*1i,2], [4-2*1i,11+4*1i,2*1i]];
13 complex[int, int] B = [[3,-1], [1,2], [3,-1]];
 
15 complex[int, int] C = A * B;
16 complex[int, int] Ct = C';
17 complex[int, int] D = C * Ct;
18 complex[int, int] vectors(D.n, D.n);
19 complex[int] values(D.n);
20 zgeev(D, values, vectors);