Computational Science

2

C ++ - Bibliothek zur nichtlinearen eingeschränkten Minimierung

Ich versuche derzeit, das nichtlineare Problem der eingeschränkten Minimierung zu lösen, wie es in der matlab-Funktion "fmincon" implementiert ist. Meine Erwartungen sind: Minimieren (fun1, x0, uB, lB, fun2), wobei x0 der Anfangszustand ist, fun1 eine Funktion ist, die minimiert werden muss, uB obere Grenzen sind, lB untere Grenzen sind und …

9 c++ nonlinear-programming

1

Hybride räumliche Schemata für CFD: Gibt es Nachteile beim Mischen gegenüber dem Umschalten?

Gibt es neben den zusätzlichen Rechenkosten, die durch die Berechnung beider Flüsse über einen bestimmten Bereich entstehen, einen Nachteil, wenn zwei Flussbewertungen für ein Hybridschema in einer Methode mit endlichem Volumen gemischt werden müssen? Die Flussbewertung würde folgendermaßen aussehen: F.i + 12= Λi + 12F.ci + 12+ ( 1 - …

9 fluid-dynamics finite-volume

4

Wenn wir Bernstein-Polynome in der Anwendung verwenden

Wenn es vorzuziehen ist, Bernstein-Polynome zu verwenden, um eine stetige Funktion zu approximieren, anstatt die einzigen folgenden vorläufigen numerischen Analysemethoden zu verwenden: "Lagrange-Polynome", "Einfache Operatoren für endliche Differenzen". Die Frage ist, diese Methode zu vergleichen.

9 finite-element finite-difference interpolation

6

Weiche Frage: Wo passt Python ins Bild?

Ich habe überlegt, ob ich Python lernen soll oder nicht. Aus dem Gespräch mit meinen Professoren scheint Matlab die gemeinsame Sprache in der angewandten Mathematik / Computational Science verwendet werden, soweit die Wissenschaft angeht; Während meiner Zeit in der Industrie haben meine Professoren (insbesondere diejenigen, die in der Industrie gearbeitet …

9 matlab python c++

1

Matlab Pde Toolbox: Zeichnen Sie die Lösung in einer Linie oder auf einer Untervielfalt

Ich verwende die Matlab pde-Toolbox, um eine bestimmte elliptische Gleichung in 2D zu lösen. Die Lösung ist in Ordnung, obwohl ich sie entlang einer bestimmten Linie zeichnen muss, dh um eine planare Schicht aus dem 3D-Netz zu schneiden, das die Lösung darstellt. Ich kann nicht herausfinden, wie die Toolbox-Funktionen intelligent …

9 pde matlab

4

Hosting-Site für eine kleine wissenschaftliche Bibliothek

Für meine Forschungsarbeit habe ich eine kleine C ++ - Bibliothek entwickelt, die die Kommunikation zwischen C ++ - Rechencodes und Octave / Matlab erleichtern soll (wenn letzterer für Nachbearbeitungszwecke verwendet wird). Ich möchte eine solche Bibliothek unter GPL auf einer kostenlosen Hosting-Site wie Sourceforge oder Github veröffentlichen. Ehrlich gesagt …

9 matlab libraries c++

2

Wie entferne ich starre Körperbewegungen bei linearer Elastizität?

Ich möchte lösen, wobei meine Steifheitsmatrix ist. Es können jedoch einige Einschränkungen fehlen und daher kann immer noch eine Bewegung des starren Körpers im System vorhanden sein (aufgrund des Eigenwerts Null). Da ich CG zum Lösen des linearen Systems verwende, ist dies nicht akzeptabel, da CG manchmal nicht bei halbpositiven …

9 finite-element iterative-method conjugate-gradient

3

Schnellster Algorithmus zur Berechnung der Bedingungsnummer einer großen Matrix in Matlab / Octave

Aus der Definition der Bedingungsnummer geht hervor, dass eine Matrixinversion erforderlich ist, um sie zu berechnen. Ich frage mich, ob für eine generische quadratische Matrix (oder besser, wenn eine symmetrische positive Bestimmtheit möglich ist) eine Matrixzerlegung ausgenutzt werden kann, um die Bedingungsnummer in a zu berechnen schneller Weg.

9 linear-algebra matrix condition-number

4

Fourier-Pseudospektralmethode und numerische Dissipation

Durchführen einer direkten numerischen Simulation isotroper Turbulenzen mit der Fourier-Pseudospektralmethode (Orzag & Patterson, PRL, 1972) unter Verwendung von FFT. In diesem Kurs finden Sie Hintergrundinformationen zu der in der Turbulenz-Community weit verbreiteten Methode: http://www.math.ualberta.ca/~bowman/m655/lab3d.pdf Unter Verwendung der so genannten - Regeln für dealiasing die besteht aus tun u ( k …

9 fluid-dynamics simulation

1

Schnelle (ungefähre) Auswertung des Chebyshev-Polynoms

Gibt es eine bevorzugte Methode zur Implementierung einer schnellen (ungefähren) Auswertung des Chebyshev-Interpolationspolynoms auf einem einheitlichen Gitter (angesichts der Funktionswerte an den Chebyshev-Knoten)? Mein Problem ist, dass die Interpolation langsam wird, wenn der Grad des Interpolationspolynoms zunimmt. Die folgenden Ideen kamen mir in den Sinn: Versuchen Sie, ungleichmäßige FFT-Techniken (NFFT) …

9 interpolation fourier-analysis fftw

1

Empfindlichkeit von BFGS gegenüber anfänglichen hessischen Näherungen

Ich versuche, die Broyden-Fletcher-Goldfarb-Shanno-Methode zu implementieren, um das Minimum einer Funktion zu finden. Ich brauche zwei anfängliche Vermutungen & und eine anfängliche hessische Matrixnäherung . Die einzigen Anforderungen, die ich für finde, sind, dass wenn der Hessische symmetrisch positiv definit ist, auch . Wenn ich auf Wikipedia schaue, sehe ich, …

9 optimization algorithms

3

Empfehlungen für einen leichten / nicht installierten C- oder C ++ - dichten linearen Algebra-Löser

Der größte Teil meiner Programmierung besteht aus einmaligen Forschungscodes in C für meinen eigenen Gebrauch. Ich habe niemals Code an andere als enge Mitarbeiter verteilt. Ich habe einen Algorithmus entwickelt, den ich in einer wissenschaftlichen Zeitschrift veröffentliche. Ich möchte den Quellcode und möglicherweise den ausführbaren Code in der Online-Ergänzung zum …

9 linear-solver libraries c++ c

2

Der schnellste Weg, um Eigenpaare einer kleinen unsymmetrischen Matrix auf einer GPU im gemeinsamen Speicher zu finden

Ich habe ein Problem, bei dem ich alle positiven (wie im Eigenwert positiv ist) Eigenpaare einer kleinen (normalerweise kleiner als 60x60) unsymmetrischen Matrix finden muss. Ich kann aufhören zu berechnen, wenn der Eigenwert kleiner als ein bestimmter Schwellenwert ist. Ich weiß, dass die Eigenwerte real sind. Irgendwelche Vorschläge zu Algorithmen, …

9 performance eigensystem gpu

3

Wie kann man intelligent versuchen, Konvexität auszuschließen?

Ich möchte eine komplizierte Zielfunktion minimieren und bin mir nicht sicher, ob sie konvex ist. Gibt es einen netten Algorithmus, der zu beweisen versucht, dass er nicht konvex ist? Natürlich könnte der Algorithmus dies nicht beweisen. In diesem Fall würde ich nicht wissen, ob er konvex ist oder nicht, und …

9 convex-optimization

2

Gleitkomma-Arithmetik mit höherer Genauigkeit in numerischer PDE

Ich habe aus sehr unterschiedlichen Quellen und Gesprächen mit Forschern den Eindruck, dass eine wachsende Nachfrage nach hochpräzisen Berechnungen in numerischen partiellen Differentialgleichungen besteht. Hohe Präzision bedeutet hier mehr Präzision als nur die standardmäßige 64-Bit-Doppelpräzision. Ich frage mich über den Stand der Technik zu diesem Thema. Zum Vergleich gibt es …

9 reference-request precision