English 中文(简体)
在NumPy计算Tp,q,r = Sp,q Sp,r的最佳方法是什么?
原标题:What s the best way to calculate Tp,q,r = Sp,q Sp,r in NumPy?
  • 时间:2023-10-04 02:19:40
  •  标签:
  • python
  • numpy

I m 试图将每一件零件(<>subm,n)的系列米投入序列转换成一个帐篷,其中tp,q,r = ssp,r

虽然我的法典确实可行,但我认为必须找到更好的解决办法。 这里是我所说的话。

# nseqs is the number of sequences
# seq_length is the sequence length
# seq is a list of sequences

output = np.empty((nseqs, seq_length, seq_length))
for n in range(nseqs):
    for i, j in enumerate(seq[n]):
        output[n, i, :] = j
    output[n, :, :] *= output[n, :, :].T

更确切地说,有办法对<编码>产出加以改动。 不要忘记,重复阶段是在一个步骤中进行的,没有这些漏洞?

问题回答

你没有给出任何数字或实例,因此我只举了一个小小小小小例子(如果你没有实际工作的话,就不得不告诉我)。

Consider the following:

import numpy as np

nseqs = 4
seq_length = 5

seq = np.arange(nseqs*seq_length).reshape(nseqs, seq_length)

从内部通道开始,我们可以通过广播实现同样的目标。

output = np.empty((nseqs, seq_length, seq_length))
for n in range(nseqs):
    output2[n] = seq[n,None]
    output[n, :, :] *= output[n, :, :].T

Now, the multiplication doesn t actually affect the subsequent loops, so we can split this into two loops.

output = np.empty((nseqs, seq_length, seq_length))
for n in range(nseqs):
    output2[n] = seq[n,None]
for n in range(nseqs):
    output[n, :, :] *= output[n, :, :].T

从第一波段播放的节目只是重复每个次马基里的seq[n]的内容,因此,我们能够利用>>np.repeat

output = np.repeat(seq[:,None], seq_length, axis=1)

第二处是代号:代号:<>output[n]乘以代号。 可使用https://numpy.org/doc/stable/vis/generated/numpy.swapaxes.html” rel=“nofollow noreferer”np.swapaxes<>, swapping the final two axes (i.e. exercise the transpose of the submatrices)。

output *= np.swapaxes(output, 1, 2)

最后,你的法典削减到:

output = np.repeat(seq[:,None], seq_length, axis=1)
output *= np.swapaxes(output, 1, 2)




相关问题
Can Django models use MySQL functions?

Is there a way to force Django models to pass a field to a MySQL function every time the model data is read or loaded? To clarify what I mean in SQL, I want the Django model to produce something like ...

An enterprise scheduler for python (like quartz)

I am looking for an enterprise tasks scheduler for python, like quartz is for Java. Requirements: Persistent: if the process restarts or the machine restarts, then all the jobs must stay there and ...

How to remove unique, then duplicate dictionaries in a list?

Given the following list that contains some duplicate and some unique dictionaries, what is the best method to remove unique dictionaries first, then reduce the duplicate dictionaries to single ...

What is suggested seed value to use with random.seed()?

Simple enough question: I m using python random module to generate random integers. I want to know what is the suggested value to use with the random.seed() function? Currently I am letting this ...

How can I make the PyDev editor selectively ignore errors?

I m using PyDev under Eclipse to write some Jython code. I ve got numerous instances where I need to do something like this: import com.work.project.component.client.Interface.ISubInterface as ...

How do I profile `paster serve` s startup time?

Python s paster serve app.ini is taking longer than I would like to be ready for the first request. I know how to profile requests with middleware, but how do I profile the initialization time? I ...

Pragmatically adding give-aways/freebies to an online store

Our business currently has an online store and recently we ve been offering free specials to our customers. Right now, we simply display the special and give the buyer a notice stating we will add the ...

Converting Dictionary to List? [duplicate]

I m trying to convert a Python dictionary into a Python list, in order to perform some calculations. #My dictionary dict = {} dict[ Capital ]="London" dict[ Food ]="Fish&Chips" dict[ 2012 ]="...

热门标签