]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Fix undefined behavior in TridiagonalMatrix 5242/head
authorDaniel Arndt <daniel.arndt@iwr.uni-heidelberg.de>
Sat, 14 Oct 2017 12:31:53 +0000 (14:31 +0200)
committerDaniel Arndt <daniel.arndt@iwr.uni-heidelberg.de>
Sat, 14 Oct 2017 12:44:43 +0000 (14:44 +0200)
source/lac/tridiagonal_matrix.cc

index 711343ecf60e29fca183049170e0e445558796bd..16937b54457d2dd0acbb95e88f8a92d832096b8c 100644 (file)
@@ -237,7 +237,7 @@ TridiagonalMatrix<double>::compute_eigenvalues()
 
   const int nn = n();
   int info;
-  stev (&N, &nn, &*diagonal.begin(), &*right.begin(), nullptr, &one, nullptr, &info);
+  stev (&N, &nn, diagonal.data(), right.data(), nullptr, &one, nullptr, &info);
   Assert(info == 0, ExcInternalError());
 
   state = LAPACKSupport::eigenvalues;

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.