wishlist

[arom]

Hi Mark,

1. could you please add Fortran interfaces for

• magmaf_dgetri
• magmaf_dstedc
• magmaf_zhegvx
• magmaf_zheevx

2. could you please fix zhegvx.cpp and comment/remove Mymagma_ztrmm routines

3. could you please add magma_dstedc function (see attachment)

[mgates3]

These seem reasonable, so hopefully we can get to them for the next release.
-mark

[arom]

Hi Mark!

I need zstedc function, which is absent in Magma.

Is it correct implementation:

Code: Select all

#include "common_magma.h"

extern "C" magma_int_t
magma_zstedc(char range, magma_int_t n, double* d, double* e, cuDoubleComplex* z, magma_int_t ldz,
             cuDoubleComplex* work, magma_int_t lwork, double *rwork, magma_int_t lrwork,
             magma_int_t* iwork, magma_int_t liwork, magma_int_t* info)
{

        double *dwork;

        if (MAGMA_SUCCESS != magma_dmalloc( &dwork, 3*n*(n/2 + 1) )) {
            *info = -15;
            return MAGMA_ERR_DEVICE_ALLOC;
        }

        char range_t = ' ';
        if (range == 'I') range_t = 'A';
        magma_zstedx(range_t, n, 0., 0., 0, 0, d, e, z, ldz,
                     rwork, lrwork, iwork, liwork, dwork, info);

        magma_free( dwork );
}

Please, add it to the next release.

BTW, when will the next release be announced?

Alexey

[brom]

Is magma_zstedx even accelerated? It looks like a complete pass through to the host LAPACK.

[Stan Tomov]

This routine becomes more compute intensive when eigenvectors are needed. In that case most of the flops are in gemm and this is what is GPU accelerated.
Stan