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
|
function checkfiltres(filein, fileout)
fout = fopen(fileout);
pdattype = fread(fout, 1, 'int32');
numlen = fread(fout, 1, 'uint32');
num = fread_array(fout, [1,numlen], pdattype);
denumlen = fread(fout, 1, 'uint32');
denum = fread_array(fout, [1,denumlen], pdattype);
datouttype = fread(fout, 1, 'int32');
nchannout = fread(fout, 1, 'uint32');
ntotout = fread(fout, 1, 'uint32');
datout = fread_array(fout, [nchannout, ntotout], datouttype);
fclose(fout);
fin = fopen(filein);
datintype = fread(fin, 1, 'int32');
nchannin = fread(fin, 1, 'uint32');
ntotin = fread(fin, 1, 'uint32');
datin = fread_array(fin, [nchannin, ntotin], datintype);
fclose(fin);
if (nchannin ~= nchannout)
error('data params differs in the 2 files');
end
niter = 100;
tic
for k=1:niter
datmatlab = filter(num,denum,datin,[],2);
end
dt = toc;
samdt = dt/(niter*nchannin*size(datin,2));
fprintf('time per sample: %f ns\n',samdt*1e9);
% for ichann=1:6:nchann
% figure
% plot(datout(ichann,:),'b')
% hold on
% plot(datmatlab(ichann,:),'r')
% % plot(datmatlab(ichann,:),'g--')
% hold off
% axis([0 600 -2 2])
% end
diffval = abs(datmatlab - datout);
errvals = max(diffval,[],2) ./ max(abs(datmatlab),[],2);
errval = max(errvals);
fprintf('Error value = %10.10g\n',errval);
return;
function A = fread_array(fid, arrsize, typeval)
type = '*float32';
if bitand(typeval,1)
type = '*float64';
end
c = 0;
if bitand(typeval,2)
c = 1;
arrsize(1) = arrsize(1)*2;
end
A = fread(fid, arrsize, type);
if c == 1
A = A(1:2:end,:) + 1i*A(2:2:end,:);
end
return
|
