Sto usando il pacchetto di risolutore non lineare SNES di PETSc per risolvere un sistema di equazioni non lineari ottenuto discretizzando un'equazione differenziale parziale. Come posso determinare perché il risolutore non converge e cosa posso fare per risolvere con successo le mie equazioni?
Il paradigma di programmazione OpenCL promette di essere uno standard aperto senza royalty per il calcolo eterogeneo. Dovremmo investire il nostro tempo nello sviluppo di software basato su OpenCL? Vantaggi svantaggi?
Quale sarebbe una discreta discretizzazione a differenza finita per la seguente equazione: ∂ρ∂t+∇⋅(ρu)=0∂ρ∂t+∇⋅(ρu)=0\frac{\partial \rho}{\partial t} + \nabla \cdot \left(\rho u\right)=0 ? Possiamo prendere il caso 1D: ∂ρ∂t+ddx(ρu)=0∂ρ∂t+ddx(ρu)=0\frac{\partial \rho}{\partial t} + \frac{d}{dx}\left(\rho u\right)=0 Per qualche ragione tutti gli schemi che riesco a trovare sono per la formulazione in coordinate lagrangiane. Ho …
Ho una funzione bidimensionale cui vorrei campionare i valori. La funzione è molto costosa da calcolare e ha una forma complessa, quindi ho bisogno di trovare un modo per ottenere il maggior numero di informazioni sulla sua forma usando il minor numero di punti campione.f(x,y)f(x,y)f(x,y) Quali buoni metodi ci sono …
Come sappiamo, i metodi di trasformazioni ortogonali (rotazioni di Givens e riflessioni di Housholder) per i sistemi di equazioni lineari sono più costosi dell'eliminazione gaussiana, ma teoricamente hanno proprietà di stabilità più belle, nel senso che non cambiano il numero di condizioni del sistema. Anche se conosco solo un esempio …
Andare a insegnare agli studenti di livello universitario un corso intitolato Introduzione alla programmazione informatica. Sono un po 'confuso. Nella fisica computazionale gli scienziati usano C / C ++ o Python o Fortran, CUDA ecc ..... questo è il momento di costruire la loro base. Cosa dovrei usare? So che …
Qual è l'approccio preferito ed efficiente per l'interpolazione di dati multidimensionali? Cose di cui sono preoccupato: prestazioni e memoria per costruzione, valutazione singolo / batch dimensioni di movimentazione da 1 a 6 lineare o di ordine superiore capacità di ottenere gradienti (se non lineari) griglia regolare vs sparsa usando come …
Sto cercando libri o articoli, o post di blog, o qualsiasi materiale pubblicato in generale, che riguardi specificamente gli usi delle funzionalità moderne del C ++ (spostare la semantica, l'STL, gli iteratori, la valutazione lenta, ecc.) Nell'informatica scientifica. Puoi suggerire qualcuno? Penso che queste nuove funzionalità renderanno più semplice la …
Inizierò con la mia esperienza personale nel nostro laboratorio. Negli ifort 9 e 10 giorni, eravamo piuttosto aggressivi con le ottimizzazioni, compilando con -O3 e flag specifici del processore (-xW -xSSE4.2 per esempio). Ma a partire da ifort 11, abbiamo iniziato a notare: 1. alcune incongruenze nei risultati (perché la …
Come si adattano le matrici Python / Numpy all'aumentare delle dimensioni dell'array? Questo si basa su alcuni comportamenti che ho notato durante il benchmarking del codice Python per questa domanda: Come esprimere questa espressione complicata usando le sezioni intorpidite Il problema riguardava principalmente l'indicizzazione per popolare un array. Ho scoperto …
Esistono diverse librerie che risolvono un sistema di equazioni sparse e lineari, tuttavia trovo difficile capire quali siano le differenze. Per quanto ne so, ci sono tre pacchetti principali: Trilinos , PETSc e Intel MKL . Sono tutti in grado di eseguire soluzioni a matrice sparsa, sono tutti veloci (per …
Sto cercando algoritmi per disegnare grafici 2D standard per funzioni che possono o meno avere singolarità. Lo scopo è scrivere un "Mini-CAS", quindi non ho una conoscenza a priori dei tipi di funzioni, gli utenti vogliono rappresentare graficamente. Questo problema è molto vecchio, quindi immagino che ci debbano essere alcuni …
Le formule standard delle differenze finite sono utilizzabili per calcolare numericamente una derivata nell'aspettativa che tu abbia valori di funzione in punti equidistanti, in modo che sia una costante. Cosa succede se ho punti spaziati in modo non uniforme, in modo che ora varia da una coppia di punti adiacenti …
Sto cercando di risolvere un problema di ottimizzazione vincolata in cui conosco i limiti di alcune delle variabili (in particolare un vincolo inscatolato). argminuf( u , x )argminuf(u,x) \arg \min_u f(u,x) soggetto a a ≤ d ( u , x ) ≤ bc ( u , x ) = 0c(u,x)=0 …
Esiste un metodo per risolvere sistemi lineari della forma dove è una matrice SPD fissa e sono matrici diagonali positive?k ( D i + A ) x i = b i A D iO(n3+n2k)O(n3+n2k)O(n^3+n^2 k)kkk(Di+A)xi=bi(Di+A)xi=bi(D_i + A) x_i = b_iAAADiDiD_i Ad esempio, se ogni è scalare, è sufficiente calcolare la …
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.