[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 csum_length;

extern int verbose;
extern int am_server;

extern int remote_version;

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) {
		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->n, (double)len);

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