FORTRAN INTERFACE

use elpa
class(elpa_t), pointer :: elpa

call elpa%skew_eigenvectors (a, ev, q, error)

With the definitions of the input and output variables:

class(elpa_t) :: elpa ! returns an instance of the ELPA object

datatype :: a

The matrix a for which the eigenvalues should be computed. The dimensions of matrix a must be set BEFORE with the methods elpa_set(3) and elpa_setup(3). The datatype of the matrix can be one of "real(kind=c_double)", "real(kind=c_float)", "complex(kind=c_double)", or "complex(kind=c_float)". The matrix has to be skew-symmetric, this is not checked by the routine.

datatype :: ev

The vector ev where the eigenvalues will be stored in ascending order. The datatype of the vector ev can be either "real(kind=c_double)", or "real(kind=c_float)", depending of the datatype of the matrix. Note that complex hermitian matrices also have real valued eigenvalues. datatype :: q The storage space for the computed eigenvectors. The dimensions of matrix a must be set BEFORE with the methods elpa_set(3) and elpa_setup(3). The datatype of the matrix can be one of "real(kind=c_double)", "real(kind=c_float)", "complex(kind=c_double)", or "complex(kind=c_float)". Note, that for a skew-symmetric matrix the eigenvectors are complex. The routines returns separately the real and imaginary parts of the complex eigenvectors. Thus, the storage space has to be of dimension q(#numer_of_rows,2*#number_of_column).

integer, optional :: error

The return error code of the function. Should be "ELPA_OK". The error code can be querried with the function elpa_strerr(3)

C INTERFACE

#include <elpa/elpa.h>
elpa_t handle;

void elpa_eigenvalues(elpa_t handle, datatype *a, datatype *ev, datatype *q, int *error);

With the definitions of the input and output variables:

elpa_t handle;

The handle to the ELPA object

datatype *a;

The matrix a for which the eigenvalues should be computed. The dimensions of the matrix must be set BEFORE with the methods elpa_set(3) and elpa_setup(3). The datatype can be one of "double", "float", "double complex", or "float complex". The matrix has to be symmetric or hermitian, this is not checked by the routine.

datatype *ev;

The storage for the computed eigenvalues. Eigenvalues will be stored in ascendig order. The datatype can be either "double" or "float". Note that the eigenvalues of complex hermitian matrices are also real.

datatype *q;

The storage space for the computed eigenvectors. The dimensions of the matrix must be set BEFORE with the methods elpa_set(3) and elpa_setup(3). The datatype can be one of "double", "float", "double complex", or "float complex". Note, that for a skew-symmetric matrix the eigenvectors are complex. The routines returns separately the real and imaginary parts of the complex eigenvectors. Thus, the storage space has to be of dimension q(#numer_of_rows,2*#number_of_column).

int *error;

The error code of the function. Should be "ELPA_OK". The error codes can be querried with elpa_strerr(3)