// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
// Copyright (C) ???? - INRIA - Farid BELAHCENE
// Copyright (C) 2008 - INRIA - Pierre MARECHAL
//
// 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 z = bitand(x, y)

    // BITAND function
    // Given x,y two positives integers this function returns the decimal
    // number whose the binary form is the AND of the binary representations
    // of x and y
    //
    // -Inputs :
    //  x, y :  scalars/vectors/matices/hypermatices of positives integers, x and y must have the same size
    // -Output :
    //  z : scalar/vector/matrix/hypermatice
    //
    // F.Belahcene

    // P. Marechal, 5 Feb 2008
    //   - fix bug 2691 and 2692
    //   - Add argument check

    // Check input arguments
    // =========================================================================

    rhs = argn(2);

    if rhs <> 2 then
        error(msprintf(gettext("%s: Wrong number of input argument(s): %d expected.\n"),"bitand",2));
    end

    if typeof(x)<>typeof(y)
        error(msprintf(gettext("%s: Wrong type for input arguments: Same types expected.\n"),"bitand"));
    end

    if or(size(x)<>size(y)) then
        error(msprintf(gettext("%s: Wrong size for input arguments: Same size expected.\n"),"bitand"));
    end

    if isempty(x) & isempty(x)
        z=[];
        return
    end

    if    (type(x)==1  & (x-floor(x)<>0 | x<0)) ..
        | (type(x)==8  & (inttype(x)<10) ) ..
        | (type(x)<>1  & type(x)<>8) then

        error(msprintf(gettext("%s: Wrong input argument #%d: Scalar/matrix of unsigned integers expected.\n"),"bitand",1));
    end

    if    (type(y)==1  & (y-floor(y)<>0 | y<0)) ..
        | (type(y)==8  & (inttype(y)<10) ) ..
        | (type(y)<>1  & type(y)<>8) then

        error(msprintf(gettext("%s: Wrong input argument #%d: Scalar/matrix of unsigned integers expected.\n"),"bitand",2));
    end

    // Algorithm
    // =========================================================================

    if type(x)==8 then
        z = x & y;
    else
        a = 2^32;
        z = double( uint32(x/a) & uint32(y/a) ) * a;
        z = z + double(uint32(x-z) & uint32(y-z));
    end

endfunction