abacus-develop/scalapack__connector_8h_source.html

#ifndef SCALAPACK_CONNECTOR_H

#define SCALAPACK_CONNECTOR_H


#include <limits>

// =========================================================

// Tolerances for LAPACK/ScaLAPACK eigenvalue routines

// =========================================================

// Huawei KML strictly validates input parameters.

// It rejects abstol=0 and orfac=-1, which are standard

// defaults in open-source ScaLAPACK to trigger internal logic.

// We must explicitly pass the mathematically equivalent defaults.

#ifdef __KML

    constexpr double SCALAPACK_ABSTOL = 2*std::numeric_limits<double>::min();          // 2*PDLAMCH('S')

    constexpr double SCALAPACK_ORFAC  = 2.0e-12;                                       // Default value

#else

    constexpr double SCALAPACK_ABSTOL = 0.0;

    constexpr double SCALAPACK_ORFAC  = -1.0;

#endif

// =========================================================


#ifdef __MPI


#include <complex>


extern "C"

{

    int numroc_( const int *n, const int *nb, const int *iproc, const int *srcproc, const int *nprocs );

    void descinit_(

        int *desc,

        const int *m, const int *n, const int *mb, const int *nb, const int *irsrc, const int *icsrc,

        const int *ictxt, const int *lld, int *info);


    void pddot_(int* n, double* dot, double* x, int* ix, int* jx, int* descx, int* incx,

        double* y, int* iy, int* jy, int* descy, int* incy);

    void pzdotc_(int* n, std::complex<double>* dot, std::complex<double>* x, int* ix, int* jx, int* descx, int* incx,

        std::complex<double>* y, int* iy, int* jy, int* descy, int* incy);


    void pdpotrf_(char *uplo, int *n, double *a, int *ia, int *ja, int *desca, int *info);

//  void pzpotrf_(char *uplo, int *n, double _Complex *a, int *ia, int *ja, int *desca, int *info);

    void pzpotrf_(char *uplo, int *n, std::complex<double> *a, int *ia, int *ja, int *desca, int *info);


    void pdtran_(const int* m, const int* n,

        const double* alpha, const double* a, const int* ia, const int* ja, const  int* desca,

        const double* beta, double* c, const int* ic, const int* jc, const int* descc);


    void pztranu_(const int *m,const int*n,

        const std::complex<double>* alpha, const std::complex<double>* a, const int* ia, const int* ja, const int* desca,

        const std::complex<double> *beta ,  std::complex<double> *c , const int *ic ,const int *jc ,const int *descc);


    void pzgemv_(

        const char *transa,

        const int *M, const int *N,

        const double *alpha,

        const std::complex<double> *A, const int *IA, const int *JA, const int *DESCA,

        const std::complex<double> *B, const int *IB, const int *JB, const int *DESCB, const int *K,

        const double *beta, std::complex<double> *C, const int *IC, const int *JC, const int *DESCC,const int *L);

    void pdgemv_(

        const char *transa,

        const int *M, const int *N,

        const double *alpha,

        const double *A, const int *IA, const int *JA, const int *DESCA,

        const double *B, const int *IB, const int *JB, const int *DESCB, const int *K,

        const double *beta, double *C, const int *IC, const int *JC, const int *DESCC,const int *L);

    // C = a * A.? * B.? + b * C

    void pdgemm_(

        const char *transa, const char *transb,

        const int *M, const int *N, const int *K,

        const double *alpha,

        const double *A, const int *IA, const int *JA, const int *DESCA,

        const double *B, const int *IB, const int *JB, const int *DESCB,

        const double *beta,

        double *C, const int *IC, const int *JC, const int *DESCC);

    void pzgemm_(

        const char *transa, const char *transb,

        const int *M, const int *N, const int *K,

        const std::complex<double> *alpha,

        const std::complex<double> *A, const int *IA, const int *JA, const int *DESCA,

        const std::complex<double> *B, const int *IB, const int *JB, const int *DESCB,

        const std::complex<double> *beta,

        std::complex<double> *C, const int *IC, const int *JC, const int *DESCC);

    void pdsymm_(char *side , char *uplo , int *m , int *n ,

        double *alpha , double *a , int *ia , int *ja , int *desca ,

        double *b , int *ib , int *jb , int *descb ,

        double *beta ,  double *c , int *ic , int *jc , int *descc );

    void pdtrmm_(char *side , char *uplo , char *transa , char *diag , int *m , int *n ,

        double *alpha , double *a , int *ia , int *ja , int *desca ,

        double *b , int *ib , int *jb , int *descb );

//  void pztrmm_(char *side , char *uplo , char *transa , char *diag , int *m , int *n ,

//      double *alpha , double _Complex *a , int *ia , int *ja , int *desca ,

//      double _Complex *b , int *ib , int *jb , int *descb );

    void pztrmm_(char *side , char *uplo , char *transa , char *diag , int *m , int *n ,

        std::complex<double> *alpha , std::complex<double> *a , int *ia , int *ja , int *desca ,

        std::complex<double> *b , int *ib , int *jb , int *descb );

    void pzhemm_(char* side , char* uplo , int* m , int* n ,

        std::complex<double>* alpha , std::complex<double>* a , int* ia , int* ja , int* desca ,

        std::complex<double>* b , int* ib , int* jb , int* descb ,

        std::complex<double>* beta ,  std::complex<double>* c , int* ic , int* jc , int* descc );

    void pzgetrf_(

        const int *M, const int *N,

        std::complex<double> *A, const int *IA, const int *JA, const int *DESCA,

        int *ipiv,  int *info);


    void pzgesv_(

            const int *n, const int *nrhs,

            const std::complex<double> *A, const int *ia, const int *ja, const int *desca,

            int *ipiv, std::complex<double>* B, const int* ib, const int* jb, const int*descb, const int *info

        );


    void pdsygvx_(const int* itype, const char* jobz, const char* range, const char* uplo,

        const int* n, double* A, const int* ia, const int* ja, const int*desca, double* B, const int* ib, const int* jb, const int*descb,

        const double* vl, const double* vu, const int* il, const int* iu,

        const double* abstol, int* m, int* nz, double* w, const double*orfac, double* Z, const int* iz, const int* jz, const int*descz,

        double* work, int* lwork, int*iwork, int*liwork, int* ifail, int*iclustr, double*gap, int* info);


    void pzhegvx_(const int* itype, const char* jobz, const char* range, const char* uplo,

        const int* n, std::complex<double>* A, const int* ia, const int* ja, const int*desca, std::complex<double>* B, const int* ib, const int* jb, const int*descb,

        const double* vl, const double* vu, const int* il, const int* iu,

        const double* abstol, int* m, int* nz, double* w, const double*orfac, std::complex<double>* Z, const int* iz, const int* jz, const int*descz,

        std::complex<double>* work, int* lwork, double* rwork, int* lrwork, int*iwork, int*liwork, int* ifail, int*iclustr, double*gap, int* info);


    void pssygvx_(const int* itype, const char* jobz, const char* range, const char* uplo,

        const int* n, float* A, const int* ia, const int* ja, const int*desca, float* B, const int* ib, const int* jb, const int*descb,

        const float* vl, const float* vu, const int* il, const int* iu,

        const float* abstol, int* m, int* nz, float* w, const float*orfac, float* Z, const int* iz, const int* jz, const int*descz,

        float* work, int* lwork, int*iwork, int*liwork, int* ifail, int*iclustr, float*gap, int* info);


    void pchegvx_(const int* itype, const char* jobz, const char* range, const char* uplo,

        const int* n, std::complex<float>* A, const int* ia, const int* ja, const int*desca, std::complex<float>* B, const int* ib, const int* jb, const int*descb,

        const float* vl, const float* vu, const int* il, const int* iu,

        const float* abstol, int* m, int* nz, float* w, const float*orfac, std::complex<float>* Z, const int* iz, const int* jz, const int*descz,

        std::complex<float>* work, int* lwork, float* rwork, int* lrwork, int*iwork, int*liwork, int* ifail, int*iclustr, float*gap, int* info);


    void pzgetri_(

        const int *n,

        const std::complex<double> *A, const int *ia, const int *ja, const int *desca,

        int *ipiv, const std::complex<double> *work, const int *lwork, const int *iwork, const int *liwork, const int *info);


    void pzgeadd_(

        const char *transa,

        const int *m, const int *n,

        const std::complex<double> *alpha,

        const std::complex<double> *a, const int *ia, const int *ja, const int *desca,

                        const std::complex<double> *beta,

        const std::complex<double> *c, const int *ic, const int *jc, const int *descc);


    void pztranc_(

        const int *M, const int *N,

        const std::complex<double> *alpha,

        const std::complex<double> *A, const int *IA, const int *JA, const int *DESCA,

        const std::complex<double> *beta,

        std::complex<double> *C, const int *IC, const int *JC, const int *DESCC);


    void pdgemr2d_(const int *M, const int *N,

        double *A, const int *IA, const int *JA, const int *DESCA,

        double *B, const int *IB, const int *JB, const int *DESCB,

        const int *ICTXT);


    void pzgemr2d_(const int *M, const int *N,

        std::complex<double> *A, const int *IA, const int *JA, const int *DESCA,

        std::complex<double> *B, const int *IB, const int *JB, const int *DESCB,

        const int *ICTXT);


    // Scalapack wrappers to copy 2D blocks of data

    // more info:

    // https://netlib.org/scalapack/explore-html/da/db5/pigemr_8c.html

    // https://netlib.org/scalapack/explore-html/dd/dcd/pdgemr_8c.html

    // https://netlib.org/scalapack/explore-html/d5/dd4/pzgemr_8c.html

    // https://netlib.org/scalapack/explore-html/d5/deb/psgemr_8c.html

    // https://netlib.org/scalapack/explore-html/d4/dad/pcgemr_8c.html

    void Cpigemr2d (int m, int n, int *ptrmyblock, int ia, int ja, int *ma, int *ptrmynewblock, int ib, int jb, int *mb, int globcontext);

    void Cpdgemr2d (int m, int n, double *ptrmyblock, int ia, int ja, int *ma, double *ptrmynewblock, int ib, int jb, int *mb, int globcontext);

    void Cpzgemr2d (int m, int n, std::complex<double> *ptrmyblock, int ia, int ja, int *ma, std::complex<double> *ptrmynewblock, int ib, int jb, int *mb, int globcontext);

    void Cpsgemr2d (int m, int n, float *ptrmyblock, int ia, int ja, int *ma, float *ptrmynewblock, int ib, int jb, int *mb, int globcontext);

    void Cpcgemr2d (int m, int n, std::complex<float> *ptrmyblock, int ia, int ja, int *ma, std::complex<float> *ptrmynewblock, int ib, int jb, int *mb, int globcontext);

}


    template <typename T>


    struct block2d_data_type

    {

        constexpr static bool value = std::is_same<T, double>::value || std::is_same<T, std::complex<double>>::value || std::is_same<T, float>::value || std::is_same<T, std::complex<float>>::value || std::is_same<T, int>::value;

    };


    template <typename T>


    typename std::enable_if<block2d_data_type<T>::value,void>::type Cpxgemr2d(int M, int N, T *A, int IA, int JA, int *DESCA, T *B, int IB, int JB, int *DESCB, int ICTXT)

    {

        if (std::is_same<T,double>::value) Cpdgemr2d(M, N, reinterpret_cast<double*>(A),IA, JA, DESCA,reinterpret_cast<double*>(B),IB,JB, DESCB,ICTXT);

        if (std::is_same<T,std::complex<double>>::value) Cpzgemr2d(M, N, reinterpret_cast<std::complex<double>*>(A),IA, JA, DESCA,reinterpret_cast<std::complex<double>*>(B),IB,JB, DESCB,ICTXT);

        if (std::is_same<T,float>::value) Cpsgemr2d(M, N, reinterpret_cast<float*>(A),IA, JA, DESCA,reinterpret_cast<float*>(B),IB,JB, DESCB,ICTXT);

        if (std::is_same<T,std::complex<float>>::value) Cpcgemr2d(M, N, reinterpret_cast<std::complex<float>*>(A),IA, JA, DESCA,reinterpret_cast<std::complex<float>*>(B),IB,JB, DESCB,ICTXT);

        if (std::is_same<T,int>::value) Cpigemr2d(M, N, reinterpret_cast<int*>(A),IA, JA, DESCA,reinterpret_cast<int*>(B),IB,JB, DESCB,ICTXT);

    };


