diff --git a/README b/README index 7b56d13..2e33bb6 100644 --- a/README +++ b/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 diff --git a/src/mcio.c b/src/mcio.c index 4441e8d..7b26af0 100644 --- a/src/mcio.c +++ b/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) {