.. SPDX-FileCopyrightText: 2019-2020 Intel Corporation .. .. SPDX-License-Identifier: CC-BY-4.0 .. _onemkl_lapack_hegvd: hegvd ===== Computes all eigenvalues and, optionally, eigenvectors of a real generalized symmetric definite eigenproblem using a divide and conquer method. .. container:: section .. rubric:: Description ``hegvd`` supports the following precisions. .. list-table:: :header-rows: 1 * - T * - ``std::complex`` * - ``std::complex`` The routine computes all the eigenvalues, and optionally, the eigenvectors of a complex generalized Hermitian positive-definite eigenproblem, of the form :math:`Ax = \lambda Bx, ABx = \lambda x`, or :math:`BAx =\lambda x`. Here :math:`A` and :math:`B` are assumed to be Hermitian and :math:`B` is also positive definite. It uses a divide and conquer algorithm. hegvd (Buffer Version) ---------------------- .. container:: section .. rubric:: Syntax .. code-block:: cpp namespace oneapi::mkl::lapack { void hegvd(cl::sycl::queue &queue, std::int64_t itype, oneapi::mkl::job jobz, oneapi::mkl::uplo upper_lower, std::int64_t n, cl::sycl::buffer &a, std::int64_t lda, cl::sycl::buffer &b, std::int64_t ldb, cl::sycl::buffer &w, cl::sycl::buffer &scratchpad, std::int64_t scratchpad_size) } .. container:: section .. rubric:: Input Parameters queue The queue where the routine should be executed. itype Must be 1 or 2 or 3. Specifies the problem type to be solved: if :math:`\text{itype} = 1`, the problem type is :math:`Ax = \lambda Bx;` if :math:`\text{itype} = 2`, the problem type is :math:`ABx = \lambda x;` if :math:`\text{itype} = 3`, the problem type is :math:`BAx = \lambda x`. jobz Must be ``job::novec`` or ``job::vec``. If ``jobz = job::novec``, then only eigenvalues are computed. If ``jobz = job::vec``, then eigenvalues and eigenvectors are computed. upper_lower Must be ``uplo::upper`` or ``uplo::lower``. If ``upper_lower = uplo::upper``, ``a`` and ``b`` store the upper triangular part of :math:`A` and :math:`B`. If ``upper_lower = uplo::lower``, ``a`` and ``b`` stores the lower triangular part of :math:`A` and :math:`B`. n The order of the matrices :math:`A` and :math:`B` (:math:`0 \le n`). a Buffer, size ``a(lda,*)`` contains the upper or lower triangle of the Hermitian matrix :math:`A`, as specified by upper_lower. The second dimension of ``a`` must be at least :math:`\max(1, n)`. lda The leading dimension of ``a``; at least :math:`\max(1,n)`. b Buffer, size ``b(ldb,*)`` contains the upper or lower triangle of the Hermitian matrix :math:`B`, as specified by upper_lower. The second dimension of ``b`` must be at least :math:`\max(1, n)`. ldb The leading dimension of ``b``; at least :math:`\max(1,n)`. scratchpad_size Size of scratchpad memory as a number of floating point elements of type ``T``. Size should not be less than the value returned by :ref:`onemkl_lapack_hegvd_scratchpad_size` function. .. container:: section .. rubric:: Output Parameters a On exit, if ``jobz = job::vec``, then if :math:`\text{info} = 0`, ``a`` contains the matrix :math:`Z` of eigenvectors. The eigenvectors are normalized as follows: if :math:`\text{itype} = 1` or :math:`\text{itype} = 2`, :math:`Z^{H}BZ = I`; if :math:`\text{itype} = 3`, :math:`Z^{H}B^{-1}Z = I`; If ``jobz = job::novec``, then on exit the upper triangle (if ``upper_lower = uplo::upper``) or the lower triangle (if ``upper_lower = uplo::lower``) of :math:`A`, including the diagonal, is destroyed. b On exit, if :math:`\text{info} \le n`, the part of ``b`` containing the matrix is overwritten by the triangular factor :math:`U` or :math:`L` from the Cholesky factorization :math:`B = U^{H}U`\ or :math:`B = LL^{H}`. w Buffer, size at least :math:`n`. If :math:`\text{info} = 0`, contains the eigenvalues of the matrix :math:`A` in ascending order. scratchpad Buffer holding scratchpad memory to be used by routine for storing intermediate results. .. container:: section .. rubric:: Throws This routine shall throw the following exceptions if the associated condition is detected. An implementation may throw additional implementation-specific exception(s) in case of error conditions not covered here. :ref:`oneapi::mkl::host_bad_alloc` :ref:`oneapi::mkl::device_bad_alloc` :ref:`oneapi::mkl::unimplemented` :ref:`oneapi::mkl::unsupported_device` :ref:`oneapi::mkl::lapack::invalid_argument` :ref:`oneapi::mkl::lapack::computation_error` Exception is thrown in case of problems during calculations. The ``info`` code of the problem can be obtained by `info()` method of exception object: If :math:`\text{info}=-i`, the :math:`i`-th parameter had an illegal value. For :math:`\text{info} \le n`: If :math:`\text{info}=i`, and ``jobz = oneapi::mkl::job::novec``, then the algorithm failed to converge; :math:`i` indicates the number of off-diagonal elements of an intermediate tridiagonal form which did not converge to zero; If :math:`\text{info}=i`, and ``jobz = oneapi::mkl::job::vec``, then the algorithm failed to compute an eigenvalue while working on the submatrix lying in rows and columns :math:`\text{info}/(n+1)`` through :math:`\text{mod}(\text{info}, n+1)`. For :math:`\text{info}>n`: If :math:`\text{info}=n+i`, for :math:`1 \le i \le n`, then the leading minor of order :math:`i` of :math:`B` is not positive-definite. The factorization of :math:`B` could not be completed and no eigenvalues or eigenvectors were computed. If ``info`` equals to value passed as scratchpad size, and `detail()` returns non zero, then passed scratchpad is of insufficient size, and required size should not be less than value return by `detail()` method of exception object. hegvd (USM Version) ---------------------- .. container:: section .. rubric:: Syntax .. code-block:: cpp namespace oneapi::mkl::lapack { cl::sycl::event hegvd(cl::sycl::queue &queue, std::int64_t itype, oneapi::mkl::job jobz, oneapi::mkl::uplo upper_lower, std::int64_t n, T *a, std::int64_t lda, T *b, std::int64_t ldb, RealT *w, T *scratchpad, std::int64_t scratchpad_size, const std::vector &events = {}) } .. container:: section .. rubric:: Input Parameters queue The queue where the routine should be executed. itype Must be 1 or 2 or 3. Specifies the problem type to be solved: if :math:`\text{itype} = 1`, the problem type is :math:`Ax = \lambda Bx;` if :math:`\text{itype} = 2`, the problem type is :math:`ABx = \lambda x;` if :math:`\text{itype} = 3`, the problem type is :math:`BAx = \lambda x`. jobz Must be ``job::novec`` or ``job::vec``. If ``jobz = job::novec``, then only eigenvalues are computed. If ``jobz = job::vec``, then eigenvalues and eigenvectors are computed. upper_lower Must be ``uplo::upper`` or ``uplo::lower``. If ``upper_lower = uplo::upper``, ``a`` and ``b`` store the upper triangular part of :math:`A` and :math:`B`. If ``upper_lower = uplo::lower``, ``a`` and ``b`` stores the lower triangular part of :math:`A` and :math:`B`. n The order of the matrices :math:`A` and :math:`B` (:math:`0 \le n`). a Pointer to array of size ``a(lda,*)`` containing the upper or lower triangle of the Hermitian matrix :math:`A`, as specified by upper_lower. The second dimension of ``a`` must be at least :math:`\max(1, n)`. lda The leading dimension of ``a``; at least :math:`\max(1,n)`. b Pointer to array of size ``b(ldb,*)`` containing the upper or lower triangle of the Hermitian matrix :math:`B`, as specified by upper_lower. The second dimension of ``b`` must be at least :math:`\max(1, n)`. ldb The leading dimension of ``b``; at least :math:`\max(1,n)`. scratchpad_size Size of scratchpad memory as a number of floating point elements of type ``T``. Size should not be less than the value returned by :ref:`onemkl_lapack_hegvd_scratchpad_size` function. events List of events to wait for before starting computation. Defaults to empty list. .. container:: section .. rubric:: Output Parameters a On exit, if ``jobz = job::vec``, then if :math:`\text{info} = 0`, ``a`` contains the matrix :math:`Z` of eigenvectors. The eigenvectors are normalized as follows: if :math:`\text{itype} = 1`` or :math:`\text{itype} = 2`, :math:`Z^{H}BZ = I`; if :math:`\text{itype} = 3`, :math:`Z^{H} B^{-1} Z = I`; If ``jobz = job::novec``, then on exit the upper triangle (if ``upper_lower = uplo::upper``) or the lower triangle (if ``upper_lower = uplo::lower``) of :math:`A`, including the diagonal, is destroyed. b On exit, if :math:`\text{info} \le n`, the part of ``b`` containing the matrix is overwritten by the triangular factor :math:`U` or :math:`L` from the Cholesky factorization :math:`B = U^{H}U`\ or :math:`B` = :math:`LL^{H}`. w Pointer to array of size at least n. If :math:`\text{info} = 0`, contains the eigenvalues of the matrix :math:`A` in ascending order. scratchpad Pointer to scratchpad memory to be used by routine for storing intermediate results. .. container:: section .. rubric:: Throws This routine shall throw the following exceptions if the associated condition is detected. An implementation may throw additional implementation-specific exception(s) in case of error conditions not covered here. :ref:`oneapi::mkl::host_bad_alloc` :ref:`oneapi::mkl::device_bad_alloc` :ref:`oneapi::mkl::unimplemented` :ref:`oneapi::mkl::unsupported_device` :ref:`oneapi::mkl::lapack::invalid_argument` :ref:`oneapi::mkl::lapack::computation_error` Exception is thrown in case of problems during calculations. The ``info`` code of the problem can be obtained by `info()` method of exception object: If :math:`\text{info}=-i`, the :math:`i`-th parameter had an illegal value. For :math:`\text{info} \le n`: If :math:`\text{info}=i`, and ``jobz = oneapi::mkl::job::novec``, then the algorithm failed to converge; :math:`i` indicates the number of off-diagonal elements of an intermediate tridiagonal form which did not converge to zero; If :math:`\text{info}=i`, and ``jobz = oneapi::mkl::job::vec``, then the algorithm failed to compute an eigenvalue while working on the submatrix lying in rows and columns :math:`\text{info}/(n+1)` through :math:`\text{mod}(\text{info},n+1)`. For :math:`\text{info}>n`: If :math:`\text{info}=n+i`, for :math:`1 \le i \le n`, then the leading minor of order :math:`i` of :math:`B` is not positive-definite. The factorization of :math:`B` could not be completed and no eigenvalues or eigenvectors were computed. If ``info`` equals to value passed as scratchpad size, and `detail()` returns non zero, then passed scratchpad is of insufficient size, and required size should not be less than value return by `detail()` method of exception object. .. container:: section .. rubric:: Return Values Output event to wait on to ensure computation is complete. **Parent topic:** :ref:`onemkl_lapack-singular-value-eigenvalue-routines`