[SciPy-user] unsubscribe

Vern Pierre Pierre.Vern at fr.thalesgroup.com
Fri Jul 4 12:26:34 EDT 2003


unsubscribe

-----Message d'origine-----
De : scipy-user-request at scipy.net [mailto:scipy-user-request at scipy.net]
Envoyé : vendredi 4 juillet 2003 19:00
À : scipy-user at scipy.net
Objet : SciPy-user digest, Vol 1 #433 - 6 msgs


Send SciPy-user mailing list submissions to
	scipy-user at scipy.net

To subscribe or unsubscribe via the World Wide Web, visit
	http://www.scipy.net/mailman/listinfo/scipy-user
or, via email, send a message with subject or body 'help' to
	scipy-user-request at scipy.net

You can reach the person managing the list at
	scipy-user-admin at scipy.net

When replying, please edit your Subject line so it is more specific
than "Re: Contents of SciPy-user digest..."


Today's Topics:

   1. Fwd: scipy.signal.lfilter: memory leak fixed (Andrew Straw)
   2. help for the math-challenged (Andrew Straw)
   3. Re: help for the math-challenged (Andrew Straw)
   4. Re: Re: SciPy-user digest, Vol 1 #430 - 4 msgs (Travis Oliphant)
   5. SuperLU, UMFPACK, ... (Nils Wagner)
   6. The Matrix Computation Toolbox (Nils Wagner)

--__--__--

Message: 1
Date: Thu, 3 Jul 2003 18:38:31 +0200
From: Andrew Straw <andrew.straw at adelaide.edu.au>
To: scipy-user at scipy.net
Subject: [SciPy-user] Fwd: scipy.signal.lfilter: memory leak fixed
Reply-To: scipy-user at scipy.net


--Apple-Mail-14--513782946
Content-Transfer-Encoding: 7bit
Content-Type: text/plain;
	charset=US-ASCII;
	format=flowed

From: Andrew Straw <astraw at insightscientific.com>
Date: Thu Jul 3, 2003  6:31:52  PM Europe/Berlin
To: scipy-user at scipy.net
Subject: scipy.signal.lfilter: memory leak fixed

Hi All,

I found and fixed a memory leak when lfilter was called with initial 
conditions for the filter delays set in zi.  Here's the patch.  Please 
apply it to CVS.

Cheers!
Andrew


--Apple-Mail-14--513782946
Content-Disposition: attachment;
	filename=lfilter_mem_leak.patch
Content-Transfer-Encoding: 7bit
Content-Type: application/octet-stream;
	x-unix-mode=0666;
	name="lfilter_mem_leak.patch"

Index: sigtoolsmodule.c
===================================================================
RCS file: /home/cvsroot/world/scipy/Lib/signal/sigtoolsmodule.c,v
retrieving revision 1.9
diff -c -r1.9 sigtoolsmodule.c
*** sigtoolsmodule.c	1 Feb 2003 00:10:01 -0000	1.9
--- sigtoolsmodule.c	3 Jul 2003 17:24:24 -0000
***************
*** 2135,2141 ****
  	}
  	else {
  	  free(vi); free(vf);
! 	  return Py_BuildValue("(OO)",arY,arVf);
   	}
  
   fail:
--- 2135,2141 ----
  	}
  	else {
  	  free(vi); free(vf);
! 	  return Py_BuildValue("(NN)",arY,arVf);
   	}
  
   fail:

--Apple-Mail-14--513782946--


--__--__--

Message: 2
Date: Thu, 3 Jul 2003 18:37:45 +0200
Subject: [SciPy-user] help for the math-challenged
From: Andrew Straw <andrew.straw at adelaide.edu.au>
To: scipy-user at scipy.net
Reply-To: scipy-user at scipy.net

John Hunter wrote:

