26 unsigned int nrow =
a.Nrow();
31 for (
unsigned int i = 0; i <
nrow; i++) {
35 s(i) = 1. / std::sqrt(
si);
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);
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));
ROOT::Detail::TRangeCast< T, true > TRangeDynCast
TRangeDynCast is an adapter class that allows the typed iteration through a TCollection.
Class describing a symmetric matrix of size n.
int mnvert(LASymMatrix &t)
Inverts a symmetric matrix.
tbb::task_arena is an alias of tbb::interface7::task_arena, which doesn't allow to forward declare tb...