MoFEM::SchurDGESV Struct Reference

Static Public Member Functions
static auto	invertMat (MatrixDouble &m, MatrixDouble &inv)

Detailed Description

Definition at line 964 of file Schur.cpp.

Member Function Documentation

invertMat()

static auto MoFEM::SchurDGESV::invertMat ( MatrixDouble &  m, MatrixDouble &  inv  )

inlinestatic

Definition at line 966 of file Schur.cpp.

                                                            {
    MoFEMFunctionBeginHot;
 
    const auto nb = m.size1();
#ifndef NDEBUG
    if (nb != m.size2()) {
      SETERRQ2(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
               "It should be square matrix %d != %d", nb, m.size2());
    }
#endif
 
    inv.resize(nb, nb, false);
    inv.swap(m);
 
    VectorInt ipiv(nb);
    int info;
 
    info = lapack_dgetrf(nb, nb, &*inv.data().begin(), nb, &*ipiv.begin());
    if (info)
      SETERRQ1(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
               "lapack error info = %d", info);
    info = lapack_dgetri(nb, &*inv.data().begin(), nb, &*ipiv.begin(),
                         &*m.data().begin(), m.data().size());
    if (info)
      SETERRQ1(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
               "lapack error info = %d", info);
 
    MoFEMFunctionReturnHot(0);
  }

---

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

• src/finite_elements/impl/Schur.cpp