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...
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...
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...
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...
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...
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...
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)...
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...
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...
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...
In einem kürzlich erschienenen Artikel der SIAM News gibt es einen langen Artikel, der eine systematische Organisation der finiten Elemente beschreibt und treffend als Periodensystem der finiten Elemente bezeichnet wird . Es ist wirklich faszinierend zu sehen, wie die Klassifizierung über die...
Für einen albernen Bildschirmschoner, den ich entwickeln möchte, möchte ich zufällig ein divergenzfreies 2D-Array von 2D-Vektoren generieren und dann ein linienintegrales Faltungsdiagramm erstellen. Ich habe gehört , dass ein Weg , dies zu tun ist , zufälliges Rauschen zu erzeugen, und dann die...
Ich habe einen Vektor VVV der in den Eigenraum des hermitischen Sparse-Operators zerlegt werden kann :MMM V=∑ivim^iV=∑ivim^iV = \sum_i v_i \hat{m}_i Gibt es eine Möglichkeit, den (den Eigenvektor selbst) zu finden, der dem größten (in der Größe) entspricht?m^im^i\hat{m}_iviviv_i Ich möchte im...
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 ....
Ich habe versucht, zu den Primärquellen zu gehen, um zu verstehen, wie man Butcher-Tabellen verwendet, um die Algebra zu vereinfachen, die ich tun muss, wenn ich beispielsweise Taylor-Reihen verwende, um die Reihenfolge der Genauigkeit eines Schemas zu ermitteln. Vielleicht aufgrund des Mangels an...
Es ist ein ziemlich bekannter Trick, um eine Division bei der Berechnung von Quadratwurzeln zu vermeiden, um die Newtonsche Methode auf das Finden von anzuwenden1 / x- -- -√1/.x1/\sqrt{x} , und wahrscheinlich besser bekannt, wenn man die Newtonsche Methode verwendet, um Kehrwerte ohne Division zu...
Ich habe eine ziemlich spezifische Frage bezüglich der Bedingungsnummer. Ich führe FEM-Simulationen mit mehreren Längenskalen durch, was zu einer großen Ungleichheit zwischen den größten Einträgen und den kleinsten Einträgen in meiner Matrix führt. Die Bedingungsnummer kann unter bestimmten...
Gibt es eine Bibliothek / Toolbox, in der inkrementelle SVD in MATLAB implementiert ist? Ich habe dieses Papier implementiert , es ist schnell, funktioniert aber nicht gut. Ich habe dies versucht, aber auch hier breitet sich der Fehler schnell aus (beim Aktualisieren von 5-10 Punkten ist der Fehler...
Die Bewegungsgleichungen für einen elastischen Festkörper sind gegeben durch ∇⋅σ+f=ρu¨σ=Cεε=12(∇u+[∇u]T)∇⋅σ+f=ρu¨σ=Cεε=12(∇u+[∇u]T)\begin{align} &\nabla \cdot \boldsymbol{\sigma} + \mathbf{f} = \rho \ddot{\mathbf{u}}\\ &\boldsymbol{\sigma} = \mathbb{C}\boldsymbol{\varepsilon}\\...
Angenommen, wir haben ein Anfangswertproblem der Form wobei genau bekannt ist (dh mit unbegrenzter Genauigkeit) und wir effizient auswerten können mit beliebiger Genauigkeit. Das heißt, wir haben eine Blackbox, die bei gegebenem Vektor und einer ganzen Zahl eine Annäherung an zurückgibt, die...