Re: [PD] tabread4~ broken interpolation algorithm - was Re:, Max Smoother Audio than Pd?

2010-04-03 Thread Matteo Sisti Sette
Charles Henry escribió: The hardest class I ever had was stochastic analysis (as recent as 4 years ago), where we solved problems like this. Fundamentally, it's not too hard, but the details of the calculus are tricky. I'd prefer to stay away unless there's a real good reason to do so :)

Re: [PD] tabread4~ broken interpolation algorithm - was Re:, Max Smoother Audio than Pd?

2010-04-02 Thread Charles Henry
You're trying to restrict the analysis to a convenient (but reasonable) class of signals, and to assume that the signal to be interpolated, x, belongs to that class. Right? Well, sort of. What works well as an interpolator for one signal may not work well for another. The point I started

Re: [PD] tabread4~ broken interpolation algorithm - was Re:, Max Smoother Audio than Pd?

2010-04-01 Thread Charles Henry
On Wed, Mar 31, 2010 at 5:12 PM, Matteo Sisti Sette matteosistise...@gmail.com wrote: It occurs to me that there exists one very obvious function for which the squared error is minimized for a 4-point interpolator.  4-point interpolator impulse functions have to be 0 outside the interval

Re: [PD] tabread4~ broken interpolation algorithm - was Re:, Max Smoother Audio than Pd?

2010-04-01 Thread Matteo Sisti Sette
Charles Henry escribió: The error depends on x the signal. Here, I want to make the *convenient* assumption that the spectrum of x is flat, since we want some kind of generality and we want to minimize average error across frequencies. This would make the problem equivalent to using just

Re: [PD] tabread4~ broken interpolation algorithm - was Re:, Max Smoother Audio than Pd?

2010-04-01 Thread Charles Henry
I get what you're saying too, and I'm at least a little skeptical myself. But as I think about it generally, my entire approach to looking at these problems has been very similar. I basically thought that when comparing interpolators, I could disregard the signals involved and just look at the

Re: [PD] tabread4~ broken interpolation algorithm - was Re:, Max Smoother Audio than Pd?

2010-04-01 Thread Matteo Sisti Sette
Charles Henry escribió: When it comes to the general class of functions with flat spectra, the only difference is in phase, right? But the error is the same in time domain as in frequency domain thanks to the isometric property of the Fourier transform. Our interpolation is the same as a

Re: [PD] tabread4~ broken interpolation algorithm - was Re: Max Smoother Audio than Pd?

2010-03-31 Thread Charles Henry
I don't know either.  We have the formulas for each, so we can calculate squared error vs. sinc(x), but there also appears to be differences in which frequencies the distortion occurs and some could be more audible. It occurs to me that there exists one very obvious function for which the

Re: [PD] tabread4~ broken interpolation algorithm - was Re:, Max Smoother Audio than Pd?

