Итак, у меня есть 2 матрицы: каждая 100x100.
Я хочу рассчитать третью матрицу, такую что: M3 [i] = sqrt (M1 [i] ^ 2 + M2 [i] ^ 2).
Я, очевидно, могу делать ForLoops, но я уверен, что есть что-то быстрее.
Я покопался в фреймворке ускорений и заблудился в мире Linpack
Любая помощь, чтобы получить меня на правильном пути ...
Спасибо
Фреймворк Accelerate — хорошая идея.
Вы можете использовать такую функцию, как vDSP_vsq для работы с одним столбцом матрицы за раз, помещая результат в соответствующий столбец M3. Возможно, вам придется создать две результирующие матрицы (одна будет M1^2, другая M2^2), а затем сложить их, чтобы сформировать правильную результирующую матрицу M3, добавив столбцы с вызовом vDSP_vadd, который, опять же, , сможет работать с одной колонкой за раз.
Пример кода (показывающий, как добавить два вектора и т. д.) здесь, на странице разработчиков Apple.
Я думаю, что самый быстрый способ - использовать
vDSP_vpythg Вектор Пифагора; одинарная точность.
Вычитает вектор C из A и возводит в квадрат разности, вычитает вектор D из B и возводит в квадрат разности, складывает два набора квадратов разностей, а затем записывает квадратные корни сумм в вектор E.
Очевидно, передайте ZERO вектор для C и D .
