Ho un ODE: u ( 0 ) = - 1u′=−1000u+sin(t)u′=−1000u+sin(t)u'=-1000u+sin(t) u(0)=−11000001u(0)=−11000001u(0)=-\frac{1}{1000001} So che questo particolare ODE è rigido, analiticamente. So anche che se usiamo un metodo esplicito (forward) di time stepping (Euler, Runge-Kutta, Adams, ecc.), Il metodo dovrebbe restituire errori molto grandi se la fase temporale è troppo grande. Quindi, …
Esiste almeno un'enciclopedia abbastanza completa di regole di quadratura che non sembra essere stata aggiornata da un po 'di tempo e ha un accesso limitato. Questa fonte si riferisce a diverse fonti classiche e moderne e generalmente è ben messa insieme. Tuttavia, si avvicina alla costruzione di regole di quadratura …
A quanto ho capito, il successivo rilassamento si basa sulla scelta di un parametro 0 ≤ ω ≤ 20≤ω≤20\leq\omega\leq2 e sull'utilizzo di una combinazione lineare di un'iterazione (quasi) di Gauss-Seidel e del valore al precedente timestep ... uk + 1= ( ω ) ugSk + 1+ ( 1 - ω …
Sto cercando di risolvere un sistema lineare di equazioni che è scarso, ma privo di qualsiasi tipo di struttura a bande. Ho sentito che esiste un modo per estendere i principi di un risolutore multigrid per schemi impliciti di differenza finita a un problema lineare generale (se non sbaglio, si …
Dopo alcuni passaggi di aggiornamento bayesiano, mi rimane una distribuzione posteriore della forma di una miscela di distribuzioni normali,Vale a dire, il parametro è tratto da una distribuzione il cui PDF è dato come una miscela ponderata di PDF normali e non è una somma di camper normali. Vorrei disegnare …
Quali situazioni di analisi numerica diventano più / meno stabili, hanno una convergenza più veloce / più lenta o sono piuttosto diverse quando si tratta di funzioni di variabili complesse anziché di variabili reali?
Mi è sempre stato detto che il diagramma di Voronoi è il doppio del problema della triangolazione di Delaunay. In che senso possono essere duali l'uno dell'altro? Ho pensato che i doppi problemi (cioè nella programmazione lineare) dovrebbero produrre la stessa risposta. Chiaramente, i due problemi non hanno la stessa …
Che viene calcolato più velocemente, oppure registra a c o b √un'Baba^blogun'clogac\log_a c ? un,becsono reali positivi conb>1.c√Bcb\sqrt[b]{c}un'aaBbbcccb > 1b>1b>1 Quali tipi di algoritmi utilizzerai nel confronto? Quali sono le loro complessità? Ad esempio, quando o c ≈ a bc ≡ aBc≡abc \equiv a^bc ≈ aBc≈abc \approx a^b Questa domanda …
È possibile ottenere la precisione del secondo ordine per gli elementi finiti esaedrici con meno di 8 punti Gauss senza introdurre modalità non fisiche? Un singolo punto Gauss centrale introduce una modalità di taglio non fisico e la disposizione simmetrica standard di 8 punti Gauss è costosa rispetto alle discretizzazioni …
In questo momento ho un codice che utilizza l'algoritmo Crank-Nicholson, ma penso che mi piacerebbe passare a un algoritmo di ordine superiore per il timestrap. So che l'algoritmo Crank-Nicholson è stabile nel dominio in cui voglio lavorare, ma temo che alcuni altri algoritmi potrebbero non esserlo. So calcolare la regione …
Vari importanti set di polinomi (Legendre, Chebyshev, ecc.) Sono ortogonali su un intervallo reale con una certa ponderazione. Esistono famiglie di polinomi che sono ortogonali rispetto ad altre curve sul piano complesso? Ad esempio, vorrei una base per i polinomi di grado n che è ortogonale sul cerchio, per esempio …
Ho usato SVD di Intel MKL ( dgesvdtramite SciPy) e ho notato che i risultati sono significativamente diversi quando cambio precisione tra float32e float64quando la mia matrice è mal condizionata / non di rango massimo. Esiste una guida sulla quantità minima di regolarizzazione da aggiungere per rendere insensibili i risultati …
Sono nuovo nella scienza computazionale e ho già imparato i metodi di base per l'integrazione, l'interpolazione, metodi come RK4, Numerov ecc. Su c ++, ma recentemente il mio professore mi ha chiesto di imparare come usare LAPACK per risolvere i problemi relativi alle matrici. Come ad esempio la ricerca di …
Sto cercando riferimenti credibili che affermino quante risorse i supercomputer spendono per coordinare rispetto a svolgere un vero lavoro relativo alle attività. Le risorse potrebbero essere disponibili potenza di elaborazione, ma anche Watts sembra un'unità valida. Credo che uno dei miei professori o libri di testo abbia detto una volta …
Ho bisogno di calcolare molte inversioni di matrici (per la decomposizione polare di iterazione di Newton), con un numero molto piccolo di casi degeneri ( ).< 0,1 %3 × 33×33\times3< 0,1 %<0.1%<0.1\% L'inverso esplicito (tramite matrici minori divise per determinante) sembra funzionare ed è circa 32 ~ 40 flop fusi …
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.