Scienza computazionale

1

Quali sono i metodi possibili per risolvere le equazioni di Eulero comprimibili

Vorrei scrivere il mio solutore per le equazioni di Eulero comprimibili e, soprattutto, voglio che funzioni in modo robusto in tutte le situazioni. Vorrei che fosse basato su FE (DG è ok). Quali sono i metodi possibili? Sono consapevole di fare DG del 0 ° ordine (volumi finiti) e che …

13 pde hyperbolic-pde finite-element

4

Conteggio FLOP per le funzioni di libreria

Quando si valuta il numero di FLOP in una semplice funzione, spesso si può semplicemente abbassare l'espressione calcolando gli operatori aritmetici di base. Tuttavia, nel caso di affermazioni matematiche che coinvolgono una divisione uniforme, non si può fare questo e aspettarsi di essere in grado di confrontare con i conteggi …

13 performance complexity floating-point

1

Come costruire volumi finiti ben bilanciati e metodi Galerkin discontinui per PDE iperbolici con termini di origine?

I termini di origine, come quelli dovuti alla batimetria nelle equazioni delle acque poco profonde, devono essere integrati in modo speciale al fine di preservare gli stati fisici stabili. Esiste un modo generale per costruire metodi ben bilanciati o richiede tecniche speciali per ogni equazione?

13 pde hyperbolic-pde

3

Matematicamente, perché funziona la matrici di massa / il vettore di carico?

So che le persone spesso sostituiscono matrici di massa coerenti con matrici diagonali raggruppate. In passato, ho anche implementato un codice in cui il vettore di carico è assemblato in modo aggregato anziché in modo coerente con FEM. Ma non ho mai esaminato il motivo per cui ci è permesso …

13 finite-element

2

Progetta modelli scrivendo software numerico in C ++

Sto cercando risorse su modelli e principi di progettazione per software numerici, potenzialmente ma non necessariamente con un focus sugli approcci orientati agli oggetti ai codici numerici. Sono a conoscenza del testo, Scrivere software scientifico: una guida per un buon stile , ma questo si occupa ancora molto di cose …

13 design-patterns

2

precondizionare un metodo krylov con un altro metodo krylov

In metodi come gmres o bicgstab potrebbe essere attraente usare un altro metodo krylov come precondizionatore. Dopotutto sono facili da implementare in modo privo di matrice e in un ambiente parallelo. Ad esempio, un coul usa alcune (diciamo ~ 5) iterazioni di bigcstab non condizionate come precondizionatore per gmres, o …

13 linear-algebra linear-solver preconditioning krylov-method gmres

3

I file HDF5 sono adatti per il controllo di revisione git?

Non ho familiarità con il formato di file utilizzato in HDF5, ma mi chiedo se i file HDF5 sono adatti per il controllo di revisione con git (o per esempio Mercurial o Subversion)? Immagino che cosa intendo: i file HDF5 sono adatti per il diff linea o si dovrà trattare …

13 data-management data-storage

1

Perché bloccare un punto per rimuovere uno spazio nullo è negativo?

Un'equazione di Poisson con tutte le condizioni al contorno di Neumann ha un unico spazio nullo dimensionale costante. Quando si risolve tramite un metodo Krylov, lo spazio nullo può essere rimosso sottraendo la media della soluzione ogni iterazione o bloccando il valore di un singolo vertice. L'inserimento di un singolo …

13 krylov-method

1

Perché l'integrazione con leapfrog è semplice e RK4 no, se quest'ultima è più accurata?

In un sistema in cui l'energia dovrebbe teoricamente essere conservata, la simulazione più accurata conserverà l'energia (oltre a fornire posizioni, velocità e così via precisi). RK4 è più preciso di leapfrog, ma leapfrog conserva energia e RK4 no. Perchè è questo?

13 computational-physics time-integration molecular-dynamics runge-kutta

1

Soluzioni forti contro deboli di PDE

La forma forte di un PDE richiede che la soluzione sconosciuta appartenga a . Ma la forma debole richiede solo che la soluzione sconosciuta appartenga a H 1 .H2H2H^2H1H1H^1 Come conciliare questo?

13 pde weak-solution

3

Adattamento di superfici implicite a set di punti orientati

Ho una domanda riguardante l'adattamento quadric a un insieme di punti e corrispondenti normali (o equivalentemente, tangenti). L'adattamento di superfici quadric a punti dati è ben esplorato. Alcuni lavori sono i seguenti: Montaggio diretto vincolato al tipo di superfici quadrate , James Andrews, Carlo H. Sequin Computer-Aided Design & Applications, …

13 computational-geometry regression geometry curve-fitting

1

L'equazione di avanzamento con velocità variabile può essere conservativa?

Sto cercando di capire un po 'meglio l'equazione di avanzamento con coefficiente di velocità variabile. In particolare non capisco come l'equazione possa essere conservativa. L' equazione di avanzamento , ∂u∂t+ ∂∂X( v u ) = 0∂u∂t+∂∂x(vu)=0 \frac{\partial u}{\partial t} + \frac{\partial}{\partial x}(\boldsymbol{v}u) = 0 Interpretiamo u ( x , t …

13 advection conservation

2

Qual è il modo più efficiente per scrivere loop 'for' in Matlab?

Ho letto che se, ad esempio, ho un doppio forciclo che corre sopra gli indici di una matrice, quindi mettere la colonna che esegue l'indice nel ciclo esterno è più efficiente. Per esempio: a=zeros(1000); for j=1:1000 for i=1:1000 a(i,j)=1; end end Qual è il modo più efficiente per codificarlo se …

12 matlab efficiency

2

Per dati disturbati o strutturati in modo preciso, esistono quadrature migliori della regola del punto medio?

Solo le prime due sezioni di questa lunga domanda sono essenziali. Gli altri sono solo a scopo illustrativo. sfondo Le quadrature avanzate come Newton – Cotes composito di grado superiore, Gauß – Legendre e Romberg sembrano essere destinate principalmente ai casi in cui è possibile campionare con precisione la funzione …

12 reference-request quadrature

2

osservazioni puntuali vs. continue nel problema inverso della PDE

Lavoro su un problema inverso per il mio dottorato. ricerca, che per semplicità diremo che sta determinando inββ\beta L ( β) u ≡ - ∇ ⋅ ( k0eβ∇ u ) = fL(β)u≡−∇⋅(k0eβ∇u)=fL(\beta)u \equiv -\nabla\cdot(k_0e^\beta\nabla u) = f da alcune osservazioni ; k 0 è una costante e f è noto. …

12 inverse-problem