Menu
▾
▴
numpy-discussion — Discussion list for all users of Numerical Python
You can subscribe to this list here.
| 2000 |
Jan
(8) |
Feb
(49) |
Mar
(48) |
Apr
(28) |
May
(37) |
Jun
(28) |
Jul
(16) |
Aug
(16) |
Sep
(44) |
Oct
(61) |
Nov
(31) |
Dec
(24) |
|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 2001 |
Jan
(56) |
Feb
(54) |
Mar
(41) |
Apr
(71) |
May
(48) |
Jun
(32) |
Jul
(53) |
Aug
(91) |
Sep
(56) |
Oct
(33) |
Nov
(81) |
Dec
(54) |
| 2002 |
Jan
(72) |
Feb
(37) |
Mar
(126) |
Apr
(62) |
May
(34) |
Jun
(124) |
Jul
(36) |
Aug
(34) |
Sep
(60) |
Oct
(37) |
Nov
(23) |
Dec
(104) |
| 2003 |
Jan
(110) |
Feb
(73) |
Mar
(42) |
Apr
(8) |
May
(76) |
Jun
(14) |
Jul
(52) |
Aug
(26) |
Sep
(108) |
Oct
(82) |
Nov
(89) |
Dec
(94) |
| 2004 |
Jan
(117) |
Feb
(86) |
Mar
(75) |
Apr
(55) |
May
(75) |
Jun
(160) |
Jul
(152) |
Aug
(86) |
Sep
(75) |
Oct
(134) |
Nov
(62) |
Dec
(60) |
| 2005 |
Jan
(187) |
Feb
(318) |
Mar
(296) |
Apr
(205) |
May
(84) |
Jun
(63) |
Jul
(122) |
Aug
(59) |
Sep
(66) |
Oct
(148) |
Nov
(120) |
Dec
(70) |
| 2006 |
Jan
(460) |
Feb
(683) |
Mar
(589) |
Apr
(559) |
May
(445) |
Jun
(712) |
Jul
(815) |
Aug
(663) |
Sep
(559) |
Oct
(930) |
Nov
(373) |
Dec
|
|
From: Bill B. <wb...@gm...> - 2006-10-26 02:52:38
|
On 10/26/06, jeremito <jer...@gm...> wrote: > > > Hmmm. I guess I'll have to find another way to get some > "standard/well-known" matrices to work with. Thanks for your help. Well if all you want is some matrices, there's nothing stopping you from grabbing the matrices in the LAPACK distribution and using them yourself. Robert's just saying they won't be included in Numpy. There's also the matrix market, whcih has a large number of (sparse-only?) example matrices. http://math.nist.gov/MatrixMarket/index.html --bb |
|
From: jeremito <jer...@gm...> - 2006-10-26 02:43:19
|
Robert Kern wrote: > Robert Kern wrote: > > jeremito wrote: > >> If I understand correctly, the LAPACK package has some standard > >> matrices included with it for testing. Does the linalg portion of > >> numpy also have access to these matrices/arrays? > > > > No. > > To clarify, these are routines provided with the LAPACK source distribution in > the TESTING/MATGEN/ directory. However, they are not part of the LAPACK library > itself, and thus we *cannot* provide wrappers for them in numpy. > > -- > Robert Kern > > "I have come to believe that the whole world is an enigma, a harmless enigma > that is made terrible by our own mad attempt to interpret it as though it had > an underlying truth." > -- Umberto Eco > Hmmm. I guess I'll have to find another way to get some "standard/well-known" matrices to work with. Thanks for your help. Jeremy |
|
From: <tr...@16...> - 2006-10-26 00:31:30
|
<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
<HTML><HEAD><TITLE>新建网页 1</TITLE>
<META http-equiv=Content-Language content=zh-cn>
<META http-equiv=Content-Type content="text/html; charset=gb2312">
<META content="MSHTML 6.00.2900.2523" name=GENERATOR></HEAD>
<BODY>
<TABLE id=table1 width="100%" border=1>
<TBODY>
<TR>
<TD>
<P style="LINE-HEIGHT: 200%" align=center><FONT color=#0000ff><FONT
size=2>主办单位： 易腾企业管理咨询有限公司</FONT><SPAN lang=zh-cn><FONT
size=6><B><BR>供应商管理与采购成本降低技巧研修班</B></FONT></SPAN></FONT></P>
<P style="LINE-HEIGHT: 150%"><B><SPAN class=noo><FONT color=#0000ff
size=3>课程背景</FONT></SPAN></B><BR><FONT lang=ZH-CN size=2> 采购人员每年在做降价工作，但企业为了控制库存，采购周期越来越短、采购批量越来越小，供应商怨声载道，加上原材料的价格不断上涨，降价的工作越来越富有挑战。通过对本课程的学习学员可以了解:<br>
现代采购管理的发展趋势,改善企业的采购组织以及采购流程的设定，完善供应商管理体系，提升采购谈判能力。从而帮助采购人员选择最佳供应商和采购策略，确保采购工作高质量、高效率及低成本执行，使企业具有最佳的供货状态，同时与供应商保持良好的战略伙伴关系。</FONT><BR><SPAN
class=r3><B><FONT color=#0000ff size=3>课程大纲</FONT></B></SPAN><BR>
<FONT
face=宋体 size=2><font color="#0000FF">一、采购管理的新型理念</font><br>
从职能管理转向过程管理<br>
从采购管理转向供应管理<br>
从企业间交易性管理转向关系性管理<br>
从零和竞争转向多赢竞争<br>
从简单的多元化经营转向核心竞争力管理<br>
<font color="#0000FF">二、采购组织设计与流程优化</font><br>
供应链管理管理模式下采购职能的新定位<br>
从采购员到采购工程师的角色转换<br>
采购部门与PMC部门的高效率运作关系<br>
采购流程设计的原则<br>
如何优化采购流程缩短事务型工作的处理周期<br>
<font color="#0000FF">三、供应商选择与管理</font><br>
如何依据公司发展目标制定供应商需求体系<br>
供应商开发与认可程序<br>
采购商业体系、质量体系的构建<br>
供应商选择与评价的考评因素<br>
供应商商业、服务、质量审核要素<br>
批量生产中供应商的日常管理<br>
供应商定期评估、等级划分与双赢模式建立<br>
◆案例分析<BR><font color="#0000FF">四、如何管理供应商</font><br>
供应商之交期管理（Delivery）<br>
供应商品质管理（Quality ）<br>
供应商成本管理（Cost）<br>
供应商服务管理（Service）<br>
如何管理主要供应商管理<br>
建立与策略性供应商的伙伴关系<br>
如何管理单一供应商<br>
如何实施供应商的早期参与<br>
◆案例研讨<br>
<font color="#0000FF">五、如何进行供应商绩效评估</font><br>
建立供应商绩效考核标准<br>
供应商绩效分析<br>
供应商绩效考评的关键指标<br>
如何奖励优秀供应商，促进其继续进步<br>
如何协助供应商改善绩效<br>
如何进行供应商发展<br>
◆案例研讨</FONT><FONT
size=2><BR><font color="#0000FF">六、采购核心价值与采购成本控制</font><br>
成本导向采购管理的目标<br>
成本导向采购管理的成功要素<br>
主动采购取代被动采购<br>
如何核算采购总成本,正确决策<br>
采购中的成本影响因素分析<br>
供应商通常依据哪些要素进行报价<br>
如何分析供应商的价格市场定位与走向<br>
如何运用价格分析工具来分析报价<br>
价值分析/价值评价（VA/VE ）<br>
如何实施有效的招标<br>
小批量、多批次和大批量、少批次的供货成本权衡<br>
如何根据物资类别建立低成本的供应合作关系<br>
◆ 案例研讨<br>
<font color="#0000FF">七、采购谈判的步骤与结构</font><br>
计划准备阶段：如何进行谈判方案设计<br>
谈判开始阶段：该阶段特点分析及对策<br>
谈判过渡阶段：该阶段特点分析及对策<br>
实质性谈判阶段：该阶段特点分析及对策<br>
交易明确阶段：该阶段特点分析及对策<br>
谈判结束阶段：谈判总过程回顾；收关；拟定协议<br>
模拟谈判（确立角色，设立目标、策略实施）<br>
<font color="#0000FF">八、采购谈判技巧</font><br>
如何掌握卖方真实的销售心理<br>
如何分析销售方的需求<br>
如何运用技术分析手段实现“不谈”的谈判<br>
如何利用买卖双方的优劣进行谈判<br>
如何利用各级别的权限进行议价<br>
买方占优势时应采用何种采购策略<br>
卖方占优势时应采用何种采购策略<br>
现货采购时的谈判策略<br>
订单式供应时的谈判策略<br>
◆ 案例演示<br>
<font color="#0000FF">九、采购绩效管理的实施与评估</font><br>
通过改善供应商绩效来提升采购绩效表现<br>
通过采购职能的整合提升采购绩效表现<br>
通过有效的成本管理方法提升采购成本竞争力<br>
通过供应链管理技术实现材料库存合理化<br>
采购绩效年度总结<br>
十、综合案例分析与实际问题解答<br>
</FONT><B><FONT color=#0000ff
size=3>导师简介</FONT></B><BR> <FONT color=#000000 size="2">jack.wu，易腾管理顾问高级顾问师、物流与供应</FONT><FONT color=#000000 size="2">专业课程讲授专家。吴老师在外资企业的供应链管理领域具有十多年的工作经验,曾担任IBM中国区采购总监,结合多年的实际操作经验和丰富的管理咨询经验，讲授或辅导过的企业
包括TDK、松下、联想、华为、中兴、汇源果汁、可口可乐、百事可乐、Cadbury（吉百利）、华硕电脑、美的空调、厦华集团、中国国际航空、ABB、万达集团、深圳威尼斯酒店、上海粤海酒店、NEC东金电子、步步高电子、太原钢铁集团、PHILIPS、大冷王运输制冷、三洋华强、TCL、荣仕达、好孩子、清华同方、正大集团、伊莱克斯等两百多家,受到广泛好评.</FONT><FONT
color=#ff0000 size=2><BR></FONT><B><FONT
size=2 color="#0000ff">时间地点:</FONT> </B><FONT face=宋体 size=2>11月4-5日</FONT><FONT
size=2> (周六日) 上海 <BR></FONT><B><FONT
color=#0000ff size=2>费用:</FONT></B><FONT size=2> 1800/人<FONT
face=宋体>(含课程费、教材、午餐等)</FONT> 四人以上参加,赠予一名名额<BR>报名/咨询:
021-51187131 张小姐 <br>
注:如您不需要此邮件,请来电说明。 </FONT></P></TD></TR></TBODY></TABLE></BODY></HTML>
|
|
From: Robert K. <rob...@gm...> - 2006-10-26 00:13:51
|
Robert Kern wrote: > jeremito wrote: >> If I understand correctly, the LAPACK package has some standard >> matrices included with it for testing. Does the linalg portion of >> numpy also have access to these matrices/arrays? > > No. To clarify, these are routines provided with the LAPACK source distribution in the TESTING/MATGEN/ directory. However, they are not part of the LAPACK library itself, and thus we *cannot* provide wrappers for them in numpy. -- Robert Kern "I have come to believe that the whole world is an enigma, a harmless enigma that is made terrible by our own mad attempt to interpret it as though it had an underlying truth." -- Umberto Eco |
|
From: Robert K. <rob...@gm...> - 2006-10-26 00:06:54
|
jeremito wrote: > If I understand correctly, the LAPACK package has some standard > matrices included with it for testing. Does the linalg portion of > numpy also have access to these matrices/arrays? No. -- Robert Kern "I have come to believe that the whole world is an enigma, a harmless enigma that is made terrible by our own mad attempt to interpret it as though it had an underlying truth." -- Umberto Eco |
|
From: jeremito <jer...@gm...> - 2006-10-26 00:04:20
|
If I understand correctly, the LAPACK package has some standard matrices included with it for testing. Does the linalg portion of numpy also have access to these matrices/arrays? Thanks, Jeremy |
|
From: Tom D. <tom...@al...> - 2006-10-25 21:32:58
|
Thanks Robert. Will do. On 10/25/06, Robert Kern <rob...@gm...> wrote: > > Tom Denniston wrote: > > Hi, I've been using the lexsort and it is really nice. Much simpler and > > faster than the implementation I had before. The only problem is it > > seems to crash on strings. > > > > Do I just need to upgrade numpy or is it a current bug? > > Upgrade numpy. I believe this bug has been fixed. > > http://projects.scipy.org/scipy/numpy/ticket/298 > > -- > Robert Kern > > "I have come to believe that the whole world is an enigma, a harmless > enigma > that is made terrible by our own mad attempt to interpret it as though > it had > an underlying truth." > -- Umberto Eco > > > ------------------------------------------------------------------------- > Using Tomcat but need to do more? Need to support web services, security? > Get stuff done quickly with pre-integrated technology to make your job > easier > Download IBM WebSphere Application Server v.1.0.1 based on Apache Geronimo > http://sel.as-us.falkag.net/sel?cmd=lnk&kid=120709&bid=263057&dat=121642 > _______________________________________________ > Numpy-discussion mailing list > Num...@li... > https://lists.sourceforge.net/lists/listinfo/numpy-discussion > |
|
From: Andrew J. <a.h...@gm...> - 2006-10-25 21:29:22
|
OS X 10.4, PPC In [5]: import numpy as N In [6]: N.__version__ Out[6]: '1.0.1.dev3399' In [7]: print N.longfloat(0).itemsize 16 In [8]: a = N.exp(N.array([1000],dtype=N.longfloat)) In [9]: print str(a) [inf] In [10]: print str(a[0]) inf In [11]: print a.itemsize 16 |
|
From: Brian G. <ell...@gm...> - 2006-10-25 20:46:19
|
Pearu, Fantastic, I will try this out. Thanks for your work on this. Brian On 10/25/06, pe...@ce... <pe...@ce...> wrote: > > On Wed, 25 Oct 2006, Brian Granger wrote: > > > HI, > > > > I have noticed that Pearu has been doing lots of work on f2py g3. I > > have heard that this will have support for derived types. What is the > > status of this work? Is is ready for public use? Are the derived > > types implemented? Thanks! > > At the moment derived types support is implemented for types that have > only scalar components. Work on the array component support is going > on. If you would like to try out f2py g3, use --3g-numpy option: > > f2py --3g-numpy ... > > See numpy/f2py/lib/test_*.py files for examples that currently work with > f2py g3. > > Pearu > > > ------------------------------------------------------------------------- > Using Tomcat but need to do more? Need to support web services, security? > Get stuff done quickly with pre-integrated technology to make your job easier > Download IBM WebSphere Application Server v.1.0.1 based on Apache Geronimo > http://sel.as-us.falkag.net/sel?cmd=lnk&kid=120709&bid=263057&dat=121642 > _______________________________________________ > Numpy-discussion mailing list > Num...@li... > https://lists.sourceforge.net/lists/listinfo/numpy-discussion > |
|
From: <pe...@ce...> - 2006-10-25 19:36:16
|
On Wed, 25 Oct 2006, Brian Granger wrote: > HI, > > I have noticed that Pearu has been doing lots of work on f2py g3. I > have heard that this will have support for derived types. What is the > status of this work? Is is ready for public use? Are the derived > types implemented? Thanks! At the moment derived types support is implemented for types that have only scalar components. Work on the array component support is going on. If you would like to try out f2py g3, use --3g-numpy option: f2py --3g-numpy ... See numpy/f2py/lib/test_*.py files for examples that currently work with f2py g3. Pearu |
|
From: Robert K. <rob...@gm...> - 2006-10-25 19:28:48
|
Tom Denniston wrote: > Hi, I've been using the lexsort and it is really nice. Much simpler and > faster than the implementation I had before. The only problem is it > seems to crash on strings. > > Do I just need to upgrade numpy or is it a current bug? Upgrade numpy. I believe this bug has been fixed. http://projects.scipy.org/scipy/numpy/ticket/298 -- Robert Kern "I have come to believe that the whole world is an enigma, a harmless enigma that is made terrible by our own mad attempt to interpret it as though it had an underlying truth." -- Umberto Eco |
|
From: Tom D. <tom...@al...> - 2006-10-25 18:53:40
|
Hi, I've been using the lexsort and it is really nice. Much simpler and
faster than the implementation I had before. The only problem is it seems
to crash on strings.
Do I just need to upgrade numpy or is it a current bug?
In [2]: import numpy
nu
In [3]: numpy.array(['a', 'b', 'c'])
Out[3]:
array([a, b, c],
dtype='|S1')
In [4]: strArray = numpy.array(['a', 'b', 'c'])
In [5]: intArray = numpy.array([1,2,3])
In [6]: numpy.__version__
Out[6]: '1.0b5'
In [7]: numpy.lexsort((intArray,))
Out[7]: array([0, 1, 2])
In [8]: numpy.lexsort((intArray,intArray))
Out[8]: array([0, 1, 2])
In [9]: numpy.lexsort((strArray,intArray))
Segmentation fault (core dumped)
|
|
From: Brian G. <ell...@gm...> - 2006-10-25 17:21:35
|
HI, I have noticed that Pearu has been doing lots of work on f2py g3. I have heard that this will have support for derived types. What is the status of this work? Is is ready for public use? Are the derived types implemented? Thanks! Brian |
|
From: David M. C. <co...@ph...> - 2006-10-25 15:45:42
|
On Oct 22, 2006, at 17:38 , Simon Burton wrote: > > Apologies if this is too off-topic for these lists, but > I hope some people here find this interesting! It looks very nifty. Interval arithmetic for rigourous bounds is something that should be used more often, I think. You should add it to the topical software list at http:// www.scipy.org/Topical_Software. [and it's on topic, btw; announcements of Python numerical software is always welcome here]. > PyDX - first public announcement > ================================ > > Overview > -------- > > PyDX is a package for working with calculus (differential > geometry), > arbitrary precision arithmetic (using gmpy), and interval arithmetic. > > It provides: > * multivariate automatic differentiation (to arbitrary order) > * Tensor objects for computations in differential geometry > * Interval scalars (based on libMPFI) for calculating rigorous > bounds > on arithmetic operations (validated numerics). > * Arbitrary order validated ODE solver. > > PyDX uses lazy computation techniques to greatly enhance > performance of > the resulting functions. > > This code grew out of a research project at The Australian National > University, > Department of Physics, which involved computing bounds on geodesics in > relativistic space-time's. -- |>|\/|< /------------------------------------------------------------------\ |David M. Cooke http://arbutus.physics.mcmaster.ca/dmc/ |co...@ph... |
|
From: David M. C. <co...@ph...> - 2006-10-25 14:38:58
|
On Oct 22, 2006, at 22:05 , Jeremy R. Fishman wrote: > Hi, I was wondering if anyone can give me some advice on how to go > about cross-compiling NumPy. I have been searching around and > can't find any support in distutils for cross compilation. Is > there some way I can still compile Numerical Python using a mipsel- > linux compiler, on say a Cygwin host? > > The end result would be to use NumPy with an embedded Python > interpreter running on a mipsel host. In specific, The OPEN-R > system by Sony for use on Aibo robotic dogs. I need to figure out > how to get setup to use the mipsel compiler, not how to set up the > compiler. I suppose you could run the config step (using the host's compiler), then go and fix up the generated config files in the build/ directory. Hmm, looks to me if you run 'python setup.py build_src', then go poke around in build/src.<something>, you should find a numpy/ __config__.py, numpy/distutils/__config__.py (which look to be the same file, for some reason), and a numpy/core/config.h I *think* that's all you you'd have to fiddle. Then you should be able to build the modules by setting CC to your cross compiler, and doing 'python build_ext build_py' (Almost; looks like one config test is done in build_ext). Hope that helps; patches to make this cleaner are welcome. -- |>|\/|< /------------------------------------------------------------------\ |David M. Cooke http://arbutus.physics.mcmaster.ca/dmc/ |co...@ph... |
|
From: Lisandro D. <da...@gm...> - 2006-10-25 13:44:53
|
FYI, I got the following warning building numpy-1.0 from sources on a
FC2 box, using distribution provided blas and lapack rpm's.
$ rpm -q blas
blas-3.0-24.1
$ rpm -q lapack
lapack-3.0-24.1
$ python setup.py build > build.log
Running from numpy source directory.
/u/dalcinl/Build/Python/numpy-1.0/numpy/distutils/system_info.py:934:
UserWarning:
*********************************************************************
Lapack library (from ATLAS) is probably incomplete:
size of /usr/lib/liblapack.so is 3617k (expected >4000k)
Follow the instructions in the KNOWN PROBLEMS section of the file
numpy/INSTALL.txt.
*********************************************************************
warnings.warn(message)
My lapack shared library is only 3.6MB, But I believe it is complete.
Am I right?
$ ls -lh /usr/lib/liblapack.so.3.0.3
-rwxr-xr-x 1 root root 3.6M Mar 2 2004 /usr/lib/liblapack.so.3.0.3
--=20
Lisandro Dalc=EDn
---------------
Centro Internacional de M=E9todos Computacionales en Ingenier=EDa (CIMEC)
Instituto de Desarrollo Tecnol=F3gico para la Industria Qu=EDmica (INTEC)
Consejo Nacional de Investigaciones Cient=EDficas y T=E9cnicas (CONICET)
PTLC - G=FCemes 3450, (3000) Santa Fe, Argentina
Tel/Fax: +54-(0)342-451.1594
|
|
From: David H. <dav...@gm...> - 2006-10-25 13:25:51
|
[snip] Vincent, I may be wrong, but it looks like you could use a weighted histogramdd (see other thread where I sent a patch to do so). You run it once to find the number of sample in each bin, then again with weights=value, normed=False to find the total weight in each bin, then you can divide each bin by the number of samples to find the average. N, d = histogramdd(cells) W, d = histogramdd(cells, weights = value) average = W/N David |
|
From: Albert S. <fu...@gm...> - 2006-10-25 13:16:40
|
Hello all
I'm trying to get SciPy compiling with MSVC and the Intel Visual Fortran
compiler, version 9.1. I'm running into problems with the version
detection code failing.
The current version_cmd in intel.py is:
'version_cmd' : [fc_exe, "-FI -V -c %(fname)s.f -o %(fname)s.o" \
% {'fname':dummy_fortran_file()}],
Even when running inside a Command Prompt with the necessary environment
variables for the compiler set up, this command causes the following
error:
LINK : fatal error LNK1104: cannot open file 'ifconsol.lib'
Changing version_cmd to to pass /help /quiet instead of trying to
compile a dummy file does the trick. (/quiet so that the compiler
doesn't warn about it being called ifort now, causing the regular
expression that detects the version to fail.)
Can we fix this? Is it important to try to compile the dummy file?
Cheers,
Albert
|
|
From: David H. <dav...@gm...> - 2006-10-25 13:13:47
|
Nicolas, thanks for the bug report, I fooled around with argument passing and should have checked every case. You'll find the histogram function that deals with weights on the numpy tra= c ticket 189, <http://projects.scipy.org/scipy/numpy/ticket/189> I'm waiting for some hints as to where the histogram function should reside (numpy.histogram, numpy.stats.histogram, ...) before submitting a patch . Salut, David 2006/10/25, Nicolas Champavert <nic...@ob...>: > > Hi, > > it would be great if you could add the weight option in the 1D > histogram too. > > Nicolas > > David Huard a =E9crit : > > Xavier, > > Here is the patch against svn. Please report any bug. I haven't had > > the time to test it extensively, something that should be done before > > commiting the patch to the repo. I'd appreciate your feedback. > > > > David > > > > 2006/10/24, David Huard <dav...@gm... > > <mailto:dav...@gm...>>: > > > > Hi Xavier, > > > > You could tweak histogram2d to do what you want, or you could give > > me a couple of days and I'll do it and let you know. If you want > > to help, you could write a test using your particular application > > and data. > > > > David > > > > > > 2006/10/24, Xavier Gnata < gn...@ob... > > <mailto:gn...@ob...>>: > > > > Hi, > > > > I have a set of 3 1D large arrays. > > The first 2 one stand for the coordinates of particules and > > the last one > > for their masses. > > I would like to be able to plot this data ie to compute a 2D > > histogram > > summing the masses in each bin. > > I cannot find a way to do that without any loop on the indices > > resulting > > too a very slow function. > > > > I'm looking for an elegant way to do that with numpy (or > > scipy??) function. > > > > For instance, scipy.histogram2d cannot do the job because it > > only counts > > the number of samples in each bin. > > There is no way to deal with weights. > > > > Xavier. > > > > > > -- > > ############################################ > > Xavier Gnata > > CRAL - Observatoire de Lyon > > 9, avenue Charles Andr=E9 > > 69561 Saint Genis Laval cedex > > Phone: +33 4 78 86 85 28 > > Fax: +33 4 78 86 83 86 > > E-mail: gn...@ob... <mailto:gn...@ob...= > > > ############################################ > > > > > > > ------------------------------------------------------------------------- > > Using Tomcat but need to do more? Need to support web > > services, security? > > Get stuff done quickly with pre-integrated technology to make > > your job easier > > Download IBM WebSphere Application Server v.1.0.1 based on > > Apache Geronimo > > > http://sel.as-us.falkag.net/sel?cmd=3Dlnk&kid=3D120709&bid=3D263057&dat= =3D121642 > > < > http://sel.as-us.falkag.net/sel?cmd=3Dlnk&kid=3D120709&bid=3D263057&dat= =3D121642> > > _______________________________________________ > > Numpy-discussion mailing list > > Num...@li... > > <mailto:Num...@li...> > > https://lists.sourceforge.net/lists/listinfo/numpy-discussion > > > > > > > > > > ------------------------------------------------------------------------- > Using Tomcat but need to do more? Need to support web services, security? > Get stuff done quickly with pre-integrated technology to make your job > easier > Download IBM WebSphere Application Server v.1.0.1 based on Apache Geronim= o > http://sel.as-us.falkag.net/sel?cmd=3Dlnk&kid=3D120709&bid=3D263057&dat= =3D121642 > _______________________________________________ > Numpy-discussion mailing list > Num...@li... > https://lists.sourceforge.net/lists/listinfo/numpy-discussion > |
|
From: Nicolas C. <nic...@ob...> - 2006-10-25 12:25:42
|
Hi,
it would be great if you could add the weight option in the 1D=20
histogram too.
Nicolas
David Huard a =E9crit :
> Xavier,
> Here is the patch against svn. Please report any bug. I haven't had=20
> the time to test it extensively, something that should be done before=20
> commiting the patch to the repo. I'd appreciate your feedback.
>
> David
>
> 2006/10/24, David Huard <dav...@gm...=20
> <mailto:dav...@gm...>>:
>
> Hi Xavier,
>
> You could tweak histogram2d to do what you want, or you could give
> me a couple of days and I'll do it and let you know. If you want
> to help, you could write a test using your particular application
> and data.
>
> David
>
>
> 2006/10/24, Xavier Gnata < gn...@ob...
> <mailto:gn...@ob...>>:
>
> Hi,
>
> I have a set of 3 1D large arrays.
> The first 2 one stand for the coordinates of particules and
> the last one
> for their masses.
> I would like to be able to plot this data ie to compute a 2D
> histogram
> summing the masses in each bin.
> I cannot find a way to do that without any loop on the indices
> resulting
> too a very slow function.
>
> I'm looking for an elegant way to do that with numpy (or
> scipy??) function.
>
> For instance, scipy.histogram2d cannot do the job because it
> only counts
> the number of samples in each bin.
> There is no way to deal with weights.
>
> Xavier.
>
>
> --
> ############################################
> Xavier Gnata
> CRAL - Observatoire de Lyon
> 9, avenue Charles Andr=E9
> 69561 Saint Genis Laval cedex
> Phone: +33 4 78 86 85 28
> Fax: +33 4 78 86 83 86
> E-mail: gn...@ob... <mailto:gn...@ob...=
>
> ############################################
>
>
> ---------------------------------------------------------------=
----------
> Using Tomcat but need to do more? Need to support web
> services, security?
> Get stuff done quickly with pre-integrated technology to make
> your job easier
> Download IBM WebSphere Application Server v.1.0.1 based on
> Apache Geronimo
> http://sel.as-us.falkag.net/sel?cmd=3Dlnk&kid=3D120709&bid=3D26=
3057&dat=3D121642
> <http://sel.as-us.falkag.net/sel?cmd=3Dlnk&kid=3D120709&bid=3D2=
63057&dat=3D121642>
> _______________________________________________
> Numpy-discussion mailing list
> Num...@li...
> <mailto:Num...@li...>
> https://lists.sourceforge.net/lists/listinfo/numpy-discussion
>
>
>
>
|
|
From: Nicolas C. <nic...@ob...> - 2006-10-25 11:56:04
|
Hi,
there is a problem when using this patch of numpy.histogram2d.
I've tried:
a=3Dpylab.rand(10)
numpy.histogram2d(a,a)
and there is an exception (exceptions.AttributeError) in=20
numpy.histogramdd function called by numpy.histogram2d.
But numpy.histogram2d(a,a,bins=3D[10,10]) works.
So I've tried to use directly numpy.histogramdd function:
numpy.histogramdd([a,a])
and it works.
I've looked at the histogram2d code:
try:
N =3D len(bins)
except TypeError:
N =3D 1
bins =3D [bins]
I think the last line (bins =3D [bins]) have to be deleted because if yo=
u=20
do so, numpy.histogram2d(a,a) works.
Or maybe the good think to do is:
try:
N =3D len(bins)
except TypeError:
N =3D 1
bins =3D 2*[bins]
Nicolas
David Huard a =E9crit :
> Xavier,
> Here is the patch against svn. Please report any bug. I haven't had=20
> the time to test it extensively, something that should be done before=20
> commiting the patch to the repo. I'd appreciate your feedback.
>
> David
>
> 2006/10/24, David Huard <dav...@gm...=20
> <mailto:dav...@gm...>>:
>
> Hi Xavier,
>
> You could tweak histogram2d to do what you want, or you could give
> me a couple of days and I'll do it and let you know. If you want
> to help, you could write a test using your particular application
> and data.
>
> David
>
>
>
>
> =20
>
>
>
>
>
>
> 2006/10/24, Xavier Gnata < gn...@ob...
> <mailto:gn...@ob...>>:
>
> Hi,
>
> I have a set of 3 1D large arrays.
> The first 2 one stand for the coordinates of particules and
> the last one
> for their masses.
> I would like to be able to plot this data ie to compute a 2D
> histogram
> summing the masses in each bin.
> I cannot find a way to do that without any loop on the indices
> resulting
> too a very slow function.
>
> I'm looking for an elegant way to do that with numpy (or
> scipy??) function.
>
> For instance, scipy.histogram2d cannot do the job because it
> only counts
> the number of samples in each bin.
> There is no way to deal with weights.
>
> Xavier.
>
>
> --
> ############################################
> Xavier Gnata
> CRAL - Observatoire de Lyon
> 9, avenue Charles Andr=E9
> 69561 Saint Genis Laval cedex
> Phone: +33 4 78 86 85 28
> Fax: +33 4 78 86 83 86
> E-mail: gn...@ob... <mailto:gn...@ob...=
>
> ############################################
>
>
> ---------------------------------------------------------------=
----------
> Using Tomcat but need to do more? Need to support web
> services, security?
> Get stuff done quickly with pre-integrated technology to make
> your job easier
> Download IBM WebSphere Application Server v.1.0.1 based on
> Apache Geronimo
> http://sel.as-us.falkag.net/sel?cmd=3Dlnk&kid=3D120709&bid=3D26=
3057&dat=3D121642
> <http://sel.as-us.falkag.net/sel?cmd=3Dlnk&kid=3D120709&bid=3D2=
63057&dat=3D121642>
> _______________________________________________
> Numpy-discussion mailing list
> Num...@li...
> <mailto:Num...@li...>
> https://lists.sourceforge.net/lists/listinfo/numpy-discussion
>
>
>
> -----------------------------------------------------------------------=
-
>
> =20
|
|
From: A. M. A. <per...@gm...> - 2006-10-25 11:44:39
|
On 25/10/06, Travis Oliphant <oli...@ie...> wrote: > It takes old "type-codes" as one big string so you say > > vectorize(f,otypes='d') > > This should be modernized to handle a list of dtype objects. > > I've fixed vectorize in SVN. Thanks! A. M. Archibald |
|
From: Travis O. <oli...@ie...> - 2006-10-25 11:38:22
|
A. M. Archibald wrote: > Hi, > > Vectorize is a very handy function, but it has at least one pitfall: > > def f(x): > if 1.3<x<2.5: > return sqrt(x) > else: > return 0 > > Now vectorize(f)(2)=1.41421356237 but vectorize(f)(array([1,2]))=array([0,1]). > > The problem is that, when given an array as input, vectorize feeds in > the first element, looks at the return type, and returns an array of > that type - and I didn't put a "." after the zero. > > This should perhaps be in the docstring of vectorize, since I can't > see any way to work around it, but it can easily lead to > difficult-to-find bugs. It may seem like an artificial example, but it > came up with a function I was working on. But it's confusing > behaviour. > > vectorize appears to support an "otypes" argument, but it doesn't take > standard numpy type objects, and it doesn't do anything obvious. > It takes old "type-codes" as one big string so you say vectorize(f,otypes='d') This should be modernized to handle a list of dtype objects. I've fixed vectorize in SVN. -Travis |
|
From: Vincent S. <sc...@sa...> - 2006-10-25 11:07:39
|
Andreas Eisele wrote: > Recently, there were several requests and discussions on this list about > how to > increment an array a in cells pointed to from a second integer array b > (optionally by > values from a third array c), such as: > > >> Yes, that'd be >> a[b] += c >> >> On 10/8/06, Daniel Mahler <dmahler@gm...> wrote: >> > Is there a 'loop free' way to do this in Numeric >> > >> > for i in arange(l): >> > a[b[i]]+=c[i] >> > >> > where l == len(b) == len(c) >> > >> > thanks >> > Daniel >> >> >> > or > >> It is clear to me that the numpy += operator in combination with the use >> of arrays of indexes, as is explained in the Tentative Numpy Tutorial >> >> (http://www.scipy.org/Tentative_NumPy_Tutorial#head-3f4d28139e045a442f78c5218c379af64c2c8c9e), >> >> the limitation being that indexes that appear more than 1 time in the >> indexes-array will get incremented only once. >> >> Does anybody know a way to work around this? >> >> I am using this to fill up a custom nd-histogram, and obviously each bin >> should be able to get incremented more than once. Looping over the >> entire array and incrementing each bin succesively takes waaay to long >> (these are pretty large arrays, like 4000x2000 items, or even larger) >> > I just came across a function that seems to provide the solution to both > requests, > which is called bincount. > > The first usecase could be written as > > a += bincount(b,c) > > (assuming a has already the right dimension, otherwise a = bincount(b,c) > would create an > array with the minimal required size), the second case is even simpler: > > counts = bincount(index) > > On my machine, this does 20M counting operations per second, which is _much_ > faster than anything that could be done in an explicit for loop. > > Hope this helps, > > Andreas > Andreas, thanks for this tip! And thanks to this, I stumbled across the related function 'digitize' which is also very useful for me. Now the only problem left is that bincount has no way to deal with nd-histograms (where you have multiple index arrays, that together point to a bin in a multi-dimensional histogram 'grid'). Anyone any ideas about that? Only thing I can think of is to create a unique number for each possible nd index combination and use that as an indermediate step in bincount... Would probably work ok and fast enough, but gets picky when you use lots of bins (so your unique numbers will need to be very large). In my case (something like 10x10x10 bins) it would however be OK I guess. Anyway thanks a lot for sharing this. VS > > > ------------------------------------------------------------------------- > Using Tomcat but need to do more? Need to support web services, security? > Get stuff done quickly with pre-integrated technology to make your job easier > Download IBM WebSphere Application Server v.1.0.1 based on Apache Geronimo > http://sel.as-us.falkag.net/sel?cmd=lnk&kid=120709&bid=263057&dat=121642 > _______________________________________________ > Numpy-discussion mailing list > Num...@li... > https://lists.sourceforge.net/lists/listinfo/numpy-discussion > |
|
From: A. M. A. <per...@gm...> - 2006-10-25 10:33:56
|
Hi,
Vectorize is a very handy function, but it has at least one pitfall:
def f(x):
if 1.3<x<2.5:
return sqrt(x)
else:
return 0
Now vectorize(f)(2)=1.41421356237 but vectorize(f)(array([1,2]))=array([0,1]).
The problem is that, when given an array as input, vectorize feeds in
the first element, looks at the return type, and returns an array of
that type - and I didn't put a "." after the zero.
This should perhaps be in the docstring of vectorize, since I can't
see any way to work around it, but it can easily lead to
difficult-to-find bugs. It may seem like an artificial example, but it
came up with a function I was working on. But it's confusing
behaviour.
vectorize appears to support an "otypes" argument, but it doesn't take
standard numpy type objects, and it doesn't do anything obvious.
Perhaps I should file a ticket, but first I'd like to understand the
correct behaviour...
A. M. Archibald
|
258 messages has been excluded from this view by a project administrator.
