Re: [PD] cross-correlation in pd ?

```
Alright here is my last version of the patches.  Even if no one
else is interested, I consider the problem solved at this point.
Whenever the calculations are reset, there is a brief transient.
digression:
That's because I can't find a way to precisely calculate a signal
which represents the sum of squares of s1, as follows:
output(k)=sum(i= -32+k:31+k; s1(i)^2) for k=-32:31
in the same block as sxcov~
(I tried to do that first, but no easy solution came up, so I chose
something simpler, if only a little less accurate)

anyway, each of the patches from previous posts have been modified and
corrected.  Hopefully, this method gives you the most flexibility for
computing sxcov~ and sxcorr~ in real time with abstractions.

Chuck

On 6/20/07, Charles Henry [EMAIL PROTECTED] wrote:

This one works within block sizes of 64, by using an [block~ 128 2] .

Now, I've got one where you add an argument and use arbitrary block
sizes.  It's much more useful, this way.
[sxcov~ 2048] works within a blocksize of 2048, and calculates the
symmetric cross covariance of two signals.

Still having aggravating time, trying to normalize by the proper
calculation of variance on each signal.  It might be a while before I
get it together.

Chuck

sxcov~.pd
Description: Binary data

sxcov~-test.pd
Description: Binary data

sxcorr~.pd
Description: Binary data

sxcorr~-test.pd
Description: Binary data

conj_mult~.pd
Description: Binary data
___
PD-list@iem.at mailing list
UNSUBSCRIBE and account-management -
http://lists.puredata.info/listinfo/pd-list

```

Re: [PD] cross-correlation in pd ?

```
This one works within block sizes of 64, by using an [block~ 128 2] .

Now, I've got one where you add an argument and use arbitrary block
sizes.  It's much more useful, this way.
[sxcov~ 2048] works within a blocksize of 2048, and calculates the
symmetric cross covariance of two signals.

Still having aggravating time, trying to normalize by the proper
calculation of variance on each signal.  It might be a while before I
get it together.

Chuck

sxcov~-test.pd
Description: Binary data

conj_mult~.pd
Description: Binary data

sxcov~.pd
Description: Binary data
___
PD-list@iem.at mailing list
UNSUBSCRIBE and account-management -
http://lists.puredata.info/listinfo/pd-list

```

Re: [PD] cross-correlation in pd ?

```
Hi,

On Sun, 2007-06-17 at 23:01 -0500, Charles Henry wrote:

Jamie Bullock has cc~ in his flib collection.  It computes cross
correlation in the time domain, and cross-covariance using the freq.
domain.  I worked on it a while back...but have been lazy about trying
to change the cross-covariance into a cross-correlation.  Sorry,
Jamie.

flib is now deprecated, and has been surpassed by libxtract, which is a
generic C library that includes as an example a PD external that wraps
it's functionality.

It can be found at http://sourceforge.net/projects/libxtract. If you
want to use the fft-based autocorrelation you will need the fftw
library, and to compile libxtract with --enable-fft.

It includes a time domain and fft-based autocorrelation, but no
arbitrary cross correlation functions.

Chuck, I must remember to add your cross-covariance code into libxtract
(with acknowledgement of course).

best,

Jamie

___
PD-list@iem.at mailing list
UNSUBSCRIBE and account-management -
http://lists.puredata.info/listinfo/pd-list

```

Re: [PD] cross-correlation in pd ?

```Hallo!

The problem is that two of the objects I have in the patch that Andy sent me
1- tab_mul
2- tab_cross_corr

these should be in the iem_tab external (included in latest pd extended)

LG
Georg

___
PD-list@iem.at mailing list
UNSUBSCRIBE and account-management -
http://lists.puredata.info/listinfo/pd-list

```

Re: [PD] cross-correlation in pd ?