class ScalapackConnector

{

public:

    static inline


    void geadd(

        const char transa,

        const int m, const int n,

        const std::complex<double> alpha,

        const std::complex<double> *a, const int ia, const int ja, const int *desca,

        const std::complex<double> beta,

        const std::complex<double> *c, const int ic, const int jc, const int *descc)

    {

        pzgeadd_(&transa, &m, &n, &alpha, a, &ia, &ja, desca, &beta, c, &ic, &jc, descc);

    }


    static inline


    void dot(int n,

        double& dot,

        double* a,

        int ia,

        int ja,

        int inca,

        double* b,

        int ib,

        int jb,

        int incb,

        int* desc)

    {

        pddot_(&n, &dot, a, &ia, &ja, desc, &inca, b, &ib, &jb, desc, &incb);

    }


    static inline


    void dot(int n,

        std::complex<double>& dotc,

        std::complex<double>* a,

        int ia,

        int ja,

        int inca,

        std::complex<double>* b,

        int ib,

        int jb,

        int incb,

        int* desc)

    {

        pzdotc_(&n, &dotc, a, &ia, &ja, desc, &inca, b, &ib, &jb, desc, &incb);

    }


    static inline


        void gemm(

            const char transa, const char transb,

            const int M, const int N, const int K,

            const double alpha,

            const double* A, const int IA, const int JA, const int* DESCA,

            const double* B, const int IB, const int JB, const int* DESCB,

            const double beta,

            double* C, const int IC, const int JC, const int* DESCC)

