Ich habe versucht, Testfälle für CG auszuführen, und ich muss Folgendes generieren: symmetrische positive definitive Matrizen mit einer Größe> 10.000 VOLL DICHTE Verwenden Sie nur Matrixindizes und gegebenenfalls 1 Vektor (wie A(i,j)=x(i)−x(j)(i+j)A(i,j)=x(i)−x(j)(i+j)A(i,j) = \dfrac{x(i) -...
Ich muss nur eine kleine Anzahl niederfrequenter Fourier-Komponenten eines komplexen zweidimensionalen Arrays berechnen. Ich werde immer wieder dieselben Fourier-Komponenten berechnen, wenn sich das Eingabearray ändert. In der Grenze, in der ich nur eine Fourier-Komponente haben möchte, wäre es...
Ich habe CG in FORTRAN implementiert, indem ich es mit Intel MKL verknüpft habe. Wenn es Aussagen gibt wie: ( Siehe Wikipedia ) p=r; x=x+alpha*p r=r-alpha*Ap; oder ähnliche in QMR (in viel größerer Menge) v_tld = r; y = v_tld; rho = norm( y ); w_tld = r; z = w_tld; xi = norm( z ); (and more) Ist es...
Ich versuche MD mit Wasserdampf zu machen. Wie ich weiß, gibt es einige Wassermodelle für flüssiges Wasser, wie SPC, SPC / E, TIP3P, aber gelten sie auch für den Dampfzustand von Wasser? Und was ist der Unterschied zwischen Simulationen von flüssigem und dampfförmigem
Es scheint, dass Matlab's eigsmir schlechte Annäherungen an die kleinsten Eigenvektoren einer Matrix gibt. Ich gehe davon aus, dass ich einige langsamere Methoden verwenden kann, die auch genauer wären ... Ich suche nach dem zweitkleinsten Eigenvektor einer Lapalcian-Matrix (bekannt als "Fiedler"...
Ich erhalte die Warnung im Betreff des Beitrags, wenn ich versuche, eine Funktion in Python mit der Funktion scipy.optimize.fmin_bfgs zu optimieren . Die komplette Ausgabe: Warnung: Gewünschter Fehler wird aufgrund von Präzisionsverlust nicht unbedingt erreicht Current function value: nan...
Soweit ich weiß, stimmt die Reihenfolge, in der nicht blockierende Punkt-zu-Punkt-MPI-Nachrichten (Isend und Irecv) empfangen werden, mit der Reihenfolge überein, in der sie gesendet werden. Gibt es Techniken, um bestimmten Nachrichten Vorrang vor anderen zu geben? Zum Beispiel haben wir einen...
Ich habe eine hermitische Matrix die von zwei Parametern abhängt, z. B. x und y . Wenn ich es an zwei engen Punkten ( x 1 , y 1 ) und ( x 2 , y 2 ) diagonalisiere, erhalte ich zwei enge Eigenwerte ( ε 1 und ε 2 ) und zwei entsprechende Eigenräume ( S 1 und S 2 ) derselben Dimension...
Wenn ich daran denke, nichtlineare Gleichungen zu lösen, denke ich im Allgemeinen daran, zuerst zu linearisieren und dann einen Vorkonditionierer auf die lineare Matrix anzuwenden. Mir kam der Gedanke, dass es möglich sein könnte, die nichtlinearen Gleichungen zu konditionieren, bevor sie...
Jeder kennt Referenzen, die die Geschichte und die Ideen hinter Supercomputing zusammenfassen, einschließlich der Erwähnung von Entwicklungen in parallelen Programmiersprachen, Anwendungen, Startup-Unternehmen (einige wurden gegründet, aber nicht überlebt), Anforderungen der Industrie...
Meine Herausforderung besteht darin, das folgende Gleichungssystem zu lösen, das die Gasverbrennung in porösen Medien beschreibt: 1) Kontinuität ε ∂ρG∂t+ ∂∂x( ρGux) =0ε∂ρg∂t+∂∂x(ρgux)=0\varepsilon \frac{\partial \rho_g}{\partial t} +\frac{\partial}{\partial x} \left(\rho_g u_x\right)=0 2)...
Was wären gute Methoden (und / oder Softwarepakete), um ein Problem zu lösen, das eine quadratische Funktion minimiert? , st 0 ≤ x i ≤ 1 und Es gibt mehr Einschränkungen, von denen einige nicht linear (und nicht differenzierbar) sind, z. B. ∑ i x i 1 x i > a < b ?f( x ) = ∑N.i = 1( xich-...
Ich habe gelernt, dass es für den BiCGStab-Solver für spärliche lineare Systeme so gut wie immer erforderlich ist, einen Vorkonditionierer zu verwenden. Inzwischen wurde mir klar, dass die Auswahl eines guten problemabhängig ist. Beim Surfen im Internet habe ich herausgefunden, dass es viele...
Verschiedene numerische lineare Algebra- und numerische Optimierungsmethoden haben unterschiedliche Größenbereiche, in denen sie zusätzlich zu ihren eigenen Eigenschaften eine „gute Idee“ sind. Beispielsweise werden für sehr große Optimierungsprobleme Gradienten-, stochastische Gradienten- und...
Ich möchte die Lane-Emden-Isothermengleichung lösen [PDF, Gl. 15.2.9] d2ψdξ2+2ξdψdξ=e−ψd2ψdξ2+2ξdψdξ=e−ψ\frac{d^2 \!\psi}{d \xi^2} + \frac{2}{\xi} \frac{d \psi}{d \xi} = e^{-\psi} mit den Anfangsbedingungen ψ(ξ=0)=0dψdξ∣∣∣ξ=0=0ψ(ξ=0)=0dψdξ|ξ=0=0\psi(\xi = 0) = 0 \quad \left.\frac{d\psi}{d...
Ich bin neu auf dem Gebiet der CFD. Wann sollte man sich für ein strukturiertes Raster entscheiden und wann sollte man sich für ein unstrukturiertes Raster entscheiden? (Ja, es hängt stark von der Geometrie des Problems ab.) Insbesondere möchte ich den Unterschied in der erforderlichen...
So ziemlich die Frage. Wie wichtig ist bei einer allgemein spärlichen, nicht symmetrischen (sowohl numerisch als auch strukturell) Matrix das Sparsity-Muster (dh die Zeilen- / Spaltenpermutation von Matrix / Vektor) für iterative Löser? Ich kann sehen, dass es für Direktlöser (LU) oder...
Ich stehe vor dem Problem, das Ellipsoid ( ist eine symmetrische positive definitive Matrix) mit maximalem Volumen innerhalb einer konvexen Menge die als eine Menge linearer Ungleichungen . Ich habe verstanden, wie es als konvexes Optimierungsproblem formalisiert wird wie in "Convex Optimization,...
Ich bin neu in der linearen Programmierung und habe ein lineares Programm (LP) mit einer Reihenfolge von Variablen und Einschränkungen formuliert , obwohl die Einschränkungsmatrix extrem spärlich ist. 10 131013101310^{13}1013101310^{13} Ich wollte wissen, ob eine LP dieser Größenordnung handhabbar...
Ich interessiere mich für lokale AMR auf unstrukturierten Netzen. Derzeit arbeite ich mit der OpenFOAM-Bibliothek - sie unterstützt völlig unstrukturierte lokale AMR: Zellverfeinerungskriterien bestimmen eine Liste der Zellen, die geschnitten werden Ausgewählte Zellen werden verfeinert: Das gesamte...