Scienza computazionale

Domande e risposte per gli scienziati che utilizzano i computer per risolvere problemi scientifici

4
Scalabilità di Fast Fourier Transform (FFT)
Per utilizzare la Fast Fourier Transform (FFT) su dati campionati uniformemente, ad es. In connessione con i solutori PDE, è noto che FFT è un algoritmo ). Quanto bene la scala FFT elaborata in parallelo per n → ∞ (cioè molto grande)?O (nlog( n )O(nlog⁡(n)\mathcal{O}(n\log(n)n → ∞n→∞n\to\infty
3
metodo del volume finito: mesh non strutturata vs adattamento dell'octree + taglio delle celle
Sto lavorando con la libreria OpenFOAM C ++ Computational Continuum Mechanics (può gestire l'interazione fluido-solido, flussi MHD ...) che utilizza mesh arbitrarie non strutturate. Ciò è stato guidato dall'idea di utilizzare il vantaggio della generazione rapida (di solito automatica) di mesh non strutturate per simulare problemi in geometrie complesse. Tuttavia, …
3
In quali casi applicativi gli schemi di precondizionamento additivo sono superiori a quelli moltiplicativi?
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 …
1
Quali discretizzazioni spaziali funzionano per un flusso incomprimibile con maglie di contorno anisotrope?
I flussi di numeri di High Reynolds producono strati limite molto sottili. Se la risoluzione a parete viene utilizzata nella simulazione a grande Eddy, le proporzioni potrebbero essere dell'ordine di 10610610^6 . Molti metodi diventano instabili in questo regime perché la costante inf-sup degrada come la radice quadrata delle proporzioni …
5
Massimizzazione globale di costose funzioni oggettive
Sono interessato a massimizzare a livello globale una funzione di molti ( ) parametri reali (risultato di una simulazione complessa). Tuttavia, la funzione in questione è relativamente costosa da valutare e richiede circa 2 giorni per ciascun set di parametri. Sto confrontando diverse opzioni e mi chiedevo se qualcuno avesse …
4
Metodo di interpolazione efficiente per le reti non strutturate?
Vorrei conoscere un buon metodo per interpolare i dati tra due griglie non strutturate, in cui una griglia è una versione più grossolana dell'altra. L'efficienza è molto importante per me poiché sto risolvendo un problema transitorio di PDE in cui ho bisogno di trasferire dati tra le reti in ogni …
2
Qual è la funzione LAPACK corrispondente dietro Matlab [Q, R, E] = qr (A)?
Io attualmente cercando di calcolare una stima a buon mercato buon grado per una matrice . Pertanto computo un columnt che ruota la decomposizione QR usandoUNAA [Q,R,E]=qr(A) in Matlab. Stimo il rango di usandoUNAA tol = size(A,n)*eps*norm(A,'fro'); r = sum(abs(diag(R))>tol) Funziona bene e una trama su tutte le voci diagonali …
12 matlab  lapack 
1
problema SVD ponderato?
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 \ …
Utilizzando il nostro sito, riconosci di aver letto e compreso le nostre Informativa sui cookie e Informativa sulla privacy.
Licensed under cc by-sa 3.0 with attribution required.