Вопросы с тегом «lapack»

LAPACK (Linear Algebra PACKage) - это широко используемая библиотека подпрограмм для задач числовой линейной алгебры, включая решения линейных наборов уравнений, линейный метод наименьших квадратов, задачи на собственные значения и разложение по сингулярным числам. Подпрограммы LAPACK могут использоваться с fortran, C и родственными языками, а также множеством других языков.

3
Какова взаимосвязь BLAS, LAPACK и других библиотек линейной алгебры?
Я искал библиотеки линейной алгебры C ++ для проекта, над которым я работал. То, что я до сих пор не понимаю, это соединение BLAS и LAPACK с другими библиотеками линейной алгебры. Просматривая эту статью о библиотеках линейной алгебры, я обнаружил, что интересно: некоторые библиотеки независимы от BLAS и LAPACK некоторые …
3
Почему не мое матрично-векторное масштабирование?
Извините за длинный пост, но я хотел включить все, что я думал, было актуально с первого взгляда. Что я хочу Я реализую параллельную версию методов подпространств Крылова для плотных матриц. В основном GMRES, QMR и CG. Я понял (после профилирования), что моя рутина DGEMV была жалкой. Поэтому я решил сосредоточиться …
4
Быстрое определение, является ли плотная матрица низкого ранга
В программном проекте, над которым я работаю, некоторые вычисления намного проще для плотных матриц низкого ранга. В некоторых проблемных случаях используются плотные матрицы низкого ранга, но они даны мне полностью, а не как факторы, поэтому мне придется проверять ранг и фактор матрицы, если я хочу воспользоваться преимуществами структуры низкого ранга …
3
Есть ли какая-то польза от компиляции LAPACK из исходного кода по сравнению с установкой готового пакета из Ubuntu?
Я знаю, что ATLAS может оптимизировать себя для машины, на которой он компилируется, и, таким образом, максимальная выгода достигается компиляцией из исходного кода. Есть ли польза от компиляции LAPACK из исходного кода? Было бы намного проще просто установить готовый пакет.
3
Понимание того, как Numpy делает SVD
Я использовал разные методы для вычисления как ранга матрицы, так и решения матричной системы уравнений. Я наткнулся на функцию linalg.svd. Сравнивая это с моими собственными усилиями по решению системы с устранением по Гауссу, она кажется более быстрой и точной. Я пытаюсь понять, как это возможно. Насколько я знаю, функция linalg.svd …
2
решить для используя LAPACK и BLAS
Я портирую существующий код из MATLAB на C ++ и имею линейную систему для решения (вместо более типичной формы )A x = bх = ЬxA=bxA=bA x = bAx=bAx=b Матрица плотная и общего вида, но не больше 1000x1000. Таким образом, в MATLAB решение находится по функции или через косую чертуAAAmrdivide(b,A)x = …
2
Какова соответствующая функция LAPACK за Matlab [Q, R, E] = qr (A)?
Я в настоящее время пытаюсь дешево вычислить оценку хорошего ранга для матрицы . Поэтому я вычисляю разложение QR-кода с помощью Columnt, используяAAA [Q,R,E]=qr(A) в Matlab. Я оцениваю ранг используяAAA tol = size(A,n)*eps*norm(A,'fro'); r = sum(abs(diag(R))>tol) Это прекрасно работает, и график по всем диагональным элементам R выглядит следующим образом: Если перенести …
12 matlab  lapack 
2
Каковы наиболее быстрые доступные реализации BLAS / LAPACK или других подпрограмм линейной алгебры в системах с графическим процессором?
Например, у nVidia есть CUBLAS, который обещает ускорение в 7-14 раз. Наивно, это далеко не теоретическая пропускная способность любой из видеокарт nVidia. Каковы проблемы в ускорении линейной алгебры на графических процессорах, и есть ли уже доступные более быстрые линейные алгебраические маршруты?
4
Матричная экспонента косоэрмитовой матрицы с фортраном 95 и LAPACK
Я просто заправлен в Фортран 95 для моделирования квантовой механики. Честно говоря, я был избалован Октавой, поэтому я воспринял возведение в матрицу как должное. Учитывая (малую, ) косоэрмитову матрицу размера n × n , каков наиболее эффективный способ использования LAPACK для решения этой проблемы? Я не использую оболочку LAPACK95, просто …
11 fortran  lapack 
1
Требуются эталонные задачи для алгоритмов переупорядочения собственных значений
Каждая действительная матрица можно свести к вещественной форме Шура Т = U T A U с помощью ортогонального преобразования similiary U . Здесь матрица T имеет квазитреугольную форму с 1 на 1 или 2 на 2 блока на главной диагонали. Каждый 1 по 1 блок соответствует реальному собственному значению А …
1
Использует ли когда-нибудь PETSc библиотеки LAPACK для математики с разреженной матрицей?
Влияет ли компиляция PETSc с внешней библиотекой BLAS / LAPACK существенным образом на производительность разреженных матриц, или он использует эти библиотеки только для математики с плотной матрицей?
3
Матричная экспонента вещественной асимметричной матрицы с Fortran 95 и LAPACK
Недавно я задал вопрос в том же духе для косоэрмитовых матриц. Вдохновленный успехом этого вопроса, и после того, как пару часов ударился головой о стену, я смотрю на экспоненциальную матрицу вещественных асимметричных матриц. Путь к поиску собственных значений и собственных векторов кажется довольно запутанным, и я боюсь, что я заблудился. …
3
Как начать использовать LAPACK в C ++?
Я новичок в вычислительной науке, и я уже изучил основные методы для интеграции, интерполяции, методы, такие как RK4, Numerov и т. Д. На C ++, но недавно мой профессор попросил меня научиться использовать LAPACK для решения проблем, связанных с матрицами. Как, например, поиск собственных значений комплексной матрицы. Я никогда не …
10 lapack 
3
По какой причине LAPACK использует
QR-код LAPACK хранит Q в качестве отражателей для домохозяев. Он масштабирует вектор отражения с , поэтому первый элемент результата становится , поэтому его не нужно сохранять. И он хранит отдельный вектор , который содержит необходимые масштабные коэффициенты. Таким образом, матрица отражателя выглядит так:vvv1 /v11/v11/v_1111ττ\tauЧАС= Я- τvvT,ЧАСзнак равноя-τvvT,H=I-\tau v v^T, где …
1
Алгоритм Матрицы Балансировки
Я пишу набор инструментов системы управления с нуля и чисто на Python3 (бесстыдный плагин:) harold. Из моих прошлых исследований я всегда жаловался на решатель Riccati care.mпо техническим / неуместным причинам. Следовательно, я писал свой собственный набор процедур. Одна вещь, которую я не могу найти, это получить высокопроизводительный алгоритм балансировки, по …
Используя наш сайт, вы подтверждаете, что прочитали и поняли нашу Политику в отношении файлов cookie и Политику конфиденциальности.
Licensed under cc by-sa 3.0 with attribution required.