442 {
446 const int nb =
m.size1();
447
449 for (int cc = 0; cc != nb; ++cc) {
451 }
452 }
453
454 inv.resize(nb, nb, false);
455 inv.clear();
456 for (int cc = 0; cc != nb; ++cc)
457 inv(cc, cc) = -1;
458 ipiv.resize(nb, false);
459 lapack_work.resize(nb * nb, false);
460 const auto info =
461 lapack_dsysv(
'L', nb, nb, &*
m.data().begin(), nb, &*ipiv.begin(),
462 &*inv.data().begin(), nb, &*lapack_work.begin(), nb * nb);
463 if (info != 0)
465 "Can not invert matrix info = %d", info);
467 };
#define MoFEMFunctionReturnHot(a)
Last executable line of each PETSc function used for error handling. Replaces return()
@ MOFEM_OPERATION_UNSUCCESSFUL
#define MoFEMFunctionBeginHot
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
FTensor::Index< 'm', SPACE_DIM > m
static __CLPK_integer lapack_dsysv(char uplo, __CLPK_integer n, __CLPK_integer nrhs, __CLPK_doublereal *a, __CLPK_integer lda, __CLPK_integer *ipiv, __CLPK_doublereal *b, __CLPK_integer ldb, __CLPK_doublereal *work, __CLPK_integer lwork)
UBlasVector< double > VectorDouble
UBlasVector< int > VectorInt