[rsync/rsync.git] / match.c

/* 
   Copyright (C) Andrew Tridgell 1996
   Copyright (C) Paul Mackerras 1996
   
   This program is free software; you can redistribute it and/or modify
   it under the terms of the GNU General Public License as published by
   the Free Software Foundation; either version 2 of the License, or
   (at your option) any later version.
   
   This program is distributed in the hope that it will be useful,
   but WITHOUT ANY WARRANTY; without even the implied warranty of
   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
   GNU General Public License for more details.
   
   You should have received a copy of the GNU General Public License
   along with this program; if not, write to the Free Software
   Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
*/

#include "rsync.h"

extern int verbose;
extern int am_server;
extern int do_progress;

typedef unsigned short tag;

#define TABLESIZE (1<<16)
#define NULL_TAG (-1)

static int false_alarms;
static int tag_hits;
static int matches;
static int64 data_transfer;

static int total_false_alarms;
static int total_tag_hits;
static int total_matches;

extern struct stats stats;

struct target {
  tag t;
  int i;
};

static struct target *targets;

static int *tag_table;

#define gettag2(s1,s2) (((s1) + (s2)) & 0xFFFF)
#define gettag(sum) gettag2((sum)&0xFFFF,(sum)>>16)

static int compare_targets(struct target *t1,struct target *t2)
{
  return((int)t1->t - (int)t2->t);
}


static void build_hash_table(struct sum_struct *s)
{
  int i;

  if (!tag_table)
    tag_table = (int *)malloc(sizeof(tag_table[0])*TABLESIZE);

  targets = (struct target *)malloc(sizeof(targets[0])*s->count);
  if (!tag_table || !targets) 
    out_of_memory("build_hash_table");

  for (i=0;i<(int) s->count;i++) {
    targets[i].i = i;
    targets[i].t = gettag(s->sums[i].sum1);
  }

  qsort(targets,s->count,sizeof(targets[0]),(int (*)())compare_targets);

  for (i=0;i<TABLESIZE;i++)
    tag_table[i] = NULL_TAG;

  for (i=s->count-1;i>=0;i--) {    
    tag_table[targets[i].t] = i;
  }
}


static OFF_T last_match;


/**
 * Transmit a literal and/or match token.
 *
 * This delightfully-named function is called either when we find a
 * match and need to transmit all the unmatched data leading up to it,
 * or when we get bored of accumulating literal data and just need to
 * transmit it.  As a result of this second case, it is called even if
 * we have not matched at all!
 *
 * @param i If >0, the number of a matched token.  If 0, indicates we
 * have only literal data.
 **/
static void matched(int f,struct sum_struct *s,struct map_struct *buf,
		    OFF_T offset,int i)
{
	OFF_T n = offset - last_match;
	OFF_T j;

	if (verbose > 2 && i >= 0)
		rprintf(FINFO,"match at %.0f last_match=%.0f j=%d len=%d n=%.0f\n",
			(double)offset,(double)last_match,i,s->sums[i].len,(double)n);

	send_token(f,i,buf,last_match,n,i<0?0:s->sums[i].len);
	data_transfer += n;

	if (i >= 0) {
		stats.matched_data += s->sums[i].len;
		n += s->sums[i].len;
	}
  
	for (j=0;j<n;j+=CHUNK_SIZE) {
		int n1 = MIN(CHUNK_SIZE,n-j);
		sum_update(map_ptr(buf,last_match+j,n1),n1);
	}


	if (i >= 0)
		last_match = offset + s->sums[i].len;
	else
		last_match = offset;

	if (buf && do_progress) {
		show_progress(last_match, buf->file_size);

		if (i == -1) end_progress(buf->file_size);
	}
}


static void hash_search(int f,struct sum_struct *s,
			struct map_struct *buf,OFF_T len)
{
	OFF_T offset, end;
	int j,k, last_i;
	char sum2[SUM_LENGTH];
	uint32 s1, s2, sum; 
	schar *map;

	/* last_i is used to encourage adjacent matches, allowing the RLL coding of the
	   output to work more efficiently */
	last_i = -1;

	if (verbose > 2)
		rprintf(FINFO,"hash search b=%ld len=%.0f\n",
			(long) s->blength, (double)len);

