0
ответов

псевдо инверсия разреженной матрицы в python

Я работаю с данными нейровизуализации, и из-за большого количества данных я хотел бы использовать разреженные матрицы для своего кода (scipy.sparse.lil_matrix или csr_matrix). В частности, мне понадобится ...
вопрос задан: 3 July 2013 15:42
0
ответов

Левое матричное деление и Numpy Solve

Я пытаюсь преобразовать код, содержащий оператор \, из Matlab (Octave) в Python. Пример кода B = [2; 4] b = [4; 4] B \ b Это работает и дает 1,2 в качестве ответа. Используя эту веб-страницу http: // ...
вопрос задан: 3 July 2013 15:19
0
ответов

Инверсия большой матрицы

Я пытаюсь взять обратную большую матрицу, общий размер которой составляет 1000 x 1000, но иногда превышает 100 000 x 100 000 (, которая в настоящее время не работает из-за времени и памяти ). Я знаю, что нормальные...
вопрос задан: 3 July 2013 15:06
0
ответов

Есть ли способ эффективно инвертировать массив матриц с помощью numpy?

Обычно я бы инвертировал массив матриц 3x3 в цикле for, как в примере ниже. К сожалению, циклы for работают медленно. Есть ли более быстрый и эффективный способ сделать это? импортировать numpy как np А = np....
вопрос задан: 3 July 2013 15:03
0
ответов

Наиболее эффективное обращение матрицы в MATLAB

При вычислении обратной для некоторой квадратной матрицы A в MATLAB с использованием Ai = inv (A) % должно быть таким же, как: Ai = A ^ -1 MATLAB обычно сообщает мне, что это не самый эффективный способ инвертирования. ...
вопрос задан: 3 July 2013 15:00
0
ответов

scipy.linalg.eig возвращает комплексные собственные значения для ковариационной матрицы?

Собственные значения ковариационной матрицы должны быть действительными и неотрицательными, потому что ковариационные матрицы симметричны и полуположительно определены. Однако взгляните на следующий эксперимент с ...
вопрос задан: 12 May 2013 21:49
0
ответов

Как сравнить векторы примерно в Eigen?

Есть ли в Eigen функция для сравнения векторов (матриц) с использованием относительной и абсолютной толерантности, иначе говоря numpy.allclose? Стандартный isApprox не работает, если один из векторов очень близок к нулю.
вопрос задан: 24 February 2013 11:39
0
ответов

numpy линейная алгебра произвольной точности