```Hi Tania,

Download the distribution from Paul#347; website: http://aubio.piem.org/pub/
tar xvzf aubio-0.3.2.tar.gz
There is no makefile because
./configure
should create one, then you type
make  make install.

do a
ls /usr/local/lib/pd/extra/aubio.pd_darwin
if you´ve installed PD from source. If you downloaded a precompiled binary then
you need to copy the files manually to the Pd-0.40-2 app folder (it is the same
as the application).
cd plugins/puredata/
cp aubio.pd_darwin /path/to/Pd-0.40-2.app/Contents/Resources/extra
cp help/* /path/to/Pd-0.40-2.app/Contents/Resources/doc/5.reference

if you are using a .pdrc file in your home directory, you need to add
-lib aubio
to it. If you are using Hans´ Pd-extended distribution, you will need to alter
the org.puredata.pd.plist file to reflect the new lib.

Best,
Ed

tania habib [EMAIL PROTECTED] wrote: Hello Ed,
Thanks for writing down all the details, everything went smoothly with a
problem that the aubio package does no have the makefile in it.
It returns with an error in terminal.you can check the attached file
Can you please send the package aubio, if you have it with you?

ps. I am not using the lpc coding and formant analysis yet but I will keep you
posted if there is some progress in that area.

regards,
tania

On 6/17/07, Ed Kelly [EMAIL PROTECTED]   wrote: Hi Tania,

aubio is quite easy, provided you have libsamplerate, libsndfile and fftw3
installed.  It is pretty essential that you instal pkgconfig as well, since
some libraries rely on that to send information as to where they all reside on
the hard disk for compilation.
So, aubio in steps:

First, URLs:
pkgconfig: http://pkg-config.freedesktop.org/wiki/
libsndfile: www.mega-nerd.com/libsndfile/
libsamplerate:  www.mega-nerd.com/SRC/
fftw:  http://www.fftw.org/
and aubio-0.3.2.tar.gz
unpack them all in the same directory. Put the attached script in the directory
all the other packages are in then...

Then you can use this script to compile them all (see attachment). You will
need to be logged in as the administrator, and know the administrator password.
You must open Applications - Utilities - Terminal

# then type

cd  /path/to/where-you-unpacked-the-files/
sudo chmod 744 instaubio.sh
password
./instaubio.sh

have a look at the script using textedit - it uses fairly common unix stuff,
and I'm sure you could work it out yourself, or you know it. Maybe I saved you
some time! ;-D but aubio needs things in the right order.

PS I know Paul is writing up   his thesis right now, so maybe he's a little
quiet!
PPS do you know much about linear predictive coding and formant analysis?

Best,
Ed

tania habib [EMAIL PROTECTED] wrote:  The problem is that two of the objects
I have in the patch that Andy sent me
1- tab_mul
2- tab_cross_corr

tare not created, I have the IEMLIB included in pd and the rest is ok.
Please tell me whether they are coming from same library or they are part of
some  other one which needs to be included with it

Can you please check them on yours pd?

On 6/17/07, Hans-Christoph Steiner  [EMAIL PROTECTED] wrote:
On Jun 17, 2007, at 12:46 PM, tania habib wrote:

1- how can i find this IEMLIB_R1.15 library, as some of the objects in the
patch you sent me can not be created.

2- I am using pd extended   0.40 on mac osx

In Pd-extended, that should all be in the iemlib library.

3- how can I include aubio in pd when I download it, a little guidance needed
here

That I don't know, I haven't played with it yet.

.hc

thanks
tania

On 6/18/07, Andy Farnell [EMAIL PROTECTED] wrote:

I found a helpfile for table cross corr IEMLIB_R1.15
Don't actually have the external though myself for some
reason.

Looking at the helpfile it would be suitable for audio
after using [tabsend~] to fill the tables.

On Sun, 17 Jun 2007 17:52:42 +0200
tania habib   [EMAIL PROTECTED] wrote:

hello list,
I was wondering if some body can tell me if there is   an object to generate
the cross correlation of  two signals.
If there is, please let me know which library has it.

regards,
tania

--
Use the source

___
PD-list@iem.at mailing list
UNSUBSCRIBE and account-management -
http://lists.puredata.info/listinfo/pd-list

There is no way to peace, peace is the way.   -  A.J. Muste

___
PD-list@iem.at mailing list
UNSUBSCRIBE and   account-management -
http://lists.puredata.info/listinfo/pd-list

Lone Shark Aviation out now on  http://www.pyramidtransmissions.com
http://www.myspace.com/sharktracksInbox cluttering up with junk? Clean
up with  Yahoo! Mail.

Lone Shark Aviation out now on   http://www.pyramidtransmissions.com
http://www.myspace.com/sharktracks

-
Web email has come of age. Don't settle for ```

