Ich habe eine Sammlung von Rechenmodellen, die als asynchrone zellulare Automaten beschrieben werden könnten. Diese Modelle ähneln dem Ising-Modell, sind jedoch etwas komplizierter. Es scheint, als würden solche Modelle davon profitieren, auf einer GPU anstatt auf einer CPU betrieben zu werden....
Ich portiere einen vorhandenen Code von MATLAB nach C ++ und habe ein lineares System zum Lösen von (anstelle der typischeren Form A x = b )xA=bxA=bxA=bAx=bAx=bAx=b Die Matrix ist dicht und von allgemeiner Form, aber nicht größer als 1000 × 1000. In MATLAB wird die Lösung also durch die Funktion...
Ich simuliere einen inkompressiblen Fluss über einen Zylinder mit einer Reynold-Zahl von 500. Ich löse die Navier-Stokes-Gleichung mithilfe der Druckkorrekturmethode. Meine Lösung wird nach einer gewissen Zeit (ca. 5s) instabil. Ich habe versucht, mein Netz in Schritten von 0,05 zu verfeinern...
Für die numerische Quadratur gibt es mehrere Methodenfamilien. Wenn ich eine bestimmte Klasse von Integranden habe, wie wähle ich die ideale Methode aus? Was sind die relevanten Fragen, die sowohl zum Integranden (z. B. ist er glatt? Hat er Singularitäten?) Als auch zum Rechenproblem (z. B....
Angenommen, ich hatte ein Randwertproblem: dud2udx2+ dvdx= f in Ωd2udx2+dvdx=f im Ω\frac{d^2u}{dx^2} + \frac{dv}{dx}=f \text{ in } \Omega u=h in ∂Ωdudx+ d2vdx2= g in Ωdudx+d2vdx2=G im Ω\frac{du}{dx} +\frac{d^2v}{dx^2} =g \text{ in } \Omega u = h in ∂Ωu=h im ∂Ωu=h \text{ in }...
Bei der FEM-Diskretisierung und Lösung eines Reaktionsdiffusionsproblems, zB mit (singuläre Störung) zeigt die Lösung des diskreten Problems typischerweise Oszillationsschichten nahe der Grenze. Mit , und linearen finiten Elementen sieht die Lösung aus- & egr; & Dgr; u + u = 1 auf Ωu =...
Mein System ist ein symmetrisches FE-Problem mit Lagrange-Multiplikatoren (z. B. inkompressibler Stokes-Fluss): ( ABBTC)(EINBTBC)\begin{pmatrix}A & B^T \\ B & C\end{pmatrix} Dabei ist der typische Fall (ich habe sogar darauf geachtet, dass die Gleichungen nummeriert sind, damit die...
Hey, und danke, dass du dir Zeit genommen hast, meine Frage anzuschauen. Dies ist eine aktualisierte Version meiner Frage, die ich zuvor unter physics.stackexchange.com gepostet habe Ich studiere derzeit ein 2D-Excitonen-Spinor-Bose-Einstein-Kondensat und bin gespannt auf den Grundzustand dieses...
Ich versuche herauszufinden, ob es einen schnelleren Weg gibt, alle Eigenwerte und Eigenvektoren einer sehr großen und spärlichen Adjazenzmatrix zu berechnen, als mit scipy.sparse.linalg.eigsh Soweit ich weiß, verwendet diese Methode nur die Spärlichkeit und Symmetrieattribute der Matrix. Eine...
In den inkompressiblen Navier-Stokes-Gleichungen Der Druckbegriff wird oft als Lagrange-Multiplikator bezeichnet, der die Inkompressibilitätsbedingung.ρ ( ut+ ( u ⋅ ∇ ) u )Weiterempfehlen ∇ ⋅ u= - ∇ p + & mgr; & Dgr; u + f= 0ρ(ut+(u⋅∇)u)=-∇p+μΔu+f∇⋅u=0\begin{align*} \rho\left(\mathbf{u}_t +...
Ich bin mit dem in HDF5 verwendeten Dateiformat nicht vertraut, frage mich jedoch, ob HDF5-Dateien für die Versionskontrolle mit Git (oder zum Beispiel Mercurial oder Subversion) geeignet sind. Ich denke, was ich meine ist: Sind HDF5-Dateien für zeilenbasiertes Vergleichen geeignet oder muss git...
Ich bin daran interessiert, eine Funktion vieler ( ) realer Parameter (ein Ergebnis einer komplexen Simulation) global zu maximieren . Die Bewertung der betreffenden Funktion ist jedoch relativ teuer und erfordert für jeden Parametersatz etwa 2 Tage. Ich vergleiche verschiedene Optionen und habe...
Wir arbeiten an einem Bayes'schen Modell für einen Raum-Zeit-Prozess und verwenden einen No-U-Turn-Sampler (NUTS), für den ein Modell für die Log-Wahrscheinlichkeit und den Gradienten in Bezug auf die Modellparameter erforderlich ist. Um es kurz zu machen, wir haben eine ziemlich komplizierte...
Sei und f ( → x ) : [ 0 , 1 ] n → C sei eine Funktion in diesen Variablen.x⃗ =(x1,x2,…,xn)∈[0,1]nx→=(x1,x2,…,xn)∈[0,1]n\vec{x} = (x_1, x_2, \dots, x_n) \in [0,1]^nf(x⃗ ):[0,1]n→Cf(x→):[0,1]n→Cf(\vec{x}): [0,1]^n \to \mathbb{C} Gibt es ein rekursives Schema für dieses iterierte...
Angenommen, ich habe langsame und schnelle Versionen eines Codes und möchte eine Beschleunigungszahl melden, die die beiden vergleicht. Ich führe die langsame Version mal und die schnelle Version m- mal aus und produziere die Zeiten ( s 1 , … , s n ) und ( f 1 , … , f m ) . Der einfachste Weg, eine...
Dies ist möglicherweise eine Frage, die sich besser für die Software-Empfehlungsseite von SE eignet. Ich glaube jedoch, dass Personen, die diesen Teil von SE häufig verwenden, diese Frage mit größerer Wahrscheinlichkeit beantworten können. Ich suche eine kostenlose (nicht nur in Freiheit)...
Wenn ich ein relativ kleines Problem lösen würde, das heißt ein Problem, das mit einer direkten Methode wie LU behandelt werden kann, hat dann die Bedingungsnummer des linearen Operators Einfluss auf die Genauigkeit der Lösung? Eines der Forschungsprobleme, an denen ich arbeite, konzentriert sich...
Ich habe gerade Kapitel 3 in "A Multigrid Tutorial" von Briggs / Henson / McCormick gelesen, Link . Der Text handelt von Multigrid-Zyklen wie V-Zyklus, Mu-Zyklus, FMG. Was mir aufgefallen ist: Bei den meisten iterativen Verfahren wird geprüft, ob die gewünschte Toleranz / Genauigkeit erreicht...
Derzeit arbeite ich mit "A Multigrid Tutorial" von Briggs et al., Kapitel 8. Der Aufbau des Interpolationsoperators ist gegeben als: Der Aufbau des Beschränkungsbetreibers und des Feinnetzbetreibers ist dann gegeben als: Nehmen wir an, wir haben drei Gitterpunkte x0, x1, x2, wobei der mittlere x1...
Ich arbeite an einem inversen Problem für meinen Ph.D. Forschung, die der Einfachheit halber in bestimmtββ\beta L ( β) u ≡ - ∇ ⋅ ( k0eβ∇ u ) = fL(β)u≡-∇⋅(k0eβ∇u)=fL(\beta)u \equiv -\nabla\cdot(k_0e^\beta\nabla u) = f nach einigen Beobachtungen ; k 0 ist eine Konstante und f ist bekannt. Dies wird...