Tag: numpy

Вставьте 0s в массив 2d

У меня есть массив x : x = [0, -1, 0, 3] и я хочу y : y = [[0, -2, 0, 2], [0, -1, 0, 3], [0, 0, 0, 4]] где первая строка – x-1 , вторая строка – x , а третья – x+1 . Все четные индексы столбцов равны нулю. Я делаю: […]

Самый быстрый способ вычислить k наибольших собственных значений и соответствующих собственных векторов с numpy

Я имею большую NxN плотную симметричную матрицу и хочу, чтобы собственные векторы соответствовали k наибольшим собственным значениям. Каков наилучший способ их поиска (желательно, используя numpy, но, возможно, вообще используя blas / atlas / lapack, если это единственный способ пойти)? В общем случае N намного больше k (скажем N> 5000, k <10). Кажется, что у Numpy […]

python, numpy boolean array: отрицание в where statement

с: import numpy as np array = get_array() Мне нужно сделать следующее: for i in range(len(array)): if random.uniform(0, 1) < prob: array[i] = not array[i] с массивом – numpy.array. Хотел бы я сделать что-то похожее: array = np.where(np.random.rand(len(array)) < prob, not array, array) но я получаю следующий результат (ссылаясь на «не массив»): Значение истинности массива […]

Numpy octuple точность float и 128 бит int. Почему и как?

Это в основном вопрос из любопытства. Я заметил, что набор тестов numpy содержит тесты для 128-битных целых чисел , а модуль numerictypes относится к int128 , float256 ( octuple precision ?) И другим типам, которые, похоже, не сопоставляются с numpy dtypes на моей машине. Моя машина 64 бит, но я могу использовать четырехместные 128-битные поплавки […]

Гаус-Лежандр с интервалами -x -> бесконечность: адаптивный алгоритм для эффективного преобразования весов и узлов

Хорошо, я знаю, что это было задано раньше с ограниченным примером для масштабирования [-1, 1] интервалов [a, b] Различные интервалы для квадратурности Гаусса-Лежандра в numpy, но никто не опубликовал, как обобщать это для [-a, Infinity] (как это делается ниже, но не очень быстро). Также это показывает, как вызвать сложную функцию (в количественном определении цены) в […]

Построение с прозрачным маркером, но непрозрачным краем

Я пытаюсь сделать заговор в matplotlib с прозрачными маркерами с фиксированным краем цвета. Тем не менее, я не могу представить маркер с прозрачным заполнением. У меня есть минимальный рабочий пример: import numpy as np import matplotlib.pyplot as plt x = np.arange(10) y1 = 2*x + 1 y2 = 3*x – 5 plt.plot(x,y1, 'o-', lw=6, ms=14) […]

Более эффективная гистограмма диаграммы matplotlib – как рассчитать нижние значения

Мне нужна помощь в создании набора диаграмм столбчатых диаграмм в python с matlibplot. Мой базовый код ниже, но мои проблемы состоят в том, как эффективно генерировать значение для дна для любого элемента за пределами второго. Я могу получить примерный граф для правильного стека (всегда a, b, c, d снизу вверх) import numpy as np import […]

Исключение с плавающей точкой с помощью Numpy и PyTables

У меня есть довольно большой файл HDF5, созданный PyTables, который я пытаюсь прочитать в кластере. У меня проблема с NumPy, когда я читаю в отдельном куске. Перейдем к примеру: Общая форма массива внутри файла HDF5, In [13]: data.shape Out[13]: (21933063, 800, 3) Каждая запись в этом массиве – np.float64 . У меня каждый узел считывает […]

Ограниченная линейная регрессия в Python

У меня есть классическая проблема линейной регрессии формы: y = X b где y – вектор отклика X – матрица входных переменных, а b – вектор подходящих параметров, которые я ищу. Python предоставляет b = numpy.linalg.lstsq( X , y ) для решения задач этой формы. Однако, когда я использую это, я имею тенденцию получать либо […]

Как заставить нулевой перехват в линейной регрессии?

Я немного извиняюсь, так что извиняюсь, если на этот вопрос уже был дан ответ, я посмотрел и не смог найти конкретно то, что искал. У меня есть некоторые более или менее линейные данные формы x = [0.1, 0.2, 0.4, 0.6, 0.8, 1.0, 2.0, 4.0, 6.0, 8.0, 10.0, 20.0, 40.0, 60.0, 80.0] y = [0.50505332505407008, 1.1207373784533172, […]

 
Interesting Posts for Van-Lav
Python - лучший язык программирования в мире.