MoFEM::SchurDSYSV Struct Reference

Static Public Member Functions
static auto	invertMat (const SchurElemMats *row_ptr, MatrixDouble &inv, double eps)

Detailed Description

Definition at line 1147 of file Schur.cpp.

Member Function Documentation

invertMat()

static auto MoFEM::SchurDSYSV::invertMat ( const SchurElemMats *  row_ptr, MatrixDouble &  inv, double  eps  )

inlinestatic

Definition at line 1148 of file Schur.cpp.

                                    {
    MoFEMFunctionBeginHot;
 
    auto &m = row_ptr->getMat();
 
    VectorInt ipiv;
    VectorDouble lapack_work;
    const int nb = m.size1();
 
    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);
    lapack_work.resize(nb * nb, false);
    const auto info =
        lapack_dsysv('L', nb, nb, &*m.data().begin(), nb, &*ipiv.begin(),
                     &*inv.data().begin(), nb, &*lapack_work.begin(), nb * 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