Python中文
首页
教程
问答
标签
搜索
登录
注册
cython:numpy数组中可能的优化
回答此问题可获得
20
贡献值,回答如果被采纳可获得
50
分。
<p>下面是我从多元正态分布中提取的Cython代码。我使用循环,因为每次我有不同的密度。(conLSigma是Cholesky因子)</p> <p>这需要很多时间,因为我要对每个循环进行逆分解和Cholesky分解。它比纯python代码更快,但我想知道是否有任何方法可以提高速度。在</p> <pre><code>from __future__ import division import numpy as np cimport numpy as np ctypedef np.float64_t dtype_t cimport cython @cython.boundscheck(False) @cython.wraparound(False) def drawMetro(np.ndarray[dtype_t, ndim = 2] beta, np.ndarray[dtype_t, ndim = 3] H, np.ndarray[dtype_t, ndim = 2] Sigma, float s): cdef int ncons = betas.shape[0] cdef int nX = betas.shape[1] cdef int con cdef np.ndarray betas_cand = np.zeros([ncons, nX], dtype = np.float64) cdef np.ndarray conLSigma = np.zeros([nX, nX], dtype = np.float64) for con in xrange(ncons): conLSigma = np.linalg.cholesky(np.linalg.inv(H[con] + Sigma)) betas_cand[con] = betas[con] + s * np.dot(conLSigma, np.random.standard_normal(size = nX)) return(betas_cand) </code></pre>
0 条评论
分类:
Python问答
请先
登录
后评论
默认排序
时间排序
1 个回答
匿名
1天前
擅长:python、mysql、java
<p>先计算cholesky分解,然后反代换求下三角矩阵。这应该比利纳格乔尔斯基(利纳格投资公司(S) )。在</p>
请先
登录
后评论
针对此问题:
更多的回答
关注
89
关注
收藏
1
收藏,
216
浏览
网友 提问于 2天前
相关Python问题
我想从用户inpu创建一个类的实例
7 回答
我想从用户导入值,为此
7 回答
我想从用户那里得到一个整数输入,然后让for循环遍历该数字,然后调用一个函数多次
9 回答
我想从用户那里收到一个列表,并在其中执行一些步骤,然后在步骤完成后将其打印回来,但它没有按照我想要的方式工作
5 回答
我想从用户那里获取输入,并将值传递给(average=dict[x]/6),然后在那里获取resu
2 回答
我想从第一个列表中展示第一个词,然后从第二个列表中展示十个词,以此类推- Python
2 回答
我想从第一个空lin开始解析文本文件
8 回答
我想从简历、简历中提取特定部分
5 回答
我想从给定字典(python)的字符串中删除\u00a9、\u201d和类似的字符。
6 回答
我想从给定的网站Lin下载许多文件扩展名相同的Wget或Python文件
9 回答
我想从网上搜集一些关于抵押贷款的数据
9 回答
我想从网站上删除电子邮件地址
10 回答
我想从网站上读取数据该网站包含可下载的文件,然后我想用python脚本把它发送给oracle如何?
10 回答
我想从网站中提取数据,然后将其显示在我的网页上
10 回答
我想从网页上提取统计数据。
6 回答
我想从网页上解析首都城市,并在用户输入国家时在终端上打印它们
4 回答
我想从色彩图中删除前n个颜色,而不丢失原始颜色数
3 回答
我想从课堂上打印字典里的键
10 回答
我想从费用表中获取学生上次支付的费用,其中学生id=id
4 回答
我想从较低的顺序对多重列表进行排序,但我无法在一行中生成结果
10 回答