Scienza computazionale

2

Quali nuove strutture di dati sono utilizzate nella FEM adattiva?

Un sacco di librerie FEM adattivi utilizzare strutture dati in rete più avanzate per gestire l'aggiunta / rimozione di nodi, bordi, triangoli, tetraedri, ecc, ad esempio, i p4est utilizza la libreria octree strutture dati per adaptive mesh refinement; non troverai spesso ocre utilizzate per i calcoli su una mesh statica. …

9 linear-algebra finite-element sparse-matrix

5

La Galerkin discontinua è davvero più parallelizzabile della Galerkin continua?

Ho sempre sentito che la semplice parallelizzazione era uno dei vantaggi dei metodi DG, ma non vedo davvero perché uno di questi motivi non si applichi anche al Galerkin continuo.

9 finite-element parallel-computing discontinuous-galerkin

1

Implementazione del metodo Jacobi-Davidson per il problema degli autovalori cubici

Ho un grosso problema di autovalori cubici: (A0+λA1+λ2A2+λ3A3)x=0.(A0+λA1+λ2A2+λ3A3)x=0.\left(\mathbf{A}_0 + \lambda\mathbf{A}_1 + \lambda^2\mathbf{A}_2 + \lambda^3\mathbf{A}_3\right)\mathbf{x} = 0. Potrei risolverlo convertendomi in un problema di autovalore lineare ma si tradurrebbe in un sistema di dimensioni:32323^2 ⎡⎣⎢−A0000I000I⎤⎦⎥⎡⎣⎢xyz⎤⎦⎥=λ⎡⎣⎢A1I0A20IA300⎤⎦⎥⎡⎣⎢xyz⎤⎦⎥,[−A0000I000I][xyz]=λ[A1A2A3I000I0][xyz],\begin{bmatrix} -\mathbf{A}_0 & 0 & 0 \\ 0 & \mathbf{I} & 0 \\ 0 & 0 & …

9 c++ eigensystem eigenvalues

10

È possibile ottimizzare questo codice di integrazione in modo che funzioni più velocemente?

double trap(double func(double), double b, double a, double N) { double j; double s; double h = (b-a)/(N-1.0); //Width of trapezia double func1 = func(a); double func2; for (s=0,j=a;j<b;j+=h){ func2 = func(j+h); s = s + 0.5*(func1+func2)*h; func1 = func2; } return s; } Quanto sopra è il mio codice …

9 c++ performance

2

Rappresenta i numeri di Eisenstein senza galleggianti

Ho un progetto in cui devo usare i campi quadratici In particolare i numeri del modulo con . a,b∈Qa + b - 3---√un'+B-3a + b \sqrt{-3}a , b ∈ Qun',B∈Qa,b \in \mathbb{Q} Ad esempio, ecco i numeri primi negli interi di Eisenstein : Non voglio usare la salvia. Vorrei scrivere …

9 linear-algebra precision

3

Metodo degli elementi finiti vs Metodo degli elementi finiti esteso (FEM vs XFEM)

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?

9 finite-element numerical-analysis adaptive-mesh-refinement numerics

2

Sforzo di calcolo degli algoritmi

Considera il problema di ottimizzazione senza strettamente convessoLascia che x_ \ text {opt} denoti i suoi minimi unici e x_0 sia una data approssimazione iniziale di x_ \ text {opt}. Chiameremo un vettore x una soluzione \ epsilon- close di \ mathcal {O} if \ begin {equation} \ frac {|| …

9 convex-optimization complexity

3

Sorta di problemi in cui SOR è più veloce di Gauss-Seidel?

Esiste una semplice regola empirica per dire se vale la pena fare SOR invece di Gauss-Seidel? (e modo possibile come stimare il parametro di realxation )ωω\omega Intendo solo guardando la matrice o la conoscenza di un problema particolare che la matrice rappresenta? Stavo leggendo la risposta a queste domande: ci …

9 linear-solver iterative-method heuristics

2

suggerimenti per la gestione delle corse di simulazione?

Questa domanda potrebbe essere un po 'fuori tema in comp-sci. se necessario, suggerisci dove si adatta. La domanda riguarda come gestire in modo efficiente tutte le simulazioni. diciamo, ad esempio, che una simulazione richiede la correzione di 2 parametri che devono essere definiti in un determinato intervallo di valori suggerito. …

9 monte-carlo data-management

3

Il mio progetto può sopravvivere senza Orientamento agli oggetti?

Sto scrivendo un piccolo pacchetto MATLAB che risolverà alcuni problemi numerici di classe. Ci sono 3 fasi dell'algoritmo e l'utente ha 5 scelte per ogni fase. Ho implementato l'intero problema usando funzioni e 3 casi di commutazione (uno per ogni fase dell'algoritmo). Funziona bene ma sto pensando di fargli fare …

9 matlab

1

Esiste una buona mailing list per `Computational Science`?

Mi chiedo se ci sia una mailing list o gruppi di Google molto validi per la scienza computazionale, dove possiamo discutere domande invece di solo porre e rispondere alle domande. In realtà, sono più interessato al calcolo parallelo e alle soluzioni numeriche dei PDE. Ma non so cosa e come …

9 parallel-computing reference-request

2

Errore peculiare nella risoluzione dell'equazione di Poisson su un metodo a volume finito mesh non uniforme (solo 1D)

Ho cercato di eseguire il debug di questo errore negli ultimi giorni, mi chiedevo se qualcuno avesse dei consigli su come procedere. Sto risolvendo l'equazione di Poisson per una distribuzione della carica a gradino (un problema comune nella fisica elettrostatica / dei semiconduttori) su una maglia a volume finito non …

9 discretization finite-volume poisson

4

Learning Fortran

Mi chiedevo qual è il modo migliore per imparare a programmare in Fortran? Ho cercato libri, gli esempi che mi sono sembrati OK erano "Fortran 95/2003 Explained" e "Fortran For Scientists And Engineers" . Non so cosa ne pensiate di questi ragazzi? Un'altra cosa che mi chiedo è quale compilatore …

9 fortran

1

Ci sono vantaggi numerici nel risolvere la matrice simmetrica rispetto alle matrici senza simmetria?

Sto applicando il metodo della differenza finita a un sistema di 3 equazioni accoppiate. Due delle equazioni non sono accoppiate, tuttavia la terza equazione si accoppia ad entrambe le altre due. Ho notato che cambiando l'ordine delle equazioni, dire da a ( x , z , y ) che la …

9 finite-difference symmetry

3

Multigrid FAS più lento della correzione lineare dei difetti?

Ho implementato un solutore multigrid V-Cycle utilizzando sia una correzione lineare dei difetti (LDC) sia uno schema di approssimazione completo (FAS). Il mio problema è il seguente: usando LDC il residuo è ridotto di un fattore di ~ 0,03 per ciclo. L'implementazione FAS converge anche con un fattore lineare, ma …

9 multigrid