I got a new problem when running it on FVC2004DB1.
The program failed to handler all the images of FVC2004DB1A. Actually, it does not work for most of the images of this database. But I have got another one that when I tried to assign block size as 32*32 to fft_enhance_cubs(img, BLKSZ). This function got a problem of matrix border.
??? Index exceeds matrix dimensions.
Error in ==> fft_enhance_cubs at 151
I started to read your code 1 week ago. Performance is good. And I have a question about general procedure of your codes. Is it run in the way like
Enhance gray image,
Binarize gray image,
I found a PDF in you .zip file. So, is it everything implemented based on that paper?
While reading your post I was reminded of a long and tedious hunt I recently did for an efficient algorithm to find the distance between two confocal ellipses. One of the things I came across was exactly this problem (point - ellipse), and I remembered I still had the code for it somewhere.
So, I brushed it a bit, and I've just posted it (I made an acknowledgement to you as well). I hope my code can help you out.