Eigen::internal::lapacke_helpers::lapacke_partial_lu< Scalar, StorageOrder > Struct Template Reference

#include <PartialPivLU_LAPACKE.h>

Static Public Member Functions
static lapack_int	blocked_lu (Index rows, Index cols, Scalar *lu_data, Index luStride, lapack_int *row_transpositions, lapack_int &nb_transpositions, lapack_int maxBlockSize=256)

Member Function Documentation

blocked_lu()

template<typename Scalar , int StorageOrder>

static lapack_int Eigen::internal::lapacke_helpers::lapacke_partial_lu< Scalar, StorageOrder >::blocked_lu ( Index  rows, Index  cols, Scalar *  lu_data, Index  luStride, lapack_int *  row_transpositions, lapack_int &  nb_transpositions, lapack_int  maxBlockSize = 256  )

inlinestatic

performs the LU decomposition in-place of the matrix represented

                                                                                             {
    EIGEN_UNUSED_VARIABLE(maxBlockSize);
    // Set up parameters for getrf
    lapack_int matrix_order = StorageOrder == RowMajor ? LAPACK_ROW_MAJOR : LAPACK_COL_MAJOR;
    lapack_int lda = to_lapack(luStride);
    Scalar* a = lu_data;
    lapack_int* ipiv = row_transpositions;
    lapack_int m = to_lapack(rows);
    lapack_int n = to_lapack(cols);
    nb_transpositions = 0;
 
    lapack_int info = getrf(matrix_order, m, n, to_lapack(a), lda, ipiv);
    eigen_assert(info >= 0);
 
    for (int i = 0; i < m; i++) {
      ipiv[i]--;
      if (ipiv[i] != i) nb_transpositions++;
    }
    lapack_int first_zero_pivot = info;
    return first_zero_pivot;
  }

References a, cols, eigen_assert, EIGEN_UNUSED_VARIABLE, i, info, LAPACK_COL_MAJOR, lapack_int, LAPACK_ROW_MAJOR, lda, m, n, Eigen::RowMajor, rows, and Eigen::internal::lapacke_helpers::to_lapack().

---

The documentation for this struct was generated from the following file:

• PartialPivLU_LAPACKE.h