MoFEM::SchurDSYSV Struct Reference

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

Detailed Description

Definition at line 1001 of file Schur.cpp.

Member Function Documentation

invertMat()

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

inlinestatic

Definition at line 1003 of file Schur.cpp.

                                                            {
    MoFEMFunctionBegin;
 
    const auto nb = m.size1();
#ifndef NDEBUG
    if (nb != m.size2()) {
      SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
              "It should be square matrix %ld != %ld", nb, m.size2());
    }
#endif
 
    inv.resize(nb, nb, false);
    inv.swap(m);
    m.resize(2 * nb, 2 * nb, false);
 
    VectorInt ipiv(nb);
    int info;
 
    // FIXME: add lapack_dsytrf and lapack_dsytri
 
    info = lapack_dgetrf(nb, nb, &*inv.data().begin(), nb, &*ipiv.begin());
    if (info)
      SETERRQ(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)
      SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
               "lapack error info = %d", info);
 
    MoFEMFunctionReturn(0);
  }

---

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

• src/finite_elements/impl/Schur.cpp