0
ответов

Как присвоить значение квадратным подматрицам в большой матрице без петель в numpy

Как я могу векторизовать этот цикл, который заполняет две квадратные подматрицы большей матрицы (также сохраняет большую матрицу симметричной), используя массивы numpy: для x в диапазоне (n): assert m[x].shape == (n,) ...
вопрос задан: 12 June 2012 00:19
0
ответов

Напишите цикл for/while с «if/ еще» более элегантным способом?

Я написал этот код: A — матрица размера nXm [nA, mA] = size(A); токВектор(нА,мА) = 0; для i = 1: нА для j = 1: мА, если A (i,j) ~= 0 currentVector(i,j) = ceil(log10(abs(...
вопрос задан: 2 June 2012 05:25
0
ответов

Генерируют ли JIT-компиляторы JVM код, использующий векторизованные инструкции с плавающей запятой?

Предположим, что узким местом моей Java-программы на самом деле являются несколько узких циклов для вычисления множества векторных скалярных произведений. Да, я профилировал, да, это узкое место, да, это существенно, да, именно так...
вопрос задан: 28 May 2012 12:48
0
ответов

ggplot2: как прозрачно затенить чередующиеся дни на графике

Что я делаю не так, пожалуйста? Я пытаюсь закрасить чередующиеся 24-часовые ежедневные прямоугольники прозрачным серым цветом. Но рисуется только последний прямоугольник из цикла for (?!?). Если я делаю что-то вручную...
вопрос задан: 12 May 2012 13:08
0
ответов

Векторизованные триггерные функции в C?

Я хочу вычислить высоко распараллеленные триггерные функции (в блоке, например, 1024), и я хотел бы воспользоваться хотя бы некоторыми из параллелизм, присущий современной архитектуре. Когда я компилирую ...
вопрос задан: 19 March 2012 17:07
0
ответов

Быстрый код в R

К вашему сведению: я значительно отредактировал это со времени своего первого издания. Это моделирование было сокращено с 14 часов до 14 минут. Я новичок в программировании, но я сделал симуляцию, которая пытается...
вопрос задан: 7 March 2012 06:19
0
ответов

Минимизация накладных расходов на векторную индексацию

У меня есть векторизованная функция, которая вычисляет расстояние до большого набора точек. Для повышения производительности я ограничиваю количество баллов, выбирая только необходимые. Поэтому вместо ...
вопрос задан: 24 February 2012 18:38
0
ответов

Загрузка данных для векторных расширений GCC

Векторные расширения GCC предлагают хороший, достаточно переносимый способ доступа к некоторым инструкциям SIMD на различных аппаратных архитектурах без использования встроенных аппаратных функций (или автоматических -...
вопрос задан: 17 February 2012 00:20
0
ответов

R: Apply FUN to kxk subsections of array

Язык - R. У меня есть матрица nxm, и я хотел бы разбить ее на 3x3 секции и вычислить среднее (или любую функцию) в каждой из них. (Если остался кусочек, который не 3x3, то используйте просто ...
вопрос задан: 7 February 2012 06:26
0
ответов

Ускорение условного заполнения огромных разреженных матриц

Мне было интересно, есть ли способ ускорить (может быть, с помощью векторизации?) Условное заполнение огромных разреженных матриц (например, ~ 1e10 x 1e10). Вот пример кода, в котором у меня есть вложенный цикл, и ...
вопрос задан: 25 January 2012 21:18
0
ответов

SSE медленнее, чем FPU?

У меня есть большой фрагмент кода, часть тела которого содержит этот фрагмент кода: result = (nx * m_Lx + ny * m_Ly + m_Lz) / sqrt (nx * nx + ny * ny + 1); который я векторизовал следующим образом (все ...
вопрос задан: 13 January 2012 08:46
0
ответов

'for' loop vs vectorization in MATLAB

Я программировал кое-что в MATLAB и, как рекомендуется, всегда стараюсь использовать векторизацию. Но в итоге программа получалась довольно медленной. Тогда я обнаружил, что в одном месте код ...
вопрос задан: 29 December 2011 00:58
0
ответов

векторизуется, если в matlab

у меня есть логический массив, назовите его флагом. У меня есть два числовых массива ifTrue, ifFalse. Все эти массивы имеют одинаковый размер. Для целей этого вопроса предположим, что каждый элемент в этих массивах уникален. Я ...
вопрос задан: 14 December 2011 14:55
0
ответов

r: операция цикла for с вложенными индексами выполняется очень медленно

У меня есть операция, которую я хотел бы выполнить для каждой строки фрейма данных, изменяя один столбец. Я работаю с apply/ddply/sqldf, но я использую циклы, когда они имеют смысл, и я думаю, что это один из таких случаев. ...
вопрос задан: 30 November 2011 18:48
0
ответов

Написание функций, которые принимают 1-D и 2-D массивы numpy?

Насколько я понимаю, 1-D массивы в numpy можно интерпретировать либо как вектор, ориентированный на столбцы, либо как вектор, ориентированный на строки. Например, одномерный массив с формой (8,) можно рассматривать как двумерный массив ...
вопрос задан: 27 November 2011 16:28
0
ответов

Добавление векторизованных массивов C #

Есть ли способ «векторизовать» добавление элементов по массивам в режиме SIMD? Например, я хотел бы включить: var a = new [] {1, 2, 3, 4}; var b = new [] {1, 2, 3, 4}; {{1} } var c = new [] {...
вопрос задан: 21 November 2011 16:11
0
ответов

Почему GCC не выполняет автоматическую векторизацию этого цикла ?

Я пытаюсь оптимизировать цикл, на который приходится большая часть времени вычислений моей программы. Но когда я включаю автоматическую векторизацию с -O3 -ffast-math -ftree-vectorizer-verbose = 6, GCC выводит, что он ...
вопрос задан: 15 November 2011 22:54
0
ответов

Сравните два вектора неравной длины, чтобы получить логический массив

Мне нужно векторизовать следующий код: a = [1 2 3 2 3 1]; b = [1 2 3]; для i = 1: длина (a) для j = 1: длина (b) r (i, j) = (a (i) == b (j)); end end Выход r должен быть логическим ...
вопрос задан: 11 November 2011 22:36
0
ответов

Массив структуры Matlab: быстрое присвоение

Есть ли способ "вектор" назначить массив структуры. В настоящее время я могу edge (1000000) = struct ('weight', 1.0); // Это действительно не присваивает значение, я проверил на 2009A. для i = 1: 1000000; Edge (i) ...
вопрос задан: 28 October 2011 15:47
0
ответов

Быстрое векторизованное преобразование из RGB в BGRA

В ответ на некоторые предыдущие вопросы о преобразовании RGB в RGBA и ARGB в BGR , Я хотел бы ускорить преобразование RGB в BGRA с помощью SSE. Предположим, 32-битная машина, и вы хотели бы использовать ...
вопрос задан: 27 September 2011 09:39
0
ответов

Применяют LM для подмножества кадра данных, определенные третьим столбцом кадра

, у меня есть кадр данных, содержащий вектор значения x, Вектор y значений y и вектора идентификаторов: x <- rep (0: 3, 3) Y <- Runif (12) ID <- C (REP («A», 4), REP («B», 4), REP («C», 4)) DF <- Data ....
вопрос задан: 14 September 2011 19:11
0
ответов

Векторизация (SIMD) операции с деревом

Каковы некоторые общие советы / указатели по операциям векторизации дерева? С точки зрения компоновки памяти, алгоритмов и т. Д. Некоторые вещи, специфичные для предметной области: у каждого родительского узла будет довольно много (20 - 200) дочерних уз
вопрос задан: 26 August 2011 23:44
0
ответов

Как избежать огромной дополнительной памяти потребление при использовании numpy vectorize?

Этот код ниже лучше всего иллюстрирует мою проблему: вывод на консоль (NB требуется ~ 8 минут для запуска даже первого теста) показывает, что выделение массива 512x512x512x16 бит занимает не более ...
вопрос задан: 16 August 2011 12:49
0
ответов

SSE-векторизация математической функции 'pow' gcc

Я пытался векторизовать цикл, содержащий использование функции 'pow' в математической библиотеке. Я знаю, что компилятор Intel поддерживает использование 'pow' для инструкций sse, но я не могу заставить его ...
вопрос задан: 3 August 2011 01:59
0
ответов

MATLAB: Как векторно умножить два массива матриц?

У меня есть два трехмерных массива, первые два измерения которых представляют матрицы, а последнее имеет значение через пространство параметров, в качестве простого примера возьмем A = repmat ([1,2; 3,4], [1 1 4]); (но ...
вопрос задан: 12 July 2011 11:51
0
ответов

Выравнивание массивов кучи в C и C ++ для упрощения векторизации компилятора (GCC)

В настоящее время я готовлю класс шаблона контейнера оболочки для std :: vector, который автоматически создает мульти- пирамида разрешения элементов в его std :: vector. Ключевой проблемой сейчас является то, что я хочу ...
вопрос задан: 6 July 2011 14:46
0
ответов

эффективный способ преобразования индексов разброса в индексы сбора?

Я пытаюсь написать сжатие потока (взять массив и избавиться от пустых элементов) с внутренними функциями SIMD. Каждая итерация цикла обрабатывает 8 элементов за раз (ширина SIMD). С помощью встроенных функций SSE ...
вопрос задан: 7 June 2011 19:36
0
ответов

Как обобщить внешние до n измерений?

Стандартное выражение R external (X, Y, f) вычисляет матрицу, чья (i, j) -я запись имеет значение f (X [i], Y [j]). Я хотел бы реализовать функцию multi.outer, n-мерное обобщение ...
вопрос задан: 31 May 2011 20:31
0
ответов

Почему векторизация выгодна для программ Matlab? То же самое для NumPy и Boost (uBLAS)?

Использование векторизации для замены циклов for может значительно увеличить скорость программ Matlab. Это потому, что векторизованные коды выполняются параллельно? Полезна ли векторизация для программы ...
вопрос задан: 16 May 2011 09:15
0
ответов

Векторизация по строкам массива

У меня есть массив X, и я хочу применить функцию f ко всем строкам X: # глупый пример X = numpy.array ([[1, 2, 3, 4, 5], [6, 7, 8, 9, 0]], 'i') def f (row): вернуть сумму (row) y =. ..
вопрос задан: 9 May 2011 19:18