quick-test: do not use lapack/create_matrix.h in scalapack test

diff --git a/tests/quick_tests/scalapack.cc b/tests/quick_tests/scalapack.cc
index 2e3876b040d763984ec09e7b2f87f4893d7e9244..f87dbbba66ebf3ce7c180919e2967fef49cfea6c 100644 (file)
--- a/tests/quick_tests/scalapack.cc
+++ b/tests/quick_tests/scalapack.cc
@@ -14,7 +14,6 @@
 // ---------------------------------------------------------------------
 
 #include "../tests.h"
-#include "../lapack/create_matrix.h"
 
 // adaptation of /scalapack/scalapack_02.cc to a quick test
 
@@ -32,6 +31,33 @@
 #include <iostream>
 
 
+
+template <typename FullMatrix>
+void create_spd (FullMatrix &A)
+{
+  const unsigned int size = A.n();
+  Assert (size == A.m(), ExcDimensionMismatch(size,A.m()));
+
+  for (unsigned int i = 0; i < size; ++i)
+    for (unsigned int j = i; j < size; ++j)
+      {
+        const double val = random_value<typename FullMatrix::value_type>();
+        Assert (val >= 0. && val <= 1.,
+                ExcInternalError());
+        if (i==j)
+          // since A(i,j) < 1 and
+          // a symmetric diagonally dominant matrix is SPD
+          A(i,j) = val + size;
+        else
+          {
+            A(i,j) = val;
+            A(j,i) = val;
+          }
+      }
+}
+
+
+
 template <typename NumberType>
 void test(const unsigned int size, const unsigned int block_size)
 {