Als «lapack» getaggte Fragen

Ich habe mir C ++ - Bibliotheken für lineare Algebra für ein Projekt angesehen, an dem ich gearbeitet habe. Was ich immer noch nicht verstehe, ist die Verbindung von BLAS und LAPACK mit anderen Bibliotheken für lineare Algebra. Als ich diesen Artikel über lineare Algebra-Bibliotheken durchgesehen...

Entschuldigung für den langen Beitrag, aber ich wollte alles, was ich für relevant hielt, gleich mit einbeziehen. Was ich möchte Ich implementiere eine parallele Version der Krylov-Subraummethoden für dichte Matrizen. Hauptsächlich GMRES, QMR und CG. Mir wurde (nach der Profilerstellung) klar,...

Ich habe verschiedene Methoden verwendet, um sowohl den Rang einer Matrix als auch die Lösung eines Matrixgleichungssystems zu berechnen. Ich bin auf die Funktion linalg.svd gestoßen. Vergleicht man dies mit meiner eigenen Anstrengung, das System mit der Gaußschen Eliminierung zu lösen, scheint es...

In einem Softwareprojekt, an dem ich arbeite, sind bestimmte Berechnungen für dichte Matrizen mit niedrigem Rang erheblich einfacher. Bei einigen Problemfällen handelt es sich um dichte Matrizen mit niedrigem Rang, die mir jedoch vollständig und nicht als Faktoren zur Verfügung stehen. Daher muss...

Ich weiß, dass ATLAS in der Lage ist, sich für die Maschine zu optimieren, auf der es kompiliert wird, und somit maximale Vorteile durch das Kompilieren aus dem Quellcode erzielt. Hat das Kompilieren von LAPACK aus dem Quellcode einen Vorteil? Es wäre viel einfacher, nur das vorgefertigte Paket zu...

Ich versuche derzeit, eine gute Rangschätzung für eine Matrix billig zu berechnen . Daher berechne ich eine columnt schwenkbare QR-Zerlegung mitAAA [Q,R,E]=qr(A) in Matlab. Ich schätze den Rang von mitAAA tol = size(A,n)*eps*norm(A,'fro'); r = sum(abs(diag(R))>tol) Dies funktioniert gut und ein...

Ich portiere einen vorhandenen Code von MATLAB nach C ++ und habe ein lineares System zum Lösen von (anstelle der typischeren Form A x = b )xA=bxA=bxA=bAx=bAx=bAx=b Die Matrix ist dicht und von allgemeiner Form, aber nicht größer als 1000 × 1000. In MATLAB wird die Lösung also durch die Funktion...

Ich stecke mich gerade für einige quantenmechanische Simulationen in fortran 95. Ehrlich gesagt wurde ich von Octave verwöhnt, daher habe ich die Matrixexponentiation als selbstverständlich angesehen. Was ist bei einer (kleinen, ) Schräg- Hermitianischen Matrix der Größe der effizienteste Weg, um...

nVidia hat zum Beispiel CUBLAS, das eine 7-14-fache Beschleunigung verspricht. Naiv ist dies bei weitem nicht der theoretische Durchsatz einer der GPU-Karten von nVidia. Was sind die Herausforderungen bei der Beschleunigung der linearen Algebra auf GPUs, und gibt es bereits schnellere lineare...

Ich habe kürzlich eine ähnliche Frage für schief-hermitische Matrizen gestellt. Inspiriert vom Erfolg dieser Frage und nachdem ich ein paar Stunden lang meinen Kopf gegen eine Wand geschlagen habe, betrachte ich die Exponentialmatrix realer asymmetrischer Matrizen. Der Weg zum Finden der Eigenwerte...

Jede reelle Matrix kann unter Verwendung einer orthogonalen Ähnlichkeitstransformation U auf die reelle Schurform T = U T A U reduziert werden . Hier ist die Matrix T quasi dreieckig mit 1 mal 1 oder 2 mal 2 Blöcken auf der Hauptdiagonale. Jeweils 1 x 1 Block entspricht einem realen Eigenwert von A...

Ich bin neu in der Computerwissenschaft und habe bereits grundlegende Methoden für Integration, Interpolation, Methoden wie RK4, Numerov usw. in C ++ gelernt. Kürzlich hat mich mein Professor gebeten, zu lernen, wie man LAPACK zur Lösung von Problemen im Zusammenhang mit Matrizen verwendet. Wie zum...

Beeinträchtigt das Kompilieren von PETSc mit einer externen BLAS / LAPACK-Bibliothek die Leistung auf dünn besetzten Matrizen erheblich oder werden diese Bibliotheken nur für die Mathematik mit dichter Matrix

Ich habe eine Steuerungssystem-Toolbox von Grund auf neu und rein in Python3 geschrieben (schamloser Plug :) harold. Aus meiner bisherigen Forschung habe ich mich immer care.maus technischen / irrelevanten Gründen über den Riccati-Löser beschwert . Daher habe ich meine eigenen Routinen geschrieben....

Ich beabsichtige, Ax = b zu lösen, wobei A eine komplexe, spärliche, unsymmetrische und stark schlecht konditionierte (Bedingungsnummer ~ 1E + 20) quadratische oder rechteckige Matrix ist. Ich konnte das System mit ZGELSS in LAPACK genau lösen. Mit zunehmenden Freiheitsgraden in meinem System...

In meinem Projekt muss ich bei jedem Zeitschritt ein paar tridiagonale Matrizen lösen, daher ist es wichtig, einen guten Löser für diese zu haben. Ich habe meine eigene Implementierung gemacht, nur die klassische Art, wie es auf Wikipedia beschrieben wird. Ich habe dann versucht, stattdessen Lapack...

Die QR-Routine von LAPACK speichert Q als Reflektoren für Haushalte. Es skaliert den Reflexionsvektor vvv mit 1 / v11/.v11/v_1 , sodass das erste Element des Ergebnisses 111 wird und nicht gespeichert werden muss. Und es speichert einen separaten ττ\tau Vektor, der die erforderlichen...