Update samplers to genereate symm pos def matrices

guimspace · May 26, 2021 · 047c7e4 · 047c7e4
1 parent 9bfaece
commit 047c7e4
Show file tree

Hide file tree

Showing 2 changed files with 10 additions and 12 deletions.
diff --git a/c/test/sampler.m b/c/test/sampler.m
@@ -1,13 +1,12 @@
 clear all
 
-n = 10;
+n = 100;
 A = rand(n);
+A = 0.5 * (A + A');
+A = A * A';
+A = A + n * eye(n);
 
-[L, U, P] = lu(A);
-clear L U
-
-A = P * A;
-b = A * [1:n]';
+b = A * [n:-1:1]';
 x = A \ b;
 
 save A.mat A;

diff --git a/python/test/sampler.m b/python/test/sampler.m
@@ -1,13 +1,12 @@
 clear all
 
-n = 10;
+n = 100;
 A = rand(n);
+A = 0.5 * (A + A');
+A = A * A';
+A = A + n * eye(n);
 
-[L, U, P] = lu(A);
-clear L U
-
-A = P * A;
-b = A * [1:n]';
+b = A * [n:-1:1]';
 x = A \ b;
 
 save A.mat A;