public void testThreeSparse() {
int K = 60;
int N = thousandconstant;
SparseMatrixLil A = SparseMatrixLil.rand(N, K);
SparseMatrixLil B = SparseMatrixLil.rand(K, N);
// System.out.println(A.toDense());
Timer timer = new Timer();
SparseMatrixLil C = B.mmul(A);
timer.printTimeCheckMilliseconds();
SparseMatrixLil C1 = B.mmul(A);
timer.printTimeCheckMilliseconds();
SparseMatrixLil C2 = B.mmul(A);
timer.printTimeCheckMilliseconds();
DenseMatrix CDense = B.toDense().mmul(A.toDense());
timer.printTimeCheckMilliseconds("for dense");
assertTrue(C.equals(CDense));
K = 400;
N = thousandconstant;
A = SparseMatrixLil.rand(N, K);
B = SparseMatrixLil.rand(K, N);
// System.out.println(A.toDense());
timer = new Timer();
C = B.mmul(A);
timer.printTimeCheckMilliseconds();
C1 = B.mmul(A);
timer.printTimeCheckMilliseconds();
C2 = B.mmul(A);
timer.printTimeCheckMilliseconds();
CDense = B.toDense().mmul(A.toDense());
timer.printTimeCheckMilliseconds("for dense");
assertTrue(C.equals(CDense));
}
public void testToCCS() {
SparseMatrixLil A = sprand(10, 10);
SparseMatrixLil B = A.toCCS().toLil();
assertTrue(A.equals(B));
A = sprand(thousandconstant, thousandconstant);
B = A.toCCS().toLil();
assertTrue(A.equals(B));
A = A.t();
B = A.toCCS().toLil();
assertTrue(A.equals(B));
}
public void testCreateSparse() {
DenseMatrix A = rand(5, 5);
SparseMatrixLil B = A.toSparseLil();
assertTrue(B.equals(B.mmul(eye(5))));
assertTrue(B.equals(B.t().mmul(eye(5)).t()));
B = spzeros(2, 2);
B.append(1, 0, 5);
B.append(0, 1, 7);
B.append(0, 0, 11);
System.out.println(B.toDense());
System.out.println(B.mmul(eye(2)));
assertTrue(B.equals(B.mmul(eye(2))));
}
public void testSliceSparse() {
SparseMatrixLil A = sprand(5, 8);
SparseMatrixLil B = sprand(5, 12);
SparseMatrixLil C = sprand(4, 8);
SparseMatrixLil AB = A.concatRight(B);
assertTrue(A.equals(AB.cols(0, 8)));
assertTrue(B.equals(AB.cols(8, 20)));
SparseMatrixLil AC = A.concatDown(C);
assertTrue(A.equals(AC.rows(0, 5)));
assertTrue(C.equals(AC.rows(5, 9)));
A = spzeros(1, 5);
A.append(0, 2, 3);
assertEquals(A.row(0).cols, 5);
}
public void testSparseMultiply() {
SparseMatrixLil A = sprand(3, 5);
System.out.println(A);
SparseMatrixLil A2 = A.mmul(speye(5));
System.out.println(A2);
assertTrue(A.equals(A2));
}
public void testRowsIndexed() {
SparseMatrixLil A = new SparseMatrixLil("1 5 3; 2 8 5; 1 9 4; 2 5 3");
DenseMatrix indexes = new DenseMatrix("2; 1; 1");
SparseMatrixLil B = A.rows(indexes);
System.out.println(B.toDense());
assertTrue(B.equals(new SparseMatrixLil("1 9 4; 2 8 5; 2 8 5")));
}
public void testConstructorByStringSparse() {
SparseMatrixLil foo = new SparseMatrixLil("1 0; 5 7; 0 5");
assertEquals(3, foo.rows);
assertEquals(2, foo.cols);
assertTrue(foo.equals(new DenseMatrix(new double[][] {{1, 0}, {5, 7}, {0, 5}})));
assertEquals(5.0, foo.toDense().get(1, 0));
assertEquals(7.0, foo.toDense().get(1, 1));
assertEquals(4, foo.size);
}
public void testOneSparse() {
SparseMatrixLil A = SparseMatrixLil.rand(3, 3).t();
System.out.println(A);
SparseMatrixLil B = SparseMatrixLil.rand(3, 3);
System.out.println(B);
SparseMatrixLil C = A.mmul(B);
System.out.println(C);
DenseMatrix CDense = A.toDense().mmul(B.toDense());
System.out.println(CDense);
assertTrue(C.equals(CDense));
}
public void testSort() {
SparseMatrixLil B = spzeros(2, 2);
B.append(1, 0, 5);
B.append(0, 1, 7);
B.append(0, 0, 11);
System.out.println(B.toDense());
System.out.println(B);
B.sort();
System.out.println(B.toDense());
System.out.println(B);
assertEquals(11.0, B.values[0]);
assertEquals(5.0, B.values[1]);
assertEquals(7.0, B.values[2]);
SparseMatrixLil A = sprand(5, 8);
B = A.add(0); // make copy
System.out.println(B);
B.sort();
for (int i = 1; i < B.size; i++) {
assertTrue(B.colIdx[i] >= B.colIdx[i - 1]);
if (B.colIdx[i] == B.colIdx[i - 1]) {
assertTrue(B.rowIdx[i] >= B.rowIdx[i - 1]);
}
}
assertTrue(A.equals(B));
System.out.println(B);
SparseMatrixLil Bt = B.t();
System.out.println(Bt);
Bt.sort();
for (int i = 1; i < Bt.size; i++) {
assertTrue(Bt.colIdx[i] >= Bt.colIdx[i - 1]);
if (Bt.colIdx[i] == Bt.colIdx[i - 1]) {
assertTrue(Bt.rowIdx[i] >= Bt.rowIdx[i - 1]);
}
}
System.out.println(Bt);
assertTrue(A.toDense().t().equals(Bt));
}