	/* cast is to make s->blength signed; it should always be reasonably
	 * small */
	k = MIN(len, (OFF_T) s->blength);
	
	map = (schar *)map_ptr(buf,0,k);
	
	sum = get_checksum1((char *)map, k);
	s1 = sum & 0xFFFF;
	s2 = sum >> 16;
	if (verbose > 3)
		rprintf(FINFO, "sum=%.8x k=%d\n", sum, k);
	
	offset = 0;
	
	end = len + 1 - s->sums[s->count-1].len;
	
	if (verbose > 3)
		rprintf(FINFO, "hash search s->blength=%ld len=%.0f count=%ld\n",
			(long) s->blength, (double) len, (long) s->count);
	
	do {
		tag t = gettag2(s1,s2);
		int done_csum2 = 0;
			
		j = tag_table[t];
		if (verbose > 4)
			rprintf(FINFO,"offset=%.0f sum=%08x\n",(double)offset,sum);
		
		if (j == NULL_TAG) {
			goto null_tag;
		}

		sum = (s1 & 0xffff) | (s2 << 16);
		tag_hits++;
		for (; j < (int) s->count && targets[j].t == t; j++) {
			int l, i = targets[j].i;
			
			if (sum != s->sums[i].sum1) continue;
			
			/* also make sure the two blocks are the same length */
			l = MIN(s->blength,len-offset);
			if (l != s->sums[i].len) continue;			

			if (verbose > 3)
				rprintf(FINFO,"potential match at %.0f target=%d %d sum=%08x\n",
					(double)offset,j,i,sum);
			
			if (!done_csum2) {
				map = (schar *)map_ptr(buf,offset,l);
				get_checksum2((char *)map,l,sum2);
				done_csum2 = 1;
			}
			
			if (memcmp(sum2,s->sums[i].sum2,s->s2length) != 0) {
				false_alarms++;
				continue;
			}

			/* we've found a match, but now check to see
                           if last_i can hint at a better match */
			for (j++; j < (int) s->count && targets[j].t == t; j++) {
				int i2 = targets[j].i;
				if (i2 == last_i + 1) {
					if (sum != s->sums[i2].sum1) break;
					if (memcmp(sum2,s->sums[i2].sum2,s->s2length) != 0) break;
					/* we've found an adjacent match - the RLL coder 
					   will be happy */
					i = i2;
					break;
				}
			}

			last_i = i;
			
			matched(f,s,buf,offset,i);
			offset += s->sums[i].len - 1;
			k = MIN((len-offset), s->blength);
			map = (schar *)map_ptr(buf,offset,k);
			sum = get_checksum1((char *)map, k);
			s1 = sum & 0xFFFF;
			s2 = sum >> 16;
			matches++;
			break;
		}
		
	null_tag:
		/* Trim off the first byte from the checksum */
		map = (schar *)map_ptr(buf,offset,k+1);
		s1 -= map[0] + CHAR_OFFSET;
		s2 -= k * (map[0]+CHAR_OFFSET);
		
		/* Add on the next byte (if there is one) to the checksum */
		if (k < (len-offset)) {
			s1 += (map[k]+CHAR_OFFSET);
			s2 += s1;
		} else {
			--k;
		}

		/* By matching early we avoid re-reading the
		   data 3 times in the case where a token
		   match comes a long way after last
		   match. The 3 reads are caused by the
		   running match, the checksum update and the
		   literal send. */
		if (offset > last_match &&
		    offset-last_match >= CHUNK_SIZE+s->blength && 
		    (end-offset > CHUNK_SIZE)) {
			matched(f,s,buf,offset - s->blength, -2);
		}
	} while (++offset < end);
	
	matched(f,s,buf,len,-1);
	map_ptr(buf,len-1,1);
}


/**
 * Scan through a origin file, looking for sections that match
 * checksums from the generator, and transmit either literal or token
 * data.
 *
 * Also calculates the MD4 checksum of the whole file, using the md
 * accumulator.  This is transmitted with the file as protection
 * against corruption on the wire.
 *
 * @param s Checksums received from the generator.  If <tt>s->count ==
 * 0</tt>, then there are actually no checksums for this file.
 *
 * @param len Length of the file to send.
 **/
