// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
// Copyright (C) 2010 - INRIA - Serge Steer <serge.steer@inria.fr>
//
// 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 out=%s_i_st(varargin)

    // - Modified by Vincent COUVERT (16/08/2004) so that insertion of an empty matrix
    //   is understood as an element deletion
    //   Only one non-colon index can be used
    // - Modified by Serge Steer INRIA (04/05/2010) to fix problems in element
    //   deletion part

    if size(varargin)>=3 & isempty(varargin($-1)) & ..
        and(type(varargin(1))<>[10 15]) then
        // st(i,:)=[] or st(:,j)=[] or st(i)=[] or st(:,j,:,:)=[]
        //remove the substruct
        out=varargin($);
        dims=double(out.dims);

        // Make the dimensions and the indices fit
        Ndims=size(dims,"*")
        nindex=size(varargin)-2
        if nindex>Ndims then
            //index in excess must be equal to 1 or to :
            for k=Ndims+1:nindex
                i=varargin(k)
                if size(i,"*")>1|(i<>1&i<>eye()) then
                    error(msprintf(_("%s: A null assignment can have only one non-colon index.\n"),"%s_i_st"));
                end
            end
            nindex=Ndims
        elseif nindex<Ndims then
            //collapse dimensions in excess
            dims=[dims(1:nindex-1) prod(dims(nindex:$))]
            Ndims=nindex;
            if size(dims,"*")==1 then dims=[dims 1],end
        end

        // Check the compatibility of the index (at most one index cannot span
        // all the elements ot the associated struct dimension)
        cj=[];
        for k=1:nindex
            ind=varargin(k)
            if or(size(ind)<>[-1 -1]) then
                if or(type(ind)==[2,129]) then // size implicit index ($ based)
                    ind=horner(ind,dims(k));
                end
                ind=floor(ind);
                //check if index is valid
                if ~isreal(ind)|or(ind<=0) then
                    error(21)
                end
                //remove indices that exceed the associated struct dimension
                ind(ind>dims(k))=[];
                //compute the complement with respect to the associated dimension of st
                ind=setdiff(1:dims(k),ind)
                if ind<>[]&cj==[] then
                    cj=ind
                    loc=k,
                else
                    error(msprintf(_("%s: A null assignment can have only one non-colon index.\n"),"%s_i_st"));
                end
            end
        end

        // Generate the result
        if cj==[] then  //st(:,:)=[]  --> empty struct
            Fout=getfield(1,out)
            Fout=Fout(3:$)
            for f=Fout
                out(f)=list()
            end
            out.dims=int32([0 0])
        else
            //replace st(:,j,:,:)=[] by st=st(:,cj,:,:) where cj is the
            //complement of j with respect to the associated dimension of st
            out.dims=int32(dims)
            varargin(loc)=cj
            out=out(varargin(1:Ndims))
        end

    elseif lstsize(varargin)==3 & type(varargin(1))==10 then // out.i=in
        i=varargin(1);
        in=varargin(2);
        out=varargin(3);
        out=generic_i_st(i,in,out)
    else
        error(msprintf(_("%s: Not yet implemented.\n"),"%s_i_st"));
    end
endfunction