2010-03-31 Thread Matteo Sisti Sette
It occurs to me that there exists one very obvious function for which the squared error is minimized for a 4-point interpolator. 4-point interpolator impulse functions have to be 0 outside the interval [-2,2]. So, E=|f(x)-sinc(x)|^2 is minimized when f(x)={sinc(x) -2x2 ,0 elsewhere

[PD] tabread4~ broken interpolation algorithm - was Re: Max Smoother Audio than Pd?

2010-03-31 Thread Mathieu Bouchard
A workaround is to up-sample everybody by a factor of 2 - this dramatically reduces error, usually by 24 dB. Every upsampling is some kind of interpolation. So what you mean, is to reduce the error of your interpolator by getting much of the job done by a different interpolator instead ?

Re: [PD] tabread4~ broken interpolation algorithm - was Re: Max Smoother Audio than Pd?

2010-03-30 Thread Roman Haefeli
On Mon, 2010-03-29 at 21:06 -0400, Matt Barber wrote: LONG, sorry. Thanks again for your time and patience. One really good way to think, then, is in terms of the continuous frequency response of the interpolator. In that long, long discussion a couple years ago, Chuck Henry made the

Re: [PD] tabread4~ broken interpolation algorithm - was Re:, Max Smoother Audio than Pd?

2010-03-30 Thread Matteo Sisti Sette
Is it really possible to express a cubic interpolator (such as Lagrange or Hermite, i.e. such as tabread4 or tabread4c) in terms of impulse response? Is it equivalent to a convolution? That is to ask: is it linear??? Or is that an approximation? -- Matteo Sisti Sette matteosistise...@gmail.com

Re: [PD] tabread4~ broken interpolation algorithm - was Re:, Max Smoother Audio than Pd?

2010-03-30 Thread Matteo Sisti Sette
So to me it still remains unclear in what aspect [tabread4~] is superior to [tabread4c~], from both a theoretical and from an empirical perspective. The answer may be here: http://www.aes.org/e-lib/browse.cfm?elib=8151 Btw does anybody have access to that article? The analysis reveals an

Re: [PD] tabread4~ broken interpolation algorithm - was Re:, Max Smoother Audio than Pd?

2010-03-30 Thread Matt Barber
Now regarding Matt's words: I have read that the Lagrange interpolators have better stopband attenuation and Hermites have flatter passband response, but I'm not sure this is true Is it possible that it is exactly viceversa? I think it probably is exactly vice-versa -- I was quoting

Re: [PD] tabread4~ broken interpolation algorithm - was Re:, Max Smoother Audio than Pd?

2010-03-30 Thread Matt Barber
Yes, as far as I know it's identical -- when you do one of these interpolations with four points, you can either think of it in terms of a cubic polynomial formula involving those four points, or in terms of the sum of four scaled basis functions - the latter seems to me intuitively equivalent to

Re: [PD] tabread4~ broken interpolation algorithm - was Re:, Max Smoother Audio than Pd?

2010-03-30 Thread Roman Haefeli
On Tue, 2010-03-30 at 14:15 +0200, Matteo Sisti Sette wrote: However, even in presence of a tradeoff that makes some sense (i.e. each of the two choices has advantages and disadvantages), it seems to me that for audio applications the generated high-frequency noise due to discontinuities

Re: [PD] tabread4~ broken interpolation algorithm - was Re:, Max Smoother Audio than Pd?

2010-03-30 Thread Matteo Sisti Sette
Roman Haefeli escribió: On Tue, 2010-03-30 at 14:15 +0200, Matteo Sisti Sette wrote: However, even in presence of a tradeoff that makes some sense (i.e. each of the two choices has advantages and disadvantages), it seems to me that for audio applications the generated high-frequency noise due

Re: [PD] tabread4~ broken interpolation algorithm - was Re:, Max Smoother Audio than Pd?

2010-03-30 Thread cyrille henry
i think this pdf can add lot's of useful information to this thread : http://www.student.oulu.fi/~oniemita/dsp/deip.pdf cyrille Matteo Sisti Sette a écrit : Roman Haefeli escribió: On Tue, 2010-03-30 at 14:15 +0200, Matteo Sisti Sette wrote: However, even in presence of a tradeoff that

Re: [PD] tabread4~ broken interpolation algorithm - was Re:, Max Smoother Audio than Pd?

2010-03-30 Thread Matteo Sisti Sette
Matt Barber escribió: Yes, as far as I know it's identical -- when you do one of these interpolations with four points, you can either think of it in terms of a cubic polynomial formula involving those four points, or in terms of the sum of four scaled basis functions - the latter seems to me

Re: [PD] tabread4~ broken interpolation algorithm - was Re: Max Smoother Audio than Pd?

2010-03-30 Thread Mathieu Bouchard
On Mon, 29 Mar 2010, Matteo Sisti Sette wrote: Mathieu Bouchard escribió: both are truly cubic interpolations. IIRC, one kind of cubic interpolation is designed to go through all four points, and the other kind is designed to be pieced with other cubic interpolations, and Miller confused the

Re: [PD] tabread4~ broken interpolation algorithm - was Re:, Max Smoother Audio than Pd?

2010-03-30 Thread Matt Barber
On Tue, Mar 30, 2010 at 9:52 AM, Matteo Sisti Sette matteosistise...@gmail.com wrote: Matt Barber escribió: Yes, as far as I know it's identical -- when you do one of these interpolations with four points, you can either think of it in terms of a cubic polynomial formula involving those four

[PD] tabread4~ broken interpolation algorithm - was Re: Max Smoother Audio than Pd?

2010-03-29 Thread Matteo Sisti Sette
Hi, Well, if you use [tabread4~] or any of the many other Pd objects that use the same broken interpolation algorithm (copy/paste programming), Broken? What's wrong in the tabread4~ interpolation? (note this is not a sarchastical question implying there's nothing wrong in it, it's a

Re: [PD] tabread4~ broken interpolation algorithm - was Re: Max Smoother Audio than Pd?

2010-03-29 Thread Roman Haefeli
On Mon, 2010-03-29 at 13:49 +0200, Matteo Sisti Sette wrote: Claude wrote: If you use [tabread4] to interpolate graphical parameters for animations, the discontinuities in the derivatives are really obvious. [] But IMHO if you're doing piecewise cubic interpolation, it's a bit

Re: [PD] tabread4~ broken interpolation algorithm - was Re: Max Smoother Audio than Pd?

2010-03-29 Thread Matteo Sisti Sette
Roman Haefeli escribió: Check this thread: http://lists.puredata.info/pipermail/pd-list/2008-06/062878.html I checked it out (not read the _whole_ thread to the end) but, In what way can the current tabread4~ interpolation, which is discontinuous even in its 1st derivative, be superior to

Re: [PD] tabread4~ broken interpolation algorithm - was Re: Max Smoother Audio than Pd?

2010-03-29 Thread Matteo Sisti Sette
By the way tabread4c~ is not in Pd Extended, is it? Roman Haefeli escribió: On Mon, 2010-03-29 at 13:49 +0200, Matteo Sisti Sette wrote: Claude wrote: If you use [tabread4] to interpolate graphical parameters for animations, the discontinuities in the derivatives are really obvious.

Re: [PD] tabread4~ broken interpolation algorithm - was Re: Max Smoother Audio than Pd?

2010-03-29 Thread Miller Puckette
Hi all- I haven't looked at Cyrille's interpolator but... tabread4~ uses true cubic interpolation (which perhaps Cyrille's object also does in some other way). The tabread4~ algorithm is to put a cubic through the 4 points surrounding the input point. However, this cubic curve does not

Re: [PD] tabread4~ broken interpolation algorithm - was Re: Max Smoother Audio than Pd?

2010-03-29 Thread cyrille henry
Matteo Sisti Sette a écrit : By the way tabread4c~ is not in Pd Extended, is it? no. it is there : http://www.chnry.net/ch/?083-Nusmuk-audio c Roman Haefeli escribió: On Mon, 2010-03-29 at 13:49 +0200, Matteo Sisti Sette wrote: Claude wrote: If you use [tabread4] to interpolate

Re: [PD] tabread4~ broken interpolation algorithm - was Re: Max Smoother Audio than Pd?

2010-03-29 Thread Matteo Sisti Sette
cyrille henry escribió: Matteo Sisti Sette a écrit : By the way tabread4c~ is not in Pd Extended, is it? no. it is there : http://www.chnry.net/ch/?083-Nusmuk-audio Hi, I downloaded the zip file but Windows tells me he can't open it. Is it something different than a normal .zip file? (it

Re: [PD] tabread4~ broken interpolation algorithm - was Re: Max Smoother Audio than Pd?

2010-03-29 Thread Matteo Sisti Sette
Miller Puckette escribió: The tabread4~ algorithm is to put a cubic through the 4 points surrounding the input point. However, this cubic curve does not necessarily match the next curve over in first derivative. Oh I see! I thought it did. I confuded that technique with natural cubic

Re: [PD] tabread4~ broken interpolation algorithm - was Re: Max Smoother Audio than Pd?

2010-03-29 Thread John Harrison
The link works and extracts fine here in WinXP SP3. I used the built-in compressed (zipped) folders tool in the explorer shell. On Mon, Mar 29, 2010 at 11:38 AM, Matteo Sisti Sette matteosistise...@gmail.com wrote: cyrille henry escribió: Matteo Sisti Sette a écrit : By the way

[PD] tabread4~ broken interpolation algorithm - was Re: Max Smoother Audio than Pd?

2010-03-29 Thread Mathieu Bouchard
I checked it out (not read the _whole_ thread to the end) but, In what way can the current tabread4~ interpolation, which is discontinuous even in its 1st derivative, be superior to true cubic interpolation? Even at transpositions near to zero, I can't see what's the advantage, nor what it

Re: [PD] tabread4~ broken interpolation algorithm - was Re: Max Smoother Audio than Pd?

2010-03-29 Thread Matteo Sisti Sette
John Harrison escribió: The link works and extracts fine here in WinXP SP3. I used the built-in compressed (zipped) folders tool in the explorer shell. Crazy. Internet Explorer breaks downloaded files whenever it takes you more than a few seconds to select the folder to download in. I

Re: [PD] tabread4~ broken interpolation algorithm - was Re: Max Smoother Audio than Pd?

2010-03-29 Thread Matteo Sisti Sette
Mathieu Bouchard escribió: both are truly cubic interpolations. IIRC, one kind of cubic interpolation is designed to go through all four points, and the other kind is designed to be pieced with other cubic interpolations, and Miller confused the two and left the bug there. According to his

Re: [PD] tabread4~ broken interpolation algorithm - was Re: Max Smoother Audio than Pd?

2010-03-29 Thread Matt Barber
I checked it out (not read the _whole_ thread to the end) but, In what way can the current tabread4~ interpolation, which is discontinuous even in its 1st derivative, be superior to true cubic interpolation? Even at transpositions near to zero, I can't see what's the advantage, nor what it

Re: [PD] tabread4~ broken interpolation algorithm - was Re: Max Smoother Audio than Pd?

2010-03-29 Thread Roman Haefeli
Hi Matt Thanks for the detailed explanation. I still have troubles getting the idea of the Lagrange interpolator in the context of [tabread4~]. You say, that it finds the cubic polynomial which hits all four points. But what is the advantage of that? As I understand [tabread4~], if the index is

Re: [PD] tabread4~ broken interpolation algorithm - was Re: Max Smoother Audio than Pd?

2010-03-29 Thread Matt Barber
LONG, sorry. On Mon, Mar 29, 2010 at 7:03 PM, Roman Haefeli reduzie...@yahoo.de wrote: Hi Matt Thanks for the detailed explanation.  I still have troubles getting the idea of the Lagrange interpolator in the context of [tabread4~]. You say, that it finds the cubic polynomial which hits all