void match_sums(int f, struct sum_struct *s, struct map_struct *buf, OFF_T len)
{
	char file_sum[MD4_SUM_LENGTH];
	extern int write_batch;  /*  dw */

	last_match = 0;
	false_alarms = 0;
	tag_hits = 0;
	matches=0;
	data_transfer=0;

	sum_init();

	if (len > 0 && s->count>0) {
		build_hash_table(s);
		
		if (verbose > 2) 
			rprintf(FINFO,"built hash table\n");
		
		hash_search(f,s,buf,len);
		
		if (verbose > 2) 
			rprintf(FINFO,"done hash search\n");
	} else {
		OFF_T j;
		/* by doing this in pieces we avoid too many seeks */
		for (j=0;j<(len-CHUNK_SIZE);j+=CHUNK_SIZE) {
			int n1 = MIN(CHUNK_SIZE,(len-CHUNK_SIZE)-j);
			matched(f,s,buf,j+n1,-2);
		}
		matched(f,s,buf,len,-1);
	}

	sum_end(file_sum);

	if (verbose > 2)
		rprintf(FINFO,"sending file_sum\n");
	write_buf(f,file_sum,MD4_SUM_LENGTH);
	if (write_batch) /* dw */
		write_batch_delta_file(file_sum, MD4_SUM_LENGTH);

	if (targets) {
		free(targets);
		targets=NULL;
	}
	
	if (verbose > 2)
		rprintf(FINFO, "false_alarms=%d tag_hits=%d matches=%d\n",
			false_alarms, tag_hits, matches);
	
	total_tag_hits += tag_hits;
	total_false_alarms += false_alarms;
	total_matches += matches;
	stats.literal_data += data_transfer;
}

void match_report(void)
{
	if (verbose <= 1)
		return;

	rprintf(FINFO,
		"total: matches=%d  tag_hits=%d  false_alarms=%d data=%.0f\n",
		total_matches,total_tag_hits,
		total_false_alarms,
		(double)stats.literal_data);
}
Commit	Line	Data
c627d613 AT	1	/*
	2	Copyright (C) Andrew Tridgell 1996
	3	Copyright (C) Paul Mackerras 1996
	4
	5	This program is free software; you can redistribute it and/or modify
	6	it under the terms of the GNU General Public License as published by
	7	the Free Software Foundation; either version 2 of the License, or
	8	(at your option) any later version.
	9
	10	This program is distributed in the hope that it will be useful,
	11	but WITHOUT ANY WARRANTY; without even the implied warranty of
	12	MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
	13	GNU General Public License for more details.
	14
	15	You should have received a copy of the GNU General Public License
	16	along with this program; if not, write to the Free Software
	17	Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
	18	*/
	19
	20	#include "rsync.h"
	21
	22	extern int verbose;
	23	extern int am_server;
16417f8b	24	extern int do_progress;
c627d613 AT	25
	26	typedef unsigned short tag;
	27
	28	#define TABLESIZE (1<<16)
100e5241	29	#define NULL_TAG (-1)
c627d613 AT	30
	31	static int false_alarms;
	32	static int tag_hits;
	33	static int matches;
a800434a	34	static int64 data_transfer;
c627d613	35
3a6a366f AT	36	static int total_false_alarms;
	37	static int total_tag_hits;
	38	static int total_matches;
c627d613	39
a800434a	40	extern struct stats stats;
c627d613 AT	41
	42	struct target {
	43	tag t;
	44	int i;
	45	};
	46
