diffreg.c

Upload User: acmefrp
Upload Date: 2010-03-06
Package Size: 23768k
Code Size: 24k
Category: OS Develop
Development Platform: C/C++
  1. #ifndef lint
  2. static char *sccsid = "@(#)diffreg.c 1.1 92/07/30 SMI; from UCB 4.16 86/03/29";
  3. #endif
  4. #include "diff.h"
  5. /*
  6.  * diff - compare two files.
  7.  */
  8. /*
  9.  * Uses an algorithm due to Harold Stone, which finds
  10.  * a pair of longest identical subsequences in the two
  11.  * files.
  12.  *
  13.  * The major goal is to generate the match vector J.
  14.  * J[i] is the index of the line in file1 corresponding
  15.  * to line i file0. J[i] = 0 if there is no
  16.  * such line in file1.
  17.  *
  18.  * Lines are hashed so as to work in core. All potential
  19.  * matches are located by sorting the lines of each file
  20.  * on the hash (called ``value''). In particular, this
  21.  * collects the equivalence classes in file1 together.
  22.  * Subroutine equiv replaces the value of each line in
  23.  * file0 by the index of the first element of its 
  24.  * matching equivalence in (the reordered) file1.
  25.  * To save space equiv squeezes file1 into a single
  26.  * array member in which the equivalence classes
  27.  * are simply concatenated, except that their first
  28.  * members are flagged by changing sign.
  29.  *
  30.  * Next the indices that point into member are unsorted into
  31.  * array class according to the original order of file0.
  32.  *
  33.  * The cleverness lies in routine stone. This marches
  34.  * through the lines of file0, developing a vector klist
  35.  * of "k-candidates". At step i a k-candidate is a matched
  36.  * pair of lines x,y (x in file0 y in file1) such that
  37.  * there is a common subsequence of length k
  38.  * between the first i lines of file0 and the first y 
  39.  * lines of file1, but there is no such subsequence for
  40.  * any smaller y. x is the earliest possible mate to y
  41.  * that occurs in such a subsequence.
  42.  *
  43.  * Whenever any of the members of the equivalence class of
  44.  * lines in file1 matable to a line in file0 has serial number 
  45.  * less than the y of some k-candidate, that k-candidate 
  46.  * with the smallest such y is replaced. The new 
  47.  * k-candidate is chained (via pred) to the current
  48.  * k-1 candidate so that the actual subsequence can
  49.  * be recovered. When a member has serial number greater
  50.  * that the y of all k-candidates, the klist is extended.
  51.  * At the end, the longest subsequence is pulled out
  52.  * and placed in the array J by unravel.
  53.  *
  54.  * With J in hand, the matches there recorded are
  55.  * check'ed against reality to assure that no spurious
  56.  * matches have crept in due to hashing. If they have,
  57.  * they are broken, and "jackpot" is recorded--a harmless
  58.  * matter except that a true match for a spuriously
  59.  * mated line may now be unnecessarily reported as a change.
  60.  *
  61.  * Much of the complexity of the program comes simply
  62.  * from trying to minimize core utilization and
  63.  * maximize the range of doable problems by dynamically
  64.  * allocating what is needed and reusing what is not.
  65.  * The core requirements for problems larger than somewhat
  66.  * are (in words) 2*length(file0) + length(file1) +
  67.  * 3*(number of k-candidates installed),  typically about
  68.  * 6n words for files of length n. 
  69.  */
  70. #define prints(s) fputs(s,stdout)
  71. FILE *input[2];
  72. FILE *fopen();
  73. struct cand {
  74. int x;
  75. int y;
  76. int pred;
  77. } cand;
  78. struct line {
  79. int serial;
  80. int value;
  81. } *file[2], line;
  82. int len[2];
  83. struct line *sfile[2]; /* shortened by pruning common prefix and suffix */
  84. int slen[2];
  85. int pref, suff; /* length of prefix and suffix */
  86. int *class; /* will be overlaid on file[0] */
  87. int *member; /* will be overlaid on file[1] */
  88. int *klist; /* will be overlaid on file[0] after class */
  89. struct cand *clist; /* merely a free storage pot for candidates */
  90. int clen = 0;
  91. int *Jvect; /* will be overlaid on class */
  92. long *ixold; /* will be overlaid on klist */
  93. long *ixnew; /* will be overlaid on file[1] */
  94. diffreg()
  95. {
  96. register int i, j;
  97. char *origfile1, *origfile2;
  98. FILE *f1, *f2;
  99. char buf1[BUFSIZ], buf2[BUFSIZ];
  100. if (hflag) {
  101. diffargv[0] = "diffh";
  102. execv(diffh, diffargv);
  103. fprintf(stderr, "diff: ");
  104. perror(diffh);
  105. status = 2; /* error */
  106. done();
  107. }
  108. origfile1 = file1;
  109. origfile2 = file2;
  110. if ((stb1.st_mode & S_IFMT) == S_IFDIR)
  111. origfile1 = file1 = splice(file1, file2);
  112. else if ((stb2.st_mode & S_IFMT) == S_IFDIR)
  113. origfile2 = file2 = splice(file2, file1);
  114. else if (!strcmp(file1, "-")) {
  115. if (!strcmp(file2, "-")) {
  116. fprintf(stderr, "diff: can't specify - -n");
  117. status = 2; /* error */
  118. done();
  119. }
  120. file1 = copytemp();
  121. } else if (!strcmp(file2, "-"))
  122. file2 = copytemp();
  123. if ((f1 = fopen(file1, "r")) == NULL) {
  124. fprintf(stderr, "diff: ");
  125. perror(file1);
  126. status = 2; /* error */
  127. done();
  128. }
  129. if ((f2 = fopen(file2, "r")) == NULL) {
  130. fprintf(stderr, "diff: ");
  131. perror(file2);
  132. fclose(f1);
  133. status = 2; /* error */
  134. done();
  135. }
  136. if (fstat(fileno(f1), &stb1) < 0) {
  137. fprintf(stderr, "diff: Can't fstat ");
  138. perror(file1);
  139. fclose(f1);
  140. fclose(f2);
  141. status = 2; /* error */
  142. done();
  143. }
  144. if (fstat(fileno(f2), &stb2) < 0) {
  145. fprintf(stderr, "diff: Can't fstat ");
  146. perror(file2);
  147. fclose(f1);
  148. fclose(f2);
  149. status = 2; /* error */
  150. done();
  151. }
  152. if (stb1.st_size != stb2.st_size)
  153. goto notsame;
  154. for (;;) {
  155. i = fread(buf1, 1, BUFSIZ, f1);
  156. j = fread(buf2, 1, BUFSIZ, f2);
  157. if (ferror(f1) || ferror(f2)) {
  158. fprintf(stderr, "diff: Error reading ");
  159. perror(ferror(f1) ? file1: file2);
  160. fclose(f1);
  161. fclose(f2);
  162. status = 2; /* error */
  163. done();
  164. }
  165. if (i != j)
  166. goto notsame;
  167. if (i == 0 && j == 0) {
  168. if (opt == D_IFDEF) {
  169. /*
  170.  * since files are the same all we need
  171.  * to do is print one of them out.
  172.  */
  173. rewind(f1);
  174. while (i = fread(buf1, 1, BUFSIZ, f1)) 
  175. fwrite(buf1, 1, i, stdout);
  176. }
  177. fclose(f1);
  178. fclose(f2);
  179. status = 0; /* files don't differ */
  180. goto same;
  181. }
  182. for (j = 0; j < i; j++)
  183. if (buf1[j] != buf2[j])
  184. goto notsame;
  185. }
  186. notsame:
  187. /*
  188.  * Files certainly differ at this point; set status accordingly
  189.  */
  190. status = 1;
  191. if (binaryfile(f1) || binaryfile(f2)) {
  192. if (ferror(f1) || ferror(f2)) {
  193. fprintf(stderr, "diff: Error reading ");
  194. perror(ferror(f1) ? file1: file2);
  195. fclose(f1);
  196. fclose(f2);
  197. status = 2; /* error */
  198. done();
  199. }
  200. printf("Binary files %s and %s differn", file1, file2);
  201. fclose(f1);
  202. fclose(f2);
  203. done();
  204. }
  205. prepare(0, f1, file1, origfile1);
  206. prepare(1, f2, file2, origfile2);
  207. fclose(f1);
  208. fclose(f2);
  209. prune();
  210. sort(sfile[0],slen[0]);
  211. sort(sfile[1],slen[1]);
  212. member = (int *)file[1];
  213. equiv(sfile[0], slen[0], sfile[1], slen[1], member);
  214. member = (int *)ralloc((char *)member,(slen[1]+2)*sizeof(int));
  215. class = (int *)file[0];
  216. unsort(sfile[0], slen[0], class);
  217. class = (int *)ralloc((char *)class,(slen[0]+2)*sizeof(int));
  218. klist = (int *)talloc((slen[0]+2)*sizeof(int));
  219. clist = (struct cand *)talloc(sizeof(cand));
  220. i = stone(class, slen[0], member, klist);
  221. free((char *)member);
  222. free((char *)class);
  223. Jvect = (int *)talloc((len[0]+2)*sizeof(int));
  224. unravel(klist[i]);
  225. free((char *)clist);
  226. free((char *)klist);
  227. ixold = (long *)talloc((len[0]+2)*sizeof(long));
  228. ixnew = (long *)talloc((len[1]+2)*sizeof(long));
  229. check();
  230. output();
  231. status = anychange;
  232. same:
  233. if (opt == D_CONTEXT && anychange == 0)
  234. printf("No differences encounteredn");
  235. done();
  236. }
  237. char *
  238. copytemp()
  239. {
  240. char buf[BUFSIZ];
  241. register int i, f;
  242. signal(SIGHUP,done);
  243. signal(SIGINT,done);
  244. signal(SIGPIPE,done);
  245. signal(SIGTERM,done);
  246. tempfile = mktemp("/tmp/dXXXXXX");
  247. f = creat(tempfile,0600);
  248. if (f < 0) {
  249. fprintf(stderr, "diff: ");
  250. perror(tempfile);
  251. done();
  252. }
  253. while ((i = read(0,buf,BUFSIZ)) > 0)
  254. if (write(f,buf,i) != i) {
  255. fprintf(stderr, "diff: Write error on ");
  256. if (i < 0)
  257. perror(tempfile);
  258. else
  259. fprintf(stderr, "%s: Premature EOF on outputn",
  260.     tempfile);
  261. done();
  262. }
  263. if (i < 0)
  264. readerr("standard input");
  265. close(f);
  266. return (tempfile);
  267. }
  268. char *
  269. splice(dir, file)
  270. char *dir, *file;
  271. {
  272. char *tail;
  273. char buf[BUFSIZ];
  274. if (!strcmp(file, "-")) {
  275. fprintf(stderr, "diff: can't specify - with other arg directoryn");
  276. done();
  277. }
  278. tail = rindex(file, '/');
  279. if (tail == 0)
  280. tail = file;
  281. else
  282. tail++;
  283. sprintf(buf, "%s/%s", dir, tail);
  284. return (savestr(buf));
  285. }
  286. prepare(i, fd, fname, origfname)
  287. int i;
  288. FILE *fd;
  289. char *fname;
  290. char *origfname;
  291. {
  292. register struct line *p;
  293. register j,h;
  294. fseek(fd, (long)0, 0);
  295. p = (struct line *)talloc(3*sizeof(line));
  296. for(j=0; h=readhash(fd,fname,origfname);) {
  297. p = (struct line *)ralloc((char *)p,(++j+3)*sizeof(line));
  298. p[j].value = h;
  299. }
  300. len[i] = j;
  301. file[i] = p;
  302. }
  303. prune()
  304. {
  305. register i,j;
  306. for(pref=0;pref<len[0]&&pref<len[1]&&
  307. file[0][pref+1].value==file[1][pref+1].value;
  308. pref++ ) ;
  309. for(suff=0;suff<len[0]-pref&&suff<len[1]-pref&&
  310. file[0][len[0]-suff].value==file[1][len[1]-suff].value;
  311. suff++) ;
  312. for(j=0;j<2;j++) {
  313. sfile[j] = file[j]+pref;
  314. slen[j] = len[j]-pref-suff;
  315. for(i=0;i<=slen[j];i++)
  316. sfile[j][i].serial = i;
  317. }
  318. }
  319. equiv(a,n,b,m,c)
  320. struct line *a, *b;
  321. int *c;
  322. {
  323. register int i, j;
  324. i = j = 1;
  325. while(i<=n && j<=m) {
  326. if(a[i].value <b[j].value)
  327. a[i++].value = 0;
  328. else if(a[i].value == b[j].value)
  329. a[i++].value = j;
  330. else
  331. j++;
  332. }
  333. while(i <= n)
  334. a[i++].value = 0;
  335. b[m+1].value = 0;
  336. j = 0;
  337. while(++j <= m) {
  338. c[j] = -b[j].serial;
  339. while(b[j+1].value == b[j].value) {
  340. j++;
  341. c[j] = b[j].serial;
  342. }
  343. }
  344. c[j] = -1;
  345. }
  346. stone(a,n,b,c)
  347. int *a;
  348. int *b;
  349. register int *c;
  350. {
  351. register int i, k,y;
  352. int j, l;
  353. int oldc, tc;
  354. int oldl;
  355. k = 0;
  356. c[0] = newcand(0,0,0);
  357. for(i=1; i<=n; i++) {
  358. j = a[i];
  359. if(j==0)
  360. continue;
  361. y = -b[j];
  362. oldl = 0;
  363. oldc = c[0];
  364. do {
  365. if(y <= clist[oldc].y)
  366. continue;
  367. if(clist[c[k]].y<y) /*quick look for typical case*/
  368. l = k+1;
  369. else
  370. l = search(c, k, y);
  371. if(l!=oldl+1)
  372. oldc = c[l-1];
  373. if(l<=k) {
  374. if(clist[c[l]].y <= y)
  375. continue;
  376. tc = c[l];
  377. c[l] = newcand(i,y,oldc);
  378. oldc = tc;
  379. oldl = l;
  380. } else {
  381. c[l] = newcand(i,y,oldc);
  382. k++;
  383. break;
  384. }
  385. } while((y=b[++j]) > 0);
  386. }
  387. return(k);
  388. }
  389. newcand(x,y,pred)
  390. {
  391. register struct cand *q;
  392. clist = (struct cand *)ralloc((char *)clist,++clen*sizeof(cand));
  393. q = clist + clen -1;
  394. q->x = x;
  395. q->y = y;
  396. q->pred = pred;
  397. return(clen-1);
  398. }
  399. search(c, k, y)
  400. register int *c;
  401. register y;
  402. {
  403. register int i, j, l;
  404. register int t;
  405. register struct cand *clistp = clist;
  406. i = 0;
  407. j = k+1;
  408. while (1) {
  409. l = i + j;
  410. if ((l >>= 1) <= i) 
  411. break;
  412. t = clistp[c[l]].y;
  413. if(t > y)
  414. j = l;
  415. else if(t < y)
  416. i = l;
  417. else
  418. return(l);
  419. }
  420. return(l+1);
  421. }
  422. unravel(p)
  423. {
  424. register int i;
  425. register struct cand *q;
  426. register *J = Jvect;
  427. for(i=0; i<=len[0]; i++)
  428. J[i] = i<=pref ? i:
  429. i>len[0]-suff ? i+len[1]-len[0]:
  430. 0;
  431. for(q=clist+p;q->y!=0;q=clist+q->pred)
  432. J[q->x+pref] = q->y+pref;
  433. }
  434. /* check does double duty:
  435. 1.  ferret out any fortuitous correspondences due
  436. to confounding by hashing (which result in "jackpot")
  437. 2.  collect random access indexes to the two files */
  438. check()
  439. {
  440. register int i, j;
  441. int jackpot;
  442. register long ctold, ctnew;
  443. register int c,d;
  444. register FILE *f0, *f1;
  445. register int *J = Jvect;
  446. if ((input[0] = fopen(file1,"r")) == NULL) {
  447. fprintf(stderr, "diff: ");
  448. perror(file1);
  449. done();
  450. }
  451. if ((input[1] = fopen(file2,"r")) == NULL) {
  452. fprintf(stderr, "diff: ");
  453. perror(file2);
  454. done();
  455. }
  456. f0 = input[0];
  457. f1 = input[1];
  458. j = 1;
  459. ixold[0] = ixnew[0] = 0;
  460. jackpot = 0;
  461. ctold = ctnew = 0;
  462. for(i=1;i<=len[0];i++) {
  463. if(J[i]==0) {
  464. ixold[i] = ctold += skipline(0);
  465. continue;
  466. }
  467. while(j<J[i]) {
  468. ixnew[j] = ctnew += skipline(1);
  469. j++;
  470. }
  471. if(bflag || wflag || iflag) {
  472. for(;;) {
  473. c = getc(f0);
  474. d = getc(f1);
  475. ctold++;
  476. ctnew++;
  477. if(bflag && isspace(c) && isspace(d)) {
  478. do {
  479. if(c=='n')
  480. break;
  481. ctold++;
  482. } while(isspace(c=getc(f0)));
  483. do {
  484. if(d=='n')
  485. break;
  486. ctnew++;
  487. } while(isspace(d=getc(f1)));
  488. } else if ( wflag ) {
  489. while( isspace(c) && c!='n' ) {
  490. c=getc(f0);
  491. ctold++;
  492. }
  493. while( isspace(d) && d!='n' ) {
  494. d=getc(f1);
  495. ctnew++;
  496. }
  497. }
  498. if(c==EOF || d==EOF) {
  499. if(ferror(f0))
  500. readerr(file1);
  501. if(ferror(f1))
  502. readerr(file2);
  503. if(c != d) {
  504. jackpot++;
  505. J[i] = 0;
  506. if(c!='n' && c!=EOF)
  507. ctold += skipline(0);
  508. if(d!='n' && d!=EOF)
  509. ctnew += skipline(1);
  510. break;
  511. }
  512. break;
  513. } else {
  514. if(chrtran[c] != chrtran[d]) {
  515. jackpot++;
  516. J[i] = 0;
  517. if(c!='n')
  518. ctold += skipline(0);
  519. if(d!='n')
  520. ctnew += skipline(1);
  521. break;
  522. }
  523. if(c=='n')
  524. break;
  525. }
  526. }
  527. } else {
  528. for(;;) {
  529. ctold++;
  530. ctnew++;
  531. if((c=getc(f0)) != (d=getc(f1))) {
  532. /* jackpot++; */
  533. J[i] = 0;
  534. if(c==EOF) {
  535. if(ferror(f0))
  536. readerr(file1);
  537. } else if(c!='n')
  538. ctold += skipline(0);
  539. if(d==EOF) {
  540. if(ferror(f1))
  541. readerr(file2);
  542. } else if(d!='n')
  543. ctnew += skipline(1);
  544. break;
  545. }
  546. if(c == EOF) {
  547. if(ferror(f0))
  548. readerr(file1);
  549. if(ferror(f1))
  550. readerr(file2);
  551. break;
  552. } else if(c=='n')
  553. break;
  554. }
  555. }
  556. ixold[i] = ctold;
  557. ixnew[j] = ctnew;
  558. j++;
  559. }
  560. for(;j<=len[1];j++) {
  561. ixnew[j] = ctnew += skipline(1);
  562. }
  563. fclose(f0);
  564. fclose(f1);
  565. /*
  566. if(jackpot)
  567. fprintf(stderr, "jackpotn");
  568. */
  569. }
  570. sort(a,n) /*shellsort CACM #201*/
  571. struct line *a;
  572. {
  573. struct line w;
  574. register int j,m;
  575. register struct line *ai;
  576. register struct line *aim;
  577. int k;
  578. if (n == 0)
  579. return;
  580. for(j=1;j<=n;j*= 2)
  581. m = 2*j - 1;
  582. for(m/=2;m!=0;m/=2) {
  583. k = n-m;
  584. for(j=1;j<=k;j++) {
  585. for(ai = &a[j]; ai > a; ai -= m) {
  586. aim = &ai[m];
  587. if(aim < ai)
  588. break; /*wraparound*/
  589. if(aim->value > ai[0].value ||
  590.    aim->value == ai[0].value &&
  591.    aim->serial > ai[0].serial)
  592. break;
  593. w.value = ai[0].value;
  594. ai[0].value = aim->value;
  595. aim->value = w.value;
  596. w.serial = ai[0].serial;
  597. ai[0].serial = aim->serial;
  598. aim->serial = w.serial;
  599. }
  600. }
  601. }
  602. }
  603. unsort(f, l, b)
  604. register struct line *f;
  605. int *b;
  606. {
  607. register int *a;
  608. register int i;
  609. a = (int *)talloc((l+1)*sizeof(int));
  610. for(i=1;i<=l;i++)
  611. a[f[i].serial] = f[i].value;
  612. for(i=1;i<=l;i++)
  613. b[i] = a[i];
  614. free((char *)a);
  615. }
  616. skipline(fnum)
  617. {
  618. register i, c;
  619. register FILE *f;
  620. f = input[fnum];
  621. for(i=1; (c=getc(f))!='n'; i++)
  622. if(c == EOF) {
  623. if(ferror(f))
  624. readerr(fnum == 0 ? file1 : file2);
  625. return(i);
  626. }
  627. return(i);
  628. }
  629. output()
  630. {
  631. int m;
  632. register int i0, i1, j1;
  633. int j0;
  634. register FILE *f;
  635. register int *J = Jvect;
  636. input[0] = fopen(file1,"r");
  637. input[1] = fopen(file2,"r");
  638. m = len[0];
  639. J[0] = 0;
  640. J[m+1] = len[1]+1;
  641. if(opt!=D_EDIT) for(i0=1;i0<=m;i0=i1+1) {
  642. while(i0<=m&&J[i0]==J[i0-1]+1) i0++;
  643. j0 = J[i0-1]+1;
  644. i1 = i0-1;
  645. while(i1<m&&J[i1+1]==0) i1++;
  646. j1 = J[i1+1]-1;
  647. J[i1] = j1;
  648. change(i0,i1,j0,j1);
  649. } else for(i0=m;i0>=1;i0=i1-1) {
  650. while(i0>=1&&J[i0]==J[i0+1]-1&&J[i0]!=0) i0--;
  651. j0 = J[i0+1]-1;
  652. i1 = i0+1;
  653. while(i1>1&&J[i1-1]==0) i1--;
  654. j1 = J[i1-1]+1;
  655. J[i1] = j1;
  656. change(i1,i0,j1,j0);
  657. }
  658. if(m==0)
  659. change(1,0,1,len[1]);
  660. if (opt==D_IFDEF) {
  661. f = input[0];
  662. for (;;) {
  663. #define c i0
  664. c = getc(f);
  665. if (c < 0)
  666. return;
  667. putchar(c);
  668. }
  669. #undef c
  670. }
  671. if (anychange && opt == D_CONTEXT)
  672. dump_context_vec();
  673. }
  674. /*
  675.  * The following struct is used to record change information when
  676.  * doing a "context" diff.  (see routine "change" to understand the
  677.  * highly mneumonic field names)
  678.  */
  679. struct context_vec {
  680. int a; /* start line in old file */
  681. int b; /* end line in old file */
  682. int c; /* start line in new file */
  683. int d; /* end line in new file */
  684. };
  685. struct context_vec *context_vec_start,
  686. *context_vec_end,
  687. *context_vec_ptr;
  688. #define MAX_CONTEXT 128
  689. /* indicate that there is a difference between lines a and b of the from file
  690.    to get to lines c to d of the to file.
  691.    If a is greater then b then there are no lines in the from file involved
  692.    and this means that there were lines appended (beginning at b).
  693.    If c is greater than d then there are lines missing from the to file.
  694. */
  695. change(a,b,c,d)
  696. {
  697. int ch;
  698. int lowa,upb,lowc,upd;
  699. struct stat stbuf;
  700. if (opt != D_IFDEF && a>b && c>d)
  701. return;
  702. if (anychange == 0) {
  703. anychange = 1;
  704. if(opt == D_CONTEXT) {
  705. printf("*** %s %s", file1, ctime(&stb1.st_mtime));
  706. printf("--- %s %s", file2, ctime(&stb2.st_mtime));
  707. context_vec_start = (struct context_vec *) 
  708. malloc(MAX_CONTEXT *
  709.    sizeof(struct context_vec));
  710. if (context_vec_start == NULL) {
  711. fprintf(stderr, "diff: ran out of memoryn");
  712. done();
  713. }
  714. context_vec_end = context_vec_start + MAX_CONTEXT;
  715. context_vec_ptr = context_vec_start - 1;
  716. }
  717. }
  718. if (a <= b && c <= d)
  719. ch = 'c';
  720. else
  721. ch = (a <= b) ? 'd' : 'a';
  722. if(opt == D_CONTEXT) {
  723. /*
  724.  * if this new change is within 'context' lines of
  725.  * the previous change, just add it to the change
  726.  * record.  If the record is full or if this
  727.  * change is more than 'context' lines from the previous
  728.  * change, dump the record, reset it & add the new change.
  729.  */
  730. if ( context_vec_ptr >= context_vec_end ||
  731.      ( context_vec_ptr >= context_vec_start &&
  732.        a > (context_vec_ptr->b + 2*context) &&
  733.        c > (context_vec_ptr->d + 2*context) ) )
  734. dump_context_vec();
  735. context_vec_ptr++;
  736. context_vec_ptr->a = a;
  737. context_vec_ptr->b = b;
  738. context_vec_ptr->c = c;
  739. context_vec_ptr->d = d;
  740. return;
  741. }
  742. switch (opt) {
  743. case D_NORMAL:
  744. case D_EDIT:
  745. range(a,b,",");
  746. putchar(a>b?'a':c>d?'d':'c');
  747. if(opt==D_NORMAL)
  748. range(c,d,",");
  749. putchar('n');
  750. break;
  751. case D_REVERSE:
  752. putchar(a>b?'a':c>d?'d':'c');
  753. range(a,b," ");
  754. putchar('n');
  755. break;
  756.         case D_NREVERSE:
  757.                 if (a>b)
  758.                         printf("a%d %dn",b,d-c+1);
  759.                 else {
  760.                         printf("d%d %dn",a,b-a+1);
  761.                         if (!(c>d))
  762.                            /* add changed lines */
  763.                            printf("a%d %dn",b, d-c+1);
  764.                 }
  765.                 break;
  766. }
  767. if(opt == D_NORMAL || opt == D_IFDEF) {
  768. fetch(ixold,a,b,input[0],file1,"< ", 1);
  769. if(a<=b&&c<=d && opt == D_NORMAL)
  770. prints("---n");
  771. }
  772. fetch(ixnew,c,d,input[1],file2,opt==D_NORMAL?"> ":"", 0);
  773. if ((opt ==D_EDIT || opt == D_REVERSE) && c<=d)
  774. prints(".n");
  775. if (inifdef) {
  776. fprintf(stdout, "#endif %sn", endifname);
  777. inifdef = 0;
  778. }
  779. }
  780. range(a,b,separator)
  781. char *separator;
  782. {
  783. printf("%d", a>b?b:a);
  784. if(a<b) {
  785. printf("%s%d", separator, b);
  786. }
  787. }
  788. fetch(f,a,b,lb,filename,s,oldfile)
  789. long *f;
  790. register FILE *lb;
  791. char *filename;
  792. char *s;
  793. {
  794. register int i, j;
  795. register int c;
  796. register int col;
  797. register int nc;
  798. int oneflag = (*ifdef1!='') != (*ifdef2!='');
  799. /*
  800.  * When doing #ifdef's, copy down to current line
  801.  * if this is the first file, so that stuff makes it to output.
  802.  */
  803. if (opt == D_IFDEF && oldfile){
  804. long curpos = ftell(lb);
  805. /* print through if append (a>b), else to (nb: 0 vs 1 orig) */
  806. nc = f[a>b? b : a-1 ] - curpos;
  807. for (i = 0; i < nc; i++) {
  808. c = getc(lb);
  809. if(c != EOF)
  810. putchar(c);
  811. else {
  812. if(ferror(lb))
  813. readerr(filename);
  814. putchar('n');
  815. break;
  816. }
  817. }
  818. }
  819. if (a > b)
  820. return;
  821. if (opt == D_IFDEF) {
  822. if (inifdef)
  823. fprintf(stdout, "#else %s%sn", oneflag && oldfile==1 ? "!" : "", ifdef2);
  824. else {
  825. if (oneflag) {
  826. /* There was only one ifdef given */
  827. endifname = ifdef2;
  828. if (oldfile)
  829. fprintf(stdout, "#ifndef %sn", endifname);
  830. else
  831. fprintf(stdout, "#ifdef %sn", endifname);
  832. }
  833. else {
  834. endifname = oldfile ? ifdef1 : ifdef2;
  835. fprintf(stdout, "#ifdef %sn", endifname);
  836. }
  837. }
  838. inifdef = 1+oldfile;
  839. }
  840. for(i=a;i<=b;i++) {
  841. fseek(lb,f[i-1],0);
  842. if (opt != D_IFDEF)
  843. prints(s);
  844. col = 0;
  845. while(c=getc(lb)) {
  846. if (c == EOF) {
  847. if (ferror(lb))
  848. readerr(filename);
  849. break;
  850. } else if (c == 'n')
  851. break;
  852. else if (c == 't' && tflag)
  853. do
  854. putchar(' ');
  855. while (++col & 7);
  856. else {
  857. putchar(c);
  858. col++;
  859. }
  860. }
  861. putchar('n');
  862. }
  863. if (inifdef && !wantelses) {
  864. fprintf(stdout, "#endif %sn", endifname);
  865. inifdef = 0;
  866. }
  867. }
  868. #define POW2 /* define only if HALFLONG is 2**n */
  869. #define HALFLONG 16
  870. #define low(x) (x&((1L<<HALFLONG)-1))
  871. #define high(x) (x>>HALFLONG)
  872. /*
  873.  * hashing has the effect of
  874.  * arranging line in 7-bit bytes and then
  875.  * summing 1-s complement in 16-bit hunks 
  876.  */
  877. readhash(f,fname,origfname)
  878. register FILE *f;
  879. char *fname;
  880. char *origfname;
  881. {
  882. register long sum;
  883. register unsigned shift;
  884. register t;
  885. register space;
  886. sum = 1;
  887. space = 0;
  888. if(!bflag && !wflag) {
  889. if(iflag)
  890. for(shift=0;(t=getc(f))!='n';shift+=7) {
  891. if(t==EOF) {
  892. if(ferror(f))
  893. readerr(fname);
  894. if(shift) {
  895. fprintf(stderr,"Warning: missing newline at end of file %sn", origfname);
  896. break;
  897. } else
  898. return(0);
  899. }
  900. sum += (long)chrtran[t] << (shift
  901. #ifdef POW2
  902.     &= HALFLONG - 1);
  903. #else
  904.     %= HALFLONG);
  905. #endif
  906. }
  907. else
  908. for(shift=0;(t=getc(f))!='n';shift+=7) {
  909. if(t==EOF) {
  910. if(ferror(f))
  911. readerr(fname);
  912. if(shift) {
  913. fprintf(stderr,"Warning: missing newline at end of file %sn", origfname);
  914. break;
  915. } else
  916. return(0);
  917. }
  918. sum += (long)t << (shift
  919. #ifdef POW2
  920.     &= HALFLONG - 1);
  921. #else
  922.     %= HALFLONG);
  923. #endif
  924. }
  925. } else {
  926. for(shift=0;;) {
  927. switch(t=getc(f)) {
  928. case EOF:
  929. if(ferror(f))
  930. readerr(fname);
  931. if(shift) {
  932. fprintf(stderr,"Warning: missing newline at end of file %sn", origfname);
  933. break;
  934. } else
  935. return(0);
  936. case 't':
  937. case ' ':
  938. space++;
  939. continue;
  940. default:
  941. if(space && !wflag) {
  942. shift += 7;
  943. space = 0;
  944. }
  945. sum += (long)chrtran[t] << (shift
  946. #ifdef POW2
  947.     &= HALFLONG - 1);
  948. #else
  949.     %= HALFLONG);
  950. #endif
  951. shift += 7;
  952. continue;
  953. case 'n':
  954. break;
  955. }
  956. break;
  957. }
  958. }
  959. sum = low(sum) + high(sum);
  960. return((short)low(sum) + (short)high(sum));
  961. }
  962. binaryfile(f)
  963. FILE *f;
  964. {
  965. char buf[BUFSIZ];
  966. int cnt;
  967. fseek(f, (long)0, 0);
  968. cnt = fread(buf, 1, BUFSIZ, f);
  969. if (ferror(f))
  970. return (1);
  971. return (isbinary(buf, cnt));
  972. }
  973. /* dump accumulated "context" diff changes */
  974. dump_context_vec()
  975. {
  976. register int a, b, c, d;
  977. register char ch;
  978. register struct context_vec *cvp = context_vec_start;
  979. register int lowa, upb, lowc, upd;
  980. register int do_output;
  981. if ( cvp > context_vec_ptr )
  982. return;
  983. lowa = max(1, cvp->a - context);
  984. upb  = min(len[0], context_vec_ptr->b + context);
  985. lowc = max(1, cvp->c - context);
  986. upd  = min(len[1], context_vec_ptr->d + context);
  987. printf("***************n*** ");
  988. range(lowa,upb,",");
  989. printf(" ****n");
  990. /*
  991.  * output changes to the "old" file.  The first loop suppresses
  992.  * output if there were no changes to the "old" file (we'll see
  993.  * the "old" lines as context in the "new" list).
  994.  */
  995. do_output = 0;
  996. for ( ; cvp <= context_vec_ptr; cvp++)
  997. if (cvp->a <= cvp->b) {
  998. cvp = context_vec_start;
  999. do_output++;
  1000. break;
  1001. }
  1002. if ( do_output ) {
  1003. while (cvp <= context_vec_ptr) {
  1004. a = cvp->a; b = cvp->b; c = cvp->c; d = cvp->d;
  1005. if (a <= b && c <= d)
  1006. ch = 'c';
  1007. else
  1008. ch = (a <= b) ? 'd' : 'a';
  1009. if (ch == 'a')
  1010. fetch(ixold,lowa,b,input[0],file1,"  ");
  1011. else {
  1012. fetch(ixold,lowa,a-1,input[0],file1,"  ");
  1013. fetch(ixold,a,b,input[0],file1,ch == 'c' ? "! " : "- ");
  1014. }
  1015. lowa = b + 1;
  1016. cvp++;
  1017. }
  1018. fetch(ixold, b+1, upb, input[0], file1, "  ");
  1019. }
  1020. /* output changes to the "new" file */
  1021. printf("--- ");
  1022. range(lowc,upd,",");
  1023. printf(" ----n");
  1024. do_output = 0;
  1025. for (cvp = context_vec_start; cvp <= context_vec_ptr; cvp++)
  1026. if (cvp->c <= cvp->d) {
  1027. cvp = context_vec_start;
  1028. do_output++;
  1029. break;
  1030. }
  1031. if (do_output) {
  1032. while (cvp <= context_vec_ptr) {
  1033. a = cvp->a; b = cvp->b; c = cvp->c; d = cvp->d;
  1034. if (a <= b && c <= d)
  1035. ch = 'c';
  1036. else
  1037. ch = (a <= b) ? 'd' : 'a';
  1038. if (ch == 'd')
  1039. fetch(ixnew,lowc,d,input[1],file2,"  ");
  1040. else {
  1041. fetch(ixnew,lowc,c-1,input[1],file2,"  ");
  1042. fetch(ixnew,c,d,input[1],file2,ch == 'c' ? "! " : "+ ");
  1043. }
  1044. lowc = d + 1;
  1045. cvp++;
  1046. }
  1047. fetch(ixnew, d+1, upd, input[1], file2, "  ");
  1048. }
  1049. context_vec_ptr = context_vec_start - 1;
  1050. }
  1051. readerr(filename)
  1052. char *filename;
  1053. {
  1054. fprintf(stderr, "diff: Error reading ");
  1055. perror(filename);
  1056. done();
  1057. }