Kürzlich habe ich diese Frage gesehen: Wie misst man den Fehler einer Finite-Differenzen-Methode? Ich bin Student der Simulationswissenschaften und leider ist für mich völlig unklar, welche Norm in welchem Kontext zu verwenden ist. Sehr oft verwenden wir die euklidische Norm oder die L2-Norm,...
Ich habe einen Teil eines Problems, das durch die Impulserhaltungsgleichung beschrieben wird: ∂ρ∂t+1sinθ∂∂θ(ρusinθ)=0∂ρ∂t+1sinθ∂∂θ(ρusinθ)=0\frac{\partial \rho}{\partial t} + \frac{1}{\sin\theta} \frac{\partial}{\partial \theta}(\rho u \sin \theta) =0 Wobei und (konstante...
Ich habe viele Veröffentlichungen in Fachzeitschriften für Computerphysik gesehen, die unterschiedliche Metriken für die Leistung ihres Codes verwenden. Insbesondere für GPGPU-Code scheint es eine Vielzahl von Timing-Ergebnissen zu geben, die von Menschen veröffentlicht werden. Insbesondere habe...
Ich habe eine Reihe von MPI-Prozessen, die Anforderungsnachrichten hin und her austauschen. Prozesse wissen nicht, welche anderen Prozesse ihnen Nachrichten senden oder wie viele. In dieser Situation möchte ich auf effiziente Weise feststellen, ob alle anderen Prozesse das Senden von Nachrichten...
Die Wärmegleichung wird mit FV (oder FEM) im Raum diskretisiert, und es wird eine halbdiskrete Gleichung erhalten (System von ODEs). Dieser Ansatz, der als Linienmethode bekannt ist , ermöglicht den einfachen Wechsel von einer zeitlichen Diskretisierung zu einer anderen, ohne dass Code dupliziert...
Ich entschuldige mich im Voraus, wenn diese Frage dumm ist. Ich muss die Wurzel von berechnen u−f(u)=0u−f(u)=0\begin{equation} u -f(u) =0 \end{equation} Wobei ein reeller Vektor ist und f ( u ) eine reelle Vektorwertfunktion ist. Ich begann mit Newtons Methode (die funktionierte), erkannte dann...
Ich habe über die Implementierung einer dichten Matrixmultiplikation gelesen, wenn die Matrix nicht in den Cache passt. Eine der Grafiken, die ich gesehen habe (Folie 9 von diesen Folien ), zeigt plötzliche Leistungseinbußen unter Verwendung des naiven Algorithmus. Diese Tropfen betragen ungefähr...
Ich habe kürzlich einen großen Teil des Legacy-Codes geerbt, der ein sehr steifes, vorübergehendes Problem löst. Ich möchte zeigen, dass die räumlichen und zeitlichen Schrittgrößen so klein sind, dass sich die qualitative Natur der berechneten Lösung nicht ändert, wenn sie verringert wird. Mit...
Ich möchte mein Python-Programm mit Cython schnell machen, aber meine innere Schleife führt immer noch langsame Python-Aufrufe an den Zufallszahlengenerator durch! Vor einigen Jahren wurde dasselbe Problem von jemandem angesprochen, der Salbei unterstützt, und es schien zu diesem Zeitpunkt keine...
Unter der Annahme, dass eine direkte numerische Simulation durchgeführt wird, was ist eine gute Methode, um den Einfluss von kleinem Maßstab auf großen Maßstab in der Fluiddynamik abzuschätzen? Ist es beispielsweise relevant, zwei Läufe mit unterschiedlicher Gittergröße oder zwei Läufe mit...
Ich muss das System zweier gekoppelter partieller Differentialgleichungen numerisch lösen. ∂x1∂t∂x2∂t= c1∇2x1+ f1( x1, x2)= c2∇2x2+ K.∂x1∂t∂x1∂t=c1∇2x1+f1(x1,x2)∂x2∂t=c2∇2x2+K∂x1∂t\begin{align} \frac{\partial x_1}{\partial t} &= c_1\nabla ^2 x_1 + f_1(x_1,x_2)\\ \frac{\partial x_2}{\partial t} &=...
Ich versuche mehr über Volume Meshing zu erfahren und habe beschlossen, einen einfachen Volume Mesher zu implementieren. Die Strategie, die ich gewählt habe, besteht darin, meinen Raum mit einem Octree zu unterteilen, der anhand einiger Kriterien verfeinert wurde. Ich habe dafür gesorgt, dass mein...
Gibt es eine Open-Source-Implementierung von konvexen Optimierungslösern in C / C ++? Einige Bibliotheken wie NLopt , Ipopt , OPT ++ erfüllen meine Anforderungen nicht. OPT ++ und Ipopt sind nicht threadsicher, und NLopt scheint keine spezifische / effiziente Möglichkeit zu haben, konvexe...
Ich habe zwei Variablen kund tals Funktionen von zwei anderen Variablen p1und p2. Ich kenne auch ihre Maximalwerte. Ich habe keinen analytischen Ausdruck dafür. Ich möchte die Werte von kund tfinden, die ihren Maximalwerten am nächsten kommen. Gibt es eine Möglichkeit, das k = f1(p1, p2)und zu...
Ich suche nach Beispielen für Helmholtz- und Biharmonische Gleichungen in kartesischen Koordinaten mit exakten Lösungen, um meine numerischen Lösungen damit zu vergleichen. Im Internet konnte ich einige Beispiele finden, bei denen das Problem mit den Randbedingungen genau definiert war. Dies waren...
Ich versuche, neuronale Netze in Python erneut zu implementieren. Ich habe die Kostenfunktion und den Backpropagation-Algorithmus korrekt implementiert. Ich habe sie überprüft, indem ich den entsprechenden Octave-Code ausgeführt habe. Wenn ich jedoch versuche, die scipy.optimize.fmin_cgFunktion zu...
Was ist der effizienteste Weg, um eine komplexe Matrix mit geringer Dichte von meinem Fortran-Code auf PETSc zu übertragen? Ich verstehe, dass dies problemabhängig ist, deshalb habe ich versucht, so viele relevante Details wie möglich unten anzugeben. Ich habe mit dem FEAST-Eigenwertlöser [1] für...
Ich habe einen Navier Stokes implementiert, der auf FEM-Diskretisierung und PETSc basiert, um das lineare Gleichungssystem zu lösen. Um ein effizientes Lösungsverfahren zu erstellen, folge ich dem Artikel "Effiziente Vorkonditionierung der linearisierten Navier-Stokes-Gleichungen für...
Wenn ich ein großes Array verwenden möchte, z. B. mpz_t A [100000], wurde während der Kompilierung "Segmentierungsfehler (Core Dumped)" angezeigt. Gibt es einen einfacheren Weg, dies zu
Ich versuche, den Rebinning-Algorithmus der Monte-Carlo-Integration von VEGAS ( Originalveröffentlichung ( Preprint von LKlevin) und Implementierungshinweise ) zu verstehen . Ich werde versuchen, zuerst zu erklären, was ich zu verstehen glaube, und dann meine Fragen zu stellen. Nehmen wir der...