Запитання з тегом «lapack»

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

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 ++ і маю лінійну систему для вирішення (а не більш типову форму )xA=bxA=bxA=bAx=bAx=bAx=b Матриця щільна і загальної форми, але не більше 1000х1000. Отже, у MATLAB рішення знайдено за функцією або позначенням вперед-косою рисоюAAAmrdivide(b,A)x = b/A; Як я можу це вирішити у своєму …
2
Яка відповідна функція LAPACK за Matlab [Q, R, E] = qr (A)?
Я в даний час намагаюся дешево обчислити оцінку хорошого рангу для матриці . Тому я обчислюю стовпчик, що розгортається QR, використовуючиААA [Q,R,E]=qr(A) в Матлабі. Я оцінюю ранг використовуючиААA tol = size(A,n)*eps*norm(A,'fro'); r = sum(abs(diag(R))>tol) Це добре працює, і сюжет над усіма діагональними записами R виглядає так: Якщо весь алгоритм переносять …
12 matlab  lapack 
2
Які найшвидші доступні реалізації BLAS / LAPACK або інших лінійних підпрограм алгебри в системах GPU?
Наприклад, nVidia має CUBLAS, який обіцяє 7-14-кратне прискорення. Наївно, це ніде не є теоретичним пропускною здатністю жодної з GPU-карт карт nVidia. Які проблеми є прискоренням лінійної алгебри на графічних процесорах, і чи є вже більш швидкі лінійні маршрутизації алгебри?
4
Експоненція матриці косо-ермітової матриці з фортраном 95 та ЛАПАК
Мене просто втягують у фортран 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 як відбивачі для дому. Він масштабує вектор відображенняvvv з 1/v11/v11/v_1, тому першим елементом результату стає 111, тому його не потрібно зберігати. І він зберігає окремоττ\tauвектор, який містить необхідні масштабні коефіцієнти. Отже, матриця відбивача така:H=I−τvvT,H=I−τvvT,H=I-\tau v v^T, де vvvне нормалізується. Тоді як у підручниках рефлекторна матриця є …
1
Алгоритм врівноваження матриці
Я писав набір інструментів системи управління з нуля і суто в Python3 (безсоромний штекер:) harold. У моєму минулому дослідженні я завжди скаржився на вирішувач Ріккаті care.mз причин, які є технічними / неактуальними. Отже, я писав власний набір процедур. Одна річ, яку я не можу знайти, - це отримати високоефективний алгоритм …
Використовуючи наш веб-сайт, ви визнаєте, що прочитали та зрозуміли наші Політику щодо файлів cookie та Політику конфіденційності.
Licensed under cc by-sa 3.0 with attribution required.