Voglio risolvere un problema non lineare con vincoli di uguaglianza non lineari e sto usando un Lagrangiano aumentato con un termine di regolarizzazione della penalità che, come noto, rovina il numero di condizione dei miei sistemi linearizzati (ad ogni iterazione di Newton intendo) . Maggiore è il termine di penalità, …
I metodi Newton-Krylov (JFNK) Jacobian-Free e i metodi di Krylov in generale possono essere molto utili perché non richiedono l'archiviazione o la costruzione esplicite di una matrice, ma solo i risultati di prodotti a matrice vettoriale. Se in realtà formi il sistema sparso, ci sono molti precondizionatori là fuori per …
C'è qualche speranza nel risolvere efficacemente il seguente sistema lineare con un metodo iterativo? A∈Rn×n,x∈Rn,b∈Rn, with n>106A∈Rn×n,x∈Rn,b∈Rn, with n>106A \in \mathbb{R}^{n \times n}, x \in \mathbb{R}^n, b \in \mathbb{R}^n \text{, with } n > 10^6 Ax=bAx=bAx=b con , dove Δ è una matrice molto sparsa con poche diagonali, derivante dalla …
Dato il sistema dove A ∈ R n × n , leggo che, nel caso in cui l'iterazione Jacobi sia utilizzata come solutore, il metodo non converge se b ha un componente diverso da zero nello spazio nullo di A . Quindi, come si può affermare formalmente che, a condizione …
Attualmente sto lavorando alla risoluzione di sistemi simmetrici molto grandi (ma non definiti positivi), generati da alcuni algoritmi. Queste matrici hanno una bella scarsità di blocchi che può essere utilizzata per la risoluzione parallela. Ma non posso decidere se dovrei usare un approccio diretto (come Multi-frontale) o iterativo (GMRES o …
in matlab, sia linsolve che mldivide sono usati per risolvere un sistema di equazioni lineari, in tutti i casi determinati, sovradeterminati e indeterminati. Leggendo i loro documenti, mi chiedevo quali differenze ci fossero tra loro? Stanno usando quasi gli stessi algoritmi di fattorizzazione e triangolarizzazione della matrice nei tre casi? …
Voglio sapere quali dei classici solutori lineari (ad es. Gauss-Seidel, Jacobi, SOR) sono sicuri di convergere per il problema dove è positivo semi definito e ovviamenteAx=bAx=bAx=bAAAb∈im(A)b∈im(A)b \in im(A) (L'avviso è semi definito e non definito)AAA
Ho già sviluppato una soluzione funzionante del metodo agli elementi finiti per risolvere i problemi di trasferimento di calore utilizzando GPU e OpenCL utilizzando il metodo del gradiente coniugato. Il principale svantaggio di questo metodo è l'elevata domanda di memoria. Inoltre, nel caso delle schede grafiche, la memoria è spesso …
Sono molto interessato all'ottimizzazione della risoluzione del sistema lineare per piccole matrici (10x10), a volte chiamate minuscole matrici. C'è una soluzione pronta per questo? La matrice può essere considerata non singolare. Questo solutore deve essere eseguito oltre 1 000 000 volte in microsecondi su una CPU Intel. Sto parlando del …
Sto cercando una soluzione esplicita veloce (oserei dire ottimale?) Il problema reale lineare 3x3, , A ∈ R 3 × 3 , b ∈ R 3 . Ax=bAx=b\mathbf{A}\mathbf{x} = \mathbf{b}A∈R3×3,b∈R3A∈R3×3,b∈R3\mathbf{A} \in \mathbf{R}^{3 \times 3}, \mathbf{b} \in \mathbf{R}^{3} La matrice è generale, ma vicina alla matrice identità con un numero di …
Esiste una semplice regola empirica per dire se vale la pena fare SOR invece di Gauss-Seidel? (e modo possibile come stimare il parametro di realxation )ωω\omega Intendo solo guardando la matrice o la conoscenza di un problema particolare che la matrice rappresenta? Stavo leggendo la risposta a queste domande: ci …
Intendo risolvere Ax = b dove A è una matrice quadrata o rettangolare complessa, sparsa, asimmetrica e altamente mal condizionata (numero condizione ~ 1E + 20). Sono stato in grado di risolvere accuratamente il sistema con ZGELSS in LAPACK. Ma man mano che i gradi di libertà nel mio sistema …
Devo risolvere Ax = b, ma mi rendo conto che anche se è scarso, la memorizzazione dei coefficienti di matrice del mio problema richiederà troppa memoria. Quindi ora sto prendendo in considerazione l'utilizzo di un metodo senza matrice, perché gli stessi coefficienti appaiono molto tempo nella matrice, quindi potrei usare …
La maggior parte della mia programmazione è costituita da codici di ricerca una tantum in C per uso personale. Non ho mai distribuito alcun codice a parte i collaboratori stretti. Ho sviluppato un algoritmo che sto pubblicando in una rivista scientifica. Voglio fornire il codice sorgente e forse il codice …
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.