function [ ind_izl, postoji_offsp ] = racunamo_offspringe_dec (ind_ul) 
% racunamo_offspringe_dec
% Ulaz: ind_ul - indeks koeficijenta matrice
% Izlazi: ind_izl - indeksi djece koeficijenta
%         postoji_offsp - 
% ind_ul - indeks elementa matrice
% ind_izl - indeksi offspringa
% postoji offspring - zastavica koja pozivajuæoj funkciji govori postoje li
%                     ili ne djeca ulaznog koeficjenta

global n_razl
global sizey
global sizex
ind_izl=[];

postoji_offsp=1;

% i,j sad idu od nula nadalje
i=ind_ul(1,1)-1;
j=ind_ul(1,2)-1;

H_sizex=sizex/(2^n_razl);
H_sizey=sizey/(2^n_razl);

if ((i<H_sizey) & (j<H_sizex))
    if (H_sizey==1) & (H_sizex==1)
        ind_izl(1,1)=1;
        ind_izl(1,2)=2;
        ind_izl(2,1)=2;
        ind_izl(2,2)=1;
        ind_izl(3,1)=2;
        ind_izl(3,2)=2;
        ind_izl(4,1)=0;
        ind_izl(4,2)=0;
    else

        if ((mod(i,2)==0) & (mod(j,2)==1)) %i paran, j neparan (1)
            ind_izl(1,1)=i+1; %sve +1 zbog matlaba
            ind_izl(1,2)=j+H_sizex-1+1;
            ind_izl(2,1)=i+1;
            ind_izl(2,2)=j+H_sizex+1;
            ind_izl(3,1)=i+1+1;
            ind_izl(3,2)=j+H_sizex-1+1;
            ind_izl(4,1)=i+1+1;
            ind_izl(4,2)=j+H_sizex+1;
        elseif ((mod(i,2)==1) & (mod(j,2)==0))
            ind_izl(1,1)=i+H_sizey-1+1; %sve +1 zbog matlaba
            ind_izl(1,2)=j+1;
            ind_izl(2,1)=i+H_sizey-1+1;
            ind_izl(2,2)=j+1+1;
            ind_izl(3,1)=i+H_sizey+1;
            ind_izl(3,2)=j+1;
            ind_izl(4,1)=i+H_sizey+1;
            ind_izl(4,2)=j+1+1;
        elseif ((mod(i,2)==1) & (mod(j,2)==1))
            ind_izl(1,1)=i+H_sizey-1+1; %sve +1 zbog matlaba
            ind_izl(1,2)=j+H_sizex-1+1;
            ind_izl(2,1)=i+H_sizey-1+1;
            ind_izl(2,2)=j+H_sizex+1;
            ind_izl(3,1)=i+H_sizey+1;
            ind_izl(3,2)=j+H_sizex-1+1;
            ind_izl(4,1)=i+H_sizey+1;
            ind_izl(4,2)=j+H_sizex+1;
        end
    end

elseif (((i>=H_sizey) | (j>=H_sizex)) & ((i<(sizey/2))  & (j<(sizex/2)))) 
    %sve +1 jer idu od 1 pa nadalje indeksi
    ind_izl(1,1)=2*i+1;
    ind_izl(1,2)=2*j+1; 
    ind_izl(2,1)=2*i+1;
    ind_izl(2,2)=2*j+1+1;
    ind_izl(3,1)=2*i+1+1;
    ind_izl(3,2)=2*j+1;
    ind_izl(4,1)=2*i+1+1;
    ind_izl(4,2)=2*j+1+1;
else
    ind_izl = zeros (4,2);
    postoji_offsp=0; %nema offspringa
end




