Mailman 3 November 2002 - SciPy-User

Re:[SciPy-user] Way to get values which are the same on two arrays
by Michael Sorich Nov. 4, 2002

Nov. 4, 2002

Hi José, I remember reading your question last week. Coincidently, I have just written a function to find the intersection of 2 sets (using a histogram function) for my own work. I'm not sure if it is the most efficient method, but it seems to work. from Numeric import * def histogram(a, bins): """counts the number of values in a that fall into each of the bins in bin""" n = searchsorted(sort(a),bins) n = concatenate([n,[len(a)]]) return asarray(n[1:]-n[:-1]) def intersection(a,b): """takes two 1D arrays and returns an array of values that are in both arrays""" if a.typecode() != 'l' or b.typecode() != 'l': raise ValueError, "arrays must be of type Int" #make bins from intersection of ranges minab = max([minimum.reduce(a),minimum.reduce(b)]) maxab = min([maximum.reduce(a),maximum.reduce(b)]) # -> If either set is empty, or their ranges don't intersect if maxab < minab or maxab < 0: return None bins = arange(minab,maxab+1) return minab+nonzero(logical_and(histogram(a,bins),histogram(b,bins))) if __name__ == '__main__': a = array([1,2,3,4,6]) b = array([4,5,6,1,7]) print intersection(a,b) Hope this helps, Michael Sorich PhD Student School of Pharmaceutical, Molecular and Biomedical Sciences University of South Australia Email: michael.sorich(a)postgrads.unisa.edu.au mike_sorich(a)hotmail.com -------------------------------- Hi, I seem to remember a Numerical function that could produce the list of values (or array positions) which are common to two arrays. In other words: a1=[1,2,3,9,10] a2=[1,4,9,10] result is [1,9,10]. Browsing through the Numerical docs doesn't refresh my memory, so any kind souls? many thanks, José -- José L Gómez Dans PhD student Tel: +44 114 222 5582 Radar & Communications Group FAX; +44 870 132 2990 Department of Electronic Engineering University of Sheffield UK --- Outgoing mail is certified Virus Free. Checked by AVG anti-virus system (http://www.grisoft.com). Version: 6.0.408 / Virus Database: 230 - Release Date: 24/10/2002

1 0

signal.iirdesign throws a TypeError on Win32, any idea why, or what I'm doing wrong ?
by John A Ferguson Nov. 3, 2002

Nov. 3, 2002

