This repository has been archived by the owner on Nov 10, 2017. It is now read-only.
-
Notifications
You must be signed in to change notification settings - Fork 104
/
ILUTTest.java
74 lines (61 loc) · 2.23 KB
/
ILUTTest.java
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
/*
* Copyright (C) 2003-2006 Bjørn-Ove Heimsund
*
* This file is part of MTJ.
*
* This library is free software; you can redistribute it and/or modify it
* under the terms of the GNU Lesser General Public License as published by the
* Free Software Foundation; either version 2.1 of the License, or (at your
* option) any later version.
*
* This library is distributed in the hope that it will be useful, but WITHOUT
* ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
* FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public License
* for more details.
*
* You should have received a copy of the GNU Lesser General Public License
* along with this library; if not, write to the Free Software Foundation,
* Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
*/
package no.uib.cipr.matrix.sparse;
import no.uib.cipr.matrix.Matrix;
import no.uib.cipr.matrix.Vector;
import org.junit.Test;
import static org.junit.Assert.assertEquals;
/**
* Test of ILU with thresholding
*/
public class ILUTTest extends IncompleteFactorizationTestAbstract {
@Override
void testFactorization(Matrix A, Vector x) {
Vector b = A.mult(x, x.copy());
ILUT ilut = new ILUT(new FlexCompRowMatrix(A));
ilut.setMatrix(A);
ilut.apply(b, x);
Vector r = A.multAdd(-1, x, b.copy());
assertEquals(0, r.norm(Vector.Norm.TwoRobust), 1e-5);
}
/**
* Test for diagInd values exceeding row data array length after dropTol is
* applied causing ArrayOutOfBounds Error
*/
@Test
public void testILUTDropOutOfBoundsError() {
int matrixSize = 100;
FlexCompRowMatrix triDiagMatix = new FlexCompRowMatrix(matrixSize,
matrixSize);
double dropTol = 1e-6;
int p = 50;
for (int i = 0; i < matrixSize; i++) {
for (int j = 0; j < matrixSize; j++) {
if (i == j) {
triDiagMatix.set(i, j, 1);
} else if (Math.abs(i - j) == 1) {
triDiagMatix.set(i, j, dropTol * 1e-2);
}
}
}
ILUT ilut = new ILUT(new FlexCompRowMatrix(triDiagMatix), dropTol, p);
ilut.setMatrix(triDiagMatix);
}
}