Lines Matching refs:ndiff

63     uint32_t *ndiff,*ncnt,ncmp, npairs;  member
128 …uint32_t ndiff,rid,pos,rand; // rand is to shuffle sites with the same ndiff from across all chrom… member
137 fprintf(stderr,"%s:%d\t%d\t",bcf_hdr_id2name(args->qry_hdr,ds->rid),ds->pos+1,ds->ndiff); in diff_sites_debug_print()
146 if ( a->ndiff < b->ndiff ) return 1; // descending order in diff_sites_cmp()
147 if ( a->ndiff > b->ndiff ) return -1; in diff_sites_cmp()
184 static inline void diff_sites_push(args_t *args, int ndiff, int rid, int pos) in diff_sites_push() argument
188 dat->ndiff = ndiff; in diff_sites_push()
195 static inline int diff_sites_shift(args_t *args, int *ndiff, int *rid, int *pos) in diff_sites_shift() argument
199 *ndiff = dat->ndiff; in diff_sites_shift()
399 …args->ndiff = (uint32_t*) calloc(args->npairs,sizeof(*args->ndiff)); // number of differing gen… in init_data()
429 free(args->ndiff); in destroy_data()
598 int ndiff = 0; in process_line() local
618 args->ndiff[i]++; in process_line()
619 if ( args->kbs_diff ) { ndiff++; kbs_insert(args->kbs_diff, i); } in process_line()
625 if ( ndiff ) diff_sites_push(args, ndiff, qry_rec->rid, qry_rec->pos); in process_line()
689 if ( !match ) args->ndiff[idx]++; in process_line()
767 int i,ndiff,rid,pos,ndiff_tot = 0, iblock = 0; in report_distinctive_sites() local
769 while ( diff_sites_shift(args,&ndiff,&rid,&pos) ) in report_distinctive_sites()
780 if ( ndiff_dbg!=ndiff ) error("Corrupted data, fixme: %d vs %d\n",ndiff_dbg,ndiff); in report_distinctive_sites()
822 if ( args->ndiff ) in report()
827 args->ndiff[i], in report()
852 if ( args->ndiff ) in report()
857 args->ndiff[idx], in report()
885 else if ( args->ndiff ) in report()
886 arr[j].val = args->ncnt[idx] ? (double)args->ndiff[idx]/args->ncnt[idx] : 0; in report()
899 if ( args->ndiff ) in report()
904 args->ndiff[idx], in report()
933 else if ( args->ndiff ) in report()
934 arr[k].val = args->ncnt[idx] ? (double)args->ndiff[idx]/args->ncnt[idx] : 0; in report()
947 else if ( args->ndiff ) in report()
948 arr[k].val = args->ncnt[idx] ? (double)args->ndiff[idx]/args->ncnt[idx] : 0; in report()
962 if ( args->ndiff ) in report()
967 args->ndiff[idx], in report()