    {

        pdgemm_(&transa, &transb, &M, &N, &K, &alpha, A, &IA, &JA, DESCA,

            B, &IB, &JB, DESCB, &beta, C, &IC, &JC, DESCC);

    }


    static inline


    void gemm(

        const char transa, const char transb,

        const int M, const int N, const int K,

        const std::complex<double> alpha,

        const std::complex<double> *A, const int IA, const int JA, const int *DESCA,

        const std::complex<double> *B, const int IB, const int JB, const int *DESCB,

        const std::complex<double> beta,

        std::complex<double> *C, const int IC, const int JC, const int *DESCC)

    {

        pzgemm_(&transa, &transb, &M, &N, &K, &alpha, A, &IA, &JA, DESCA,

            B, &IB, &JB, DESCB, &beta, C, &IC, &JC, DESCC);

    }


    static inline


    void gemm(char transa, char transb, int M, int N, int K,

        double alpha,

        double* A,

        double* B,

        double beta,

        double* C,

        int* DESC)

    {

        int isrc = 1;

        pdgemm_(&transa,

                &transb,

                &M,

                &N,

                &K,

                &alpha,

                A,

                &isrc,

                &isrc,

                DESC,

                B,

                &isrc,

                &isrc,

                DESC,

                &beta,

                C,

                &isrc,

                &isrc,

                DESC);

    }


