Computational Science

Fragen und Antworten für Wissenschaftler, die Computer verwenden, um wissenschaftliche Probleme zu lösen

1
Verschachtelte Dissektion im regulären Raster
Bei der Lösung spärlicher linearer Systeme mit direkten Faktorisierungsmethoden wirkt sich die verwendete Ordnungsstrategie erheblich auf den Füllfaktor von Nicht-Null-Elementen in den Faktoren aus. Eine solche Bestellstrategie ist die verschachtelte Dissektion. Ich frage mich, ob es möglich ist, die verschachtelte Dissektionsreihenfolge vorzeitig zu erstellen, wenn nur die Gitterparameter gegeben sind …
1
Was ist die genaueste Interpolationsmethode für ein 3D-Flussfeld in einem strukturierten Gitter?
Ich löse komprimierbare Navier-Stokes-Gleichungen für mehrere Arten in einem strukturierten 3D-Gitter. Ich habe eine Lösung für ein bestimmtes Gitter erhalten (sagen wir eine relativ grobe). Ich möchte jetzt mein Raster verfeinern und meine vorherige Lösung auf meinem neuen Raster interpolieren, bevor ich meine Simulation neu starte. Derzeit haben wir ein …
2
Kopplung von FEM DG-Methoden an Riemann-Löser
Gibt es gute Papiere und / oder Codes, die diskontinuierliche Galerkin-Finite-Elemente-Löser mit Riemann-Lösern koppeln? Ich muss die elliptischen und hyperbolischen Kopplungsprobleme untersuchen, aber die meisten Aufteilungsmethoden sind bestenfalls ad hoc. Da ich eine große Menge FEniCS-Code habe, möchte ich nur den Riemann-Solver damit koppeln. Während ein einfacher Roe-Löser ein Anfang …
3
Numerik: Wie normalisiere ich die folgende ODE?
In dieser Frage geht es mehr darum, wie ein Problem numerisch angegangen werden kann. In einem kleinen Projekt wollte ich die Coorbitalbewegung von Janus und Epimetheus simulieren. Dies ist im Grunde ein Drei-Körper-Problem. Ich wähle Saturn als Fix am Ursprung, sei und r 2 die Ortsvektoren von Janus bzw. Epimetheus. …
9 ode 
2
numerische Integration mit möglicher Division durch 'Null'
Ich versuche mich zu integrieren ∫10t2 n + 2exp(α r0t) dt∫01t2n+2exp⁡(αr0t)dt\int^1_0 t^{2n+2}\exp\left({\frac{\alpha r_0}{t}}\right)dt Das ist eine einfache Transformation von ∫∞1x2 nexp( - α r0x ) dx∫1∞x2nexp⁡(- -αr0x)dx\int^{\infty}_1 x^{2n}\exp(-\alpha r_0 x)dx mit weil es schwierig ist, falsche Integrale numerisch zu approximieren. Dies führt jedoch zu dem Problem der Bewertung des neuen …
3
Faustregel für spärliche oder dichte Matrixspeicherung
Angenommen, ich kenne die erwartete Sparsamkeit einer Matrix (dh die Anzahl der Nicht-Nullen / die insgesamt mögliche Anzahl der Nicht-Nullen). Gibt es eine Faustregel (möglicherweise eine ungefähre Angabe) für die Entscheidung, ob ein spärlicher Matrixspeicher (insbesondere ein komprimierter Zeilenspeicher) verwendet oder als dichte Matrix gespeichert werden soll? Geschwindigkeit ist in …
6
Super C ++ Optimierung der Matrixmultiplikation mit Armadillo
Ich verwende Armadillo, um sehr intensive Matrixmultiplikationen mit Seitenlängen von , wobei bis zu 20 oder sogar mehr kann. Ich verwende Armadillo mit OpenBLAS für die Matrixmultiplikation, was in parallelen Kernen sehr gute Arbeit zu leisten scheint, außer dass ich ein Problem mit dem Formalismus der Multiplikation in Armadillo zur …
5
Motivation hinter der Galerkin-Methode
Ich habe eine Frage zur Galerkin-Methode. Ich verstehe nicht, warum die Galerkin-Methode den Rest durch die Formfunktionen gewichtet und gleich Null setzt. Ich möchte wissen, was der Grund dafür ist. Warum müssen wir Gewichtungsrestfunktionen gleich Null setzen?
3
Wie verwende ich die Polylogarithmusfunktion in c ++?
Gibt es Präprozessoranweisungen, die zur Verwendung der Polylog-Funktion verwendet werden könnten? Oder ist es in cmath enthalten? Wenn ja, nennen Sie es Li oder Polylog? EDIT: Was ich wirklich versuche, ist einen analytischen Wert für das unbestimmte Integral der Funktion anzugeben x3ex- 1x3ex- -1 \frac{x^3}{e^{x} -1} welches Polylogarithmusfunktionen beinhaltet. Aber …
4
Reihenfolge der MATLAB FFT-Frequenzen
Dieses Wikibook gibt an, dass die Ausgabe von FFTMATLABs den Wellenzahlen entspricht, die wie folgt bestellt sind: k = { 0 , 1 , . . . , n2, - n2+ 1 , - n2+ 2 , . . . , - 1 }k={0,1,...,n2,−n2+1,−n2+2,...,−1}k=\left\{0,1,...,\frac{n}{2},-\frac{n}{2}+1,-\frac{n}{2}+2,...,-1\right\} In den Beispielcodes auf derselben Seite …
2
Wie kann man inhomogene Dirichlet-Randbedingungen mit FEM richtig anwenden?
Im Allgemeinen werden Dirichlet-Randbedingungen für FEM für inhomogene Randbedingungen nicht genau erfüllt. Die FEM-Codes, die ich gesehen habe, legen die Freiheitsgrade für die Interpolation der Dirichlet-Randbedingung fest, aber ich habe keine mathematische Rechtfertigung dafür gefunden. Es scheint mir, dass das Festlegen wesentlicher Randbedingungen wahrscheinlich einige Funktionen des Fehlers minimieren sollte …
Durch die Nutzung unserer Website bestätigen Sie, dass Sie unsere Cookie-Richtlinie und Datenschutzrichtlinie gelesen und verstanden haben.
Licensed under cc by-sa 3.0 with attribution required.