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);