    static inline


    void gemm(char transa, char transb, int M, int N, int K,

        std::complex<double> alpha,

        std::complex<double>* A,

        std::complex<double>* B,

        std::complex<double> beta,

        std::complex<double>* C,

        int* DESC)

    {


        int isrc = 1;

        pzgemm_(&transa,

            &transb,

            &M,

            &N,

            &K,

            &alpha,

            A,

            &isrc,

            &isrc,

            DESC,

            B,

            &isrc,

            &isrc,

            DESC,

            &beta,

            C,

            &isrc,

            &isrc,

            DESC);

    }


    static inline


    void symm(char side,

        char uplo,

        int m,

        int n,

        double alpha,

        double* a,

        double* b,

        double beta,

        double* c,

        int* desc)

    {

        int isrc = 1;

        pdsymm_(&side, &uplo, &m, &n, &alpha, a, &isrc, &isrc, desc, b, &isrc, &isrc, desc, &beta, c, &isrc, &isrc, desc);

    }


    static inline


    void getrf(

        const int M, const int N,

        std::complex<double> *A, const int IA, const int JA, const int *DESCA,

        int *ipiv,  int *info) //fix a bug: info is output and we must use int*

    {

        pzgetrf_(&M, &N, A, &IA, &JA, DESCA, ipiv, info);

    }


    static inline


