Ich habe ein Programm / eine Bibliothek geschrieben, mit der ich in einem Artikel Ergebnisse erzielt habe. (Hier ist , aber meine Frage ist allgemein.) Ich habe Tests , dass ich regelmäßig laufen mit ctest(es dauert ein paar Minuten zu laufen). Um einige Tabellen oder Abbildungen im Artikel zu...
Es gibt ein System linearer Bedingungen . Ich möchte einen streng positiven Vektor , der diese Bedingungen erfüllt. Das heißt, wird für jede Komponente von . Wie kann ich einen LP-Solver verwenden, um einen solchen streng positiven Vektor (oder um zu bestätigen, dass kein existiert)? Ich kann...
Ich habe anekdotisch gehört, dass, wenn man versucht, ein Integral der Form numerisch zu machen ∫∞0f(x)J0(x)dx∫0∞f(x)J0(x)dx\int_0^\infty f(x) J_0(x)\,\mathrm{d}x mit glatt und gutmütig (z. B. nicht selbst stark oszillierend, nicht singulär usw.), dann hilft es Genauigkeit, es als...
Ich habe eine sehr grundlegende Vorstellung von Intervallarithmetik (IA), aber es scheint sowohl theoretisch als auch praktisch ein sehr interessanter Zweig der Computerwissenschaft zu sein. Es ist klar, dass es sich bei den offensichtlichen Anwendungen um verifizierte Computer- und unsachgemäße...
Ich habe das auf Math SE gelesen . Die Grundfrage ist: Angenommen, jemand möchte etwas Fortgeschrittenes lernen. Eine Möglichkeit, dies zu tun, besteht darin, von den Grundlagen auszugehen und aufzubauen. Dabei kann jedoch das "Gesamtbild" verloren gehen. Eine weitere Methode (die ich lieber als...
In der Zeitung stieß ich auf eine verwirrende Bemerkung PJ van der Houwen, Die Entwicklung von Runge-Kutta-Methoden für partielle Differentialgleichungen, Appl. Num. Mathematik. 20: 261,1996 In den Zeilen 8ff auf Seite 264 schreibt van der Houwen: "Für die Taylor - Polynome impliziert dies , dass...
Ich bin sehr beeindruckt von der Serienleistung von mehrstufigen inversen ILU-Vorkonditionierern , insbesondere für heterogenes Helmholtz , aber ich bin überrascht, keine Open-Source-Implementierungen finden zu können. Insbesondere stellt ILUPACK den Wissenschaftlern Binärdateien frei zur...
Ich habe ein lineares Gleichungssystem der Größe mxm, wobei m groß ist. Die Variablen, die mich interessieren, sind jedoch nur die ersten n Variablen (n ist klein im Vergleich zu m). Gibt es eine Möglichkeit, die Lösung für die ersten m-Werte zu approximieren, ohne das gesamte System lösen zu...
Ausgehend von meiner vorherigen Frage versuche ich, Randbedingungen auf dieses ungleichmäßige endliche Volumen-Netz anzuwenden. Ich möchte eine Robin-Typ-Randbedingung auf die lhs der Domäne anwenden ( x=xL)x=xL)x=x_L) , so dass σL=(dux+au)∣∣∣x=xLσL=(dux+au)|x=xL \sigma_L = \left( d u_x + a u...
Lassen Sie x, ywerden zwei Gleitkommazahlen. Wie kann man den Mittelwert richtig berechnen? Der naive Weg (x+y)/2kann zu Überläufen führen, wenn xund ysind zu groß. Ich denke 0.5 * x + 0.5 * yvielleicht besser, aber es geht um zwei Multiplikationen (was vielleicht ineffizient ist), und ich bin...
Was sind die neuesten Methoden zur numerischen Lösung von ODEs mit diskontinuierlicher rechter Seite? Mich interessieren vor allem stückweise glatte rechte Seitenfunktionen, zB sign. Ich versuche die Gleichung von folgendem Typ zu lösen: x˙v˙= v= { ( | Fextern| - | FReibung| )unterschreiben(...
Ich habe einen Einblick in die Numerische Analyse (hauptsächlich numerische Methoden wie das Finden von Wurzeln, quadratische Gleichungen und andere vorläufige Dinge) in meiner Calculus-Klasse erhalten, aber jetzt möchte ich mehr Raffinesse in meiner Arbeit. Gibt es ein gutes Buch, das mir hilft,...
Ich unterrichte eine Umfrageklasse zur numerischen Analyse und suche nach Motivation für die BFGS-Methode für Studenten mit begrenztem Hintergrund / Intuition in der Optimierung! Ich habe zwar keine Zeit, konsequent zu beweisen, dass alles konvergiert, aber ich möchte eine angemessene Motivation...
Ich benutze HDF5 seit Jahren, aber mit zunehmender Größe des Datensatzes treten die gleichen Probleme auf, die hier aufgeführt sind http://cyrille.rossant.net/moving-away-hdf5/ Können Sie mich auf ein Format wie HDF5 verweisen mit - Gute Unterstützung für paralleles Schreiben - Unterstützung für...
Gibt es einen (effizienten) Algorithmus, um eine Teilmenge von Punkten aus einer Menge von Punkten ( ) so auszuwählen, dass sie den größten Bereich "abdecken" (über alle möglichen Teilmengen der Größe )?MMMNNNM< NM<NM < NMMM Ich gehe davon aus, dass die Punkte in der 2D-Ebene liegen. Der...
Nachdem ich einige mathematische Untersuchungen zur Stabilität von Elementen im 3D-Stokes-Problem durchgeführt hatte, stellte ich leicht schockiert fest, dass für ein beliebiges Tetraedernetz nicht stabil ist. Genauer gesagt, wenn Sie ein Element haben, bei dem alle Knoten und drei von vier...
Manchmal, während ich den Code optimiere, ist es erforderlich, bestimmte Teile des Codes zeitlich zu bestimmen. Ich verwende das Folgende seit Jahren, habe mich aber gefragt, ob es einen einfacheren / besseren Weg gibt, dies zu tun? call system_clock(count_rate=clock_rate) !Find the time rate call...
Ich erstelle einen einfachen Astronomiesimulator, der die Newtonsche Physik verwenden soll, um die Bewegung von Planeten in einem System (oder anderen Objekten) zu simulieren. Alle Körper sind Kreise in einer euklidischen Ebene mit Eigenschaften wie Position, Geschwindigkeit, Masse, Radius und der...
Entschuldigung für den langen Beitrag, aber ich wollte alles, was ich für relevant hielt, gleich mit einbeziehen. Was ich möchte Ich implementiere eine parallele Version der Krylov-Subraummethoden für dichte Matrizen. Hauptsächlich GMRES, QMR und CG. Mir wurde (nach der Profilerstellung) klar,...
Ich versuche einige Ergebnisse zu verstehen und würde mich über einige allgemeine Kommentare zum Umgang mit nichtlinearen Problemen freuen. Fisher-Gleichung (eine nichtlineare Reaktions-Diffusions-PDE), ut= dux x+ βu ( 1 - u ) = F( u )ut=duxx+βu(1-u)=F(u) u_t = du_{xx} + \beta u (1 - u) = F(u) in...