/*
* 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 no.uib.cipr.matrix.sparse.CompRowMatrix;
import no.uib.cipr.matrix.sparse.ILU;
/**
* Test of ILU with thresholding
*/
public class ILUTTest extends IncompleteFactorizationTestAbstract {
@Override
void testFactorization(Matrix A, Vector x) {
Vector b = A.mult(x, x.copy());
ILU ilut = new ILU(new CompRowMatrix(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);
}
}