    void getri(

        const int n,

        const std::complex<double> *A, const int ia, const int ja, const int *desca, int *ipiv,

        const std::complex<double> *work, const int *lwork, const int *iwork, const int *liwork, int *info)

    {

        pzgetri_(&n, A, &ia, &ja, desca, ipiv, work, lwork, iwork, liwork, info);

    }


    static inline


    void gesv(

        const int n, const int nrhs,

        const std::complex<double> *A, const int ia, const int ja, const int *desca,

        int *ipiv, std::complex<double>* B, const int ib, const int jb, const int*descb, int *info)

    {

        pzgesv_(&n, &nrhs, A, &ia, &ja, desca, ipiv, B, &ib, &jb, descb, info);

    }


    static inline


    void tranu(

        const int m, const int n,

        const std::complex<double> alpha , std::complex<double> *a , const int ia , const int ja , const int *desca,

        const std::complex<double> beta ,  std::complex<double> *c , const int ic , const int jc , const int *descc)

    {

        pztranu_(&m, &n, &alpha, a, &ia, &ja, desca, &beta, c, &ic, &jc, descc);

    }


    static inline


    int potrf(char uplo, int na, double* U, int* desc)

    {

        int isrc = 1;

        int info;

        pdpotrf_(&uplo, &na, U, &isrc, &isrc, desc, &info);

        return info;

    }


    static inline


    int potrf(char uplo, int na, std::complex<double>* U, int* desc)

    {

        int isrc = 1;

        int info;

        pzpotrf_(&uplo, &na, U, &isrc, &isrc, desc, &info);

        return info;

    }


    static inline


    void trmm(char side,

        char uplo,

        char trans,

        char diag,

        int m,

        int n,

        double alpha,

        double* a,

        double* b,

        int* desc)

    {

        int isrc = 1;

        pdtrmm_(&side, &uplo, &trans, &diag, &m, &n, &alpha, a, &isrc, &isrc, desc, b, &isrc, &isrc, desc);

    }


    static inline


    void trmm(char side,

        char uplo,

        char trans,

        char diag,

        int m,

        int n,

        std::complex<double> alpha,

        std::complex<double>* a,

        std::complex<double>* b,

        int* desc)

    {

        int isrc = 1;

        pztrmm_(&side, &uplo, &trans, &diag, &m, &n, &alpha, a, &isrc, &isrc, desc, b, &isrc, &isrc, desc);

    }


    static inline


    void hemm(char side,

        char uplo,

        int na,

        std::complex<double> alpha,

        std::complex<double>* a,

        std::complex<double>* b,

        std::complex<double> beta,

        std::complex<double>* c,

        int* desc)

    {

        int isrc = 1;

        pzhemm_(&side,

        &uplo,

        &na,

        &na,

        &alpha,

        a,

        &isrc,

        &isrc,

        desc,

        b,

        &isrc,

        &isrc,

        desc,

        &beta,

        c,

        &isrc,

        &isrc,

        desc);

    }


};


#endif // __MPI


#endif

ScalapackConnector
Definition scalapack_connector.h:214

ScalapackConnector::gemm
static void gemm(char transa, char transb, int M, int N, int K, double alpha, double *A, double *B, double beta, double *C, int *DESC)
Definition scalapack_connector.h:289

ScalapackConnector::getri
static void getri(const int n, const std::complex< double > *A, const int ia, const int ja, const int *desca, int *ipiv, const std::complex< double > *work, const int *lwork, const int *iwork, const int *liwork, int *info)
Definition scalapack_connector.h:377

ScalapackConnector::gemm
static void gemm(const char transa, const char transb, const int M, const int N, const int K, const std::complex< double > alpha, const std::complex< double > *A, const int IA, const int JA, const int *DESCA, const std::complex< double > *B, const int IB, const int JB, const int *DESCB, const std::complex< double > beta, std::complex< double > *C, const int IC, const int JC, const int *DESCC)
Definition scalapack_connector.h:275

ScalapackConnector::trmm
static void trmm(char side, char uplo, char trans, char diag, int m, int n, std::complex< double > alpha, std::complex< double > *a, std::complex< double > *b, int *desc)
Definition scalapack_connector.h:438

ScalapackConnector::symm
static void symm(char side, char uplo, int m, int n, double alpha, double *a, double *b, double beta, double *c, int *desc)
Definition scalapack_connector.h:352

