26 unsigned int nrow =
a.Nrow();
31 for (
unsigned int i = 0; i < nrow; i++) {
38 for (
unsigned int i = 0; i < nrow; i++)
39 for (
unsigned int j = i; j < nrow; j++)
40 a(i, j) *= (
s(i) *
s(j));
42 for (
unsigned i = 0; i < nrow; i++) {
49 unsigned int kp1 = k + 1;
51 for (
unsigned int j = 0; j < k; j++) {
53 q(j) =
a(j, k) *
q(k);
58 for (
unsigned int j = kp1; j < nrow; j++) {
60 q(j) = -
a(k, j) *
q(k);
64 for (
unsigned int j = 0; j < nrow; j++)
65 for (k = j; k < nrow; k++)
66 a(j, k) += (pp(j) *
q(k));
69 for (
unsigned int j = 0; j < nrow; j++)
70 for (
unsigned int k = j; k < nrow; k++)
71 a(j, k) *= (
s(j) *
s(k));
Class describing a symmetric matrix of size n.
VecExpr< UnaryOp< Sqrt< T >, VecExpr< A, T, D >, T >, T, D > sqrt(const VecExpr< A, T, D > &rhs)
int mnvert(LASymMatrix &t)
Inverts a symmetric matrix.
This file contains a specialised ROOT message handler to test for diagnostic in unit tests.
static constexpr double s