Page Menu
Home
c4science
Search
Configure Global Search
Log In
Files
F88553717
filter_MEX.c
No One
Temporary
Actions
Download File
Edit File
Delete File
View Transforms
Subscribe
Mute Notifications
Award Token
Subscribers
None
File Metadata
Details
File Info
Storage
Attached
Created
Sat, Oct 19, 11:03
Size
1 KB
Mime Type
text/x-c
Expires
Mon, Oct 21, 11:03 (1 d, 21 h)
Engine
blob
Format
Raw Data
Handle
21789291
Attached To
R10025 iCAPs public version 2
filter_MEX.c
View Options
/*
Author : E. Orliac, SCITAS, EPFL
Date : 27.10.2017
Purpose:
Remarks:
*/
/*
num : numerator
lnum : length of numerator
den1 : causal
den2 : non-causal
ld1 : length causal
ld2 : length non-causal
condition: 0 = 'normal'
1 = 'transpose'
*/
#include <stdio.h>
#include <stdlib.h>
#include "filter_MEX.h"
void
filter_MEX
(
double
*
out
,
double
*
in
,
int
lin
,
double
*
num
,
int
lnum
,
double
*
den
,
int
lden
)
{
int
i
,
j
;
double
acc
;
for
(
i
=
0
;
i
<
lin
;
i
++
)
out
[
i
]
=
0.0
;
if
(
lden
==
1
&&
den
[
0
]
==
1.
d
)
{
for
(
i
=
0
;
i
<
lin
;
i
++
)
{
acc
=
0.0
d
;
for
(
j
=
0
;
j
<
lnum
;
j
++
)
{
if
(
j
<=
i
)
{
acc
+=
num
[
j
]
*
in
[
i
-
j
];
}
}
out
[
i
]
=
acc
;
}
}
else
if
(
lnum
==
1
&&
num
[
0
]
==
1.
d
)
{
for
(
i
=
0
;
i
<
lin
;
i
++
)
{
acc
=
in
[
i
];
for
(
j
=
0
;
j
<
lden
;
j
++
)
{
if
(
j
<=
i
)
{
acc
-=
den
[
j
]
*
out
[
i
-
j
];
}
}
out
[
i
]
=
acc
;
}
}
else
{
fprintf
(
stderr
,
"Fatal. Unknown case in filter_MEX.
\n
"
);
exit
(
EXIT_FAILURE
);
}
}
Event Timeline
Log In to Comment