Ho iniziato a usare Python come linguaggio di programmazione per svolgere tutti i miei compiti in CFD. Ho poca esperienza nella programmazione. Sono un esperto in ingegneria meccanica e sto perseguendo un'istruzione superiore in ingegneria aerospaziale. a volte l'aspetto informatico del CFD diventa più noioso che manipolare le equazioni o …
Sono rimasto molto colpito dalle prestazioni seriali dei precondizionatori ILU multilivello basati su inverse , in particolare per l' eterogenea Helmholtz , ma sono sorpreso di non riuscire a trovare implementazioni open source. In particolare, ILUPACK rende i binari liberamente disponibili per gli accademici, ma non sembra che rilascino il …
Sto risolvendo un sistema di due PDE accoppiati in due dimensioni spaziali e nel tempo computazionalmente. Poiché le valutazioni delle funzioni sono costose, vorrei utilizzare un metodo a più fasi (inizializzato utilizzando Runge-Kutta 4-5). Il metodo Adams-Bashforth che utilizza cinque precedenti valutazioni di funzioni ha un errore globale di (questo …
Ho un set di dati in esecuzione in milioni di punti dati in 3D. Per il calcolo che sto facendo, ho bisogno di calcolare il vicino (ricerca intervallo) per ogni punto dati in un raggio, provare ad adattare una funzione, calcolare l'errore per l'adattamento, ripetere questo per il punto dati …
Quando si applica la formula classica per l'angolo tra due vettori: α=arccosv1⋅v2∥v1∥∥v2∥α=arccosv1⋅v2‖v1‖‖v2‖\alpha = \arccos \frac{\mathbf{v_1} \cdot \mathbf{v_2}}{\|\mathbf{v_1}\| \|\mathbf{v_2}\|} si rileva che, per angoli molto piccoli / acuti, si verifica una perdita di precisione e il risultato non è accurato. Come spiegato in questa risposta Stack Overflow , una soluzione è …
Molti algoritmi utilizzati nell'informatica scientifica hanno una struttura intrinseca diversa rispetto agli algoritmi comunemente considerati nelle forme meno ingegneristiche di ingegneria del software. In particolare, i singoli algoritmi matematici tendono ad essere altamente complessi, spesso coinvolgono centinaia o migliaia di righe di codice, tuttavia non coinvolgono alcuno stato (cioè non …
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 …
Esiste un algoritmo che enumera i grafici che corrispondono ad alcune tessellazioni di punti Delaunay in 3D? In tal caso, esiste un'efficace parametrizzazione delle geometrie che corrispondono a qualsiasi "grafico di Delaunay"? Sto cercando di elencare sistematicamente tutte le geometrie stabili delle molecole di una composizione specificata senza alcuna conoscenza …
Sto provando a invertire numericamente la seguente trasformazione integrale: F(y) = ∫∞0yexp[ - 12( y2+ x2) ] I0( x y) f( x )d xF(y)=∫0∞yexp[−12(y2+x2)]I0(xy)f(x)dxF(y) = \int_{0}^{\infty} y\exp{\left[-\frac{1}{2}(y^2 + x^2)\right]} I_0\left(xy\right)f(x)\;\mathrm{d}x Quindi per una data ho bisogno di approssimare f ( x ) dove:F( y)F(y)F(y)f( x )f(x)f(x) e F ( y …
Così gli stati decomposizione di Cholesky teorema che che qualsiasi reale simmetrica definita positiva matrice ha una decomposizione di Cholesky M = L L ⊤ dove L è una matrice triangolare inferiore.MMMM=LL⊤M=LL⊤M= LL^\topLLL Dato , sappiamo già ci sono algoritmi veloci per calcolare il suo fattore di Cholesky L .MMMLLL …
Ciao comunità scicomp, Ho lavorato nell'area degli algoritmi grafici usando framework come NetworkX (Python), JUNG e YFiles (Java). Ora sto entrando nell'area dell'informatica parallela e ad alte prestazioni. Per un nuovo progetto, sto cercando una libreria di grafici C ++ con le seguenti funzionalità: ha un'interfaccia intuitiva che consente lo …
Questo è un semplice problema di enumerazione della simmetria. Fornisco qui il quadro completo, ma non è necessaria alcuna conoscenza della chimica quantistica. L'integrale di due particelle è: ⟨ i j | k l ⟩ = ∫ ψ * i ( x ) ψ * j ( x ' ) …
Sto lavorando a un problema che può essere parallelizzato utilizzando una singola operazione mpi_allgather o un'operazione mpi_scatter e un'operazione mpi_gather. Queste operazioni vengono chiamate in un ciclo while, quindi possono essere chiamate più volte. Nell'implementazione con uno schema MPI_allgather, sto raccogliendo un vettore distribuito su tutti i processi per la …
Sto facendo una diagonalizzazione Lanczos di una matrice sparsa di grandi dimensioni (~ 2 milioni di elementi). Quasi tutti i passaggi dell'algoritmo Lanzcos vengono eseguiti in parallelo sulla GPU, ad eccezione della diagonalizzazione della matrice di Lanczos per verificare la convergenza. Per questo, ho usato l'algoritmo TQLI da Ricette numeriche. …
Domanda Come classificheresti una nuvola di punti rispetto a una maglia non strutturata di celle esaedriche? Ogni cella ha un centro e un'etichetta unica per rappresentarla. Fondamentalmente ci sono due punti di nuvola (nuvola di punti originale e una nuvola di punti dei centri di celle), ma le informazioni sulla …
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.