1/*
2 *  Rate conversion Plug-In
3 *  Copyright (c) 1999 by Jaroslav Kysela <perex@perex.cz>
4 *
5 *
6 *   This library is free software; you can redistribute it and/or modify
7 *   it under the terms of the GNU Library General Public License as
8 *   published by the Free Software Foundation; either version 2 of
9 *   the License, or (at your option) any later version.
10 *
11 *   This program is distributed in the hope that it will be useful,
12 *   but WITHOUT ANY WARRANTY; without even the implied warranty of
13 *   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
14 *   GNU Library General Public License for more details.
15 *
16 *   You should have received a copy of the GNU Library General Public
17 *   License along with this library; if not, write to the Free Software
18 *   Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
19 *
20 */
21
22#include <linux/time.h>
23#include <sound/core.h>
24#include <sound/pcm.h>
25#include "pcm_plugin.h"
26
27#define SHIFT	11
28#define BITS	(1<<SHIFT)
29#define R_MASK	(BITS-1)
30
31/*
32 *  Basic rate conversion plugin
33 */
34
35struct rate_channel {
36	signed short last_S1;
37	signed short last_S2;
38};
39
40typedef void (*rate_f)(struct snd_pcm_plugin *plugin,
41		       const struct snd_pcm_plugin_channel *src_channels,
42		       struct snd_pcm_plugin_channel *dst_channels,
43		       int src_frames, int dst_frames);
44
45struct rate_priv {
46	unsigned int pitch;
47	unsigned int pos;
48	rate_f func;
49	snd_pcm_sframes_t old_src_frames, old_dst_frames;
50	struct rate_channel channels[];
51};
52
53static void rate_init(struct snd_pcm_plugin *plugin)
54{
55	unsigned int channel;
56	struct rate_priv *data = (struct rate_priv *)plugin->extra_data;
57	data->pos = 0;
58	for (channel = 0; channel < plugin->src_format.channels; channel++) {
59		data->channels[channel].last_S1 = 0;
60		data->channels[channel].last_S2 = 0;
61	}
62}
63
64static void resample_expand(struct snd_pcm_plugin *plugin,
65			    const struct snd_pcm_plugin_channel *src_channels,
66			    struct snd_pcm_plugin_channel *dst_channels,
67			    int src_frames, int dst_frames)
68{
69	unsigned int pos = 0;
70	signed int val;
71	signed short S1, S2;
72	signed short *src, *dst;
73	unsigned int channel;
74	int src_step, dst_step;
75	int src_frames1, dst_frames1;
76	struct rate_priv *data = (struct rate_priv *)plugin->extra_data;
77	struct rate_channel *rchannels = data->channels;
78
79	for (channel = 0; channel < plugin->src_format.channels; channel++) {
80		pos = data->pos;
81		S1 = rchannels->last_S1;
82		S2 = rchannels->last_S2;
83		if (!src_channels[channel].enabled) {
84			if (dst_channels[channel].wanted)
85				snd_pcm_area_silence(&dst_channels[channel].area, 0, dst_frames, plugin->dst_format.format);
86			dst_channels[channel].enabled = 0;
87			continue;
88		}
89		dst_channels[channel].enabled = 1;
90		src = (signed short *)src_channels[channel].area.addr +
91			src_channels[channel].area.first / 8 / 2;
92		dst = (signed short *)dst_channels[channel].area.addr +
93			dst_channels[channel].area.first / 8 / 2;
94		src_step = src_channels[channel].area.step / 8 / 2;
95		dst_step = dst_channels[channel].area.step / 8 / 2;
96		src_frames1 = src_frames;
97		dst_frames1 = dst_frames;
98		while (dst_frames1-- > 0) {
99			if (pos & ~R_MASK) {
100				pos &= R_MASK;
101				S1 = S2;
102				if (src_frames1-- > 0) {
103					S2 = *src;
104					src += src_step;
105				}
106			}
107			val = S1 + ((S2 - S1) * (signed int)pos) / BITS;
108			if (val < -32768)
109				val = -32768;
110			else if (val > 32767)
111				val = 32767;
112			*dst = val;
113			dst += dst_step;
114			pos += data->pitch;
115		}
116		rchannels->last_S1 = S1;
117		rchannels->last_S2 = S2;
118		rchannels++;
119	}
120	data->pos = pos;
121}
122
123static void resample_shrink(struct snd_pcm_plugin *plugin,
124			    const struct snd_pcm_plugin_channel *src_channels,
125			    struct snd_pcm_plugin_channel *dst_channels,
126			    int src_frames, int dst_frames)
127{
128	unsigned int pos = 0;
129	signed int val;
130	signed short S1, S2;
131	signed short *src, *dst;
132	unsigned int channel;
133	int src_step, dst_step;
134	int src_frames1, dst_frames1;
135	struct rate_priv *data = (struct rate_priv *)plugin->extra_data;
136	struct rate_channel *rchannels = data->channels;
137
138	for (channel = 0; channel < plugin->src_format.channels; ++channel) {
139		pos = data->pos;
140		S1 = rchannels->last_S1;
141		S2 = rchannels->last_S2;
142		if (!src_channels[channel].enabled) {
143			if (dst_channels[channel].wanted)
144				snd_pcm_area_silence(&dst_channels[channel].area, 0, dst_frames, plugin->dst_format.format);
145			dst_channels[channel].enabled = 0;
146			continue;
147		}
148		dst_channels[channel].enabled = 1;
149		src = (signed short *)src_channels[channel].area.addr +
150			src_channels[channel].area.first / 8 / 2;
151		dst = (signed short *)dst_channels[channel].area.addr +
152			dst_channels[channel].area.first / 8 / 2;
153		src_step = src_channels[channel].area.step / 8 / 2;
154		dst_step = dst_channels[channel].area.step / 8 / 2;
155		src_frames1 = src_frames;
156		dst_frames1 = dst_frames;
157		while (dst_frames1 > 0) {
158			S1 = S2;
159			if (src_frames1-- > 0) {
160				S2 = *src;
161				src += src_step;
162			}
163			if (pos & ~R_MASK) {
164				pos &= R_MASK;
165				val = S1 + ((S2 - S1) * (signed int)pos) / BITS;
166				if (val < -32768)
167					val = -32768;
168				else if (val > 32767)
169					val = 32767;
170				*dst = val;
171				dst += dst_step;
172				dst_frames1--;
173			}
174			pos += data->pitch;
175		}
176		rchannels->last_S1 = S1;
177		rchannels->last_S2 = S2;
178		rchannels++;
179	}
180	data->pos = pos;
181}
182
183static snd_pcm_sframes_t rate_src_frames(struct snd_pcm_plugin *plugin, snd_pcm_uframes_t frames)
184{
185	struct rate_priv *data;
186	snd_pcm_sframes_t res;
187
188	if (snd_BUG_ON(!plugin))
189		return -ENXIO;
190	if (frames == 0)
191		return 0;
192	data = (struct rate_priv *)plugin->extra_data;
193	if (plugin->src_format.rate < plugin->dst_format.rate) {
194		res = (((frames * data->pitch) + (BITS/2)) >> SHIFT);
195	} else {
196		res = DIV_ROUND_CLOSEST(frames << SHIFT, data->pitch);
197	}
198	if (data->old_src_frames > 0) {
199		snd_pcm_sframes_t frames1 = frames, res1 = data->old_dst_frames;
200		while (data->old_src_frames < frames1) {
201			frames1 >>= 1;
202			res1 <<= 1;
203		}
204		while (data->old_src_frames > frames1) {
205			frames1 <<= 1;
206			res1 >>= 1;
207		}
208		if (data->old_src_frames == frames1)
209			return res1;
210	}
211	data->old_src_frames = frames;
212	data->old_dst_frames = res;
213	return res;
214}
215
216static snd_pcm_sframes_t rate_dst_frames(struct snd_pcm_plugin *plugin, snd_pcm_uframes_t frames)
217{
218	struct rate_priv *data;
219	snd_pcm_sframes_t res;
220
221	if (snd_BUG_ON(!plugin))
222		return -ENXIO;
223	if (frames == 0)
224		return 0;
225	data = (struct rate_priv *)plugin->extra_data;
226	if (plugin->src_format.rate < plugin->dst_format.rate) {
227		res = DIV_ROUND_CLOSEST(frames << SHIFT, data->pitch);
228	} else {
229		res = (((frames * data->pitch) + (BITS/2)) >> SHIFT);
230	}
231	if (data->old_dst_frames > 0) {
232		snd_pcm_sframes_t frames1 = frames, res1 = data->old_src_frames;
233		while (data->old_dst_frames < frames1) {
234			frames1 >>= 1;
235			res1 <<= 1;
236		}
237		while (data->old_dst_frames > frames1) {
238			frames1 <<= 1;
239			res1 >>= 1;
240		}
241		if (data->old_dst_frames == frames1)
242			return res1;
243	}
244	data->old_dst_frames = frames;
245	data->old_src_frames = res;
246	return res;
247}
248
249static snd_pcm_sframes_t rate_transfer(struct snd_pcm_plugin *plugin,
250			     const struct snd_pcm_plugin_channel *src_channels,
251			     struct snd_pcm_plugin_channel *dst_channels,
252			     snd_pcm_uframes_t frames)
253{
254	snd_pcm_uframes_t dst_frames;
255	struct rate_priv *data;
256
257	if (snd_BUG_ON(!plugin || !src_channels || !dst_channels))
258		return -ENXIO;
259	if (frames == 0)
260		return 0;
261#ifdef CONFIG_SND_DEBUG
262	{
263		unsigned int channel;
264		for (channel = 0; channel < plugin->src_format.channels; channel++) {
265			if (snd_BUG_ON(src_channels[channel].area.first % 8 ||
266				       src_channels[channel].area.step % 8))
267				return -ENXIO;
268			if (snd_BUG_ON(dst_channels[channel].area.first % 8 ||
269				       dst_channels[channel].area.step % 8))
270				return -ENXIO;
271		}
272	}
273#endif
274
275	dst_frames = rate_dst_frames(plugin, frames);
276	if (dst_frames > dst_channels[0].frames)
277		dst_frames = dst_channels[0].frames;
278	data = (struct rate_priv *)plugin->extra_data;
279	data->func(plugin, src_channels, dst_channels, frames, dst_frames);
280	return dst_frames;
281}
282
283static int rate_action(struct snd_pcm_plugin *plugin,
284		       enum snd_pcm_plugin_action action,
285		       unsigned long udata)
286{
287	if (snd_BUG_ON(!plugin))
288		return -ENXIO;
289	switch (action) {
290	case INIT:
291	case PREPARE:
292		rate_init(plugin);
293		break;
294	default:
295		break;
296	}
297	return 0;	/* silenty ignore other actions */
298}
299
300int snd_pcm_plugin_build_rate(struct snd_pcm_substream *plug,
301			      struct snd_pcm_plugin_format *src_format,
302			      struct snd_pcm_plugin_format *dst_format,
303			      struct snd_pcm_plugin **r_plugin)
304{
305	int err;
306	struct rate_priv *data;
307	struct snd_pcm_plugin *plugin;
308
309	if (snd_BUG_ON(!r_plugin))
310		return -ENXIO;
311	*r_plugin = NULL;
312
313	if (snd_BUG_ON(src_format->channels != dst_format->channels))
314		return -ENXIO;
315	if (snd_BUG_ON(src_format->channels <= 0))
316		return -ENXIO;
317	if (snd_BUG_ON(src_format->format != SNDRV_PCM_FORMAT_S16))
318		return -ENXIO;
319	if (snd_BUG_ON(dst_format->format != SNDRV_PCM_FORMAT_S16))
320		return -ENXIO;
321	if (snd_BUG_ON(src_format->rate == dst_format->rate))
322		return -ENXIO;
323
324	err = snd_pcm_plugin_build(plug, "rate conversion",
325				   src_format, dst_format,
326				   struct_size(data, channels,
327					       src_format->channels),
328				   &plugin);
329	if (err < 0)
330		return err;
331	data = (struct rate_priv *)plugin->extra_data;
332	if (src_format->rate < dst_format->rate) {
333		data->pitch = ((src_format->rate << SHIFT) + (dst_format->rate >> 1)) / dst_format->rate;
334		data->func = resample_expand;
335	} else {
336		data->pitch = ((dst_format->rate << SHIFT) + (src_format->rate >> 1)) / src_format->rate;
337		data->func = resample_shrink;
338	}
339	data->pos = 0;
340	rate_init(plugin);
341	data->old_src_frames = data->old_dst_frames = 0;
342	plugin->transfer = rate_transfer;
343	plugin->src_frames = rate_src_frames;
344	plugin->dst_frames = rate_dst_frames;
345	plugin->action = rate_action;
346	*r_plugin = plugin;
347	return 0;
348}
349