Re: [PD] cross-correlation in pd ?

```
Here's an abstraction for computing a symmetric cross-covariance, this
way.  It outputs the cross-covariance, xcov(k)=sum(i=-32,...,31;
s1(i+k)*s2(i))
It's still pretty ugly, and the details of the math confuse me a bit.
I'm still working on the one-sided cross covariance function for
delays, like from [0,63] instead of [-32,31]

and here's the same thing except it's a forward cross-covariance,
called fxcov~.pd
It's pretty much the same thing (up to a certain point), but as I've
discovered, using xcov to calculate PSD and CSD for deconvolution can
be tricky (even when it's a full rank system).  There's a big diff
between symmetric and forward cross-covariance in this respect.  None
of my reverberation measurement patches ever worked completely due to
audio dropouts and badly formatted problems (I'm still working on it).
The fxcov~ patch contains some math for your comparison.

Chuck

fxcov~.pd
Description: Binary data

fxcov~-test.pd
Description: Binary data

conj_mult~.pd
Description: Binary data
___
PD-list@iem.at mailing list
UNSUBSCRIBE and account-management -
http://lists.puredata.info/listinfo/pd-list

```

[PD] cross-correlation in pd ?

```
hello list,
I was wondering if some body can tell me if there is an object to generate
the cross correlation of  two signals.
If there is, please let me know which library has it.

regards,
tania
___
PD-list@iem.at mailing list
UNSUBSCRIBE and account-management -
http://lists.puredata.info/listinfo/pd-list

```

Re: [PD] cross-correlation in pd ?

```

tania habib a écrit :
hello list,
I was wondering if some body can tell me if there is an object to generate
the cross correlation of  two signals.
mapping/correlation
it's data correlation (not audio)

cyrille

If there is, please let me know which library has it.

regards,
tania

___
PD-list@iem.at mailing list
UNSUBSCRIBE and account-management -
http://lists.puredata.info/listinfo/pd-list

___
PD-list@iem.at mailing list
UNSUBSCRIBE and account-management -
http://lists.puredata.info/listinfo/pd-list

```

Re: [PD] cross-correlation in pd ?

```

Hey,

I am not quite sure what exactly what you are looking for, but there
are a couple libraries that are meant for getting data from sound,
like aubio:

http://aubio.piem.org/

.hc

On Jun 17, 2007, at 11:52 AM, tania habib wrote:

hello list,
I was wondering if some body can tell me if there is an object to
generate the cross correlation of  two signals.

If there is, please let me know which library has it.

regards,
tania
___
PD-list@iem.at mailing list
UNSUBSCRIBE and account-management - http://lists.puredata.info/
listinfo/pd-list

I have the audacity to believe that peoples everywhere can have three
meals a day for their bodies, education and culture for their minds,
and dignity, equality and freedom for their spirits.  - Martin
Luther King, Jr.

___
PD-list@iem.at mailing list
UNSUBSCRIBE and account-management -
http://lists.puredata.info/listinfo/pd-list

```

Re: [PD] cross-correlation in pd ?

```

I found a helpfile for table cross corr IEMLIB_R1.15
Don't actually have the external though myself for some
reason.

Looking at the helpfile it would be suitable for audio
after using [tabsend~] to fill the tables.

On Sun, 17 Jun 2007 17:52:42 +0200
tania habib [EMAIL PROTECTED] wrote:

hello list,
I was wondering if some body can tell me if there is an object to generate
the cross correlation of  two signals.
If there is, please let me know which library has it.

regards,
tania

--
Use the source

help-tab_cross_corr.pd
Description: Binary data
___
PD-list@iem.at mailing list
UNSUBSCRIBE and account-management -
http://lists.puredata.info/listinfo/pd-list

```

