1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
|
/*
* Copyright (c) 2009-2017, Peter Abeles. All Rights Reserved.
*
* This file is part of Efficient Java Matrix Library (EJML).
*
* Licensed under the Apache License, Version 2.0 (the "License");
* you may not use this file except in compliance with the License.
* You may obtain a copy of the License at
*
* http://www.apache.org/licenses/LICENSE-2.0
*
* Unless required by applicable law or agreed to in writing, software
* distributed under the License is distributed on an "AS IS" BASIS,
* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
* See the License for the specific language governing permissions and
* limitations under the License.
*/
package org.ejml.sparse.csc;
import org.ejml.UtilEjml;
import org.ejml.data.DMatrixSparseCSC;
import org.junit.Test;
import java.util.Random;
import static org.junit.Assert.*;
/**
* @author Peter Abeles
*/
public class TestRandomMatrices_DSCC {
Random rand = new Random(324);
int numRows = 6;
int numCols = 7;
@Test
public void uniform() {
DMatrixSparseCSC a = RandomMatrices_DSCC.rectangle(numRows,numCols,10,-1,1,rand);
assertEquals(numRows,a.numRows);
assertEquals(numCols,a.numCols);
assertEquals(10,a.nz_length);
assertTrue(CommonOps_DSCC.checkStructure(a));
int count = 0;
for (int row = 0; row < numRows; row++) {
for (int col = 0; col < numCols; col++) {
double value = a.get(row,col);
if( value == 0 )
continue;
if( value > 1 || value < -1 )
fail("Out of expected range");
count++;
}
}
assertEquals(10, count);
assertTrue(CommonOps_DSCC.checkSortedFlag(a));
}
@Test
public void createLowerTriangular() {
DMatrixSparseCSC L;
for (int trial = 0; trial < 20; trial++) {
for( int length : new int[]{0,2,6,12,20} ) {
L = RandomMatrices_DSCC.triangleLower(6, 0, length,-1,1, rand);
assertEquals(Math.max(6,length),L.nz_length);
assertTrue(CommonOps_DSCC.checkStructure(L));
assertTrue(MatrixFeatures_DSCC.isLowerTriangle(L,0, UtilEjml.TEST_F64));
L = RandomMatrices_DSCC.triangleLower(6, 1, length,-1,1, rand);
// L.print();
assertEquals(Math.max(5,length),L.nz_length);
assertTrue(CommonOps_DSCC.checkStructure(L));
assertTrue(MatrixFeatures_DSCC.isLowerTriangle(L,1, UtilEjml.TEST_F64));
assertFalse(CommonOps_DSCC.checkDuplicateElements(L));
}
}
}
@Test
public void symmetric() {
for (int N = 1; N <= 10; N++) {
for (int mc = 0; mc < 30; mc++) {
int nz = (int)(N*N*0.5*(rand.nextDouble()*0.5+0.1)+0.5);
DMatrixSparseCSC A = RandomMatrices_DSCC.symmetric(N, nz,-1,1, rand);
assertTrue(CommonOps_DSCC.checkStructure(A));
// Check the matrix properties
assertTrue(MatrixFeatures_DSCC.isSymmetric(A,UtilEjml.TEST_F64));
}
}
}
@Test
public void symmetricPosDef() {
for (int N = 1; N <= 10; N++) {
for (int mc = 0; mc < 30; mc++) {
int nz = (int)(N*N*0.5*(rand.nextDouble()*0.5+0.1)+0.5);
DMatrixSparseCSC A = RandomMatrices_DSCC.symmetricPosDef(N, nz,rand);
assertTrue(CommonOps_DSCC.checkStructure(A));
assertTrue(MatrixFeatures_DSCC.isPositiveDefinite(A));
}
}
}
}
|
