// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
// Copyright (C) INRIA -
//
// This file must be used under the terms of the CeCILL.
// This source file is licensed as described in the file COPYING, which
// you should have received as part of this distribution.  The terms
// are also available at
// http://www.cecill.info/licences/Licence_CeCILL_V2.1-en.txt


function [ar]=armac(a,b,d,ny,nu,sig)
    // just build a tlist for storing armacx coefficients
    //   A(z^-1)y= B(z^-1)u + D(z^-1)sig*e(t)
    //   a=<Id,a1,..,a_r>; matrix (ny,r*ny)
    //   b=<b0,.....,b_s>; matrix (ny,(s+1)*nu)
    //   d=<Id,d1,..,d_p>; matrix (ny,p*ny);
    //   ny : dim of observation y
    //   nu : dim of control  u
    //   sig : standard deviation  (ny,ny);
    //
    //!

    [na,la]=size(a);
    if na<>ny then
        error(msprintf(gettext("%s: Wrong size for input argument #%d: row dimension must be equal to %d.\n"),..
        "armac",1, ny));

    end
    [nb,lb]=size(b);
    if nb<>0 & nb<>ny then
        error(msprintf(gettext("%s: Wrong size for input argument #%d: row dimension must be equal to %d.\n"),..
        "armac",2, ny));
    end;
    if lb<>0 & nu<>0 then
        if modulo(lb,nu)<>0 then
            error(msprintf(gettext("%s: Wrong size of input argument #%d: Number of columns are incompatible with %s.\n"),..
            "armac",2,"nu"));
        end;
    end
    [nd,ld]=size(d);
    if nd<>ny then
        error(msprintf(gettext("%s: Wrong size for input argument #%d: row dimension must be equal to %d.\n"),..
        "armac",3, ny));

    end
    ar=tlist(["ar","a","b","d","ny","nu","sig"],a,b,d,ny,nu,sig);
endfunction