Als «lapack» getaggte Fragen

LAPACK (Linear Algebra PACKage) ist eine häufig verwendete Bibliothek von Unterprogrammen für numerische lineare Algebra-Aufgaben, einschließlich Lösungen für lineare Gleichungssätze, lineare kleinste Quadrate, Eigenwertprobleme und Singularwertzerlegung. LAPACK-Routinen können mit fortran, C und Verwandten sowie einer Vielzahl anderer Sprachen verwendet werden.

3
Warum wird meine Matrix-Vektor-Multiplikationsskalierung nicht durchgeführt?
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, dass meine DGEMV-Routine erbärmlich war. Also beschloss ich, mich …
3
Verstehen, wie Numpy SVD macht
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 sowohl schneller als auch präziser zu sein. …
2
Löse
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 …
2
Was ist die entsprechende LAPACK-Funktion hinter Matlab [Q, R, E] = qr (A)?
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 Plot über alle diagonalen Einträge von R sieht aus …
12 matlab  lapack 
4
Exponentialmatrix einer Skew-Hermitian-Matrix mit Fortran 95 und LAPACK
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 dieses Problem mit LAPACK zu lösen? Ich verwende nicht …
11 fortran  lapack 
3
Wie starte ich LAPACK in c ++?
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 Beispiel das Finden von Eigenwerten …
10 lapack 
3
Was ist der Grund, warum 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 Skalierungsfaktoren enthält. Eine Reflektormatrix ist also wie …
1
Matrix-Balancing-Algorithmus
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. Eine Sache, die ich nicht umgehen kann, ist, …
Durch die Nutzung unserer Website bestätigen Sie, dass Sie unsere Cookie-Richtlinie und Datenschutzrichtlinie gelesen und verstanden haben.
Licensed under cc by-sa 3.0 with attribution required.