function maximumsauxsignal=searchlocalmax(threshold, auxsignal, maximumsauxsignal,\
beginingsegment, endsegment, windowbig, windowsmall)
#this function 'inserts' local maximums in the interval [beginingsegment, endsegment]
#to the set maximumsauxsignal
#auxsignal is the signal where we search local maximums
#maximumsauxsignal contains the indices of local maximums
#windowbig is the maximum distance between two local maximums
#windowsmall is the minimum distance between two local maximums
#beginingsegment should be the next indice after the local maximum
#endsegment should be the last indice before the local maximum
#endsegment-beginingsegment must be >=windowbig
#printf("searchlocmax");fflush(stdout);
#extract subsignal
signalsegment=auxsignal(1,beginingsegment:endsegment);
#debug
#assert(length(signalsegment)>=windowbig);#verify this together with findlocmaxpeaks
#aux=findlocmaxpeaks(signalsegment, threshold,windowbig, windowsmall);
aux=findlocmaxpeaksnew(signalsegment, threshold,windowsmall);
if(!isempty(aux))
positions=find(aux>0);#positions of local maximums in the sub segment
#insert new local maximums in the long signal: if it is near a local maximum
# keep the biggest
#first elemente from positions
if(beginingsegment>1)
if(positions(1)<windowsmall)
if(signalsegment(positions(1))>auxsignal(beginingsegment-1))
p=find(maximumsauxsignal==beginingsegment-1);
if(p==1)
maximumsauxsignal=maximumsauxsignal(1,2:end);
elseif(p==length(maximumsauxsignal))
maximumsauxsignal=maximumsauxsignal(1,1:end-1);
else
maximumsauxsignal=[maximumsauxsignal(1,1:p-1),maximumsauxsignal(1,p+1:end)];
endif
else
positions=positions(1,2:end);
endif
endif
endif
#last elemente from positions
if(endsegment<length(auxsignal))
if(length(signalsegment)-positions(end)<windowsmall)
if(signalsegment(positions(end))>auxsignal(endsegment+1))
p=find(maximumsauxsignal==endsegment+1);
if(p==length(maximumsauxsignal))
maximumsauxsignal=maximumsauxsignal(1,1:end-1);
elseif(p==1)
maximumsauxsignal=maximumsauxsignal(1,2:end);
else
maximumsauxsignal=[maximumsauxsignal(1,1:p-1),maximumsauxsignal(1,p+1:end)];
endif
else
positions=positions(1,1:end-1);
endif
endif
endif
#debug
#assert(min(diff(positions))>windowsmall);
# insert positions into maximumsauxsignal
positions=beginingsegment-1+positions;
if(beginingsegment==1)
maximumsauxsignal=[positions, maximumsauxsignal];
elseif(endsegment==length(auxsignal))
maximumsauxsignal=[maximumsauxsignal, positions];
else
p=find(maximumsauxsignal==beginingsegment-1);
maximumsauxsignal=[maximumsauxsignal(1,1:p),positions,maximumsauxsignal(1,p+1:end)];
endif
#assert(diff(maximumsauxsignal)>windowsmall);
# else
# maximumsauxsignal=maximumsauxsignal;
endif
endfunction