// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
// Copyright (C) 2005 - INRIA - Farid Belahcene
// Copyright (C) 2012 - Michael Baudin
// 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  bar(varargin)
    // bar(x,y,width,style,color)
    // This function  ...
    //
    // Input :
    // x : a real scalar or a vector
    // y : a real sclar, or a vector
    // width : a double, the bar width, it's the percentage (0<width<1) of the width max of one bar which is wanted (default: width=0.8)
    // style : a string 'grouped' or 'stacked' (default: style='grouped')

    if size(varargin)<1 | size(varargin)>5  then
        error(msprintf(gettext("%s: Wrong number of input argument(s): %d to %d expected.\n"), "bar", 1, 5));
    end

    styletab=["grouped","stacked"]
    COLORBOOL=%f
    STYLE="grouped"

    //Check RHS argument
    ListArg = varargin;

    //detect and set the current axes now:
    if type(ListArg(1)) == 9
        hdle = ListArg(1);
        if (hdle.type == "Axes")
            sca(ListArg(1));
            ListArg(1) = null(); // remove this parameter from the list
        else
            warning(msprintf(gettext("%s: Wrong type for input argument #%d: Axes handle expected.\n"),"bar",1));
            return;
        end
    end

    if size(ListArg) == 4 then
        COLOR=ListArg(4);
        if type(COLOR) <> 10 then
            error(msprintf(gettext("%s: Wrong type for input arguments #%d: A string expected.\n"),"bar",4));
        end
    end
    if size(ListArg) == 5 then
        STYLE=ListArg(5);
        if type(STYLE) <> 10 then
            error(msprintf(gettext("%s: Wrong type for input arguments #%d: A string expected.\n"),"bar",5));
        end
    end
    nv = size(ListArg)

    T=[];

    // Number of inputs arguments < 6
    if  size(ListArg)>5 then
        error(msprintf(gettext("%s: Wrong number of input arguments: %d to %d expected.\n"),"bar",1,5));
    end

    for k=1:nv
        T(k) = type(ListArg(k))
    end

    argdb=find(T==1)
    argstr=find(T==10)

    if size(argdb,"*")<> argdb($) then
        error(msprintf(gettext("%s: Wrong type for input arguments: Matrix expected for %s, %s and %s.\n"),"bar", "x", "y", "width"));
    end

    if size(argstr,"*") <> nv-argdb($) then
        error(msprintf(gettext("%s: Wrong type for input arguments: String expected for %s and %s.\n"),"bar", "color", "style"));
    end

    //set the double argument : x,y,width
    // bar(y,...)
    if size(argdb,"*")==1
        Y=ListArg(1)
        WIDTH=0.8
        if or(size(Y)==1) then
            Y=Y(:)
        end
        X=1:size(Y,1)
    end

    if size(argdb,"*")==2
        if size(ListArg(2),"*")==1 then
            // bar(x,y,...)
            if size(ListArg(1),"*")==1 then
                WIDTH=0.8
                X=ListArg(1)
                Y=ListArg(2)
            else
                //bar(y,width,...)
                WIDTH=ListArg(2)
                Y=ListArg(1)
                if or(size(Y)==1) then
                    Y=Y(:)
                end
                X=1:size(Y,1)
            end
        else
            // bar(x,y,...)
            X=ListArg(1)
            Y=ListArg(2)
            if or(size(X)==1) then
                if size(X,"*")<>1 then // X is a vector
                    if or(size(Y)==1) then // Y is a vector
                        Y=Y(:)
                    end
                    if size(X,"*")<>size(Y,1) // Y is a matrix
                        error(msprintf(gettext("%s: Wrong size for input arguments #%d and #%d: The number of rows of argument #%d must be equal to the size of argument #%d.\n"),"bar",1, 2, 2, 1));
                    end
                elseif size(Y,1)>1 then
                    error(msprintf(gettext("%s: Wrong size for input argument #%d: A scalar or a column vector expected.\n"),"bar",2));
                end
            else
                error(msprintf(gettext("%s: Wrong type for input argument #%d: A scalar or a vector expected.\n"),"bar",1));
            end
            WIDTH=0.8
        end
    end

    // bar(x,y,width,...)
    if size(argdb,"*")==3
        X=ListArg(1)
        Y=ListArg(2)
        WIDTH=ListArg(3)
        if size(WIDTH,"*")<>1 then
            error(msprintf(gettext("%s: Wrong type for input argument #%d: A scalar expected.\n"),"bar",3));
        elseif or(size(X)==1) then
            if size(X,"*")<>1 then // X is a vector
                if or(size(Y)==1) then // Y is a vector
                    Y=Y(:)
                end
                if size(X,"*")<>size(Y,1)
                    error(msprintf(gettext("%s: Wrong size for input arguments #%d and #%d: The number of rows of argument #%d must be equal to the size of argument #%d.\n"),"bar",1, 2, 2, 1))
                end
            elseif size(Y,1)>1 then
                error(msprintf(gettext("%s: Wrong size for input arguments #%d: A scalar or a column vector expected.\n"),"bar",2));
            end
        else
            error(msprintf(gettext("%s: Wrong type for input argument #%d: A scalar or a vector expected.\n"),"bar",1));
        end
    end
    X=X(:)

    // set the string argument
    for i=1:size(argstr,"*")
        // bar(...,style)
        if or(ListArg(argstr(i))==styletab) then
            STYLE=ListArg(argstr(i))
        else
            COLOR=ListArg(argstr(i))
            COLORBOOL=%t
        end
    end

    // drawlater
    curFig = gcf();
    immediate_drawing = curFig.immediate_drawing;

    wmode = warning("query");
    warning("off"); // See bug #13579 (some bar() calling sequences will lead to a plot() warning)
    if COLORBOOL
        plot(X,Y,COLOR); // plot manages immediate_drawing property itself to avoid flickering
    else
        plot(X,Y); // plot manages immediate_drawing property itself to avoid flickering
    end
    warning(wmode);

    curFig.immediate_drawing = "off";

    bar_number=size(Y,2)
    if size(X,"*")>1 then
        Xtemp=gsort(X,"r","i")
        inter=Xtemp(2)-Xtemp(1)
        for i=2:size(Xtemp,"*")-1
            inter=min(Xtemp(i+1)-Xtemp(i),inter)
        end
        if bar_number>1
            inter=inter*0.9
        end
    else
        Xtemp=X
        inter=1
    end

    wmax=inter/bar_number

    y_shift=zeros(size(X,"*"),1)

    bar_number= bar_number

    e=gce()
    a=gca()

    a.sub_ticks(1) = 0; // bar (barh => a.sub_ticks(2) = 0;)

    for i=bar_number:-1:1

        ei = e.children(i);

        // Perform x_shift
        if modulo(bar_number,2)==0 then
            x_shift=(-i+bar_number/2)*wmax+wmax/2
        elseif modulo(bar_number,2)==1 then
            x_shift=(-i+1+floor(bar_number/2))*wmax
        end

        // Perform y_shift
        if i==bar_number then
            y_shift=zeros(size(X,"*"),1)
        else
            y_shift=Y(:,bar_number-i)+y_shift
        end

        // Udate the axes data bounds
        if STYLE=="grouped"
            xmin=min(a.data_bounds(1,1),min(X)+x_shift-0.4*wmax)
            ymin=min(a.data_bounds(1,2),0,min(y_shift+Y(:,bar_number-i+1)))
            xmax=max(a.data_bounds(2,1),max(X)+x_shift+0.4*wmax)
            ymax=max(a.data_bounds(2,2),0)
            ei.x_shift=x_shift*ones(size(X,"*"),1)
        else
            wmax=inter
            xmin=min(a.data_bounds(1,1),min(X)-0.4*wmax)
            ymin=min(a.data_bounds(1,2),min(y_shift+Y(:,bar_number-i+1)))
            xmax=max(a.data_bounds(2,1),max(X)+0.4*wmax)
            ymax=max(a.data_bounds(2,2),max(y_shift+Y(:,bar_number-i+1)))
            ei.y_shift=y_shift
        end
        a.data_bounds=[xmin ymin; xmax ymax]

        a.x_ticks=tlist("ticks",Xtemp,string(Xtemp))

        w=WIDTH*wmax

        ei.bar_width=w
        ei.background=ei.foreground
        ei.polyline_style=6; // bar type
        ei.background=ei.foreground
        ei.foreground = -1; // black by default
        ei.line_mode="off";
    end

    // drawnow
    curFig.immediate_drawing = immediate_drawing;

endfunction