Re: [PD] cross-correlation in pd ?

2007-06-21 Thread Charles Henry

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 ?

2007-06-20 Thread Charles Henry

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 ?

2007-06-19 Thread Jamie Bullock

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 ?

2007-06-18 Thread Georg Holzmann
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 ?

2007-06-18 Thread Ed Kelly
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 ?

2007-06-18 Thread Charles Henry

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


Re: [PD] cross-correlation in pd ?

2007-06-17 Thread cyrille henry


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 ?

2007-06-17 Thread Andy Farnell


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 ?

2007-06-17 Thread tania habib

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 ?

2007-06-17 Thread Hans-Christoph Steiner


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 ?

2007-06-17 Thread tania habib

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 ?

2007-06-17 Thread Ed Kelly
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 ?

2007-06-17 Thread Andy Farnell

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


Re: [PD] cross-correlation in pd ?

2007-06-17 Thread Charles Henry

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 ?

2007-06-17 Thread Mathieu Bouchard

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