// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
// Copyright (C) 2010 -  INRIA - Serge Steer
//
// 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 r=%spb_sum(varargin)
    a = varargin(1)
    dims = size(a);

    nargs = size(varargin);
    select nargs
    case 1
        d = "*";
        typ = "double";
    case 2
        if or(varargin(2) == ["native", "double"]) then
            d = "*";
            typ = varargin(2);
        else
            d = varargin(2);
            typ = "double";
        end
    case 3
        d = varargin(2);
        typ = varargin(3);
    else
        error(msprintf(_("%s: Wrong number of input argument(s): %d to %d expected.\n"),"sum", 1, 3));
    end

    // Check second argument : d
    select type(d)
    case 1
        if size(d,'*') <> 1 then
            error(msprintf(_("%s: Wrong size for input argument #%d: A scalar expected.\n"),"sum", 2))
        end
        if int(d) <> d | d < 0 then
            error(msprintf(_("%s: Wrong value for input argument #%d: Integer >= %d expected.\n"),"sum", 2, 1))
        end
    case 10 
        if size(d,'*') <> 1 then
            error(msprintf(_("%s: Wrong size for input argument #%d: A string expected.\n"),"sum",2))
        end
        if and(d<>["r","c","*","m"]) then
            error(msprintf(_("%s: Wrong value for input argument #%d: Must be in the set {%s}.\n"),..
            "sum",2,"""*"",""r"",""c"",""m"""))
        end
        pos = [1,2,0,find(dims>1,1)];
        d=pos(find(d==["r","c","*","m"]))
    else
        error(msprintf(_("%s: Wrong type for input argument #%d: A string or scalar expected.\n"),"sum",2))
    end

    // Check third argument
    if type(typ)<>10 then
        error(msprintf(_("%s: Wrong type for input argument #%d: A string expected.\n"),"sum",3))
    end

    if size(typ,"*")<>1 then
        error(msprintf(_("%s: Wrong size for input argument #%d: A string expected.\n"),"sum",3))
    end

    if and(typ <> ["native" "double"]) then
        error(msprintf(_("%s: Wrong value for input argument #%d: ""%s"" or ""%s"" expected.\n"),"sum", 3, "native", "double"));
    end
    
//    rhs=argn(2)
//    if rhs==1 then
//        d=0 //"*"
//        typ="double"
//    else
//        // call sum(a, d) or sum(a, d, typ)
//        // d must be a string or scalar -> check type and size
//        if and(type(d)<> [1, 10]) then
//            error(msprintf(_("%s: Wrong type for input argument #%d: A string or scalar expected.\n"),"sum",2))
//        end
//
//        if size(d,"*")<>1 then
//            if type(d)==10 then
//                error(msprintf(_("%s: Wrong size for input argument #%d: A string expected.\n"),"sum",2))
//            else
//                error(msprintf(_("%s: Wrong size for input argument #%d: A scalar expected.\n"),"sum",2))
//            end
//        end
//
//        // call sum(a, d) with d = "native" or "double"
//        if rhs == 2 & or(d==["native","double"]) then
//            typ=d
//            d=0 //"*"
//        else
//            // If d is a string, d = "m", "r", "c" or "*"
//            // Else d is an integer > 0
//            if type(d)==10 then
//                d=find(d==["m","*","r","c"])
//                if d==[] then
//                    error(msprintf(_("%s: Wrong value for input argument #%d: Must be in the set {%s}.\n"),..
//                    "sum",2,"""*"",""r"",""c"",""m"",1:"+string(ndims(a))))
//                end
//                d=d-2
//            else
//                if d<0 then
//                    error(msprintf(_("%s: Wrong value for input argument #%d: Must be in the set {%s}.\n"),..
//                    "sum",2,"""*"",""r"",""c"",""m"",1:"+string(ndims(a))))
//                end
//            end
//
//            // call sum(a, d, typ)
//            // typ must be a string and typ must be equal to "native" or "double"
//            // check typ or size
//            if rhs == 3  then
//                if type(typ)<>10 then
//                    error(msprintf(_("%s: Wrong type for input argument #%d: A string expected.\n"),"sum",3))
//                end
//
//                if size(typ,"*")<>1 then
//                    error(msprintf(_("%s: Wrong size for input argument #%d: A string expected.\n"),"sum",3))
//                end
//
//                if and(typ <> ["native", "double"])  then
//                    error(msprintf(_("%s: Wrong value for input argument #%d: ""%s"" or ""%s"" expected.\n"),"sum", 3, "native", "double"));
//                end
//            else
//                typ="double"
//            end
//        end
//    end
//
//    dims=size(a)
//
//    if d==-1 then
//        //sum(x,'m'), determine the summation direction
//        d=find(dims>1,1)
//        if d==[] then d=0,end
//    end

    ij=spget(a);
    select d
    case 0 then
        //sum of all elements
        if typ=="native" then
            r=size(ij,1)<>0
        else
            r=size(ij,1)
        end
    case 1 then
        if typ=="native" then
            r(1,dims(2))=%f
            for k=1:dims(2)
                r(k)=or(ij(:,2)==k)
            end
        else
            r=zeros(1,dims(2));
            for k=1:dims(2)
                r(k)=size(find(ij(:,2)==k),"*")
            end
        end
        r=sparse(r)
    case 2 then
        if typ=="native" then
            r(dims(1),1)=%f
            for k=1:dims(1)
                r(k)=or(ij(:,1)==k)
            end
        else
            r=zeros(ndims(1),1);
            for k=1:dims(1)
                r(k)=size(find(ij(:,1)==k),"*")
            end
        end
        r=sparse(r)
    else
        if typ=="native" then
            r=a
        else
            r=sparse(ij,ones(size(ij,1),1),dims)
        end
    end

endfunction