diff options
author | Jonathan McCrohan <jmccrohan@gmail.com> | 2012-10-30 21:29:30 +0000 |
---|---|---|
committer | Jonathan McCrohan <jmccrohan@gmail.com> | 2012-10-30 21:29:30 +0000 |
commit | 62f81e5cafbacfb90ac0f86e52e8e3176aa5ba8f (patch) | |
tree | 8f8f6028fa880b3296b1e42bc8d87cead73aaf78 /correlation.c | |
download | dhex-62f81e5cafbacfb90ac0f86e52e8e3176aa5ba8f.tar.gz |
Imported Upstream version 0.65upstream/0.65
Diffstat (limited to 'correlation.c')
-rw-r--r-- | correlation.c | 272 |
1 files changed, 272 insertions, 0 deletions
diff --git a/correlation.c b/correlation.c new file mode 100644 index 0000000..2c0e68e --- /dev/null +++ b/correlation.c @@ -0,0 +1,272 @@ +#include "correlation.h" +#include "buffers.h" +#include <strings.h> +#include <string.h> +void find_longestmatch(tOutput* output,tCorrelation* correlation,tBuffer* buf1,tBuffer* buf2,tInt64* cursorpos1,tInt64* cursorpos2) +{ + tInt64 idx1; + tInt64 idx2; + tInt64 match; + tInt64 sumpos; + tInt64 pos1; + tInt64 pos2; + tInt16 offsx; + tInt16 offsy; + WINDOW* smallwin=NULL; + tInt64 diff; + tBool found=0; + + + if (output) + { + offsx=COLS/2-22; + offsy=LINES/2-3; + drawcenterframe(output,6,44,"correlating..."); + setcolor(output,COLOR_TEXT); + smallwin=subwin(output->win,6,19,offsy,offsx+24); + mvwprintw(output->win,offsy+1,offsx+1,"pos1:"); + mvwprintw(output->win,offsy+4,offsx+1,"longest match:"); + mvwprintw(output->win,offsy+5,offsx+1,"at"); + if (smallwin) mvwprintw(smallwin,4,1,"%18lli",0); + wrefresh(output->win); + } + + correlation->bestmatch=0; + correlation->mindiff=MAX(buf1->bufsize,buf2->bufsize); + for (pos1=-(tInt64)buf2->bufsize;pos1<(tInt64)buf1->bufsize;pos1++) + { + match=0; + diff=0; + if (smallwin) mvwprintw(smallwin,1,3,"%16lli",((tInt64)buf1->bufsize-pos1)); + else fprintf(stderr,"%16lli\r",((tInt64)buf1->bufsize-pos1)); + if (smallwin) wrefresh(smallwin); + for (pos2=0;pos2<buf2->bufsize;pos2++) + { + sumpos=pos1+pos2; + if (sumpos>=0 && sumpos<buf1->bufsize) + { + idx1=getbufferidx(buf1,sumpos); + idx2=getbufferidx(buf2,pos2); + + if (buf1->data[idx1]==buf2->data[idx2]) match++; + else { + diff++; + if (match>correlation->bestmatch) + { + correlation->bestmatch=match; + correlation->bestmatchpos=pos1; + correlation->mindiff=diff; + if (smallwin) mvwprintw(smallwin,4,1,"%18lli",match); + if (smallwin) mvwprintw(smallwin,5,3,"%16llx",pos1); + if (smallwin) wrefresh(smallwin); + found=1; + } + match=0; + } + if (match>correlation->bestmatch) + { + correlation->bestmatch=match; + correlation->bestmatchpos=pos1; + correlation->mindiff=diff; + if (smallwin) mvwprintw(smallwin,4,1,"%18lli",match); + if (smallwin) mvwprintw(smallwin,5,3,"%16llx",pos1); + if (smallwin) wrefresh(smallwin); + found=1; + } + } + } + if (smallwin) wrefresh(smallwin); + } + if (found) + { + if (correlation->bestmatch>=0) + { + *cursorpos1=correlation->bestmatchpos; + *cursorpos2=0; + } else { + *cursorpos1=0; + *cursorpos2=-correlation->bestmatchpos; + } + } + if (smallwin) delwin(smallwin); +} +void find_bestmatch(tOutput* output,tCorrelation* correlation,tBuffer* buf1,tBuffer* buf2,tInt64* cursorpos1,tInt64* cursorpos2) +{ + tInt64 idx1; + tInt64 idx2; + tInt64 match; + tInt64 sumpos; + tInt64 pos1; + tInt64 pos2; + tInt16 offsx; + tInt16 offsy; + WINDOW* smallwin=NULL; + tInt64 diff; + tBool found=0; + + + if (output) + { + offsx=COLS/2-22; + offsy=LINES/2-3; + drawcenterframe(output,6,44,"correlating..."); + setcolor(output,COLOR_TEXT); + smallwin=subwin(output->win,6,19,offsy,offsx+24); + mvwprintw(output->win,offsy+1,offsx+1,"pos1:"); + mvwprintw(output->win,offsy+4,offsx+1,"best match:"); + mvwprintw(output->win,offsy+5,offsx+1,"at"); + if (smallwin) mvwprintw(smallwin,4,1,"%18lli",0); + wrefresh(output->win); + } + + correlation->bestmatch=0; + correlation->mindiff=MAX(buf1->bufsize,buf2->bufsize); + for (pos1=-(tInt64)buf2->bufsize;pos1<(tInt64)buf1->bufsize;pos1++) + { + match=0; + diff=0; + if (output) mvwprintw(smallwin,1,3,"%16lli",((tInt64)buf1->bufsize-pos1)); + else fprintf(stderr,"%16lli\r",((tInt64)buf1->bufsize-pos1)); + for (pos2=0;pos2<buf2->bufsize;pos2++) + { + sumpos=pos1+pos2; + if (sumpos>=0 && sumpos<buf1->bufsize) + { + idx1=getbufferidx(buf1,sumpos); + idx2=getbufferidx(buf2,pos2); + + match+=(buf1->data[idx1]==buf2->data[idx2]); + diff+=(buf1->data[idx1]!=buf2->data[idx2]); + } else { + diff++; + } + } + if (match>correlation->bestmatch) + { + correlation->bestmatch=match; + correlation->bestmatchpos=pos1; + correlation->mindiff=diff; + if (smallwin) mvwprintw(smallwin,4,1,"%18lli",match); + if (smallwin) mvwprintw(smallwin,5,3,"%16llx",pos1); + found=1; + } + if (output) wrefresh(smallwin); + } + if (found) + { + if (correlation->bestmatch>=0) + { + *cursorpos1=correlation->bestmatchpos; + *cursorpos2=0; + } else { + *cursorpos1=0; + *cursorpos2=-correlation->bestmatchpos; + } + } + if (smallwin) delwin(smallwin); +} +void find_mindiff(tOutput* output,tCorrelation* correlation,tBuffer* buf1,tBuffer* buf2,tInt64* cursorpos1,tInt64* cursorpos2) +{ + tInt64 idx1; + tInt64 idx2; + tInt64 diff; + tInt64 sumpos; + tInt64 pos1; + tInt64 pos2; + tInt16 offsx; + tInt16 offsy; + WINDOW* smallwin=NULL; + tBuffer* smallbuf; + tBuffer* bigbuf; + tInt64 match; + tBool found; + + + + if (output) + { + offsx=COLS/2-22; + offsy=LINES/2-3; + drawcenterframe(output,6,44,"correlating..."); + setcolor(output,COLOR_TEXT); + smallwin=subwin(output->win,6,19,offsy,offsx+24); + mvwprintw(output->win,offsy+1,offsx+1,"pos1:"); + mvwprintw(output->win,offsy+4,offsx+1,"mindiff:"); + mvwprintw(output->win,offsy+5,offsx+1,"at"); + mvwprintw(smallwin,4,1,"%18lli",correlation->start_mindiff); + wrefresh(output->win); + } + + smallbuf=(buf1->bufsize<buf2->bufsize)?buf1:buf2; + bigbuf =(buf1->bufsize<buf2->bufsize)?buf2:buf1; + + correlation->mindiff=correlation->start_mindiff; + correlation->bestmatch=0; + if (!correlation->mindiff) correlation->mindiff=bigbuf->bufsize; + + for (pos1=-(tInt64)smallbuf->bufsize;pos1<(tInt64)bigbuf->bufsize && correlation->mindiff;pos1++) + { + diff=0; + match=0; + if (smallwin) mvwprintw(smallwin,1,3,"%16lli",((tInt64)bigbuf->bufsize-pos1)); + else fprintf(stderr,"%16lli\r",((tInt64)buf1->bufsize-pos1)); + + for (pos2=0;pos2<smallbuf->bufsize && diff<correlation->mindiff;pos2++) + { + sumpos=pos1+pos2; + if (sumpos>=0 && sumpos<bigbuf->bufsize) + { + idx1=getbufferidx(bigbuf,sumpos); + idx2=getbufferidx(smallbuf,pos2); + + diff+=(bigbuf->data[idx1]!=smallbuf->data[idx2]); + } else { + diff++; + } + } + if (diff<correlation->mindiff) + { + correlation->mindiffpos=pos1; + correlation->mindiff=diff; + if (smallwin) mvwprintw(smallwin,4,1,"%18lli",diff); + if (smallwin) mvwprintw(smallwin,5,3,"%16llx",pos1); + found=1; + } + if (smallwin) wrefresh(smallwin); + } + if (found) + { + if (correlation->mindiffpos>=0) + { + *cursorpos1=(bigbuf==buf1)?correlation->mindiffpos:0; + *cursorpos2=(bigbuf==buf1)?0:correlation->mindiffpos; + } else { + *cursorpos1=(bigbuf==buf1)?0:-correlation->mindiffpos; + *cursorpos2=(bigbuf==buf1)?-correlation->mindiffpos:0; + } + } + if (smallwin) delwin(smallwin); +} +void find_correlation(tOutput* output,tCorrelation* correlation,tBuffer* buf1,tBuffer* buf2,tInt64* cursorpos1,tInt64* cursorpos2) +{ + tInt16 ch=0; + correlation->correlated=1; + switch(correlation->algorithm) + { + case CORR_BEST_MATCH: find_bestmatch(output,correlation,buf1,buf2,cursorpos1,cursorpos2);break; + case CORR_LONGEST_MATCH: find_longestmatch(output,correlation,buf1,buf2,cursorpos1,cursorpos2);break; + case CORR_MIN_DIFF: find_mindiff(output,correlation,buf1,buf2,cursorpos1,cursorpos2);break; + } + if (output) + { + while (ch!=ERR) + { + ch=getch(); // remove all the keys which have been pressed from the buffer + } + } +} +void clear_correlation(tCorrelation* correlation) +{ + memset(correlation,0,sizeof(tCorrelation)); + correlation->algorithm=CORR_BEST_MATCH; +} |