HPC Programming on Intel Many-Integrated-Core Hardware with MAGMA Port to Xeon Phi

Algorithm 2

Cholesky factorization in MIC MAGMA.

(1) for do

(2)   = min(, );

(3)  magma_zherk_mic(MagmaUpper, MagmaConjTrans, , , _one, dA, ldda, one, dA, ldda, queue);

(4)  magma_zgetmatrix_async_mic(, , dA, ldda, work, , , queue, &event);

(5)  if then

(6)     magma_zgemm_mic(MagmaConjTrans, MagmaNoTrans,

       , , , _one, , ldd, , ldda, queue);

(7)  magma_event_sync_mic(event);

(8)  lapackf77_zpotrf(MagmaUpperStr, &, work, &, info);

(9)  if info then

(10)    info ;

(11)  magma_zsetmatrix_async_mic(, , work, , , dA(), ldda, queue, &event);

(12)  if then

(13)   magma_event_sync_mic(event);

(14)   magma_ztrsm_mic(MagmaLeft, MagmaUpper, MagmaConjTrans, MagmaNonUnit,

    , , _one, dA(), ldda, dA(, ), ldda, queue);