MoFEM::SCHUR_DGESV Struct Reference

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

Detailed Description

Definition at line 540 of file Schur.cpp.

Member Function Documentation

invertMat()

static auto MoFEM::SCHUR_DGESV::invertMat ( MatrixDouble &  m, MatrixDouble &  inv, double  eps  )

inlinestatic

Definition at line 541 of file Schur.cpp.

                                                                        {
    MoFEMFunctionBeginHot;
    VectorInt ipiv;
    const auto nb = m.size1();
#ifndef NDEBUG
    if (nb != m.size2()) {
      SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
              "It should be square matrix");
    }
#endif
 
    if (eps) {
      for (int cc = 0; cc != nb; ++cc) {
        m(cc, cc) += eps;
      }
    }
 
    inv.resize(nb, nb, false);
    inv.clear();
    auto ptr = &*inv.data().begin();
    for (int c = 0; c != nb; ++c, ptr += nb + 1)
      *ptr = -1;
    ipiv.resize(nb, false);
    const auto info = lapack_dgesv(nb, nb, &*m.data().begin(), nb,
                                   &*ipiv.begin(), &*inv.data().begin(), nb);
    if (info != 0)
      SETERRQ1(PETSC_COMM_SELF, MOFEM_OPERATION_UNSUCCESSFUL,
               "Can not invert matrix info = %d", info);
    MoFEMFunctionReturnHot(0);
  };

---

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

• src/finite_elements/impl/Schur.cpp