function HC=hough_circle_add_point(HC,x,y)

for i=1:HC.rNb
    nb=ceil(2*pi*HC.r(i)/sqrt(HC.xcDelta^2+HC.ycDelta^2)*2);
    nxc=v_to_no(x+HC.xcCircle{i},HC.xcMin,HC.xcMax,HC.xcNb);
    nyc=v_to_no(y+HC.ycCircle{i},HC.ycMin,HC.ycMax,HC.ycNb);
    pos=find((nxc>=0).*(nxc<HC.xcNb).*(nyc>=0).*(nyc<HC.ycNb));
    nxc=nxc(pos);
    nyc=nyc(pos);
    pos=nyc+1+nxc*HC.ycNb+(i-1)*HC.ycNb*HC.xcNb;
    pos=unique(pos);
    HC.stack(pos)=HC.stack(pos)+1;
end
end


function no=v_to_no(v,vmin,vmax,nb)

dv=(vmax-vmin)/(nb-1);
no=floor((v-vmin)/dv-0.5);
end




