forked from chhylp123/hifiasm
-
Notifications
You must be signed in to change notification settings - Fork 8
/
Copy pathmeta_util.h
22 lines (17 loc) · 872 Bytes
/
meta_util.h
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
#ifndef __META_UTILITY__
#define __META_UTILITY__
#include <stdint.h>
#include "Process_Read.h"
typedef struct hamt_peaks_s{
uint64_t n, max, min, mean;
} hamt_peaks_t;
double meanl(const uint16_t *counts, uint32_t l);
double stdl(const uint16_t *counts, uint32_t l, double mean);
uint64_t medianl(const uint16_t *counts, uint32_t l);
uint16_t quantile(uint16_t *d, uint32_t l, float q);
void moving_averagel(uint16_t *counts, uint32_t l, uint16_t w, uint64_t *buf);
// hamt_peaks_t hamt_find_peaks(uint32_t n_cnt, const int64_t *cnt, uint16_t l_smooth);
// uint8_t decide_category(double mean, double std, uint16_t *buf, uint32_t l);
uint8_t decide_drop(double mean, double std, uint16_t runtime_median, int round, uint8_t initmark);
void hamt_del_ovlp_by_coverage(All_reads *rs, hifiasm_opt_t asm_opt, int diff_abs, double diff_fold);
#endif //__META_UTILITY__