28 #if defined(TEMPLATE_RESAMPLE_DBL)
30 # define RENAME(N) N ## _double
31 # define FILTER_SHIFT 0
34 # define FELEM2 double
35 # define OUT(d, v) d = v
37 #elif defined(TEMPLATE_RESAMPLE_FLT)
39 # define RENAME(N) N ## _float
40 # define FILTER_SHIFT 0
44 # define OUT(d, v) d = v
46 #elif defined(TEMPLATE_RESAMPLE_S32)
48 # define RENAME(N) N ## _int32
49 # define FILTER_SHIFT 30
50 # define DELEM int32_t
51 # define FELEM int32_t
52 # define FELEM2 int64_t
53 # define FELEM_MAX INT32_MAX
54 # define FELEM_MIN INT32_MIN
55 # define OUT(d, v) (v) = ((v) + (1<<(FILTER_SHIFT-1)))>>FILTER_SHIFT;\
56 (d) = (uint64_t)((v) + 0x80000000) > 0xFFFFFFFF ? ((v)>>63) ^ 0x7FFFFFFF : (v)
58 #elif defined(TEMPLATE_RESAMPLE_S16)
60 # define RENAME(N) N ## _int16
61 # define FILTER_SHIFT 15
62 # define DELEM int16_t
63 # define FELEM int16_t
64 # define FELEM2 int32_t
65 # define FELEML int64_t
66 # define FELEM_MAX INT16_MAX
67 # define FELEM_MIN INT16_MIN
68 # define OUT(d, v) (v) = ((v) + (1<<(FILTER_SHIFT-1)))>>FILTER_SHIFT;\
69 (d) = (unsigned)((v) + 32768) > 65535 ? ((v)>>31) ^ 32767 : (v)
74 int dst_size, int64_t index2, int64_t incr)
77 const DELEM *
src = source;
80 for (dst_index = 0; dst_index < dst_size; dst_index++) {
81 dst[dst_index] = src[index2 >> 32];
87 void *dest,
const void *source,
88 int n,
int update_ctx)
91 const DELEM *
src = source;
95 int sample_index = index >>
c->phase_shift;
97 index &=
c->phase_mask;
98 for (dst_index = 0; dst_index <
n; dst_index++) {
103 for (i = 0; i <
c->filter_length; i++) {
104 val += src[sample_index + i] * (
FELEM2)filter[i];
106 OUT(dst[dst_index], val);
108 frac +=
c->dst_incr_mod;
109 index +=
c->dst_incr_div;
110 if (frac >=
c->src_incr) {
114 sample_index += index >>
c->phase_shift;
115 index &=
c->phase_mask;
127 void *dest,
const void *source,
128 int n,
int update_ctx)
131 const DELEM *
src = source;
135 int sample_index = index >>
c->phase_shift;
136 #if FILTER_SHIFT == 0
137 double inv_src_incr = 1.0 /
c->src_incr;
140 index &=
c->phase_mask;
141 for (dst_index = 0; dst_index <
n; dst_index++) {
146 for (i = 0; i <
c->filter_length; i++) {
147 val += src[sample_index + i] * (
FELEM2)filter[i];
148 v2 += src[sample_index + i] * (
FELEM2)filter[i +
c->filter_alloc];
151 val += (v2 -
val) * (
FELEML) frac /
c->src_incr;
153 # if FILTER_SHIFT == 0
154 val += (v2 -
val) * inv_src_incr * frac;
156 val += (v2 -
val) /
c->src_incr * frac;
159 OUT(dst[dst_index], val);
161 frac +=
c->dst_incr_mod;
162 index +=
c->dst_incr_div;
163 if (frac >=
c->src_incr) {
167 sample_index += index >>
c->phase_shift;
168 index &=
c->phase_mask;