[Lapack] - Come vengono calcolati gli autovalori
Hola gente.
Per la tesi sto utizzando la libreria Lapack per risolvere problemi agli autovalori. Siccome sono interessato a confrontare vari metodi mi piacerebbe sapere quali algoritmi (QR, Lanczos, altri) vengono utilizzati dalle seguenti subroutine di Lapack:
- DSYGV
- DSBGV
- DSTEV
Purtroppo sulla documentazioni non c'è scritto niente (o almeno, io non ho trovato niente).
Ho provato anche a scrivere su un forum specifico sul Lapack, ma non ho ricevuto risposte (probabilmente poco attivo il forum).
Mi chiedevo se qualcuno di voi conosce che algoritmi vengono utilizzati o dove posso reperire informazioni a riguardo.
Grazie!
Per la tesi sto utizzando la libreria Lapack per risolvere problemi agli autovalori. Siccome sono interessato a confrontare vari metodi mi piacerebbe sapere quali algoritmi (QR, Lanczos, altri) vengono utilizzati dalle seguenti subroutine di Lapack:
- DSYGV
- DSBGV
- DSTEV
Purtroppo sulla documentazioni non c'è scritto niente (o almeno, io non ho trovato niente).
Ho provato anche a scrivere su un forum specifico sul Lapack, ma non ho ricevuto risposte (probabilmente poco attivo il forum).
Mi chiedevo se qualcuno di voi conosce che algoritmi vengono utilizzati o dove posso reperire informazioni a riguardo.
Grazie!
Risposte
Eh, se non c'è scritto nella documentazione è dura!
L'unica che mi viene in mente è di scorrere il codice e leggere i commenti, per cercare di capire che cosa stia facendo e quindi dedurre l'algoritmo.
In generale, comunque, sui libri di testo trovi quali sono gli algoritmi migliori per i vari casi [matrice qualunque, matrice simmetrica, matrice hermitiana...] e questo può aiutarti ad intuire cosa faccia LAPACK.
L'unica che mi viene in mente è di scorrere il codice e leggere i commenti, per cercare di capire che cosa stia facendo e quindi dedurre l'algoritmo.
In generale, comunque, sui libri di testo trovi quali sono gli algoritmi migliori per i vari casi [matrice qualunque, matrice simmetrica, matrice hermitiana...] e questo può aiutarti ad intuire cosa faccia LAPACK.
E' quello che mi era venuto in mente di fare. Poi causa febbre non ho più riportato quando ho scoperto. Ho notato che in fondo alla documentazione disponibile online ci sono i codici ma non sono commentati. Però ho dedotto che DSYGV
e DSBGV utilizzano la riflessione di Householder per ricondursi a una forma tridiagonale e poi su questa viene usata la fattorizzazione QR. DSTEV usa anche lui la fattorizzazione QR. Però sono informazioni da prendere con le pinze....la certezza al 100% non ce l'ho.
e DSBGV utilizzano la riflessione di Householder per ricondursi a una forma tridiagonale e poi su questa viene usata la fattorizzazione QR. DSTEV usa anche lui la fattorizzazione QR. Però sono informazioni da prendere con le pinze....la certezza al 100% non ce l'ho.
Tieni anche conto che il nome di ciascuna funzione è costruito in modo da indicare che cosa sia, per quali matrici funzioni....
Magari questo può esserti d'aiuto.
Magari questo può esserti d'aiuto.