Commit 409103ac authored by Carsten Gräser's avatar Carsten Gräser

Fix bounds in SymmetricMatrix::frobenius_norm2()

parent c0345f8d
Pipeline #16272 passed with stage
in 59 minutes and 52 seconds
......@@ -96,8 +96,8 @@ public:
T frobenius_norm2() const
{
T result = 0;
for (size_t i=0; i<N; i++)
for (size_t j=0; j<=N; j++)
for (size_t i=0; i<N; ++i)
for (size_t j=0; j<N; ++j)
result += operator()(i,j)*operator()(i,j);
return result;
}
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment