Computational Science

2

Diagonalisierung der Matrix - Weglassen kleiner Matrixelemente

Ich habe mich gefragt, ob es einen Satz gibt, der es mir ermöglicht, eine Obergrenze für den Fehler festzulegen, der durch das Weglassen kleiner Matrixelemente aus einer Matrix vor der Diagonalisierung entsteht. Nehmen wir an, wir haben eine große Matrix, deren Matrixelemente zwischen 111 und . Wenn ich vor der …

8 eigensystem error-estimation

2

Methode zur Quantifizierung der geometrischen Differenz zweier unterschiedlicher Netze

Ich suche nach einer Methode oder einem Algorithmus, um einen Wert zu erzeugen, der beschreibt, wie unterschiedlich zwei Maschen geometrisch sind, aber unterschiedliche Topologien haben. Ein Beispiel wären einige CAD-Daten, bei denen ein FE-Netz angewendet wurde, und das Teil, nachdem es hergestellt und dann in ein dreieckiges 3D-Netz gescannt wurde. …

8 optimization finite-element

6

Ich suche eine kostenlose symbolische Regressionssoftware

Kennt Formulize / Eureqa, nachdem er damit begonnen hat, 2500 US-Dollar pro Jahr für die Verwendung zu berechnen und die Testversion verkrüppelt zu haben, irgendwelche Ersetzungen, die ähnliche Dinge tun können, wie eine Gleichung mit bestimmten Datenpunkten zu finden? Es war kostenlos, liegt aber jetzt außerhalb der Preisspanne von Personen, …

8 machine-learning regression symbolic-computation

1

Intuition hinter der Methode der Multiplikatoren mit alternierender Richtung

