Scienza computazionale

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

3
Domanda di approssimazione dei minimi quadrati
Sto seguendo un corso di calcolo scientifico e abbiamo appena superato l'approssimazione dei minimi quadrati. La mia domanda riguarda in particolare l'approssimazione con i polinomi. Capisco che se hai n + 1 punti dati, puoi trovare un polinomio unico di grado n che descrive tutti questi punti. Ma posso anche …
1
Metodi numerici per invertire le trasformazioni integrali?
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 …
2
Come posso tracciare la superficie di un diagramma 4D?
Sto cercando di tracciare la funzione d'onda per una particella in una scatola 3D. Questo mi richiede di tracciare 4 variabili: assi x, y, z e la funzione di densità di probabilità. La funzione di densità di probabilità è: abs((np.sin((p*np.pi*X)/a))*(np.sin((q*np.pi*Y)/b))*(np.sin((r*np.pi*Z)/c)))**2 Sto usando np.arange()per X, Y e Z. Ho letto che …
6
Esiste un'implementazione a livello di riferimento di BLAS in C / C ++?
L'implementazione BLAS di netlib è un riferimento eccellente, essendo per lo più non ottimizzata e ben documentata (ad es. Zgemm ). Tuttavia, è in Fortran 77, rendendolo in qualche modo inaccessibile a quelli con una formazione di programmazione più moderna. Esiste un'implementazione a livello di riferimento di BLAS, come netlib, …
11 blas  education 
3
Dato un sistema lineare tridiagonale SPD, possiamo precompilare in modo che tre o tre indici possano essere collegati nel tempo O (1)?
Si consideri un simmetrica definita positiva sistema lineare tridiagonale dove e . Dati tre indici , se assumiamo solo righe di equazioni strettamente tra e hold, possiamo eliminare le variabili intermedie per ottenere un'equazione della forma dove . Questa equazione mette in relazione il valore di con indipendentemente dall'influenza "esterna" …
3
Meshing dei dati di superficie 3D in Python
Ho un set di dati di punti tridimensionali per i quali mi piacerebbe costruire una mesh, usando Python. Tutto il software che ho visto richiede di fornire i bordi. Esiste un programma in Python che accetta come input una serie di punti in 3D e produce una mesh? Se possibile, …
3
Verifica se due matrici 12x12 hanno lo stesso determinante
12×1212×1212 \times 12QQQdet(Q)=det(12I−Q−J)(1)det(Q)=det(12I−Q−J)(1)\det(Q) = \det(12I-Q-J) \; \; (1)JJJ Attualmente lo sto facendo con la libreria di armadillo ma risulta essere troppo lento. Il fatto è che devo farlo per un trilione di matrici e si scopre che il calcolo dei due determinanti è il collo di bottiglia del mio programma. …
3
Le migliori metodologie per la gestione di una mesh nel calcolo di elementi finiti paralleli?
Attualmente sto sviluppando un metodo di decomposizione del dominio per la soluzione del problema di scattering. Fondamentalmente sto risolvendo un sistema di Helmholtz BVP in modo iterativo. Discretizzo le equazioni usando il metodo degli elementi finiti su maglie triangolari o tetraedriche. Sto sviluppando il codice per la mia tesi di …
3
Come dovrebbero essere trattati i coefficienti non costanti con lo schema di controvento del primo ordine a volume finito?
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à, …
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.