* Bug #8824 fixed - taucs_chfact returned a segfault (not the case in mode nwni).
[scilab.git] / scilab / modules / umfpack / tests / nonreg_tests / bug_8824.dia.ref
1 // =============================================================================
2 // Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
3 // Copyright (C) 2013 - Scilab Enterprises - Charlotte HECQUET
4 //
5 //  This file is distributed under the same license as the Scilab package.
6 // =============================================================================
7 //
8 // <-- CLI SHELL MODE -->
9 //
10 // <-- Non-regression test for bug 8824 -->
11 //
12 // <-- Bugzilla URL -->
13 // http://bugzilla.scilab.org/show_bug.cgi?id=8824
14 //
15 // <-- Short Description -->
16 // taucs_chfact returns a segfault (not the case in mode nwni)
17 stacksize("max");
18 N=512;
19 h = 1/(N+1);
20 x = h:h:1-h;
21 y = x;
22 [X,Y] = ndgrid(x,y);
23 F = (-2*(%pi^2))*(cos((2*%pi)*X) .*(sin(%pi*Y) .^2) + (sin(%pi*X).^2).*cos((2*%pi)*Y));
24 b = (h^2)*F(:);
25 s = [-1*ones(1,N-1) 2*ones(1,N) -1*ones(1,N-1)]';
26 i = [2:N 1:N 1:N-1]';
27 j = [1:N-1 1:N 2:N ]';
28 ij = [i j];
29 T = sparse(ij,s);
30 I = speye(N,N);
31 A = I .*. T + T .*. I;
32 assert_checktrue(execstr("hchol=taucs_chfact(A)","errcatch")==0);
33 assert_checktrue(execstr("taucs_chsolve(hchol,b)","errcatch")==0);
34 assert_checktrue(execstr("taucs_chdel(hchol)","errcatch")==0);