Scienza computazionale

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

3
È normale non avere controlli di convergenza in Multigrid?
Ho appena letto il capitolo 3 in "A Multigrid Tutorial" di Briggs / Henson / McCormick, link . Il testo riguarda i cicli multigrid come V-cycle, mu-cycle, FMG. Cosa ha attirato la mia attenzione: nella maggior parte delle procedure iterative si verifica se si è convertito alla tolleranza / accuratezza …
12 multigrid 
3
Alternative a Comsol Multiphysics
Questa potrebbe essere una domanda più adatta al lato delle raccomandazioni software di SE, tuttavia credo che le persone che frequentano questa parte di SE abbiano maggiori probabilità di poter rispondere a questa domanda. Sto cercando un'alternativa gratuita (non solo come in libertà) a Comsol Multiphysics . Ecco la parte …
2
Per quanto riguarda la differenziazione automatica, la trasformazione del codice sorgente (STC) è più efficiente del sovraccarico dell'operatore (OO)?
Stiamo lavorando su un modello bayesiano per un processo spazio-temporale e stiamo utilizzando un campionatore No-U-Turn (NUTS) che richiede un modello per la probabilità di log e il suo gradiente rispetto ai parametri del modello. Più sinteticamente, abbiamo una funzione di probabilità logaritmica piuttosto complicata , che coinvolge distribuzioni statistiche, …
2
integrazione numerica in molte variabili
Let e essere una funzione in queste variabili.f( → x ):[0,1]n→CX⃗ = ( x1, x2, ... , xn) ∈ [ 0 , 1 ]nx→=(x1,x2,…,xn)∈[0,1]n\vec{x} = (x_1, x_2, \dots, x_n) \in [0,1]^nf( x⃗ ) : [ 0 , 1 ]n→ Cf(x→):[0,1]n→Cf(\vec{x}): [0,1]^n \to \mathbb{C} Esiste uno schema ricorsivo per questo integrale …
1
La pressione come moltiplicatore di Lagrange
Nelle incomprimibili equazioni di Navier-Stokes, il termine pressione viene spesso indicato come un moltiplicatore di Lagrange che impone il condizione di incomprimibilità.ρ(ut+(u⋅∇)u)∇⋅u=−∇p+μΔu+f=0ρ(ut+(u⋅∇)u)=−∇p+μΔu+f∇⋅u=0\begin{align*} \rho\left(\mathbf{u}_t + (\mathbf{u} \cdot \nabla)\mathbf{u}\right) &= - \nabla p + \mu\Delta\mathbf{u} + \mathbf{f}\\ \nabla\cdot\mathbf{u} &= 0 \end{align*} In che senso è vero? Esiste una formulazione delle equazioni incomprimibili …
2
Qual è il modo più veloce per calcolare tutti gli autovalori di una matrice di adiacenza molto grande e sparsa in Python?
Sto cercando di capire se esiste un modo più veloce per calcolare tutti gli autovalori e gli autovettori di una matrice di adiacenza molto grande e sparsa rispetto all'uso di scipy.sparse.linalg.eigsh Per quanto ne so, questo metodo utilizza solo la scarsità e attributi di simmetria della matrice. Una matrice di …
3
Scelta della dimensione del passo usando ODE in matlab
Ciao e grazie per aver dato il tempo di guardare la mia domanda. Questa è una versione aggiornata della mia domanda che ho pubblicato in precedenza in physics.stackexchange.com Attualmente sto studiando uno spinatore di eccitoni 2D Bose-Einstein Condensate e sono curioso di sapere lo stato fondamentale di questo sistema. Il …
1
Selezione del metodo per la quadratura numerica
Esistono diverse famiglie di metodi per la quadratura numerica. Se ho una classe specifica di integrandi come posso selezionare il metodo ideale? Quali sono le domande pertinenti da porre sia sull'integrando (ad es. È regolare? Ha delle singolarità?) Sia sul problema computazionale (ad es. Tolleranza d'errore, budget computazionale)? In che …
12 quadrature 
2
risolvi
Sto trasferendo un codice esistente da MATLAB a C ++ e ho un sistema lineare per risolvere (piuttosto che la forma più tipica A x = b )xA=bxA=bxA=bAx=bAx=bAx=b La matrice è densa e di forma generale, ma non è più grande di 1000x1000. Quindi in MATLAB, la soluzione è trovata …
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.