fix constants in currently unused snpsortit

fix overflow bug when reading transpose packed data files > int size

README

@@ -1,4 +1,4 @@
-ADMIXTOOLS version 8.0.1  8/15/24 (for Linux and Mac)
+ADMIXTOOLS version 8.0.2  9/06/24 (for Linux and Mac)
 
 See README.INSTALL for installation info. 
 
@@ -108,6 +108,7 @@ See ./transpose_info and ./convertf/README
 The key code for this was written by Matthew Mah 
 
 Version 8.0.1 is a minor bugfix release.
+Version 8.0.2 is a minor bugfix release, affecting transpose packed files > 2GB.
 
 ----------------------------
 SOFTWARE COPYRIGHT NOTICE AGREEMENT

src/mcio.c

@@ -126,8 +126,8 @@ snpsortit (int **spos, int *indx, int n)
   int i, base[3];
 
   base[0] = 1;
-  base[1] = 10 ^ 8;
-  base[2] = 10 ^ 9;
+  base[1] = pow(10.0, 8) ;
+  base[2] = pow(10.0, 9) ;
 
   ZALLOC (lkode, n, long);
   for (i = 0; i < n; i++) {
@@ -2954,7 +2954,7 @@ inpackt (char *gname, SNP ** snpm, Indiv ** indiv, int numsnps, int numind)
   int fdes, ret;
   char *packit, *pbuff;
   int trlen ; // length of input record
-  int offset ;
+  long offset ;
   int *indmap, x, numoutind ; 
   int numread = 0 ;
 
@@ -3031,7 +3031,7 @@ inpackt (char *gname, SNP ** snpm, Indiv ** indiv, int numsnps, int numind)
 // this is key main program must set Indiv array appropriately 
 // jnum = indmap[j] ;    
    jnum = j ; 
-   offset = j*trlen + GENO_HEADER_SIZE ; 
+   offset = ((long) j) * trlen + GENO_HEADER_SIZE ; 
    lseek(fdes, offset,  SEEK_SET);
    t = read (fdes, buff, trlen);
    if (t != trlen) {