Costruzione e analisi di algoritmi per calcolare soluzioni discrete approssimative di problemi continui. Un esempio canonico è l'approssimazione di derivati tramite quozienti di differenza.
Quando si discretizza FEM e si risolve un problema di reazione-diffusione, ad es. con 0 < ε ≪ 1 (perturbazione singolare), la soluzione del problema discreto in genere esibirà strati oscillatori vicino al confine. Con Ω = ( 0 , 1 ) , ε = 10 - 5 e lineari …
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 …
Come posso sostituire il metodo Euler con il Runge-Kutta del 4 ° ordine per determinare il movimento di caduta libera in magnitudo gravitazionale non costante (es. Caduta libera da 10 000 km dal suolo)? Finora ho scritto una semplice integrazione con il metodo Euler: while() { v += getMagnitude(x) * …
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 …
A partire dall'equazione di avanzamento in forma di conservazione. ut=(a(x)u)xut=(a(x)u)x u_t = (a(x)u)_x dove è una velocità che dipende dallo spazio e è una concentrazione di una specie che viene conservata.a(x)a(x)a(x)uuu Discretizzare il flusso (dove il flusso , è definito sui bordi delle celle tra i punti di mesh) dà, …
Ho un sistema di equazioni non lineari che voglio risolvere numericamente:nnn f = ( f 1 , … , f n )f( x ) = af(X)=un'\mathbf{f}(\mathbf{x})=\mathbf{a} f= ( f1, ... , fn)x =( x1, ... , xn)f=(f1,...,fn)X=(X1,...,Xn)\mathbf{f}=(f_1,\dots,f_n)\quad\mathbf{x}=(x_1,\dots,x_n) Questo sistema ha una serie di caratteristiche che lo rendono particolarmente difficile da …
Qualche tempo fa avevo scritto un codice che tentava di calcolare il senza usare le funzioni di libreria. Ieri stavo rivedendo il vecchio codice e ho cercato di renderlo il più veloce possibile (e corretto). Ecco il mio tentativo finora:l o g( x )log(x)log(x) const double ee = exp(1); double …
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 …
Sto leggendo un documento [1] in cui risolvono la seguente equazione non lineare usando metodi a differenza finita. Analizzano anche la stabilità degli schemi usando l'analisi di stabilità di Von Neumann. Tuttavia, come gli autori realizzano, questo è applicabile solo ai PDE lineari. Così gli autori lavorano intorno a questo …
Non riesco a immaginare di essere il primo a pensare al seguente problema, quindi sarò soddisfatto di un riferimento (ma una risposta completa e dettagliata è sempre apprezzata): Supponi di avere un simmetrico positivo definito . n è considerato molto grande, quindi è impossibile tenere Σ in memoria. È possibile, …
Sono interessato a calcolare la soluzione di un sistema di lage di ODE usando un metodo krylov come in [1]. Tale metodo prevede funzioni correlate all'esponenziale (le cosiddette funzioni ). Consiste essenzialmente nel calcolare l'azione della funzione matrice costruendo un sottospazio di Krylov usando l'iterazione di Arnoldi e proiettando la …
Sto eseguendo una ricerca di linea come parte di un algoritmo BFGS quasi-Newton. In un passaggio della ricerca della linea utilizzo un'interpolazione cubica per avvicinarmi al minimizer locale. Sia la funzione di interesse. Voglio trovare un tale che .x ∗ f ′ ( x ∗ ) ≈ 0f:R→R,f∈C1f:R→R,f∈C1f : R …
Quali sono le principali differenze tra FEM e XFEM? Quando dovremmo (non) usare XFEM al posto di FEM e viceversa? In altre parole, quando incontro un nuovo problema, come posso sapere quale di questi?
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.