File: random_read_openmp.c

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 (50 lines) | stat: -rw-r--r-- 1,403 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

// an example programme that uses slow5lib to randomly access records in a SLOW5 file using multiple threads (openMP)
#include <omp.h>
#include <stdio.h>
#include <stdlib.h>
#include <slow5/slow5.h>

#define FILE_PATH "examples/example.slow5"
#define READ_LIST_SIZE 4

char * read_id_list[READ_LIST_SIZE] = {"r4", "r1", "r3", "r4"};

int main(){

    //open the SLOW5 file
    slow5_file_t *sp = slow5_open(FILE_PATH,"r");
    if(sp==NULL){
       fprintf(stderr,"Error in opening file\n");
       exit(EXIT_FAILURE);
    }

    int ret=0; //for return value

    //load the SLOW5 index (will be built if not present)
    ret = slow5_idx_load(sp);
    if (ret < 0) {
        fprintf(stderr, "Error in loading index\n");
        exit(EXIT_FAILURE);
    }

    #pragma omp parallel for
    for(int32_t i=0; i<READ_LIST_SIZE; i++) {
        slow5_rec_t *rec = NULL; //slow5 record to be read
        ret = slow5_get(read_id_list[i], &rec, sp); //fetch the read
        if (ret < 0) {
            fprintf(stderr, "Error in when fetching the read %s\n",read_id_list[i]);
        } else {
            fprintf(stderr, "Successfully fetched the read %s with %ld raw signal samples\n", rec->read_id, rec->len_raw_signal);
        }
        slow5_rec_free(rec); //free the SLOW5 record
    }

    //free the SLOW5 index
    slow5_idx_unload(sp);

    //close the SLOW5 file
    slow5_close(sp);

    return 0;
}