- android - 多次调用 OnPrimaryClipChangedListener
- android - 无法更新 RecyclerView 中的 TextView 字段
- android.database.CursorIndexOutOfBoundsException : Index 0 requested, 光标大小为 0
- android - 使用 AppCompat 时,我们是否需要明确指定其 UI 组件(Spinner、EditText)颜色
在我的工作(计量经济学/统计学)中,我经常需要将不同大小的矩阵相乘,然后对生成的矩阵执行其他操作。我一直依靠 bsxfun()
来矢量化代码,总的来说我发现它比 repmat()
更有效。但我不明白的是,为什么有时 bsxfun()
的性能在沿不同维度扩展矩阵时会大不相同。
考虑这个具体的例子:
x = ones(j, k, m);
beta = rand(k, m, s);
exp_xBeta = zeros(j, m, s);
for im = 1 : m
for is = 1 : s
xBeta = x(:, :, im) * beta(:, im, is);
exp_xBeta(:, im, is) = exp(xBeta);
end
end
y = mean(exp_xBeta, 3);
上下文:
我们有来自 m 个市场的数据,在每个市场中我们要计算 exp(X * beta) 的期望值,其中 X 是一个 j x k 矩阵,beta 是一个 k x 1 随机向量.我们通过蒙特卡洛积分计算这个期望 - 绘制 s 个 beta,为每个绘制计算 exp(X * beta),然后取平均值。通常我们使用 m > k > j 获取数据,并且我们使用非常大的 s。在这个例子中,我只是让 X 成为一个矩阵。
我使用 bsxfun()
做了 3 个版本的向量化,它们的区别在于 X 和 beta 的形状:
矢量化 1
x1 = x; % size [ j k m 1 ]
beta1 = permute(beta, [4 1 2 3]); % size [ 1 k m s ]
tic
xBeta = bsxfun(@times, x1, beta1);
exp_xBeta = exp(sum(xBeta, 2));
y1 = permute(mean(exp_xBeta, 4), [1 3 2 4]); % size [ j m ]
time1 = toc;
矢量化 2
x2 = permute(x, [4 1 2 3]); % size [ 1 j k m ]
beta2 = permute(beta, [3 4 1 2]); % size [ s 1 k m ]
tic
xBeta = bsxfun(@times, x2, beta2);
exp_xBeta = exp(sum(xBeta, 3));
y2 = permute(mean(exp_xBeta, 1), [2 4 1 3]); % size [ j m ]
time2 = toc;
矢量化 3
x3 = permute(x, [2 1 3 4]); % size [ k j m 1 ]
beta3 = permute(beta, [1 4 2 3]); % size [ k 1 m s ]
tic
xBeta = bsxfun(@times, x3, beta3);
exp_xBeta = exp(sum(xBeta, 1));
y3 = permute(mean(exp_xBeta, 4), [2 3 1 4]); % size [ j m ]
time3 = toc;
这就是它们的表现(通常我们使用 m > k > j 获取数据,并且我们使用了非常大的 s):
j = 5,k = 15,m = 100,s = 2000:
For-loop version took 0.7286 seconds.
Vectorized version 1 took 0.0735 seconds.
Vectorized version 2 took 0.0369 seconds.
Vectorized version 3 took 0.0503 seconds.
j = 10,k = 15,m = 150,s = 5000:
For-loop version took 2.7815 seconds.
Vectorized version 1 took 0.3565 seconds.
Vectorized version 2 took 0.2657 seconds.
Vectorized version 3 took 0.3433 seconds.
j = 15,k = 35,m = 150,s = 5000:
For-loop version took 3.4881 seconds.
Vectorized version 1 took 1.0687 seconds.
Vectorized version 2 took 0.8465 seconds.
Vectorized version 3 took 0.9414 seconds.
为什么版本 2 始终是最快的版本?最初,我认为性能优势是因为 s 设置为维度 1,Matlab 可能能够更快地计算,因为它以列优先顺序存储数据。但是 Matlab 的分析器告诉我,计算该平均值所花费的时间相当微不足道,并且在所有 3 个版本中都或多或少相同。 Matlab 大部分时间都在用 bsxfun()
评估该行,这也是三个版本中运行时差异最大的地方。
有没有想过为什么版本 1 总是最慢而版本 2 总是最快的?
我在这里更新了我的测试代码: Code
编辑:这篇文章的早期版本不正确。 beta
的大小应为 (k, m, s)
。
最佳答案
bsxfun
当然是向量化事物的好工具之一,但如果你能以某种方式引入 matrix-multiplication
那将是最好的方法,如 matrix multiplications are really fast on MATLAB
.
在这里你似乎可以使用matrix-multiplication
来得到exp_xBeta
-
[m1,n1,r1] = size(x);
n2 = size(beta,2);
exp_xBeta_matmult = reshape(exp(reshape(permute(x,[1 3 2]),[],n1)*beta),m1,r1,n2)
或者直接获取y
如下图——
y_matmult = reshape(mean(exp(reshape(permute(x,[1 3 2]),[],n1)*beta),2),m1,r1)
解释
为了更详细地解释它,我们将尺寸设置为 -
x : (j, k, m)
beta : (k, s)
我们的最终目标是使用矩阵乘法
从x
和beta
中“消除”k。因此,我们可以使用 permute
将 x
中的 k
“推”到末尾,并 reshape 为保持 k
的二维> 作为行,即 ( j * m , k ),然后用 beta
( k , s ) 执行矩阵乘法,得到 ( j * m , s )。然后可以将乘积重新整形为 3D 数组(j、m、s)并执行元素指数,这将是 exp_xBeta
。
现在,如果最终目标是 y
,即沿 exp_xBeta
的第三维获取平均值,则相当于计算沿行的平均值矩阵乘积 (j * m, s ),然后 reshape 为 ( j , m ) 以直接得到我们 y
。
关于performance - Matlab bsxfun() - 如何解释沿不同维度扩展时的性能差异?,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/30408795/
SQLite、Content provider 和 Shared Preference 之间的所有已知区别。 但我想知道什么时候需要根据情况使用 SQLite 或 Content Provider 或
警告:我正在使用一个我无法完全控制的后端,所以我正在努力解决 Backbone 中的一些注意事项,这些注意事项可能在其他地方更好地解决......不幸的是,我别无选择,只能在这里处理它们! 所以,我的
我一整天都在挣扎。我的预输入搜索表达式与远程 json 数据完美配合。但是当我尝试使用相同的 json 数据作为预取数据时,建议为空。点击第一个标志后,我收到预定义消息“无法找到任何内容...”,结果
我正在制作一个模拟 NHL 选秀彩票的程序,其中屏幕右侧应该有一个 JTextField,并且在左侧绘制弹跳的选秀球。我创建了一个名为 Ball 的类,它实现了 Runnable,并在我的主 Draf
这个问题已经有答案了: How can I calculate a time span in Java and format the output? (18 个回答) 已关闭 9 年前。 这是我的代码
我有一个 ASP.NET Web API 应用程序在我的本地 IIS 实例上运行。 Web 应用程序配置有 CORS。我调用的 Web API 方法类似于: [POST("/API/{foo}/{ba
我将用户输入的时间和日期作为: DatePicker dp = (DatePicker) findViewById(R.id.datePicker); TimePicker tp = (TimePic
放宽“邻居”的标准是否足够,或者是否有其他标准行动可以采取? 最佳答案 如果所有相邻解决方案都是 Tabu,则听起来您的 Tabu 列表的大小太长或您的释放策略太严格。一个好的 Tabu 列表长度是
我正在阅读来自 cppreference 的代码示例: #include #include #include #include template void print_queue(T& q)
我快疯了,我试图理解工具提示的行为,但没有成功。 1. 第一个问题是当我尝试通过插件(按钮 1)在点击事件中使用它时 -> 如果您转到 Fiddle,您会在“内容”内看到该函数' 每次点击都会调用该属
我在功能组件中有以下代码: const [ folder, setFolder ] = useState([]); const folderData = useContext(FolderContex
我在使用预签名网址和 AFNetworking 3.0 从 S3 获取图像时遇到问题。我可以使用 NSMutableURLRequest 和 NSURLSession 获取图像,但是当我使用 AFHT
我正在使用 Oracle ojdbc 12 和 Java 8 处理 Oracle UCP 管理器的问题。当 UCP 池启动失败时,我希望关闭它创建的连接。 当池初始化期间遇到 ORA-02391:超过
关闭。此题需要details or clarity 。目前不接受答案。 想要改进这个问题吗?通过 editing this post 添加详细信息并澄清问题. 已关闭 9 年前。 Improve
引用这个plunker: https://plnkr.co/edit/GWsbdDWVvBYNMqyxzlLY?p=preview 我在 styles.css 文件和 src/app.ts 文件中指定
为什么我的条形这么细?我尝试将宽度设置为 1,它们变得非常厚。我不知道还能尝试什么。默认厚度为 0.8,这是应该的样子吗? import matplotlib.pyplot as plt import
当我编写时,查询按预期执行: SELECT id, day2.count - day1.count AS diff FROM day1 NATURAL JOIN day2; 但我真正想要的是右连接。当
我有以下时间数据: 0 08/01/16 13:07:46,335437 1 18/02/16 08:40:40,565575 2 14/01/16 22:2
一些背景知识 -我的 NodeJS 服务器在端口 3001 上运行,我的 React 应用程序在端口 3000 上运行。我在 React 应用程序 package.json 中设置了一个代理来代理对端
我面临着一个愚蠢的问题。我试图在我的 Angular 应用程序中延迟加载我的图像,我已经尝试过这个2: 但是他们都设置了 src attr 而不是 data-src,我在这里遗漏了什么吗?保留 d
我是一名优秀的程序员,十分优秀!