ScalapackConnector::potrf
static int potrf(char uplo, int na, double *U, int *desc)
Definition scalapack_connector.h:404

ScalapackConnector::gemm
static void gemm(const char transa, const char transb, const int M, const int N, const int K, const double alpha, const double *A, const int IA, const int JA, const int *DESCA, const double *B, const int IB, const int JB, const int *DESCB, const double beta, double *C, const int IC, const int JC, const int *DESCC)
Definition scalapack_connector.h:261

ScalapackConnector::geadd
static void geadd(const char transa, const int m, const int n, const std::complex< double > alpha, const std::complex< double > *a, const int ia, const int ja, const int *desca, const std::complex< double > beta, const std::complex< double > *c, const int ic, const int jc, const int *descc)
Definition scalapack_connector.h:217

ScalapackConnector::dot
static void dot(int n, std::complex< double > &dotc, std::complex< double > *a, int ia, int ja, int inca, std::complex< double > *b, int ib, int jb, int incb, int *desc)
Definition scalapack_connector.h:245

ScalapackConnector::gemm
static void gemm(char transa, char transb, int M, int N, int K, std::complex< double > alpha, std::complex< double > *A, std::complex< double > *B, std::complex< double > beta, std::complex< double > *C, int *DESC)
Definition scalapack_connector.h:320

ScalapackConnector::tranu
static void tranu(const int m, const int n, const std::complex< double > alpha, std::complex< double > *a, const int ia, const int ja, const int *desca, const std::complex< double > beta, std::complex< double > *c, const int ic, const int jc, const int *descc)
Definition scalapack_connector.h:395

ScalapackConnector::trmm
static void trmm(char side, char uplo, char trans, char diag, int m, int n, double alpha, double *a, double *b, int *desc)
Definition scalapack_connector.h:422

ScalapackConnector::potrf
static int potrf(char uplo, int na, std::complex< double > *U, int *desc)
Definition scalapack_connector.h:413

ScalapackConnector::dot
static void dot(int n, double &dot, double *a, int ia, int ja, int inca, double *b, int ib, int jb, int incb, int *desc)
Definition scalapack_connector.h:229

ScalapackConnector::gesv
static void gesv(const int n, const int nrhs, const std::complex< double > *A, const int ia, const int ja, const int *desca, int *ipiv, std::complex< double > *B, const int ib, const int jb, const int *descb, int *info)
Definition scalapack_connector.h:386

ScalapackConnector::hemm
static void hemm(char side, char uplo, int na, std::complex< double > alpha, std::complex< double > *a, std::complex< double > *b, std::complex< double > beta, std::complex< double > *c, int *desc)
Definition scalapack_connector.h:454

ScalapackConnector::getrf
static void getrf(const int M, const int N, std::complex< double > *A, const int IA, const int JA, const int *DESCA, int *ipiv, int *info)
Definition scalapack_connector.h:368

N
#define N
Definition exp.cpp:24

T
#define T
Definition exp.cpp:237

Cpigemr2d
void Cpigemr2d(int m, int n, int *ptrmyblock, int ia, int ja, int *ma, int *ptrmynewblock, int ib, int jb, int *mb, int globcontext)

pztrmm_
void pztrmm_(char *side, char *uplo, char *transa, char *diag, int *m, int *n, std::complex< double > *alpha, std::complex< double > *a, int *ia, int *ja, int *desca, std::complex< double > *b, int *ib, int *jb, int *descb)

SCALAPACK_ABSTOL
constexpr double SCALAPACK_ABSTOL
Definition scalapack_connector.h:16

Cpzgemr2d
void Cpzgemr2d(int m, int n, std::complex< double > *ptrmyblock, int ia, int ja, int *ma, std::complex< double > *ptrmynewblock, int ib, int jb, int *mb, int globcontext)

pchegvx_
void pchegvx_(const int *itype, const char *jobz, const char *range, const char *uplo, const int *n, std::complex< float > *A, const int *ia, const int *ja, const int *desca, std::complex< float > *B, const int *ib, const int *jb, const int *descb, const float *vl, const float *vu, const int *il, const int *iu, const float *abstol, int *m, int *nz, float *w, const float *orfac, std::complex< float > *Z, const int *iz, const int *jz, const int *descz, std::complex< float > *work, int *lwork, float *rwork, int *lrwork, int *iwork, int *liwork, int *ifail, int *iclustr, float *gap, int *info)

