// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
// Copyright (C) 2011 - DIGITEO
//
// Copyright (C) 2012 - 2016 - Scilab Enterprises
//
// This file is hereby licensed under the terms of the GNU GPL v2.0,
// pursuant to article 5.3.4 of the CeCILL v.2.1.
// This file was originally licensed under the terms of the CeCILL v2.1,
// and continues to be available under such terms.
// For more information, see the COPYING file which you should have received
// along with this program.
test_path = SCI + "/modules/signal_processing/tests/unit_tests/";
[fact,b2,b1,b0,c1,c0,zzeros,zpoles]=syredi(1,4,[2*%pi/10,4*%pi/10,0,0],0.02,0.001);
ref_fact = 0.005976;
assert_checkfalse(norm(fact-ref_fact) > 0.0000037);
assert_checkequal(b2, [1 1 1]);
assert_checkequal(size(b1), [1 3]);
assert_checkequal(b1(3), 1);
assert_checkequal(b0, [1 1 0]);
assert_checkequal(size(c1), [1 3]);
assert_checkequal(size(c0), [1 3]);
assert_checkequal(c0(3),  0);
assert_checkequal(size(zzeros), [1 5]);
assert_checkequal(size(zpoles), [1 5]);
load(test_path+"syredi.h5");
[fact1, b21, b11, b01, c11, c01, zzeros1, zpoles1] = syredi(1, 4, [1, 2, 0, 0], 0.02, 0.001);
assert_checkequal(fact1, fact1_ref);
assert_checkequal(b21, b21_ref);
assert_checkalmostequal(b11, b11_ref);
assert_checkequal(b01, b01_ref);
assert_checkalmostequal(c11, c11_ref);
assert_checkalmostequal(c01, c01_ref, [], %eps);
assert_checkalmostequal(zzeros1, zzeros1_ref);
assert_checkalmostequal(zpoles1, zpoles1_ref);
[fact2, b22, b12, b02, c12, c02, zzeros2, zpoles2] = syredi(3, 4, [1, 2, 3, 3.01], 0.1, 0.001);
assert_checkalmostequal(fact2, fact2_ref);
assert_checkequal(b22, b22_ref);
assert_checkalmostequal(b12, b12_ref);
assert_checkequal(b02, b02_ref);
assert_checkalmostequal(c12, c12_ref);
assert_checkalmostequal(c02, c02_ref);
assert_checkalmostequal(zzeros2, zzeros2_ref);
assert_checkalmostequal(zpoles2, zpoles2_ref);