3a6a366f	47	static struct target *targets;
c627d613	48
100e5241	49	static int *tag_table;
c627d613 AT	50
	51	#define gettag2(s1,s2) (((s1) + (s2)) & 0xFFFF)
	52	#define gettag(sum) gettag2((sum)&0xFFFF,(sum)>>16)
	53
	54	static int compare_targets(struct target t1,struct target t2)
	55	{
a16bbc39	56	return((int)t1->t - (int)t2->t);
c627d613 AT	57	}
	58
	59
	60	static void build_hash_table(struct sum_struct *s)
	61	{
	62	int i;
	63
	64	if (!tag_table)
100e5241	65	tag_table = (int )malloc(sizeof(tag_table[0])TABLESIZE);
c627d613 AT	66
	67	targets = (struct target )malloc(sizeof(targets[0])s->count);
	68	if (!tag_table \|\| !targets)
	69	out_of_memory("build_hash_table");
	70
a261989c	71	for (i=0;i<(int) s->count;i++) {
c627d613 AT	72	targets[i].i = i;
	73	targets[i].t = gettag(s->sums[i].sum1);
	74	}
	75
	76	qsort(targets,s->count,sizeof(targets[0]),(int (*)())compare_targets);
	77
	78	for (i=0;i<TABLESIZE;i++)
	79	tag_table[i] = NULL_TAG;
	80
	81	for (i=s->count-1;i>=0;i--) {
	82	tag_table[targets[i].t] = i;
	83	}
	84	}
	85
	86
bcacc18b	87	static OFF_T last_match;
c627d613 AT	88
c627d613 AT	89
f5f95a38 MP	90	/**
	91	* Transmit a literal and/or match token.
	92	*
	93	* This delightfully-named function is called either when we find a
	94	* match and need to transmit all the unmatched data leading up to it,
	95	* or when we get bored of accumulating literal data and just need to
	96	* transmit it. As a result of this second case, it is called even if
	97	* we have not matched at all!
	98	*
	99	* @param i If >0, the number of a matched token. If 0, indicates we
	100	* have only literal data.
	101	**/
c6e7fcb4	102	static void matched(int f,struct sum_struct s,struct map_struct buf,
bcacc18b	103	OFF_T offset,int i)
c627d613	104	{
bcacc18b	105	OFF_T n = offset - last_match;
1309d90d	106	OFF_T j;
9e31c482	107
0d0e2e93	108	if (verbose > 2 && i >= 0)
5f808dfb AT	109	rprintf(FINFO,"match at %.0f last_match=%.0f j=%d len=%d n=%.0f\n",
5f808dfb AT	110	(double)offset,(double)last_match,i,s->sums[i].len,(double)n);
c627d613	111
45f133b9	112	send_token(f,i,buf,last_match,n,i<0?0:s->sums[i].len);
e7ebc36c	113	data_transfer += n;
70d794dc	114
a800434a AT	115	if (i >= 0) {
a800434a AT	116	stats.matched_data += s->sums[i].len;
e7ebc36c	117	n += s->sums[i].len;
a800434a	118	}
70d794dc	119
e7ebc36c AT	120	for (j=0;j<n;j+=CHUNK_SIZE) {
	121	int n1 = MIN(CHUNK_SIZE,n-j);
	122	sum_update(map_ptr(buf,last_match+j,n1),n1);
	123	}
9e31c482	124
9e31c482	125
0d0e2e93	126	if (i >= 0)
e7ebc36c	127	last_match = offset + s->sums[i].len;
0d0e2e93 AT	128	else
0d0e2e93 AT	129	last_match = offset;
eb86d661	130
16417f8b	131	if (buf && do_progress) {
4440b8aa	132	show_progress(last_match, buf->file_size);
eb86d661	133
5e12ce11 AT	134	if (i == -1) end_progress(buf->file_size);
5e12ce11 AT	135	}
9e31c482 AT	136	}
	137
	138
