Ich arbeite an einer 2D-Flüssigkeitssimulation mit Wirbelpartikeln / "Wirbeln", wie in Flüssigkeitssimulation für Videospiele beschrieben . Was ich denke, ist das gleiche wie die "diskrete Wirbelmethode". Grundsätzlich stellen Sie die Flüssigkeit mit einer Ansammlung von Partikeln mit definierter Vorticity dar und berechnen die Fluidgeschwindigkeit an einem Punkt, indem Sie die …
Ich arbeite derzeit an der Lösung sehr großer symmetrischer (aber nicht positiv definierter) Systeme, die von bestimmten Algorithmen generiert werden. Diese Matrizen haben eine schöne Blocksparsity, die zum parallelen Lösen verwendet werden kann. Ich kann mich jedoch nicht entscheiden, ob ich einen direkten Ansatz (wie Multi-Frontal) oder einen iterativen Ansatz …
Als ich heute über die Theorie der Differentialformen las, war ich beeindruckt, wie sehr sie mich an die Finite-Volumen-Methode zweiter Ordnung (FVM) erinnerte. Ich kämpfe darum herauszufinden, ob ich auf diese Weise nur trivial denke oder ob es eine tiefere Verbindung gibt. Nun, Differentialformen dienen dazu, einige Konzepte zu verallgemeinern, …
Ich sehe nur sehr wenige Nicht-Gleitkomma-Computer-Bibliotheken / -Pakete. Angesichts der verschiedenen Ungenauigkeiten der Gleitkommadarstellung stellt sich die Frage, warum es nicht zumindest einige Bereiche gibt, in denen diese erhöhte Genauigkeit die Komplexität der Arbeit mit Festkommawerten wert sein könnte. Gibt es größere Schwierigkeiten bei der Verwendung beispielsweise eines Festpunkt-Eigenwertlösers? Wie …
Ich bin neu im PETSc-Paket. Ich habe eine ~ 4000x4000 Matrix A im Matrix-Market-Format und möchte PETSc dazu bringen, dies mit mehreren Prozessoren zu lösen. Ich weiß, wie man das System auf einem einzelnen Prozessor löst, aber ich weiß nicht, wie man die Matrix und die Vektoren auf verschiedene Prozessoren …
Ich suche nach der wichtigsten Literatur und Folienreferenzen zur Modellierung der aktuellen und zukünftigen Energiekosten von Gleitkommaoperationen und Datenübertragungen über CPU, Speicher, Netzwerk und Speicher. Ich habe diese Frage als Community-Wiki markiert und möchte, dass Sie jede Antwort auf das folgende Format beschränken: Titel, Autoren, Ort / Konferenz / Zeitschrift, …
Wenn wir ein numerisches Modell haben, das ein reales physikalisches System darstellt und Chaos aufweist (z. B. Modelle der Fluiddynamik, Klimamodelle), wie können wir dann wissen, dass das Modell so funktioniert, wie es sollte? Wir können zwei Sätze von Modellausgaben nicht direkt vergleichen, da selbst kleine Änderungen der Anfangsbedingungen die …
Ich bin daran interessiert, eine Funktion maximieren , wobei θ ∈ R p ist .f( θ )f(θ)f(\mathbf \theta)θ∈Rpθ∈Rp\theta \in \mathbb R^p Das Problem ist, dass ich die analytische Form der Funktion oder ihrer Ableitungen nicht kenne. Das einzige , was ich tun kann , ist die Funktion punktweise zu bewerten, …
Wie werden Spezialfunktionen mit doppelter Genauigkeit auf dem neuesten Stand der Technik implementiert? Ich brauche die folgende Integral: fürm=0,1,2,. . . undt>0, was als niedrigere unvollständige Gammafunktion geschrieben werden kann. Hier ist meine Fortran- und C-Implementierung:F.m( t ) = ∫10u2 me- t u2du = γ( m + 12, t )2 …
Ich habe eine numerische Funktion, f(x, y)die eine doppelte Gleitkommazahl zurückgibt, die eine Formel implementiert, und ich möchte überprüfen, ob sie gegen analytische Ausdrücke für alle Kombinationen der Parameter korrekt ist xund ymich interessiert. Was ist der richtige Weg, um die berechneten und zu vergleichen? analytische Gleitkommazahlen? Nehmen wir an, …
In den letzten Jahren sind mehrere Bibliotheken / Softwareprojekte erschienen, die die eine oder andere Form der datengesteuerten Parallelität für gemeinsam genutzten Speicher für allgemeine Zwecke bieten. Die Hauptidee ist, dass Programmierer ihre Algorithmen nicht als explizit getunten Code schreiben, sondern als voneinander abhängige Aufgaben implementieren, die dann von einer …
Ich versuche, eine MPI-Anwendung mit einem stark asynchronen Kommunikationsmuster zu optimieren. Jeder Rang hat eine Liste der zu berechnenden Dinge und sendet nach Bedarf Nachrichten, wenn sich die Ein- oder Ausgänge auf einem anderen Rang befinden. Zusätzlich wird jeder Rang mit einem Thread versehen (derzeit mit einem Kommunikationsthread und 5 …
Gibt es eine offene C-Implementierung für die Lösung von Quartalsgleichungen: ax⁴+bx³+cx²+dx+e=0ax⁴+bx³+cx²+dx+e=0ax⁴+bx³+cx²+dx+e=0 Ich denke an eine Implementierung der Ferrari-Lösung. Auf Wikipedia habe ich gelesen, dass die Lösung nur für einige der möglichen Vorzeichenkombinationen der Koeffizienten rechenstabil ist. Aber vielleicht habe ich Glück ... Ich habe eine pragmatische Lösung erhalten, indem ich …
Betrachten Sie das folgende Problem bei dem der Forcierterm von ( Formulierung siehe Edit 1 unten) und und seinen ersten Ableitungen abhängen kann . Dies ist eine 1 + 1-dimensionale Wellengleichung. Wir haben Anfangsdaten, die bei .u , v W { u + v = 0 }W.u v= F.Wuv=F W_{uv} …
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.