Re: [PD] cross-correlation in pd ?

```
1- how can i find this IEMLIB_R1.15 library, as some of the objects in the
patch you sent me can not be created.

2- I am using pd extended 0.40 on mac osx

3- how can I include aubio in pd when I download it, a little guidance
needed here
thanks
tania

On 6/18/07, Andy Farnell [EMAIL PROTECTED] wrote:

I found a helpfile for table cross corr IEMLIB_R1.15
Don't actually have the external though myself for some
reason.

Looking at the helpfile it would be suitable for audio
after using [tabsend~] to fill the tables.

On Sun, 17 Jun 2007 17:52:42 +0200
tania habib [EMAIL PROTECTED] wrote:

hello list,
I was wondering if some body can tell me if there is an object to
generate
the cross correlation of  two signals.
If there is, please let me know which library has it.

regards,
tania

--
Use the source

___
PD-list@iem.at mailing list
UNSUBSCRIBE and account-management -
http://lists.puredata.info/listinfo/pd-list

```

Re: [PD] cross-correlation in pd ?

```

On Jun 17, 2007, at 12:46 PM, tania habib wrote:

1- how can i find this IEMLIB_R1.15 library, as some of the objects
in the patch you sent me can not be created.

2- I am using pd extended 0.40 on mac osx

In Pd-extended, that should all be in the iemlib library.

3- how can I include aubio in pd when I download it, a little
guidance needed here

That I don't know, I haven't played with it yet.

.hc

thanks
tania

On 6/18/07, Andy Farnell [EMAIL PROTECTED] wrote:

I found a helpfile for table cross corr IEMLIB_R1.15
Don't actually have the external though myself for some
reason.

Looking at the helpfile it would be suitable for audio
after using [tabsend~] to fill the tables.

On Sun, 17 Jun 2007 17:52:42 +0200
tania habib [EMAIL PROTECTED] wrote:

hello list,
I was wondering if some body can tell me if there is an object to
generate

the cross correlation of  two signals.
If there is, please let me know which library has it.

regards,
tania

--
Use the source

___
PD-list@iem.at mailing list
UNSUBSCRIBE and account-management - http://lists.puredata.info/
listinfo/pd-list

There is no way to peace, peace is the way.   -A.J. Muste

___
PD-list@iem.at mailing list
UNSUBSCRIBE and account-management -
http://lists.puredata.info/listinfo/pd-list

```

Re: [PD] cross-correlation in pd ?

```
The problem is that two of the objects I have in the patch that Andy sent me
1- tab_mul
2- tab_cross_corr

tare not created, I have the IEMLIB included in pd and the rest is ok.
Please tell me whether they are coming from same library or they are part of
some other one which needs to be included with it

Can you please check them on yours pd?

On 6/17/07, Hans-Christoph Steiner [EMAIL PROTECTED] wrote:

On Jun 17, 2007, at 12:46 PM, tania habib wrote:

1- how can i find this IEMLIB_R1.15 library, as some of the objects in the
patch you sent me can not be created.

2- I am using pd extended 0.40 on mac osx

In Pd-extended, that should all be in the iemlib library.

3- how can I include aubio in pd when I download it, a little guidance
needed here

That I don't know, I haven't played with it yet.

.hc

thanks
tania

On 6/18/07, Andy Farnell [EMAIL PROTECTED] wrote:

I found a helpfile for table cross corr IEMLIB_R1.15
Don't actually have the external though myself for some
reason.

Looking at the helpfile it would be suitable for audio
after using [tabsend~] to fill the tables.

On Sun, 17 Jun 2007 17:52:42 +0200
tania habib [EMAIL PROTECTED] wrote:

hello list,
I was wondering if some body can tell me if there is an object to
generate
the cross correlation of  two signals.
If there is, please let me know which library has it.

regards,
tania

--
Use the source

___
PD-list@iem.at mailing list
UNSUBSCRIBE and account-management -
http://lists.puredata.info/listinfo/pd-list

There is no way to peace, peace is the way.   -A.J. Muste

___
PD-list@iem.at mailing list
UNSUBSCRIBE and account-management -
http://lists.puredata.info/listinfo/pd-list

```

