PDL 成对行比较

Question

我创建了一个 PDL 矩阵。我需要在每一行之间进行成对比较。目前我正在使用“where”和“cov”命令返回两个切片的成对比较（在 perl 循环中生成）。

我的问题：如何使用“范围”和“切片”以成对的方式循环行？如何返回我的索引位置？我已经使用 perl 循环了矩阵。我读过，使用 perl 循环确实削弱了 PDL 的功能。

所需的输出：

indexA indexB Value
pos1   pos5   1
pos1   pos6   5
pos1   pos0   7 

需要明确的是，我只想使用 PDL 功能。

这是一些伪代码（希望）可以更好地说明我的观点。

p $b

[
 [1 0 3 0]
 [0 1 0 1]
 [1 3 1 3]   <- example piddle y
 [0 1 0 1]   <- example piddle z
]

my concept function{


slice $b (grab row z) - works fine
slice $b (grab row y) - works fine


($a, $b) = where($a,$b, $a < 3 && $b < 3 ) - works fine

p $a [1 1] 
p $b [0  0] 

cov($a $b) - works just fine.

}

我只需要一种在所有行中成对执行的方法。我需要进行阶乘（n 行）比较。

Answer 1

PDL 线程是您在这里寻找的概念。沿维度循环的一般技术是在适当的位置添加虚拟维度，以便计算生成所需的隐式线程循环。对于多维问题，可以有多种不同的方法来添加暗淡，从而创建线程循环。

对于成对行计算，您可以在切片索引上选择两个嵌套循环，其中对两个索引计数进行 Perl 循环，并将沿行生成 PDL 线程。您可以仅对索引使用一个 Perl 循环，但利用隐式线程循环一次性计算所有行。

完全 PDL 线程循环计算是为每个参数的行循环添加一个虚拟维度，以便您可以一次计算整个 N**2 行计算。下面是形状 [4,3] 数组的示例，其计算为 == 运算符：

pdl> $b = floor(random(4,3)*5)

pdl> p $b

[
 [0 4 3 3]
 [3 3 4 2]
 [4 0 1 4]
]

pdl> p $b(,*3)==$b(,,*3)

[
 [
  [1 1 1 1]
  [0 0 0 0]
  [0 0 0 0]
 ]
 [
  [0 0 0 0]
  [1 1 1 1]
  [0 0 0 0]
 ]
 [
  [0 0 0 0]
  [0 0 0 0]
  [1 1 1 1]
 ]
]

结果是形状 [4,3,3] piddle，其第 0 维对应于成对计算产生的行，并且第一个和第二个暗淡对应于 == 操作中涉及的行索引。

如果您需要这些线程循环计算之一的索引值，请使用 xvals、yvals、zvals 或 axisvals< /em> 生成一个 piddle，其索引值对应于该数组轴。

pdl> p $b->xvals

[
 [0 1 2 3]
 [0 1 2 3]
 [0 1 2 3]
]


pdl> p $b->yvals

[
 [0 0 0 0]
 [1 1 1 1]
 [2 2 2 2]
]

有很多与 PDL 线程的实现相关的细节（与 perl 线程或 posix 线程不同）。我推荐 perldl 邮件列表供其他 PDL 用户和开发人员参考和讨论。另请参阅 PDL Book 的第一个在线草稿，其中包含更全面地覆盖 PDL 计算和线程。

Answer 2

我认为您正在寻找一种方法来查找数组中所有不同的行对，然后使用 cov 处理每对行？如果这是正确的，那么我还没有听说过 cov，并且快速搜索文档也没有帮助。不过我可以说一些可能有帮助的事情。

我认为您对于从 PDL 中删除到 Perl 代码中过于谨慎，如果您所做的只是循环遍历所有行对的索引并使用 切片。下面的一些示例代码显示了这一点。

另外，您不能像 $a $a 这样调用 where 。 3 等本身就是 piddle，布尔运算符不会对它们执行您想要的操作。请改用 & 运算符，并添加一些括号以确保表达式以正确的顺序执行。

除此之外，除非您纠正我对问题的理解或指导我查看 cov 子例程的一些文档，否则我无法提供帮助。

use strict;
use warnings;

use PDL;

my $dat = pdl <<END;
[
 [1 0 3 0]
 [0 1 0 1]
 [1 3 1 3]
 [0 1 0 1]
]
END

my $max2 = $dat->dim(1) - 1;

for my $i (0 .. $max2 - 1) {
  for my $j ($i + 1 .. $max2) {

    my $row1 = $dat->slice(",($i)");
    my $row2 = $dat->slice(",($j)");

    ($row1, $row2) = where($row1, $row2, ($row1 < 3) & ($row2 < 3));

    cov($row1, $row2);
  }
}