Domande sugli aspetti algoritmici / computazionali dell'algebra lineare, compresa la soluzione di sistemi lineari, problemi dei minimi quadrati, problemi di autovalori e altre questioni simili.
Attualmente sto lavorando con "A Multigrid Tutorial" di Briggs et al, capitolo 8. La costruzione dell'operatore di interpolazione è data come: Quindi la costruzione dell'operatore di restrizione e dell'operatore di griglia fine viene data come: Supponiamo di avere tre punti della griglia x0, x1, x2 con quello centrale x1 va …
Sto cercando di capire se esiste un modo più veloce per calcolare tutti gli autovalori e gli autovettori di una matrice di adiacenza molto grande e sparsa rispetto all'uso di scipy.sparse.linalg.eigsh Per quanto ne so, questo metodo utilizza solo la scarsità e attributi di simmetria della matrice. Una matrice di …
Devo risolvere lo stesso sistema lineare rado (da 300x300 a 1000x1000) con molti lati di destra (da 300 a 1000). Oltre a questo primo problema, vorrei anche risolvere sistemi diversi, ma con gli stessi elementi diversi da zero (solo valori diversi), ovvero molti sistemi sparsi con un modello di sparsità …
Sto usando MATLAB per risolvere un problema che comporta la risoluzione di in ogni momento, dove b cambia nel tempo. In questo momento, sto realizzando questo usando MATLAB :A x = bAx=b\mathbf{A} \mathbf{x}=\mathbf{b}Bb\mathbf{b}mldivide x = A\b Ho la flessibilità di effettuare tutti i precomputamenti necessari, quindi mi chiedo se esiste …
Sto cercando una libreria che esegua operazioni con matrici su grandi matrici sparse senza sacrificare la stabilità numerica. Le matrici saranno 1000+ per 1000+ e i valori della matrice saranno compresi tra 0 e 1000. Eseguirò l'algoritmo di calcolo dell'indice, quindi genererò i vettori (sparsi) di riga della matrice in …
Mi chiedo: qual è l'algoritmo migliore per risolvere dudt= A ududt=Au\begin{equation} \frac{du}{dt} = Au \end{equation} DoveUNAAè unamatricen × nn×nn\times nreale. A non è esplicitamente dipendente dal tempo, di solito scarso ma non necessariamente legato. I suoi autovalori hanno parti reali non positive. A è anche diagonale ma potrebbe essere troppo …
L'articolo "Modelli di espressione rivisitati: un'analisi delle prestazioni delle attuali metodologie" nel SIAM Journal of Scientific Computing fa riferimento alla libreria di algebra lineare "Blaze". Non ne ho mai sentito parlare prima e non riesco a trovare riferimenti online. (Le ovvie ricerche su Google stanno restituendo il documento sopra.) Cos'è …
Sto lavorando su una libreria di matrici solo intestazione per fornire un ragionevole grado di capacità di algebra lineare nel più semplice pacchetto possibile, e sto cercando di esaminare quale sia l'attuale stato dell'arte: calcolare l'SVD di un matrice complessa. Sto facendo una decomposizione in due fasi, bidiagonalizzazione seguita da …
In entrambi i metodi di decomposizione del dominio (DD) e multigrid (MG), si può comporre l'applicazione degli aggiornamenti di blocco o correzioni grossolane come additivi o moltiplicativi . Per i risolutori puntuali, questa è la differenza tra le iterazioni Jacobi e Gauss-Seidel. Il valore più uniforme per agisce come viene …
Sto risolvendo un problema fisico usando uno schema numerico implicito. Questo mi porta a risolvere un'equazione lineare con matrice tridiagonale. Ho codificato questo algoritmo da Wikipedia. Mi chiedo se esiste una libreria efficiente che consenta di risolvere questo tipo di equazione in modo ottimizzato. Una nota importante è che la …
Dati due matrici e , mi piacerebbe trovare vettori ed , in modo tale che, In forma di matrice, sto cercando di ridurre al minimo la norma di Frobenius di A - \ mbox {diag} (x) \ cdot B \ cdot \ mbox {diag} (y) = A - B \ …
Supponiamo che A ∈ Rn × nUN∈Rn×nA\in\mathbb{R}^{n\times n} sia una matrice definita simmetrica e positiva. UNUNA è abbastanza grande da essere costoso per risolvere direttamente A x = bUNX=BAx=b . Esiste un algoritmo iterativo per trovare l'autovalore più piccolo di UNUNA che non comporta l'inversione di UNUNA in ogni iterazione? …
Nel libro di Nocedal & Wright sull'ottimizzazione numerica, c'è un'affermazione nella sezione 2.2 (pagina 27), "In generale, è più facile preservare l'invarianza della scala per gli algoritmi di ricerca di linea piuttosto che per gli algoritmi della regione di fiducia". In quella stessa sezione, parlano di avere nuove variabili che …
Il complesso prodotto interno ha due diverse definizioni decise dalle convenzioni: o . In BLAS, ho trovato le routine cdotu, zdotu e cdotc, zdotc. Le prime due routine in realtà calcolano (un prodotto interno falso!) E le ultime due routine coniugano il primo vettore nel prodotto interno. Inoltre, da una …
Ho un elenco di matrici simmetriche di cui ho bisogno per verificare la semi-definitività positiva (ovvero i loro autovalori non sono negativi).LL{\cal L} Il commento sopra implica che si potrebbe farlo calcolando i rispettivi autovalori e controllando se non sono negativi (forse dovendo occuparsi degli errori di arrotondamento). Il calcolo …
We use cookies and other tracking technologies to improve your browsing experience on our website,
to show you personalized content and targeted ads, to analyze our website traffic,
and to understand where our visitors are coming from.
By continuing, you consent to our use of cookies and other tracking technologies and
affirm you're at least 16 years old or have consent from a parent or guardian.