Ich habe in letzter Zeit viele Artikel über ADMM gelesen und auch versucht, einige Probleme damit zu lösen, bei denen es sehr effektiv war. Im Gegensatz zu anderen Optimierungsmethoden kann ich mir nicht vorstellen, wie und warum diese Methode so effektiv ist (natürlich habe ich in einigen Fällen eine Konvergenzanalyse …

8 optimization admm

1

Gibt es eine verbesserte Methode zur Berechnung des folgenden Ausdrucks?

gegeben eine symmetrische Matrix und eine beliebige Matrix und einen Vektor , ist es möglich, den folgenden Ausdruck in -Zeit zu berechnen ?Y∈Rn×nY∈Rn×nY \in \mathbb{R}^{n \times n}X∈Rn×nX∈Rn×nX \in \mathbb{R}^{n \times n}v∈Rn×1v∈Rn×1v \in \mathbb{R}^{n \times 1}O(n2)O(n2)O(n^2) diag(XTYX)⋅vdiag(XTYX)⋅vdiag(X^TYX) \cdot v wobei eine Matrix zurückgibt, deren Hauptdiagonalelemente gleich denen von und nicht diagonalen …

8 optimization algorithms performance matrix complexity

2

Ist die Aufnahmemethode die einzige allgemeine numerische Methode zur Lösung nichtlinearer Randwert-ODEs?

Während meiner Wanderung in Mathematica.se bemerkte ich allmählich, dass eine bestimmte Art von Problem der Lösung von Differentialgleichungen uns ständig "beunruhigt", dh das Randwertproblem (BVP) nichtlinearer gewöhnlicher Differentialgleichungen (ODEs). Die Aufnahmemethode , die von der Mathematica- Funktion verwendet wirdNDSolve , scheint die einzige Methode zu sein, die Benutzer von Mathematica.SE …

8 ode nonlinear-equations

2

Soll ich MPI_Init Befehlszeilenargumente übergeben oder nicht?

Beim Schreiben von MPI 3.0 - Code, soll ich weitergeben argcund argvauf den MPI_InitAnruf oder nicht, und warum? EDIT: Da die Antwort von Bill Barth einige Fragen aufwirft, möchte ich einige Anmerkungen machen: Bestehen argc/ argvist seit MPI 1.1 nicht erforderlich. Die Frage ist speziell, warum sollte / sollte man …

8 c++ c mpi

1

Finite-Elemente-Konvergenzraten für gemischte Probleme

Ich habe ein Stokes Flow-Problem mit finiten Elementen codiert und bin dabei zu überprüfen, ob es funktioniert. Ich bin mir nur nicht sicher, welche Konvergenzrate ich erwarten sollte, wenn ich das Netz global verfeinere. Ich weiß, dass ich für skalare Probleme mit linearen Basisfunktionen eine Konvergenz der Ordnung h2h2h^2 ( …

8 pde finite-element convergence

1

Veröffentlichung eines Softwarepakets zusammen mit einem Papier: Lizenzierung

Ich möchte das Softwarepaket, das ich für meine Abschlussarbeit geschrieben habe, zusammen mit einem Dokument veröffentlichen, das das Paket beschreibt. Bisher habe ich über liberale Lizenzen wie BSD nachgedacht. Mir ist jedoch jetzt bewusst, dass die Veröffentlichung Ihres Codes eine Vereinbarung mit den Lizenzbedingungen des Journals (cmp. ZB diese Frage …

8 publications

2

Planen und Einbinden von Projektmanagement-Tools in den Forschungscode

Ich gehe in mein drittes Jahr meines Doktorandenprogramms und bis jetzt besteht mein Forschungscode (numerische inverse Probleme / medizinische Bildgebung / Bildverarbeitung / usw.) hauptsächlich aus unorganisierten MATLAB-Skripten und -Funktionen mit einigen MEX'd C ++ und CUDA Routinen, die hier und da für zusätzliche Geschwindigkeit sorgen. Die Arbeit mit MATLAB …

8 open-source

2

Periodische Greensche Funktionen in Integralgleichungsmethoden in verschiedenen Frequenzbereichen

Ich frage nach der Lösung der Helmholtz-Gleichung in einem periodischen Bereich mit stückweise konstanter Wellengeschwindigkeit in verschiedenen Frequenzbereichen. Ein möglicher Ansatz zur Lösung dieses Problems besteht darin, Integralgleichungen auf die Grenzflächen in Bezug auf die Greensche Funktion des Systems aufzuschreiben. Da die Domäne periodisch ist, ist dies eine periodische Greensche …

8 integral-equations helmholtz-equation

1

Effiziente Implementierung von Gravitationsfeldern

Ich habe eine ähnliche Frage zu physics.stackexchange gestellt , da ich diese Website nicht kenne. Ich suche grundsätzlich nach einem effizienten Weg, um Gravitationsfelder zu implementieren. Ich habe einen riesigen 2D-Raum mit Tausenden von Objekten. Ich muss dann simulieren, wie diese Objekte durch die Schwerkraft des anderen beeinflusst werden. Ich …

8 optimization simulation approximation

4

Sollte ich eine Datenbank verwenden, um große Mengen an Ergebnissen zu verarbeiten?

Hintergrund: Ich führe derzeit eine große Anzahl von Parametervariationsexperimenten durch. Sie werden in Python 2.6+ mit numpy ausgeführt. Diese Experimente werden ungefähr 2 Wochen dauern. Ich variiere ungefähr 3 Parameter (unabhängige Variablen) über einen Wertebereich. Ich behebe 6 weitere unabhängige Variablen (vorerst) Ich berichte über 4 abhängige Variablen. Einer der …

8 software visualization data-analysis

1

Radiale Integration der teuren Funktion mit Bessel-Gewichten

Ich muss das Integral I=∫R0f(r)Jn(znmrR)rdrI=∫0Rf(r)Jn(znmrR)rdrI = \int_0^R f(r)J_n\left(\frac{z_{nm}r}{R}\right)rdr wobei die Bessel-Funktionen der ersten Art der Ordnung , ihre Null ist und eine reelle Funktion ist, die etwas ähnlich ist zu (aber nicht dasselbe, es ist ziemlich kompliziert und beinhaltet normalerweise Begriffe mit und manchmal ).n t h z n m …

8 quadrature integration special-functions numerics

1

GUT Pseudozufallszahlengenerationen

Ich habe MT19937in einem Test-Harness einheitliche (vorzeichenlose) 32-Bit- Werte [0, - 1] basierend auf der ursprünglichen mt19937.c- Implementierung der Autoren generiert , um ein (im Wesentlichen unerschöpfliches) Angebot an zu generieren statistisch zufälliger Strom von Bitoktetten. Dies ersetzt ein CSPRNG, das für diese speziellen Tests nicht erforderlich ist. Ich habe …

8 performance c random-number-generation