Add the <-- JVM NOT MANDATORY --> TAG
[scilab.git] / scilab / modules / umfpack / tests / unit_tests / res_with_prec.tst
1 // ============================================================================
2 // Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
3 // Copyright (C) 2007-2008 - Bruno PINCON
4 //
5 //  This file is distributed under the same license as the Scilab package.
6 // ============================================================================
7
8 // <-- JVM NOT MANDATORY -->
9
10 [A] = ReadHBSparse(SCI+"/modules/umfpack/examples/bcsstk24.rsa");
11 C_ptr = taucs_chfact(A);
12 b = rand(size(A,1),1);
13 x0 = taucs_chsolve(C_ptr, b);
14 firstNorm=norm(A*x0 - b);
15 secondNorm=norm(res_with_prec(A, x0, b));
16 if firstNorm <> secondNorm then pause, end