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...
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...
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...
max f( x ) vorbehaltlich A x = bmaxf(x) subject to Ax=b\max f(\mathbf{x}) \text{ subject to } \mathbf{Ax} = \mathbf{b} Dabei ist , und x=[x1,x2,. . . ,xN]T≤RN×1A≤RM×N.f( x ) = ∑N.i = 11 + x4ich( ∑N.i = 1x2ich)2- -- -- -- -- -- -- -- -- -- -√f(x)=∑i=1N1+xi4(∑i=1Nxi2)2f(\mathbf{x}) =...
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...
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...
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...
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...
In Matlab werden sowohl Linsolve als auch Mldivide zum Lösen eines linearen Gleichungssystems in allen bestimmten, überbestimmten und unterbestimmten Fällen verwendet. Als ich ihre Dokumente las, fragte ich mich, welche Unterschiede zwischen ihnen bestehen. Verwenden sie in den drei Fällen fast...
Geschlossen. Diese Frage ist nicht zum Thema . Derzeit werden keine Antworten akzeptiert. Möchten Sie diese Frage verbessern? Aktualisieren Sie die Frage so dass es beim Thema für Computational Science Stapel Börse. Geschlossen vor 2 Jahren . Wie der Titel schon sagt, möchte ich in der Lage sein,...
C ++ 11 führt eine Verschiebungssemantik ein, die beispielsweise die Codeleistung in Situationen verbessern kann, in denen C ++ 03 eine Kopierkonstruktion oder eine Kopierzuweisung durchführen müsste. Dieser Artikel berichtet, dass der folgende Code beim Kompilieren mit C + 11 eine 5-fache...
Wir alle kennen die vielen Berechnungsmethoden zur Lösung des linearen Standardsystems Ax=b.Ax=b. Ax=b. Ich bin jedoch gespannt, ob es "Standard" -Berechnungsmethoden zum Lösen eines allgemeineren (endlichdimensionalen) linearen Systems der Form gibt wobei beispielsweise A eine m 1 × n 1 -Matrix...
Wie alle wissen, gibt es in der Realität keinen inkompressiblen Fluss. Dies ist eine Annahme, die eingeführt wurde, um die Regelungsgleichungen zu vereinfachen. Wir können diese Annahme nicht einfach anwenden. Im Allgemeinen sind die Machzahl (M <0,3 für inkompressiblen Fluss), die...
Bei der Methode der hergestellten Lösungen (MMS) postuliert man eine exakte Lösung, setzt sie in die Gleichungen ein und berechnet den entsprechenden Quellterm. Die Lösung wird dann zur Codeüberprüfung verwendet. Für inkompressible Navier-Stokes-Gleichungen führt MMS leicht zu einem Quellterm...
Ich verwende das Crank-Nicolson-Finite-Differenzen-Schema, um eine 1D-Wärmegleichung zu lösen. Ich frage mich, ob das Maximum / Minimum-Prinzip der Wärmegleichung (dh das Maximum / Minimum tritt im Anfangszustand oder an den Grenzen auf) auch für die diskretisierte Lösung gilt. Dies wird...
Ich möchte einen Code portieren, der einen Satz partieller Differentialgleichungen (PDE) nach der Methode des endlichen Volumens in IMPLICIT-Form (für die Zeitdiskretisierung) auflöst. Als Ergebnis gibt es ein tridiagonales Gleichungssystem in x-, y- und z-Richtung, das vom ADI / TDMA-Schema...
Ich bin neu im wissenschaftlichen Rechnen. Ich suche eine Fortran (vorzugsweise f90) Implementierung eines Octree. Mein Problem erfordert einen Octree, der meine Domäne teilt, bis nicht mehr als einige N Partikel (oder Quellen, bei denen ich den Dichtewert kenne, der in eine
Ich versuche, die 2D-Poisson-Gleichung mit der diskontinuierlichen Galerkin-Methode (DG) und der folgenden Diskretisierung zu lösen (ich habe eine PNG-Datei, darf sie aber leider nicht hochladen): Gleichung: ∇⋅(κ∇T)+f=0∇⋅(κ∇T)+f=0\nabla \cdot( \kappa \nabla T) + f = 0 Neue Gleichungen:...
Ich verwende die RBF-Kernelfunktion, um einen kernelbasierten Algorithmus für maschinelles Lernen (KLPP) zu implementieren. Die resultierende Kernelmatrix ist extrem schlecht konditioniert. Die Bedingungsnummer der L2-Norm lautetKKK 1017-1064K(i,j)=exp(−(xi−xj)2σ2m)K(i,j)=exp(−(xi−xj)2σm2)K(i,j)=...
Zur Modellreduktion möchte ich die linken Singularvektoren berechnen, die den - sagen wir 20 - größten Singularwerten einer Matrix , wobei und . Leider wird meine Matrix ohne Struktur dicht sein. N ≤ 10 6 k ≤ 10 3 A.A∈RN,kA∈RN,kA \in \mathbb R^{N,k}N≈106N≈106N\approx 10^6k≈103k≈103k\approx 10^3AAA...