c6e7fcb4	139	static void hash_search(int f,struct sum_struct *s,
bcacc18b	140	struct map_struct *buf,OFF_T len)
c627d613	141	{
5f808dfb	142	OFF_T offset, end;
923fa978	143	int j,k, last_i;
e7ebc36c AT	144	char sum2[SUM_LENGTH];
e7ebc36c AT	145	uint32 s1, s2, sum;
debb4505	146	schar *map;
e7ebc36c	147
923fa978 AT	148	/* last_i is used to encourage adjacent matches, allowing the RLL coding of the
	149	output to work more efficiently */
	150	last_i = -1;
	151
e7ebc36c	152	if (verbose > 2)
e35080ce	153	rprintf(FINFO,"hash search b=%ld len=%.0f\n",
fc0257c9	154	(long) s->blength, (double)len);
e7ebc36c	155
fc0257c9	156	/* cast is to make s->blength signed; it should always be reasonably
e35080ce	157	* small */
fc0257c9	158	k = MIN(len, (OFF_T) s->blength);
e7ebc36c	159
debb4505	160	map = (schar *)map_ptr(buf,0,k);
e7ebc36c AT	161
	162	sum = get_checksum1((char *)map, k);
	163	s1 = sum & 0xFFFF;
	164	s2 = sum >> 16;
	165	if (verbose > 3)
9486289c	166	rprintf(FINFO, "sum=%.8x k=%d\n", sum, k);
e7ebc36c AT	167
	168	offset = 0;
	169
	170	end = len + 1 - s->sums[s->count-1].len;
	171
	172	if (verbose > 3)
fc0257c9 S	173	rprintf(FINFO, "hash search s->blength=%ld len=%.0f count=%ld\n",
fc0257c9 S	174	(long) s->blength, (double) len, (long) s->count);
e7ebc36c AT	175
	176	do {
	177	tag t = gettag2(s1,s2);
	178	int done_csum2 = 0;
	179
	180	j = tag_table[t];
	181	if (verbose > 4)
5f808dfb	182	rprintf(FINFO,"offset=%.0f sum=%08x\n",(double)offset,sum);
e7ebc36c AT	183
	184	if (j == NULL_TAG) {
	185	goto null_tag;
	186	}
	187
	188	sum = (s1 & 0xffff) \| (s2 << 16);
	189	tag_hits++;
a261989c	190	for (; j < (int) s->count && targets[j].t == t; j++) {
536541d5	191	int l, i = targets[j].i;
e7ebc36c	192
536541d5	193	if (sum != s->sums[i].sum1) continue;
e7ebc36c	194
536541d5	195	/* also make sure the two blocks are the same length */
fc0257c9	196	l = MIN(s->blength,len-offset);
536541d5 AT	197	if (l != s->sums[i].len) continue;
536541d5 AT	198
e7ebc36c	199	if (verbose > 3)
5f808dfb AT	200	rprintf(FINFO,"potential match at %.0f target=%d %d sum=%08x\n",
5f808dfb AT	201	(double)offset,j,i,sum);
e7ebc36c AT	202
e7ebc36c AT	203	if (!done_csum2) {
debb4505	204	map = (schar *)map_ptr(buf,offset,l);
e7ebc36c AT	205	get_checksum2((char *)map,l,sum2);
	206	done_csum2 = 1;
	207	}
	208
fc0257c9	209	if (memcmp(sum2,s->sums[i].sum2,s->s2length) != 0) {
e7ebc36c AT	210	false_alarms++;
	211	continue;
	212	}
923fa978 AT	213
	214	/* we've found a match, but now check to see
	215	if last_i can hint at a better match */
a261989c	216	for (j++; j < (int) s->count && targets[j].t == t; j++) {
34d3eed4	217	int i2 = targets[j].i;
923fa978 AT	218	if (i2 == last_i + 1) {
923fa978 AT	219	if (sum != s->sums[i2].sum1) break;
fc0257c9	220	if (memcmp(sum2,s->sums[i2].sum2,s->s2length) != 0) break;
923fa978 AT	221	/* we've found an adjacent match - the RLL coder
	222	will be happy */
	223	i = i2;
	224	break;
	225	}
	226	}
	227
	228	last_i = i;
e7ebc36c AT	229
	230	matched(f,s,buf,offset,i);
	231	offset += s->sums[i].len - 1;
fc0257c9	232	k = MIN((len-offset), s->blength);
debb4505	233	map = (schar *)map_ptr(buf,offset,k);
e7ebc36c AT	234	sum = get_checksum1((char *)map, k);
	235	s1 = sum & 0xFFFF;
	236	s2 = sum >> 16;
	237	matches++;
	238	break;
	239	}
	240
	241	null_tag:
	242	/* Trim off the first byte from the checksum */
debb4505	243	map = (schar *)map_ptr(buf,offset,k+1);
e7ebc36c AT	244	s1 -= map[0] + CHAR_OFFSET;
	245	s2 -= k * (map[0]+CHAR_OFFSET);
	246
	247	/* Add on the next byte (if there is one) to the checksum */
	248	if (k < (len-offset)) {
	249	s1 += (map[k]+CHAR_OFFSET);
	250	s2 += s1;
	251	} else {
	252	--k;
	253	}
45f133b9	254
5914bf15 PM	255	/* By matching early we avoid re-reading the
	256	data 3 times in the case where a token
	257	match comes a long way after last
	258	match. The 3 reads are caused by the
	259	running match, the checksum update and the
	260	literal send. */
0e948031	261	if (offset > last_match &&
fc0257c9	262	offset-last_match >= CHUNK_SIZE+s->blength &&
5914bf15	263	(end-offset > CHUNK_SIZE)) {
fc0257c9	264	matched(f,s,buf,offset - s->blength, -2);
45f133b9	265	}
e7ebc36c AT	266	} while (++offset < end);
	267
	268	matched(f,s,buf,len,-1);
	269	map_ptr(buf,len-1,1);
c627d613 AT	270	}
	271
	272