My Setup is binary installs as follows onto Win2K SP3: Python-2.2.2.exe win32all-148.exe wxPythonWIN32-2.3.3.1-Py22.exe Numeric-22.0.win32-py2.2.exe SciPy-0.2.0_alpha_145.4398.win32-py2.2.exe Python 2.2.2 (#37, Oct 14 2002, 17:02:34) [MSC 32 bit (Intel)] on win32 Type "help", "copyright", "credits" or "license" for more information. >>> from scipy import signal >>> signal.iirdesign(0.2, 0.3, 3, 40) Traceback (most recent call last): File "<stdin>", line 1, in ? File "C:\Python22\Lib\site-packages\scipy\signal\filter_design.py", line 378, in iirdesign return iirfilter(N, Wn, rp=gpass, rs=gstop, analog=analog, btype=btype, ftyp e=ftype, output=output) File "C:\Python22\Lib\site-packages\scipy\signal\filter_design.py", line 448, in iirfilter z, p, k = typefunc(N, rp, rs) File "C:\Python22\Lib\site-packages\scipy\signal\filter_design.py", line 1052, in ellipap disp=0) File "C:\Python22\Lib\site-packages\scipy\optimize\optimize.py", line 145, in fmin fsim[0] = apply(func,(x0,)+args) File "C:\Python22\Lib\site-packages\scipy\signal\filter_design.py", line 1016, in kratio k = special.ellipk([m,1-m]) TypeError: bad argument type for built-in operation >>>

1 0

Vector plot using gplt
by Neilen Nov. 3, 2002

Nov. 3, 2002

Hi I'd like to do a vector plot of a 2D field. If I remember correctly, the matlab quiver() command would have made this easy. Looking at the gnuplot manual, I see it has a "set arrow" command, which looks like it might help me. Is there any way I can talk to gnuplot on the level where I can do set command from scipy? Or are there any other (easier?) possibilities? Thanks Neilen -- all we are waiting for is something worth waiting for --- KMFDM

1 0

Including LaTeX commands in the output of a gplt plot
by Neilen Nov. 3, 2002

Nov. 3, 2002

Hi. I have, since my last mail, also found out that I can get output straight into LaTeX. Sweet. Now, it would be even more fantastic if I can include \whatever command in the titles. It seems gplt strips off the \. Is there any way this may be done? Thanks Neilen -- all we are waiting for is something worth waiting for --- KMFDM

1 0

Getting postscript out of gplt
by Neilen Nov. 3, 2002

Nov. 3, 2002

Hi. I'm trying to get some postscript output of a gplt plot. Is there somewhere I can find more documentation on the gplt.output function? basically, I want the plot to be black and white, with a white background. I tried replacing the png example from the plotting tutorial with ps, eps and postscript, to no avail. Thanks Neilen -- all we are waiting for is something worth waiting for --- KMFDM

1 1

Urgent Assistance
by BUNDI KYARI Nov. 2, 2002

Nov. 2, 2002

FROM:BUNDI KYARI. TEL: E-mail: kyaribundi(a)netscape.net AMSTERDAM THE NETHERLANDS. Dear sir, compliment of the day, I am BUNDI KYARI, The son of late General Kubwa Kyari of the Democratic Republic of Congo. My father was a General in the Congolese Army. In his position (My father) with the office of the presidentcy during the regime of Laurent Kabila, he was assigned on a secret mission to source and acquire arms internationally in order to strengthen the Government forces against the rebels, which already had the support of Rwandan and Uganda Army. Meanwhile, he was still negotiating for the purchase of the arms, he received on the 16th January 2001 news of the assassination of Laurent Kabila which force him to call off the assignment and deposited the sum of US$12.5M, Packed in a diplomatic case in a private security company in the Amsterdam, the Netherlands, though he registered the content as precious stones while the real content is (US12.5M) meant for the purchase of arms for the Congolese Army. My father went home for the funeral of the late president, but on his arrival he was arrested, detained and tortured, unfortunately my father suffer cardiac arrest and died on the 17th of March 2001. However, on one of our numerous visits, my mother and I paid him while in prison, my father was able to reveal this secret to me and advice that i should proceed to the Netherlands to claim the money, he handed me all the relevant documents that will enable me claim the box from the security company. Already, I have made my first visit to the security company and the availability of this box have been confirmed. On our arrival in the Netherlands few months ago, we sought for political asylum; which was granted. My mother and I are making frantic effort on the best way to handle this money. We sought advice from an attorney who advised that we must seek for a trustworthy foreign business partner whom this money could be transferred into his/her (company"s) account This we view as the best option because our refugee status dose not permit us to operate a bank account, hence we seek your assistance and hope you could be trusted. I got your contact from the commercial section of the congolese embassy in Belgium.Meanwhile, I sincerely ask for your assistance to get this money through your account, Your share for assisting us will be 25% of the total sum, 5% will be use for upsetting all the expenses incurred in the course of concluding this venture and the remaining 70% that will be for me and my family. Also you stand to gain from any investment you might introduce us into after the conclusion of the transfer. Please keep this confidential until we finalize and get this money into your account for security reasons. This is my e-mail address you can reach me: kyaribundi(a)netscape.net Thanks and GOD bless MR, BUNDI KYARI.

1 0

Compiling with weave and g++ 3.2
by Leon Widdershoven Nov. 1, 2002

Nov. 1, 2002

Hi, I wrote some code using weave.inline, and this code would not compile due to random_access_iterator: no such class. I read that this has been an issue, but could not find an answer. I have in the mean time found a way to compile the code: in Config.hxx I changed STANDARD_LIBRARY_HAS_ITERATOR_TRAITS 1 to STANDARD_LIBRARY_HAS_ITERATOR_TRAITS 0 I Don't know what it does, but at least my code runs again. Regards, Leon Widdershoven

1 0