Ich habe einen Datensatz, der sich langsam ändert, und ich muss die Eigenvektoren / Eigenwerte seiner Kovarianzmatrix verfolgen. Ich habe es benutzt scipy.linalg.eigh, aber es ist zu teuer und es nutzt nicht die Tatsache, dass ich bereits eine Zerlegung habe, die nur geringfügig falsch ist. Kann...
Ich habe die SVD von Intel MKL ( dgesvdüber SciPy) verwendet und festgestellt, dass die Ergebnisse erheblich voneinander abweichen, wenn ich die Genauigkeit zwischen float32und float64wenn meine Matrix schlecht konditioniert ist / nicht den vollen Rang hat. Gibt es einen Leitfaden zum Mindestmaß an...
In der Veröffentlichung http://www.sciencedirect.com/science/article/pii/S0045782509003521 wird auf Seite 584 Gleichung (4) eine HDG-Element-lokale Gleichung beschrieben, wobei eine der Gleichungen die folgende Form annimmt −(uh,∇q)K=−⟨u^h⋅n,q−q¯⟩∂K−(uh,∇q)K=−⟨u^h⋅n,q−q¯⟩∂K-(u_h,\nabla q)_K =...
Bei der Erforschung inverser Probleme ist es üblich, aus einem bekannten Parametersatz einen synthetischen Datensatz zu erstellen und dann zu testen, ob die Inversionstechnik diese Parameter rekonstruieren kann. Dabei ist es wichtig, den synthetischen Daten ein angemessenes Maß an zufälligem...
Sei reelle, quadratische, dichte Matrizen. G und Q sind symmetrisch. LassenA , G , Q.A,G,QA, G, QGGGQ.QQ H.= [ A.- Q.- G.- A.T.]]H=[A−G−Q−AT]H = \begin{bmatrix} A & -G \\ -Q &-A^T \end{bmatrix} sei eine Hamiltonsche Matrix. Ich möchte die Matrix Exponential von berechnen . Ich brauche das...
Jede reelle Matrix kann unter Verwendung einer orthogonalen Ähnlichkeitstransformation U auf die reelle Schurform T = U T A U reduziert werden . Hier ist die Matrix T quasi dreieckig mit 1 mal 1 oder 2 mal 2 Blöcken auf der Hauptdiagonale. Jeweils 1 x 1 Block entspricht einem realen Eigenwert von A...
Bestehende Algorithmen zum Lösen von ODEs verarbeiten Funktionen , wobei . In vielen physikalischen Systemen ist die Differentialgleichung jedoch autonom, so dass , , wobei das weggelassen wird. Welche Verbesserungen können mit dieser vereinfachenden Annahme bei bestehenden numerischen Methoden...
Ich habe festgestellt, dass ein Großteil meiner rechnergestützten Programmierung Testanforderungen hat, die nicht durch Standard-Test-Frameworks abgedeckt sind: Berechnung der Rechenzeit Um sicherzustellen, dass Algorithmen nicht langsamer werden. Ich könnte so etwas tun,
Ich löse komprimierbare Navier-Stokes-Gleichungen für mehrere Arten in einem strukturierten 3D-Gitter. Ich habe eine Lösung für ein bestimmtes Gitter erhalten (sagen wir eine relativ grobe). Ich möchte jetzt mein Raster verfeinern und meine vorherige Lösung auf meinem neuen Raster interpolieren,...
Gibt es gute Papiere und / oder Codes, die diskontinuierliche Galerkin-Finite-Elemente-Löser mit Riemann-Lösern koppeln? Ich muss die elliptischen und hyperbolischen Kopplungsprobleme untersuchen, aber die meisten Aufteilungsmethoden sind bestenfalls ad hoc. Da ich eine große Menge FEniCS-Code...
Wann immer man eine Quasi-Monte-Carlo-Methode zur Kubatur oder Optimierung verwendet, scheint es eine Vielzahl von Sequenzen mit geringer Diskrepanz zu geben, die mit den Namen van der Corput, Halton, Hammersley, Faure, Niederreiter, Sobol ', assoziiert sind. und andere Namen, an die ich mich nicht...
Bei der Lösung spärlicher linearer Systeme mit direkten Faktorisierungsmethoden wirkt sich die verwendete Ordnungsstrategie erheblich auf den Füllfaktor von Nicht-Null-Elementen in den Faktoren aus. Eine solche Bestellstrategie ist die verschachtelte Dissektion. Ich frage mich, ob es möglich ist,...
Ich versuche mich zu integrieren ∫10t2 n + 2exp(α r0t) dt∫01t2n+2exp(αr0t)dt\int^1_0 t^{2n+2}\exp\left({\frac{\alpha r_0}{t}}\right)dt Das ist eine einfache Transformation von ∫∞1x2 nexp( - α r0x ) dx∫1∞x2nexp(- -αr0x)dx\int^{\infty}_1 x^{2n}\exp(-\alpha r_0 x)dx mit weil es schwierig ist,...
Zum Beispiel möchte ich die -Norm von u = 1 numerisch berechnenL2L2L^2 ingewissen Domänedie Null enthält, habe ich versuchtGauss Quadratur und es scheitert, ist es irgendwie weit von der realenL2-Norm auf der Einheitskugel Kugelkoordinaten verwendet, zu integrieren, Gibt es einen guten Weg, dies zu...
Diese Frage ergibt sich aus der Notwendigkeit, eine Lektion über die Grenzen der Dichtefunktionaltheorie als rechnerischen Ansatz vorzubereiten. Ich möchte nicht nur die Einschränkungen kennen, sondern auch Referenztexte, mit denen ich eine ansprechendere und vollständigere Lektion vorbereiten...
Unbestimmte Matrizensysteme treten beispielsweise bei der Diskretisierung von Sattelpunktproblemen durch gemischte finite Elemente auf. Die Systemmatrix kann dann in das Formular eingefügt werden ( A.B.B.tC.)(ABtBC)\begin{pmatrix} A & B^t \\ B & C\end{pmatrix} wobei negativ (semi) -definit ist,...
Ich habe gehört, dass einige Zeitschriften höher bewertet werden als andere. Ist das wahr? Und wenn ja, nach welchen Kriterien kann der Wert eines Peer-Review-Journals gegenüber einem anderen beurteilt werden? Wie finde ich die Bewertung heraus? Wird meine Veröffentlichung weniger "wert" sein, wenn...
Ich habe gehört, dass eine schnelle Fourier-Transformation verwendet werden kann, um das Poisson-Problem zu lösen, wenn die Randbedingungen alle ein Typ sind ... Sinusreihen für Dirichlet, Cosinus für Neumann und beide für periodische. Angenommen, zwei gegenüberliegende Seiten haben periodische...
Ich habe einige Fragen zu folgenden Themen: Ich versuche, die Schrödinger-Gleichung in 1D unter Verwendung der Kurbel-Nicolson-Diskretisierung zu lösen, gefolgt von der Invertierung der resultierenden tridiagonalen Matrix. Mein Problem hat sich nun zu einem Problem mit periodischen Randbedingungen...
Der Titel sagt das meiste davon. Ich suche eine leichte und benutzerfreundliche Bibliothek, die ich für Android (NDK) -Projekte verwenden kann. Für dichtes Material verwende ich gerne Eigen, aber ich habe nicht viele umfassende (und dokumentierte!) Bibliotheken für spärliches Material gefunden, die...