f5f95a38 MP	273	/**
	274	* Scan through a origin file, looking for sections that match
	275	* checksums from the generator, and transmit either literal or token
	276	* data.
	277	*
538ba24f MP	278	* Also calculates the MD4 checksum of the whole file, using the md
	279	* accumulator. This is transmitted with the file as protection
	280	* against corruption on the wire.
	281	*
f5f95a38 MP	282	* @param s Checksums received from the generator. If <tt>s->count ==
	283	* 0</tt>, then there are actually no checksums for this file.
	284	*
	285	* @param len Length of the file to send.
	286	**/
	287	void match_sums(int f, struct sum_struct s, struct map_struct buf, OFF_T len)
c627d613	288	{
e7ebc36c	289	char file_sum[MD4_SUM_LENGTH];
6902ed17	290	extern int write_batch; /* dw */
e7ebc36c AT	291
	292	last_match = 0;
	293	false_alarms = 0;
	294	tag_hits = 0;
	295	matches=0;
	296	data_transfer=0;
	297
	298	sum_init();
	299
	300	if (len > 0 && s->count>0) {
	301	build_hash_table(s);
	302
	303	if (verbose > 2)
9486289c	304	rprintf(FINFO,"built hash table\n");
e7ebc36c AT	305
	306	hash_search(f,s,buf,len);
	307
	308	if (verbose > 2)
9486289c	309	rprintf(FINFO,"done hash search\n");
e7ebc36c	310	} else {
eb86d661 AT	311	OFF_T j;
	312	/* by doing this in pieces we avoid too many seeks */
	313	for (j=0;j<(len-CHUNK_SIZE);j+=CHUNK_SIZE) {
	314	int n1 = MIN(CHUNK_SIZE,(len-CHUNK_SIZE)-j);
	315	matched(f,s,buf,j+n1,-2);
	316	}
e7ebc36c AT	317	matched(f,s,buf,len,-1);
e7ebc36c AT	318	}
9e31c482	319
e7ebc36c	320	sum_end(file_sum);
c627d613	321
bc63ae3f S	322	if (verbose > 2)
	323	rprintf(FINFO,"sending file_sum\n");
	324	write_buf(f,file_sum,MD4_SUM_LENGTH);
	325	if (write_batch) /* dw */
	326	write_batch_delta_file(file_sum, MD4_SUM_LENGTH);
c627d613	327
e7ebc36c AT	328	if (targets) {
	329	free(targets);
	330	targets=NULL;
	331	}
	332
	333	if (verbose > 2)
9486289c	334	rprintf(FINFO, "false_alarms=%d tag_hits=%d matches=%d\n",
e7ebc36c AT	335	false_alarms, tag_hits, matches);
	336
	337	total_tag_hits += tag_hits;
	338	total_false_alarms += false_alarms;
	339	total_matches += matches;
a800434a	340	stats.literal_data += data_transfer;
c627d613 AT	341	}
	342
	343	void match_report(void)
	344	{
e7ebc36c AT	345	if (verbose <= 1)
e7ebc36c AT	346	return;
c627d613	347
9486289c	348	rprintf(FINFO,
a800434a	349	"total: matches=%d tag_hits=%d false_alarms=%d data=%.0f\n",
e7ebc36c	350	total_matches,total_tag_hits,
a800434a AT	351	total_false_alarms,
a800434a AT	352	(double)stats.literal_data);
c627d613	353	}