File: sload.c

package info (click to toggle)
biosig4c++ 1.9.3-2
  • links: PTS, VCS
  • area: main
  • in suites: buster
  • size: 6,220 kB
  • sloc: ansic: 41,218; cpp: 8,946; sh: 4,365; makefile: 1,758; python: 87; awk: 73; php: 40; perl: 36; java: 14; ruby: 7
file content (63 lines) | stat: -rw-r--r-- 1,821 bytes parent folder | download
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
/*
    Copyright (C) 2016 Alois Schloegl <alois.schloegl@gmail.com>
    This file is part of the "BioSig for C/C++" repository
    (biosig4c++/libbiosig) at http://biosig.sf.net/

    This program is free software; you can redistribute it and/or
    modify it under the terms of the GNU General Public License
    as published by the Free Software Foundation; either version 3
    of the License, or (at your option) any later version.

    This program is distributed in the hope that it will be useful,
    but WITHOUT ANY WARRANTY; without even the implied warranty of
    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
    GNU General Public License for more details.

    You should have received a copy of the GNU General Public License
    along with this program.  If not, see <http://www.gnu.org/licenses/>.
 */

#include <R.h>
#include <Rdefines.h>
#include <R_ext/Rdynload.h>

#include "../biosig.h"

HDRTYPE *hdr=NULL;


SEXP sload(SEXP filename, SEXP channels) {

    SEXP result = R_NilValue;

    // sanity check of input 
    if(!isString(filename) || length(filename) != 1)
        error("filename is not a single string");
        
    const char* fn = CHAR(asChar(filename));

    // Open file 
    HDRTYPE *hdr = sopen(fn, "r", NULL);
    if (serror2(hdr)) return R_NilValue;

    long NS = biosig_get_number_of_channels(hdr);	
  
    // allocate memory for results
    result=PROTECT(allocMatrix(REALSXP, hdr->SPR*hdr->NRec, NS));

    // read data from file and write into result
    int status = sread(REAL(result), 0, hdr->SPR*hdr->NRec, hdr); 
    if (serror2(hdr)) {
        destructHDR(hdr);
	Free(result);
	UNPROTECT(1);
        return R_NilValue;
    }

    // close file, and cleanup memory to avoid any leaks 
    destructHDR(hdr);  

    UNPROTECT(1);
    return result;
}