Re: [PD] cross-correlation in pd ?

```Hi Tania,

aubio is quite easy, provided you have libsamplerate, libsndfile and fftw3
installed.  It is pretty essential that you instal pkgconfig as well, since
some libraries rely on that to send information as to where they all reside on
the hard disk for compilation.
So, aubio in steps:

First, URLs:
pkgconfig: http://pkg-config.freedesktop.org/wiki/
libsndfile: www.mega-nerd.com/libsndfile/
libsamplerate: www.mega-nerd.com/SRC/
fftw: http://www.fftw.org/
and aubio-0.3.2.tar.gz
unpack them all in the same directory. Put the attached script in the directory
all the other packages are in then...

Then you can use this script to compile them all (see attachment). You will
need to be logged in as the administrator, and know the administrator password.
You must open Applications - Utilities - Terminal

# then type

cd /path/to/where-you-unpacked-the-files/
sudo chmod 744 instaubio.sh
password
./instaubio.sh

have a look at the script using textedit - it uses fairly common unix stuff,
and I'm sure you could work it out yourself, or you know it. Maybe I saved you
some time! ;-D but aubio needs things in the right order.

PS I know Paul is writing up his thesis right now, so maybe he's a little quiet!
PPS do you know much about linear predictive coding and formant analysis?

Best,
Ed

tania habib [EMAIL PROTECTED] wrote: The problem is that two of the objects I
have in the patch that Andy sent me
1- tab_mul
2- tab_cross_corr

tare not created, I have the IEMLIB included in pd and the rest is ok.
Please tell me whether they are coming from same library or they are part of
some other one which needs to be included with it

Can you please check them on yours pd?

On 6/17/07, Hans-Christoph Steiner [EMAIL PROTECTED] wrote:
On Jun 17, 2007, at 12:46 PM, tania habib wrote:

1- how can i find this IEMLIB_R1.15 library, as some of the objects in the
patch you sent me can not be created.

2- I am using pd extended  0.40 on mac osx

In Pd-extended, that should all be in the iemlib library.

3- how can I include aubio in pd when I download it, a little guidance needed
here

That I don't know, I haven't played with it yet.

.hc

thanks
tania

On 6/18/07, Andy Farnell [EMAIL PROTECTED] wrote:

I found a helpfile for table cross corr IEMLIB_R1.15
Don't actually have the external though myself for some
reason.

Looking at the helpfile it would be suitable for audio
after using [tabsend~] to fill the tables.

On Sun, 17 Jun 2007 17:52:42 +0200
tania habib  [EMAIL PROTECTED] wrote:

hello list,
I was wondering if some body can tell me if there is an object to generate
the cross correlation of  two signals.
If there is, please let me know which library has it.

regards,
tania

--
Use the source

___
PD-list@iem.at mailing list
UNSUBSCRIBE and account-management -
http://lists.puredata.info/listinfo/pd-list

There is no way to peace, peace is the way.   - A.J. Muste

___
PD-list@iem.at mailing list
UNSUBSCRIBE and account-management -
http://lists.puredata.info/listinfo/pd-list

Lone Shark Aviation out now on http://www.pyramidtransmissions.com
http://www.myspace.com/sharktracks
Inbox cluttering up with junk? Clean up with Yahoo! Mail.

instaubio.sh
Description: 1279115023-instaubio.sh
___
PD-list@iem.at mailing list
UNSUBSCRIBE and account-management -
http://lists.puredata.info/listinfo/pd-list

```

Re: [PD] cross-correlation in pd ?

