- html - 出于某种原因,IE8 对我的 Sass 文件中继承的 html5 CSS 不友好?
- JMeter 在响应断言中使用 span 标签的问题
- html - 在 :hover and :active? 上具有不同效果的 CSS 动画
- html - 相对于居中的 html 内容固定的 CSS 重复背景?
假设我有一个 numpy 数组
X = np.array([[1,2,3],
[4,5,6],
[7,8,9]])
我想通过添加(在左侧)通过将所有可能的列对相乘而得到的列来扩展此矩阵。在这个例子中它将变成
X = np.array([[1, 2, 3, 2, 6],
[4,5,6,20,24,30],
[7,8,9,56,63,72]])
其中第四列是 X[:,0]
和 X[:,1]
的乘积,第五列是 X 的乘积[:,0]
和 X[:,2]
第六列是 X[:,1]
和 X[ 的乘积:,2]
.
我的尝试
我想为此使用np.hstack
。不过我也知道使用循环会减慢一切速度,但我不知道如何在没有循环的情况下正确完成它。
for i in range(matrix.shape[1]-1):
for j in range(matrix.shape[1])[i:]:
matrix2 = np.hstack((matrix, (matrix[:,i]*matrix[:,j]).reshape(-1,1))).copy()
这样做的问题是它很慢,而且我必须使用不同的矩阵,否则它会继续添加列..有更好的主意吗?
最佳答案
方法#1
使用np.triu_indices
获取成对列索引。使用它们来选择通过输入数组的列索引获得的两组 block 。使用这些 block 执行逐元素乘法,最后使用 np.concatenate
将它们作为新列与输入数组堆叠在一起。
因此,实现 -
n = X.shape[1]
r,c = np.triu_indices(n,1)
out0 = X[:,r] * X[:,c]
out = np.concatenate(( X, out0), axis=1)
方法#2
对于内存效率和性能,尤其是大型数组,更多灵感来自 this post
通过循环配对组 -
m,n = X.shape
N = n*(n-1)//2
idx = np.concatenate(( [0], np.arange(n-1,0,-1).cumsum() ))+n
start, stop = idx[:-1], idx[1:]
out = np.empty((m,n+N),dtype=X.dtype)
out[:,:n] = X
for j,i in enumerate(range(n-1)):
out[:, start[j]:stop[j]] = X[:,i,None]*X[:,i+1:]
In [403]: X = np.random.randint(0,9,(10,100))
In [404]: %timeit app1(X)
...: %timeit app2(X)
...:
1000 loops, best of 3: 277 µs per loop
1000 loops, best of 3: 350 µs per loop
In [405]: X = np.random.randint(0,9,(10,1000))
In [406]: %timeit app1(X)
...: %timeit app2(X)
...:
10 loops, best of 3: 68.6 ms per loop
100 loops, best of 3: 12.5 ms per loop
In [407]: X = np.random.randint(0,9,(10,2000))
In [408]: %timeit app1(X)
...: %timeit app2(X)
...:
1 loop, best of 3: 311 ms per loop
10 loops, best of 3: 44.8 ms per loop
关于Python:添加包含列元素成对乘积的列的最快方法,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/46220362/
我正在 csv 上使用 hadoop 来分析一些数据。我使用sql/mysql(不确定)来分析数据,现在陷入了僵局。 我花了好几个小时在谷歌上搜索,却没有找到任何相关的东西。我需要一个查询,在该查询中
我正在为 Bootstrap 网格布局的“简单”任务而苦苦挣扎。我希望在大视口(viewport)上有 4 列,然后在中型设备上有 2 列,最后在较小的设备上只有 1 列。 当我测试我的代码片段时,似
对于这个令人困惑的标题,我深表歉意,我想不出这个问题的正确措辞。相反,我只会给你背景信息和目标: 这是在一个表中,一个人可能有也可能没有多行数据,这些行可能包含相同的 activity_id 值,也可
具有 3 列的数据库表 - A int , B int , C int 我的问题是: 如何使用 Sequelize 结果找到 A > B + C const countTasks = await Ta
我在通过以下功能编写此查询时遇到问题: 首先按第 2 列 DESC 排序,然后从“不同的第 1 列”中选择 只有 Column1 是 DISTINCT 此查询没有帮助,因为它首先从第 1 列中进行选择
使用 Bootstrap 非常有趣和有帮助,目前我在创建以下需求时遇到问题。 “使用 bootstrap 在桌面上有 4 列,在平板电脑上有 2 列,在移动设备上有 1 列”谁能告诉我正确的结构 最佳
我是 R 新手,正在问一个非常基本的问题。当然,我在尝试从所提供的示例中获取指导的同时做了功课here和 here ,但无法在我的案例中实现这个想法,即可能是由于我的问题中的比较维度更大。 我的实
通常我会使用 R 并执行 merge.by,但这个文件似乎太大了,部门中的任何一台计算机都无法处理它! (任何从事遗传学工作的人的附加信息)本质上,插补似乎删除了 snp ID 的 rs 数字,我只剩
我有一个 df , delta1 delta2 0 -1 2 0 -1 0 0 0 我想知道如何分配 delt
您好,我想知道是否可以执行以下操作。显然,我已经尝试在 phpMyAdmin 中运行它,但出现错误。也许还有另一种方式来编写此查询。 SELECT * FROM eat_eat_restaurants
我有 2 个列表(标题和数据值)。我想要将数据值列 1 匹配并替换为头文件列 1,以获得与 dataValue 列 1 和标题值列 2 匹配的值 头文件 TotalLoad,M0001001 Hois
我有两个不同长度的文件,file2 是一个很大的引用文件,我从中提取文件 1 的数据。 我有一行 awk,我通常会对其进行调整以在我的文件中进行查找和替换,但它总是在同一列中进行查找和替换。 所以对于
假设我有两个表,如下所示。 create table contract( c_ID number(1) primary key, c_name varchar2(50) not
我有一个带有 varchar 列的 H2 表,其检查约束定义如下: CONSTRAINT my_constraint CHECK (varchar_field <> '') 以下插入语句失败,但当我删
这是最少量的代码,可以清楚地说明我的问题: One Two Three 前 2 个 div 应该是 2 个左列。第三个应该占据页面的其余部分。最后,我将添加选项来隐藏和
在 Azure 中的 Log Analytics 中,我为 VM Heartbeat 选择一个预定义查询,我在编辑器中运行查询正常,但当我去创建警报时,我不断收到警报“查询未返回 TimeGenera
在 Azure 中的 Log Analytics 中,我为 VM Heartbeat 选择一个预定义查询,我在编辑器中运行查询正常,但当我去创建警报时,我不断收到警报“查询未返回 TimeGenera
今天我开始使用 JexcelApi 并遇到了这个:当您尝试从特定位置获取元素时,不是像您通常期望的那样使用sheet.getCell(row,col),而是使用sheet.getCell(col,ro
我有一个包含 28 列的数据库。第一列是代码,第二列是名称,其余是值。 public void displayData() { con.Open(); MySqlDataAdapter
我很沮丧:每当我缩小这个网页时,一切都变得一团糟。我如何将网页居中,以便我可以缩小并且元素不会被错误定位。 (它应该是 2 列,但所有内容都合并为 1)我试过 但由于某种原因,这不起作用。 www.o
我是一名优秀的程序员,十分优秀!