У меня есть массив numpy 2d [средний / большой размер - скажем, 500x500]. Я хочу найти собственные значения его поэлементной экспоненты. Проблема в том, что некоторые значения довольно отрицательны (-800, -1000, ...
вопрос задан: 19 February 2013 15:50
0
ответов

Объединение библиотеки линейной алгебры с Boost :: Units

Я много занимаюсь научным программированием и получил очень хороший опыт работы с Boost.Units, которые обеспечивают размерный анализ величин во время компиляции (т. Е. Теги количества with units ...
вопрос задан: 25 September 2012 16:26
0
ответов

Определите, существует ли A + B = C в массиве из n целых чисел

Это проблема, которую мой друг получил в качестве домашнего задания (в классе алгоритмов и структур данных). Он спросил меня об этом. Однако я могу ' Я решил ее и обдумывал ее некоторое время ...
вопрос задан: 15 September 2012 22:55
0
ответов

Библиотека матриц линейной алгебры Haskell для произвольных типов элементов [закрыто]

Я ищу библиотеку линейной алгебры Haskell, которая имеет следующие особенности: Матричное умножение Матричное сложение Матричное транспонирование Расчет ранга Матричная инверсия является плюсом и имеет ...
вопрос задан: 8 September 2012 15:35
0
ответов

Алгоритм поиска основы набора битовых строк?

Это для утилиты diff, которую я пишу на C++. У меня есть список из n наборов символов -{"a", "abc", "abcde", "bcd", "de"} (, взятых из алфавита из k=5 разных букв ). Мне нужен способ наблюдать, что...
вопрос задан: 30 August 2012 20:26
0
ответов

Надлежащее умножение матриц для поворота/перемещения

Чтобы вращать/перемещать объект (вращение только вокруг оси z -и перемещение только в плоскости xy )не только относительно глобального центра (центра устройства ), но и относительно других произвольных точек, я создал ан...
вопрос задан: 17 August 2012 18:35
0
ответов

Является ли руководство Google по Android OpenGL неправильным обучением линейной алгебре?

После того, как я помог другому пользователю ответить на вопрос, касающийся учебника «Ответ на сенсорные события Android», я скачал исходный код и был весьма озадачен увиденным. Туториал вроде бы не...
вопрос задан: 12 August 2012 21:37
0
ответов

Какую линейную алгебру использовать для OpenGL в Haskell?

Я пытаюсь сделать некоторое программирование OpenGL в haskell. Но меня смущает текущее состояние библиотек. OpenGL использует пакет Tensor, который определяет только несколько векторных типов (, но не делает этого...
вопрос задан: 12 August 2012 20:19
0
ответов

Расстояние Махаланобиса, инвертирующее ковариационную матрицу

Я пишу функцию для определения расстояния Махаланобиса между двумя векторами. Насколько я понимаю, это достигается с помощью уравнения а' *C^ -1 *b, где a и b — векторы, а C — ковариационная матрица....
вопрос задан: 3 August 2012 16:16
0
ответов

Решение малых симметричных положительно определенных Ax = b только на графическом процессоре

Я пытаюсь оптимизировать приложение для 3D-моделирования в реальном времени. Вычислительная часть приложения почти полностью выполняется на GPU в CUDA. Приложение требует решения небольшого (6x6 )...
вопрос задан: 29 July 2012 17:06
0
ответов

Гиперплоскость, определяемая n точками

У меня есть следующая задача :Имея n точек в пространстве, я ищу гиперплоскость, проходящую через них. Самый простой пример такой задачи — две точки (x _1=0,x _2=0 )и (1, -1 )и я бы...
вопрос задан: 21 July 2012 16:38
0
ответов

Расчет угла наклона с учетом движения корпуса

У меня есть сеанс записи для моего приложения. Когда пользователь начал сеанс записи, я начинаю собирать данные из объекта CMMotionManager устройства и сохранять их в CoreData для обработки и представления позже. В...
вопрос задан: 20 July 2012 09:20
0
ответов

Вычисление собственных значений/собственных векторов сотен малых матриц с использованием CUDA

У меня есть вопрос о собственной -декомпозиции сотен небольших матриц с использованием CUDA. Мне нужно вычислить собственные значения и собственные векторы сотен (, например. 500 )малых (64 -по -64 )действительных симметричных...
вопрос задан: 9 July 2012 18:53
0
ответов

Библиотека линейной алгебры для разреженных матриц [закрыто]

Я заинтересован в переносе моей реализации Matlab на C++ для повышения скорости. Пробовал броненосец. Он идеально подходит для переноса кода из Matlab, потому что имена/синтаксис библиотечных функций Armadillo...
вопрос задан: 5 July 2012 09:04
0
ответов

Как я могу векторизовать этот тройной -цикл над двумерными массивами в numpy?

Могу ли я исключить все циклы Python в этом вычислении :result[i,j,k] = (x[i] *y[j] *z[k] ).sum ()где x[i], y[j], z[k] — векторы длины N, а x,y,z имеют первые измерения длины A,B,C s.t....
вопрос задан: 29 June 2012 16:23
0
ответов

Ищу одну библиотеку линейной алгебры для встраиваемых систем (без malloc и бесплатно) [закрыто]

Использую для работы с микроконтроллерами. ОСРВ, которые я использую в своих приложениях, не имеют функций free и malloc (и других вызовов, таких как assert), иногда они могут быть доступны, но я предпочитаю иметь...
вопрос задан: 10 June 2012 10:03
0
ответов

Вычисление вектора Фидлера в Python

Как найти вектор Филдера лапласиана (L) в Python? Я могу получить собственные значения и собственные векторы, используя: eigenvalues, eigenvectors = linalg.eig(L) Я предполагаю, что python не возвращает ...
вопрос задан: 7 June 2012 02:45
0
ответов

Складывание листа бумаги (Computer Vision)

Я работаю над задачей, пытаясь преобразовать 2D-эскиз со складками в полное 3D-представление. Красными линиями будут долины, а синими — горы/вершины. Я хотел бы рассчитать преобразованный/...
вопрос задан: 6 June 2012 20:53
0
ответов

Умножение матриц Ладермана 3x3 всего с 23 умножениями, стоит ли оно того?

Возьмем произведение двух матриц 3x3 A*B=C. Наивно для этого требуется 27 умножений по стандартному алгоритму. Если бы кто-то был умным, вы могли бы сделать это, используя только 23 умножения, результат найден ...
вопрос задан: 31 May 2012 22:20
0
ответов

Как преобразовать матрицу в произведение ядерных матриц?

Постановка задачи. Допустим, у нас есть набор ядерных квадратных матриц = {K1, K2, .., Kn}. Для заданной матрицы A найдите произведение, включающее наименьшее количество матричных умножений, которое дает: A = Ki * ...
вопрос задан: 29 May 2012 07:45
0
ответов

Есть ли способ вычислить CGAffineTransform, необходимый для преобразования вида из кадра A в кадр B?

Скажем, у меня есть два CGRect, CGRect A и CGRect B. Кадр моего UIView такой же, как CGRect B, но я хочу создать анимацию, показывающую переход UIView от кадра A к B. Я пытаюсь сделать это ...
вопрос задан: 23 May 2012 13:06
0
ответов

Библиотека линейной алгебры конечных полей для Haskell

Я ищу библиотеку линейной алгебры конечных полей для Haskell. Что-то вроде FFLAS-FFPACK для Haskell было бы здорово :-). Конечно, я проверил hmatrix, кажется, есть некоторая поддержка...
вопрос задан: 19 May 2012 14:07
0
ответов

оптимизация quadprog

Вот интересная головоломка. Ниже приведен фрагмент R, который идентифицирует точку касания квадратичной функции по отношению к линии, проведенной из точки (0, rf) на оси Y. Для знакомых...
вопрос задан: 10 May 2012 04:18