pzgetri_
void pzgetri_(const int *n, const std::complex< double > *A, const int *ia, const int *ja, const int *desca, int *ipiv, const std::complex< double > *work, const int *lwork, const int *iwork, const int *liwork, const int *info)

pzgemr2d_
void pzgemr2d_(const int *M, const int *N, std::complex< double > *A, const int *IA, const int *JA, const int *DESCA, std::complex< double > *B, const int *IB, const int *JB, const int *DESCB, const int *ICTXT)

Cpcgemr2d
void Cpcgemr2d(int m, int n, std::complex< float > *ptrmyblock, int ia, int ja, int *ma, std::complex< float > *ptrmynewblock, int ib, int jb, int *mb, int globcontext)

pdpotrf_
void pdpotrf_(char *uplo, int *n, double *a, int *ia, int *ja, int *desca, int *info)

numroc_
int numroc_(const int *n, const int *nb, const int *iproc, const int *srcproc, const int *nprocs)

Cpdgemr2d
void Cpdgemr2d(int m, int n, double *ptrmyblock, int ia, int ja, int *ma, double *ptrmynewblock, int ib, int jb, int *mb, int globcontext)

pdgemv_
void pdgemv_(const char *transa, const int *M, const int *N, const double *alpha, const double *A, const int *IA, const int *JA, const int *DESCA, const double *B, const int *IB, const int *JB, const int *DESCB, const int *K, const double *beta, double *C, const int *IC, const int *JC, const int *DESCC, const int *L)

pztranc_
void pztranc_(const int *M, const int *N, const std::complex< double > *alpha, const std::complex< double > *A, const int *IA, const int *JA, const int *DESCA, const std::complex< double > *beta, std::complex< double > *C, const int *IC, const int *JC, const int *DESCC)

pdsygvx_
void pdsygvx_(const int *itype, const char *jobz, const char *range, const char *uplo, const int *n, double *A, const int *ia, const int *ja, const int *desca, double *B, const int *ib, const int *jb, const int *descb, const double *vl, const double *vu, const int *il, const int *iu, const double *abstol, int *m, int *nz, double *w, const double *orfac, double *Z, const int *iz, const int *jz, const int *descz, double *work, int *lwork, int *iwork, int *liwork, int *ifail, int *iclustr, double *gap, int *info)

pzhegvx_
void pzhegvx_(const int *itype, const char *jobz, const char *range, const char *uplo, const int *n, std::complex< double > *A, const int *ia, const int *ja, const int *desca, std::complex< double > *B, const int *ib, const int *jb, const int *descb, const double *vl, const double *vu, const int *il, const int *iu, const double *abstol, int *m, int *nz, double *w, const double *orfac, std::complex< double > *Z, const int *iz, const int *jz, const int *descz, std::complex< double > *work, int *lwork, double *rwork, int *lrwork, int *iwork, int *liwork, int *ifail, int *iclustr, double *gap, int *info)

pdtran_
void pdtran_(const int *m, const int *n, const double *alpha, const double *a, const int *ia, const int *ja, const int *desca, const double *beta, double *c, const int *ic, const int *jc, const int *descc)

pzpotrf_
void pzpotrf_(char *uplo, int *n, std::complex< double > *a, int *ia, int *ja, int *desca, int *info)

pddot_
void pddot_(int *n, double *dot, double *x, int *ix, int *jx, int *descx, int *incx, double *y, int *iy, int *jy, int *descy, int *incy)

pzdotc_
void pzdotc_(int *n, std::complex< double > *dot, std::complex< double > *x, int *ix, int *jx, int *descx, int *incx, std::complex< double > *y, int *iy, int *jy, int *descy, int *incy)

pdgemm_
void pdgemm_(const char *transa, const char *transb, const int *M, const int *N, const int *K, const double *alpha, const double *A, const int *IA, const int *JA, const int *DESCA, const double *B, const int *IB, const int *JB, const int *DESCB, const double *beta, double *C, const int *IC, const int *JC, const int *DESCC)

