File: slow5_get_next.md

package info (click to toggle)
libslow5lib 0.7.0%2Bdfsg-3
  • links: PTS, VCS
  • area: main
  • in suites: sid, trixie
  • size: 25,092 kB
  • sloc: ansic: 11,825; python: 1,179; sh: 547; makefile: 91; cpp: 40
file content (90 lines) | stat: -rw-r--r-- 3,158 bytes parent folder | download | duplicates (2)
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
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
# slow5_get_next

## NAME
slow5_get_next - fetches the record at the current file pointer of a slow5 file

## SYNOPSYS
`int slow5_get_next(slow5_rec_t **read, slow5_file_t *s5p)`

## DESCRIPTION
`slow5_get_next()` fetches a record from a SLOW5 file *s5p* at the current file pointer into a *slow5_rec_t* and stores the address of the *slow5_rec_t* in **read*.

If **read* is set to NULL before the call, then `slow5_get_next()` will allocate a *slow5_rec_t* for storing the record.
This *slow5_rec_t* should be freed by the user program using `slow5_rec_free()`.
Alternatively, before calling `slow5_get_next()`, **read* can contain a pointer to an allocated *slow5_rec_t* from a previous `slow5_get_next()` call.
If the allocated *slow5_rec_t* is not large enough to hold the record, `slow5_get_next()` will resize it internally.

The argument *s5p* points to a *slow5_file_t* opened using `slow5_open()`.

## RETURN VALUE
Upon successful completion, `slow5_get_next()` returns a non negative integer (>=0). Otherwise (including the end of file), a negative value is returned that indicates the error and `slow5_errno` is set to indicate the error.

## ERRORS

* `SLOW5_ERR_EOF`      
        End of file (EOF) (if blow5, EOF marker was properly found at the end).
* `SLOW5_ERR_ARG`      
        Invalid argument - Read or s5p is NULL.
* `SLOW5_ERR_RECPARSE`  
        Record parsing error.
* `SLOW5_ERR_IO`       
        I/O error (except EOF) when reading the slow5 file, for instance, `getline()` or `fread()` failed.
* `SLOW5_ERR_MEM`        
        Memory allocation error.
* `SLOW5_ERR_PRESS`      
       Record decompression error.
* `SLOW5_ERR_TRUNC`      
       Truncated blow5 file - EOF reached without seeing the blow5 EOF marker. 
*  `SLOW5_ERR_UNK`        
       Slow5 format is unknown.

## NOTES
 `slow5_get_next` does not require an index to be pre-loaded, as opposed to `slow5_get()` which so.

## 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(FILE_PATH,"r");
    if(sp==NULL){
       fprintf(stderr,"Error in opening file\n");
       exit(EXIT_FAILURE);
    }
    slow5_rec_t *rec = NULL;
    int ret=0;

    while((ret = slow5_get_next(&rec,sp)) >= 0){
        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");
    }

    if(ret != SLOW5_ERR_EOF){  //check if proper end of file has been reached
        fprintf(stderr,"Error in slow5_get_next. Error code %d\n",ret);
        exit(EXIT_FAILURE);
    }

    slow5_rec_free(rec);

    slow5_close(sp);

}

```

## SEE ALSO

[slow5_get()](slow5_open.md)