1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
|
// =============================================================================
// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
// Copyright (C) 2007-2008 - INRIA - Serge STEER <serge.steer@inria.fr>
//
// This file is distributed under the same license as the Scilab package.
// =============================================================================
// <-- Non-regression test for bug 3310 -->
//
// <-- Bugzilla URL -->
// http://bugzilla.scilab.org/show_bug.cgi?id=3310
//
// <-- Short Description -->
// Wrong implementation of input arguments
A=[10,1;1,10];
// Suppose that the rhs is in myrhs instead of b
myrhs=[11;11];
[xcomputed, flag, err, iter, res]=pcg(A,myrhs);
xexpected=[1;1];
if norm(xcomputed-xexpected)>%eps then bugmes();quit;end
|
