1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
|
# slow5lib
## NAME
slow5_get_next_mem - Gets next slow5 record from current file pointer as it is stored in s5p with length stored at pointer n.
## SYNOPSYS
`void *slow5_get_next_mem(size_t *n, const struct slow5_file *s5p)`
## DESCRIPTION
This function returns the next slow5 record as it is in a void buffer.
Returned buffer should be freed using `free()`
## RETURN VALUE
Upon successful completion, `slow5_get_next_mem()` returns a *void* pointer. Otherwise, NULL is returned, e.g., if *s5p* is NULL.
## NOTES
On error following errors are set
* slow5_errno errors:
* SLOW5_ERR_ARG
* SLOW5_ERR_IO
* SLOW5_ERR_EOF end of file reached (and blow5 eof marker found)
* slow5_is_eof errors:
* SLOW5_ERR_TRUNC
* SLOW5_ERR_MEM
* SLOW5_ERR_UNK
## EXAMPLES
```
#include <stdio.h>
#include <stdlib.h>
#include <slow5/slow5.h>
#define FILE_PATH "examples/example.slow5"
#define TO_PICOAMPS(RAW_VAL,DIGITISATION,OFFSET,RANGE) (((RAW_VAL)+(OFFSET))*((RANGE)/(DIGITISATION)))
int main(){
slow5_file_t *sp = slow5_open_with(FILE_PATH,"r",SLOW5_FORMAT_ASCII);
if(sp==NULL){
fprintf(stderr,"Error in opening file\n");
exit(EXIT_FAILURE);
}
slow5_rec_t *rec = NULL;
size_t bytes;
char *mem;
if (!(mem = (char *) slow5_get_next_mem(&bytes, from))) {
return EXIT_FAILURE;
}
if (slow5_rec_depress_parse(&mem, &bytes, NULL, &rec, &from) != 0) {
exit(EXIT_FAILURE);
} else {
free(mem);
}
printf("%s\t",rec->read_id);
uint64_t len_raw_signal = rec->len_raw_signal;
for(uint64_t i=0;i<len_raw_signal;i++){
double pA = TO_PICOAMPS(rec->raw_signal[i],rec->digitisation,rec->offset,rec->range);
printf("%f ",pA);
}
printf("\n");
slow5_rec_free(rec);
slow5_close(sp);
}
```
|