pzgemm_
void pzgemm_(const char *transa, const char *transb, const int *M, const int *N, const int *K, const std::complex< double > *alpha, const std::complex< double > *A, const int *IA, const int *JA, const int *DESCA, const std::complex< double > *B, const int *IB, const int *JB, const int *DESCB, const std::complex< double > *beta, std::complex< double > *C, const int *IC, const int *JC, const int *DESCC)

pzgemv_
void pzgemv_(const char *transa, const int *M, const int *N, const double *alpha, const std::complex< double > *A, const int *IA, const int *JA, const int *DESCA, const std::complex< double > *B, const int *IB, const int *JB, const int *DESCB, const int *K, const double *beta, std::complex< double > *C, const int *IC, const int *JC, const int *DESCC, const int *L)

pdtrmm_
void pdtrmm_(char *side, char *uplo, char *transa, char *diag, int *m, int *n, double *alpha, double *a, int *ia, int *ja, int *desca, double *b, int *ib, int *jb, int *descb)

descinit_
void descinit_(int *desc, const int *m, const int *n, const int *mb, const int *nb, const int *irsrc, const int *icsrc, const int *ictxt, const int *lld, int *info)

pzgesv_
void pzgesv_(const int *n, const int *nrhs, const std::complex< double > *A, const int *ia, const int *ja, const int *desca, int *ipiv, std::complex< double > *B, const int *ib, const int *jb, const int *descb, const int *info)

Cpxgemr2d
std::enable_if< block2d_data_type< T >::value, void >::type Cpxgemr2d(int M, int N, T *A, int IA, int JA, int *DESCA, T *B, int IB, int JB, int *DESCB, int ICTXT)
Definition scalapack_connector.h:203

pztranu_
void pztranu_(const int *m, const int *n, const std::complex< double > *alpha, const std::complex< double > *a, const int *ia, const int *ja, const int *desca, const std::complex< double > *beta, std::complex< double > *c, const int *ic, const int *jc, const int *descc)

pdgemr2d_
void pdgemr2d_(const int *M, const int *N, double *A, const int *IA, const int *JA, const int *DESCA, double *B, const int *IB, const int *JB, const int *DESCB, const int *ICTXT)

pdsymm_
void pdsymm_(char *side, char *uplo, int *m, int *n, double *alpha, double *a, int *ia, int *ja, int *desca, double *b, int *ib, int *jb, int *descb, double *beta, double *c, int *ic, int *jc, int *descc)

pzgeadd_
void pzgeadd_(const char *transa, const int *m, const int *n, const std::complex< double > *alpha, const std::complex< double > *a, const int *ia, const int *ja, const int *desca, const std::complex< double > *beta, const std::complex< double > *c, const int *ic, const int *jc, const int *descc)

pssygvx_
void pssygvx_(const int *itype, const char *jobz, const char *range, const char *uplo, const int *n, float *A, const int *ia, const int *ja, const int *desca, float *B, const int *ib, const int *jb, const int *descb, const float *vl, const float *vu, const int *il, const int *iu, const float *abstol, int *m, int *nz, float *w, const float *orfac, float *Z, const int *iz, const int *jz, const int *descz, float *work, int *lwork, int *iwork, int *liwork, int *ifail, int *iclustr, float *gap, int *info)

SCALAPACK_ORFAC
constexpr double SCALAPACK_ORFAC
Definition scalapack_connector.h:17

Cpsgemr2d
void Cpsgemr2d(int m, int n, float *ptrmyblock, int ia, int ja, int *ma, float *ptrmynewblock, int ib, int jb, int *mb, int globcontext)

pzhemm_
void pzhemm_(char *side, char *uplo, int *m, int *n, std::complex< double > *alpha, std::complex< double > *a, int *ia, int *ja, int *desca, std::complex< double > *b, int *ib, int *jb, int *descb, std::complex< double > *beta, std::complex< double > *c, int *ic, int *jc, int *descc)

pzgetrf_
void pzgetrf_(const int *M, const int *N, std::complex< double > *A, const int *IA, const int *JA, const int *DESCA, int *ipiv, int *info)

block2d_data_type
Definition scalapack_connector.h:180

block2d_data_type::value
static constexpr bool value
Definition scalapack_connector.h:181

ictxt
int ictxt
Definition tddft_test.cpp:14