```
No more news on where you might find the externs im afraid,
but in case you do get the iem table one working [tabsend~]
I suggested would be useless because two of them would
seem to be asynchronous and make nonsense of your
correlation (is that right?). Two triggered [tabwrites~]
as per helpfile would be the way.

On Mon, 18 Jun 2007 05:41:41 +0100
Andy Farnell [EMAIL PROTECTED] wrote:

I found a helpfile for table cross corr IEMLIB_R1.15
Don't actually have the external though myself for some
reason.

Looking at the helpfile it would be suitable for audio
after using [tabsend~] to fill the tables.

On Sun, 17 Jun 2007 17:52:42 +0200
tania habib [EMAIL PROTECTED] wrote:

hello list,
I was wondering if some body can tell me if there is an object to generate
the cross correlation of  two signals.
If there is, please let me know which library has it.

regards,
tania

--
Use the source

--
Use the source

___
PD-list@iem.at mailing list
UNSUBSCRIBE and account-management -
http://lists.puredata.info/listinfo/pd-list

```

[PD] cross-correlation in pd ?

```
On Sun, 17 Jun 2007, Tania Habib wrote:

I was wondering if some body can tell me if there is an object to
generate the cross correlation of two signals. If there is, please let
me know which library has it.

Do it with [fft~] just like for signal convolution, except that in
computing the [*~]s you replace the [-~] by a [+~].

_ _ __ ___ _  _ _ ...
| Mathieu Bouchard - tél:+1.514.383.3801, Montréal QC Canada___
PD-list@iem.at mailing list
UNSUBSCRIBE and account-management -
http://lists.puredata.info/listinfo/pd-list

```

Re: [PD] cross-correlation in pd ?

```
Mathieu's right here...  It can be done with abstractions.  I'm not
sure that mine are all correct, but I've used this technique a couple
times without making complete abstractions for it...
Here's an abstraction for computing a symmetric cross-covariance, this
way.  It outputs the cross-covariance, xcov(k)=sum(i=-32,...,31;
s1(i+k)*s2(i))
It's still pretty ugly, and the details of the math confuse me a bit.
I'm still working on the one-sided cross covariance function for
delays, like from [0,63] instead of [-32,31]

This one works within block sizes of 64, by using an [block~ 128 2] .
It is formatted very tricky using lrshift, but I think you will see
that it's continuous.  The test patch uses zexy for z~.  If you want
to compute xcov on larger block sizes,  just change all the numbers
according to the size you need.

Do you think I've got it right?  Should the output be reversed?  This
abstraction also could be improved to make running calculations of
cross-correlation, which add up the variance in s1 and s2 and
normalize the output.

Jamie Bullock has cc~ in his flib collection.  It computes cross
correlation in the time domain, and cross-covariance using the freq.
domain.  I worked on it a while back...but have been lazy about trying
to change the cross-covariance into a cross-correlation.  Sorry,
Jamie.

Chuck

On 6/17/07, Mathieu Bouchard [EMAIL PROTECTED] wrote:

On Sun, 17 Jun 2007, Tania Habib wrote:

I was wondering if some body can tell me if there is an object to
generate the cross correlation of two signals. If there is, please let
me know which library has it.

Do it with [fft~] just like for signal convolution, except that in
computing the [*~]s you replace the [-~] by a [+~].

_ _ __ ___ _  _ _ ...
| Mathieu Bouchard - tél:+1.514.383.3801, Montréal QC Canada
___
PD-list@iem.at mailing list
UNSUBSCRIBE and account-management -
http://lists.puredata.info/listinfo/pd-list

sxcov~.pd
Description: Binary data

sxcov~-test.pd
Description: Binary data

conj_mult~.pd
Description: Binary data
___
PD-list@iem.at mailing list
UNSUBSCRIBE and account-management -
http://lists.puredata.info/listinfo/pd-list

```

Re: [PD] cross-correlation in pd ?

```
On Sun, 17 Jun 2007, Mathieu Bouchard wrote:

On Sun, 17 Jun 2007, Tania Habib wrote:

I was wondering if some body can tell me if there is an object to generate
the cross correlation of two signals. If there is, please let me know which
library has it.
Do it with [fft~] just like for signal convolution, except that in computing
the [*~]s you replace the [-~] by a [+~].

... and you replace the [+~] by a [-~].

_ _ __ ___ _  _ _ ...
| Mathieu Bouchard - tél:+1.514.383.3801, Montréal QC Canada___
PD-list@iem.at mailing list
UNSUBSCRIBE and account-management -
http://lists.puredata.info/listinfo/pd-list

```