>>>>>> "Andrew" == Andrew Straw <astraw at insightscientific.com> writes:
>
>     Andrew> What is the "critical frequency", as mentioned in the help
>     Andrew> for scipy.signal.iirfilter?  ("Wn -- a scalar or length-2
>     Andrew> sequence giving the critical frequencies.")  More
>     Andrew> specifically, for a 1st order filter, how is the classical
>     Andrew> time constant related to the critical frequency?  By trial
>     Andrew> and error, I find 1.0/(3.19*tau*hz) works.  In other
>     Andrew> words, this is the code I use:
>
> The critical frequency is the frequency at which there is 3dB of
> attenuation of the signal amplitude.  For a first order RC filter,
> this is 1/RC = 1/tau (Hz).

That's what I originally thought (and hence the question).  However, my 
code (which you snipped out) works, which is a state I acheived only 
after I actually playing around with the output.  Therefore, my 
question is not "is my code wrong?", but "why is my code right?".  Then 
again I suppose there could be a bug in scipy and my code actually was 
right! Nah!

Cheers!
Andrew


--__--__--

Message: 3
Date: Thu, 3 Jul 2003 18:39:29 +0200
Subject: [SciPy-user] Re: help for the math-challenged
From: Andrew Straw <andrew.straw at adelaide.edu.au>
To: scipy-user at scipy.net
Reply-To: scipy-user at scipy.net

John Hunter wrote:

>>>>>> "Andrew" == Andrew Straw <astraw at insightscientific.com> writes:
>
>     Andrew> What is the "critical frequency", as mentioned in the help
>     Andrew> for scipy.signal.iirfilter?  ("Wn -- a scalar or length-2
>     Andrew> sequence giving the critical frequencies.")  More
>     Andrew> specifically, for a 1st order filter, how is the classical
>     Andrew> time constant related to the critical frequency?  By trial
>     Andrew> and error, I find 1.0/(3.19*tau*hz) works.  In other
>     Andrew> words, this is the code I use:
>
> The critical frequency is the frequency at which there is 3dB of
> attenuation of the signal amplitude.  For a first order RC filter,
> this is 1/RC = 1/tau (Hz).

That's what I originally thought (and hence the question).  However, my 
code (which you snipped out) works, which is a state I acheived only 
after I actually playing around with the output.  Therefore, my 
question is not "is my code wrong?", but "why is my code right?".  Then 
again I suppose there could be a bug in scipy and my code actually was 
right! Nah!

Cheers!
Andrew


--__--__--

Message: 4
Date: Thu, 03 Jul 2003 13:01:44 -0600
From: Travis Oliphant <"oliphant."@ee.byu.edu>
To: scipy-user at scipy.net, camartin at snet.net
Subject: Re: [SciPy-user] Re: SciPy-user digest, Vol 1 #430 - 4 msgs
Reply-To: scipy-user at scipy.net

It looks like this line will read it

 >>> a = io.read_array('AV56tst1.INT',lines=((0,4928),),atype='l')

The problem is the last line is too short. Notice I just use the lines 
option to read in all but the last line.  (I used wc to find out how 
many lines there were).  Note the atype argument specifies long integer 
array output (instead of the default double).

We should probably change the behavior here to just insert zeros in the 
missing columns.  I will look into this too see if it makes sense with 
the other features.   The error message is a bit cryptic.



> Travis,
> 
> Attached is the file. Re-looking at the file , because it is different 
> rows than columns the last line has only three elements vs. the other 
> lines each having 10.  This is an interferometry file in ASCII. I've 
> stripped the headers because I know how to handle them and what I'm 
> going to do. Thanks for the help.
> 
> Cliff Martin


--__--__--

Message: 5
Date: Fri, 04 Jul 2003 09:54:26 +0200
From: Nils Wagner <nwagner at mecha.uni-stuttgart.de>
To: "scipy-user at scipy.net" <scipy-user at scipy.net>
Subject: [SciPy-user] SuperLU, UMFPACK, ...
Reply-To: scipy-user at scipy.net

Hi all,

How do I invoke SuperLU and UMFPACK in scipy ?
A small example illustrating the main features would be appreciated.

Thanks in advance.

Nils


--__--__--

Message: 6
Date: Fri, 04 Jul 2003 10:03:18 +0200
From: Nils Wagner <nwagner at mecha.uni-stuttgart.de>
To: scipy-dev at scipy.net,
	"scipy-user at scipy.net" <scipy-user at scipy.net>
Subject: [SciPy-user] The Matrix Computation Toolbox
Reply-To: scipy-user at scipy.net

Hi all,

There exists an interesting matrix computation toolbox that has been
developed by Nicholas J. Higham.
http://www.maths.man.ac.uk/~higham/mctoolbox/

The toolbox is distributed under the terms of the GNU General Public
License (version 2 of the License, or any later version) 
as published by the Free Software Foundation. 

Maybe the features of scipy can benefit from this toolbox. However, it
is a collection of Matlab m-files.

Nils



--__--__--

_______________________________________________
SciPy-user mailing list
SciPy-user at scipy.net
http://www.scipy.net/mailman/listinfo/scipy-user


End of SciPy-user Digest




More information about the SciPy-User mailing list