擅长:python、mysql、java
<p>作为沃伦解决方案的一个替代方案,我认为这是最好的,这里有一个未记录的<code>inner1d</code>:</p>
<pre><code>>>> from numpy.core.umath_tests import inner1d
>>> a = inner1d(h, c)
>>> np.allclose(a, ans)
True
</code></pre>
<p>从其docstring:</p>
<blockquote>
<p><strong>inner1d(x1, x2[, out])</strong></p>
<p>inner on the last dimension and broadcast on the rest: (i),(i)->()</p>
</blockquote>
<p>对于这个特殊情况,在我的系统中,<code>inner1d</code>比<code>np.einsum</code>稍快:</p>
^{pr2}$