Описание тега numpy-einsum

Функция einsum в NumPy реализует соглашение Эйнштейна о суммировании для многомерных объектов массива. Используйте этот тег, чтобы задать вопросы о том, как применить ʻeinsum` к конкретной проблеме в NumPy, или задать дополнительные вопросы о том, как работает функция.
1 ответ

Python Pandas SUMPRODUCT и L Матрица расчета

Я должен столбцы в формате панд данных и хочу вывод в столбец C D, как показано ниже A B C D 1 2 1*2 1*2 3 4 (1+3)*4 (1*2)+(3*4) 5 6 (1+3+5)*6 (1*2)+(3*4)+(5*6) 7 8 (1+3+5+7)*8 (1*2)+(3*4)+(5*6)+(7*8) 9 10 .... ..... Здесь я пытаюсь закодировать фор…
03 янв '18 в 09:02
1 ответ

Умножающие тензоры, содержащие изображения в numy

У меня есть следующие тензоры 3-го порядка. Оба тензора являются матрицами: первый тензор содержит 100 матриц 10х9, а второй - 100 матриц 3х10 (которые я только что заполнил единицами для этого примера). Моя цель состоит в том, чтобы умножить матриц…
1 ответ

Векторизация тройного произведения NumPy на двумерном массиве

Я пытаюсь векторизовать следующую тройную операцию продукта на N x N массив называется p ниже: for j in range(len(p)): for k in range(len(p)): for l in range(len(p)): h[j, k, l] = p[j, k] * p[k, l] * p[l, j] - p[j, l] * p[l, k] * p[k, j] я думал num…
18 июл '16 в 01:54
1 ответ

Обработка верхних треугольных элементов только с помощью NumPy einsum

Я использую numy einsum для вычисления точечных произведений массива векторов столбцов pts, формы (3,N), с самим собой, в результате чего получается матричная точка, формы (N,N) со всеми точечными произведениями. Это код, который я использую: dotps …
13 апр '16 в 13:20
1 ответ

numpy.einsum 'ij,kl->ik' как это сделать с помощью numpy.tensordot

У меня две матрицы, 5х4 и 3х2. Я хочу получить матрицу 5х3 от них. >>>theta_ic = np.random.randint(5,size=(5,4)) >>>psi_tr = np.random.randint(5,size=(3,2)) Я могу сделать это >>>np.einsum('ij,kl->ik',theta_ic,psi_tr).s…
06 май '17 в 17:12
0 ответов

4d Array Processing (используя einsum?)

У меня есть матричная проблема, которую, я думаю, можно решить (вычислительно дешево) в одной строке кода, используя numpy (возможно, einsum?), Но я не могу найти решение. Интересно, кто-нибудь может сделать какие-либо предложения, пожалуйста? Пробл…
2 ответа

Numpy - Найти трехмерное расстояние до контрольной точки для всех точек сетки на трехмерной сетке.

Я пробовал np.hypot () и np.linalg.norm (), но у них обоих есть некоторые проблемы (по крайней мере, как я использую thm). Я почти уверен, что np.hypot может рассчитать только 2-е расстояние. Если у меня есть контрольная точка P (1,1,1) и точка сетк…
12 июн '18 в 13:59
1 ответ

Python: сумма всех перестановок внешних произведений числовых массивов массивов

У меня есть массив массивов Ai, и я хочу, чтобы каждый внешний продукт (np.outer(Ai[i],Ai[j])) суммировался с масштабным множителем для получения H. Я могу пройти и сделать их затем тензордот их с матрицей коэффициентов масштабирования. Я думаю, что…
1 ответ

Numpy Einsum ведет себя плохо. На что обращать внимание?

Что обычно терпит неудачу, когда numpy einsum выдает ошибку: Traceback (most recent call last): File "rmse_iter.py", line 30, in <module> rmse_out = np.sqrt(np.einsum('ij,ij->i',diffs,diffs)/3.0) TypeError: invalid data type for einsum Nump…
27 сен '17 в 05:19
1 ответ

NumPy состав Einsums?

Предположим, что у меня есть np.einsum который выполняет некоторые вычисления, а затем закачивает это непосредственно в еще один np.einsum сделать что-то другое. Могу ли я вообще составить эти два einsumв один einsum? Мой конкретный вариант использо…
06 ноя '18 в 18:05
0 ответов

Как einsum взаимодействует с бесшумной трансляцией?

Рассмотрим ndarrays x0=np.ones((3,3)) а также y0, у которого есть y0.shape либо (3,3), либо (1,3). Мне нужна одна команда einsum, которая вычисляет точечные произведения строк этих двух массивов; в том случае, если y0.shape является (1,3)Хочу трансл…
17 янв '18 в 14:08
1 ответ

Точечный продукт B^tDB не возвращает симметричный массив

Я пытаюсь сделать точечное произведение выражения, и оно должно было быть симметричным. Оказывается, это просто не так. B - это четырехмерный массив, который я должен перенести в два последних измерения, чтобы он стал B^ t. D - это двумерный массив.…
17 апр '18 в 23:04
2 ответа

Вычислить сумму попарных сумм двух столбцов массива

Я ищу способ избежать вложенных циклов в следующем фрагменте, где A а также B двумерные массивы, каждый из формы (m, n) с m, n быть произвольными положительными целыми числами: import numpy as np m, n = 5, 2 a = randint(0, 10, (m, n)) b = randint(0,…
21 июн '18 в 09:55
1 ответ

Numpy Einsum вещание

Может кто-нибудь объяснить, как вещание (многоточие) работает в функции numpy.einsum()? Будем весьма благодарны за некоторые примеры, показывающие, как и когда это можно использовать. Я проверил следующую официальную страницу документации, но есть т…
08 июн '16 в 23:15
1 ответ

Тензор потока Эйнсум против Матмула против Тензордо

В тензорном потоке функции tf.einsum, tf.matmul, а также tf.tensordot все могут быть использованы для тех же задач. (Я понимаю, что tf.einsum а также tf.tensordot иметь более общие определения; Я также понимаю, что tf.matmul имеет пакетную функциона…
29 мар '17 в 17:51
1 ответ

Используя более 26 индексов в tf.einsum?

Можно ли использовать более 26 строчных букв в качестве индексов в tf.einsum? Numpy позволяет также использовать заглавные буквы, т.е. np.einsum('zA,AB->zB',M1,M2) в то время как tf.einsum возвращает ошибку. Это было бы очень полезно для контракт…
18 окт '17 в 10:00
2 ответа

Использование numy einsum для вычисления внутреннего произведения векторов-столбцов матрицы

Предположим, у меня есть такая матрица: [[ 1 2 3] [ 10 100 1000]] Я хотел бы вычислить внутренний продукт каждого столбца с самим собой, поэтому результат будет: [1*1 + 10*10 2*2 + 100*100 3*3 + 1000*1000] == [101, 10004, 1000009] Я хотел бы знать, …
1 ответ

Как понять результат этого np.einsum('kij',A)?

Например, A = np.arange(24).reshape((2, 3, 4)) print np.einsum('ijk', A) это все еще A без проблем. Но если я сделаю print np.einsum('kij', A) форма (3, 4, 2), Не должно ли это быть (4, 2, 3)? Результат print np.einsum('cab', A) форма (4, 2, 3) без …
20 окт '15 в 08:42
1 ответ

NumPy: избавиться от цикла за счет трансляции

Я пытаюсь реализовать алгоритм максимизации ожидания для гауссовой модели смеси в Python. У меня есть следующая строка для вычисления гауссовской вероятности p моих данных X с учетом среднего значения mu и ковариационной сигмы гауссовского распредел…
06 фев '17 в 21:51
0 ответов

Чистая NumPy 2D средняя производная от свертки входного изображения

Я имею b 2d m x n полутоновые изображения, которые я сворачиваю с p x q фильтр, а затем делать скупой на. С чисто NumPy я хотел бы вычислить производную входного изображения и фильтра, но у меня возникли проблемы с вычислением производной входного и…
09 авг '17 в 00:51