在python numpy中，如果我用10^6长度随机生成的list生成numpy array,那么生成耗时0.1s, 但是得到这个array的mean只需要init的2%的时间。 而我自己implement的array得到mean需要十几秒。
所以numpy的本文来源[email protected]搞@^&代*@码2网array十分黑科技是应为:
1)用底层代码太厉害？
2)init的时候partially compute了某一些中间量？(应为求mean的时间比access慢，比O(n)快 )
如果是2的话能否讲一下大概思路(不需要用python O(n)就能得mean)？
感激不禁！
回复内容：
numpy的许多函数不仅是用C实现了，还使用了BLAS（一般Windows下link到MKL的，Linux下link到OpenBLAS）。基本上那些BLAS实现在每种操作上都进行了高度优化，例如使用AVX向量指令集，甚至能比你自己用C实现快上许多，更不要说和用Python实现的比。。你用blas试试numpy底层使用BLAS做向量，矩阵运算。像求平均值这种vector operation，很容易使用multi-threading或者vectorization来加速。比如MKL就有很多优化。
<code class="language-python"><span class="n">a</span><span class="o">=</span><span class="p">[];</span><span class="n">s</span><span class="o">=</span><span class="mi">0</span><span class="p">;</span><span class="n">n</span><span class="o">=</span><span class="mi">1000000</span><span class="kn">from</span> <span class="nn">time</span> <span class="kn">import</span><span class="o">*</span><span class="kn">from</span> <span class="nn">math</span> <span class="kn">import</span><span class="o">*</span><span class="kn">from</span> <span class="nn">random</span> <span class="kn">import</span><span class="o">*</span><span class="n">st</span><span class="o">=</span><span class="n">clock</span><span class="p">()</span><span class="k">for</span> <span class="n">i</span> <span class="ow">in</span> <span class="nb">range</span><span class="p">(</span><span class="n">n</span><span class="p">):</span> <span class="n">a</span><span class="o">.</span><span class="n">append</span><span class="p">(</span><span class="n">random</span><span class="p">())</span><span class="k">for</span> <span class="n">i</span> <span class="ow">in</span> <span class="n">a</span><span class="p">:</span><span class="n">s</span><span class="o">=</span><span class="n">s</span><span class="o">+</span><span class="n">i</span><span class="n">et</span><span class="o">=</span><span class="n">clock</span><span class="p">()</span><span class="k">print</span> <span class="s">"mean="</span><span class="p">,</span><span class="n">s</span><span class="o">/</span><span class="n">n</span><span class="p">,</span><span class="s">"time="</span><span class="p">,</span><span class="n">et</span><span class="o">-</span><span class="n">st</span><span class="p">,</span><span class="s">"seconds"</span></code>
搞代码网（gaodaima.com）提供的所有资源部分来自互联网，如果有侵犯您的版权或其他权益，请说明详细缘由并提供版权或权益证明然后发送到邮箱[email protected]，我们会在看到邮件的第一时间内为您处理，或直接联系QQ：872152909。本网站采用BY-NC-SA协议进行授权
转载请注明原文链接：为什么numpy的array那么快？
转载请注明原文链接：为什么numpy的array那么快？
