Computational Science

Fragen und Antworten für Wissenschaftler, die Computer verwenden, um wissenschaftliche Probleme zu lösen

4
Derivate höherer Ordnung von Splines mit SciPy
Ich habe einen Spline erstellt, um meine Daten in Python anzupassen, indem ich: spline=scipy.interpolate.UnivariateSpline(energy, fpp, k=4) Die Gleichung, die ich verwenden möchte, beinhaltet eine Summation zwischen n = 2 und n = unendlich, wobei n die Ordnung des Differentials an einem Punkt Eo ist. Verwenden Sie jedoch; UnivariateSpline.__call__(spline, e0, nu=n) …
1
Was ist der richtige Weg, um Vektoren im Gleitkomma zu vergleichen?
Ich weiß, dass ich eine Toleranz für den Vergleich von Gleitkommazahlen verwenden sollte. Aber zum Vergleichen von Vektoren kann ich mir 3 mögliche Lösungen vorstellen, die unterschiedlichen Abstandsmetriken entsprechen: Vergleichen Sie die Komponenten jedes Vektors einzeln: Die Vektoren sind gleich, wenn alle 3 innerhalb der Toleranz liegen. Diese Option würde …
4
Wie kann ich einen parallelen linearen kongruenten Pseudozufallszahlengenerator für die maximale Periode setzen?
Normalerweise verwende ich den Aufruf, wenn ich einen sequentiellen Zufallszahlengenerator in C setze srand(time(NULL)) dann benutze rand() mod N um eine zufällige ganze Zahl zwischen 0 und N-1 zu erhalten. Wenn ich dies jedoch parallel mache, sind die Aufrufe zur Zeit (NULL) so nahe beieinander, dass sie genau dieselbe Nummer …
2
Wie berechne ich den parallelen Overhead eines parallelen Codes, der auf einem einzelnen Prozessor ausgeführt wird, wenn kein sequentieller Code verfügbar ist?
Ich stelle die Leistung der linearen Löser von PETSc vor. So wie ich es verstehe, speedup=Sequential TimeParallel Time.speedup=Sequential TimeParallel Time.\text{speedup}=\frac{\text{Sequential Time}}{\text{Parallel Time}}. Ich weiß, dass das Ausführen des parallelen Codes auf einem Prozessor als Proxy für die sequentielle Leistung verwendet werden kann. Ich denke jedoch nicht, dass dies aufgrund des …
2
Einseitige nichtlineare kleinste Quadrate mit linearen Einschränkungen
Ich versuche, ein einseitiges nichtlineares Problem der kleinsten Quadrate mit linearen Einschränkungen zu lösen, dh das Problem: minx∑mi=1ri(x) s.t Ax≤bminx∑i=1mri(x) s.t Ax≤b\min_{\mathbf{x}} \quad \sum^m_{i=1} \mathbf{r}_i(\mathbf{x}) \qquad \text{ s.t } \quad A\mathbf{x} \leq \mathbf{b} wo ri(x)=fi(x)2ri(x)=fi(x)2r_i(\mathbf{x})=f_i(\mathbf{x})^2 wenn fi(x)>0fi(x)>0f_i(\mathbf{x})>0 und ri(x)=0ri(x)=0r_i(\mathbf{x})=0 sonst. Mit anderen Worten, dies kann als Problem der kleinsten Quadrate …
2
Filtern eines Datensatzes, um eine gleichmäßigere Verteilung für das Training neuronaler Netze zu erhalten
Ich möchte künstliche neuronale Netze (ANN) verwenden, um die Reaktionsgeschwindigkeiten in meiner Flüssigkeit vorherzusagen, anstatt das gesamte System steifer ODEs zu lösen. Einige Leute aus meinem Labor haben bereits daran gearbeitet, damit ich nicht bei Null anfange, sondern Probleme mit meinen Anwendungen habe. Eine davon bezieht sich meiner Meinung nach …
1
Was ist ein gutes Stoppkriterium, wenn eine iterative Methode zum Ermitteln von Eigenwerten verwendet wird?
Ich habe diese Antwort gelesen und festgestellt, dass ich den Unterschied zwischen aufeinanderfolgenden Iterationen verwendet habe, um ein Stoppkriterium für eine iterative Methode zum Finden von Eigenwerten / Vektoren zu definieren. Was sind gute Stoppkriterien für iterative Methoden, die zu Eigenwerten und Eigenvektoren konvergieren?
3
Am besten lösen Sie ein nicht symmetrisches, nicht diagonal dominantes, spärliches System
Ich erinnere mich schwach an meine frühen "numerischen" Vorlesungen, dass iterative lineare Löser für oft erfordern, wenn zerlegt wird alsAx=bEINx=bAx=bAEINA A=D+MEIN=D.+M.A=D + M Wenn D eine Diagonalmatrix ist und eine Diagonale von Null hat, sollten die Elemente von über den Einträgen in dominieren, damit iterative Löser eine gute Leistung erbringen.MM.MDD.DMM.M …
1
Den Fixpunkt eines Operators finden
Welche numerischen Methoden stehen zur Verfügung, um den Fixpunkt eines Operators , der auf die Funktionen einwirkt ? Ich suche die Funktion für die .f : [ a , b ] → [ a , b ] f A f = fEINEINAf: [ a , b ] → [ a …
2
Was ist ein guter Weg, um die Gesamtstruktur einer Codebasis zu verstehen?
Manchmal ist es in meiner Arbeit nützlich, den Open-Source-Code eines anderen zu ändern oder herauszufinden, wie Sie bestimmte Dinge für Ihre eigene Anwendung entwickeln können. Allerdings verfügen nicht alle Softwareprogramme über eine gute Dokumentation. Was ist ein guter Weg, um die Gesamtstruktur einer Codebasis zu verstehen? Welche Routinen rufen beispielsweise …
2
Gibt es einen allgemeinen Ansatz, um Projektionsmethoden für verschiedene Probleme zu erstellen?
Meine Frage wird wahrscheinlich zu allgemein sein, um sie mit ein paar Worten zu beantworten. Könnten Sie bitte in diesem Fall eine gute Lektüre vorschlagen? Projektionsmethoden werden verwendet, um die Größe des Lösungsraums für die Probleme zu reduzieren. Und es gibt mindestens zwei sehr interessante Anwendungen (aus meiner Sicht). Das …
3
Was sollten die Kriterien sein, um einzelne Werte zu akzeptieren / abzulehnen?
Ich löse ein System mit Singularwertzerlegung. Die Singularwerte (vor der Skalierung) sind: 1.82277e+29 1.95011e+27 1.15033e+23 1.45291e+21 4.79336e+17 7.48116e+15 8.31087e+12 1.71838e+11 5.63232e+08 2.17863e+08 9.02783e+07 1.72345e+07 1.73889e+05 8.09382e+02 2.16644e+00 Ich habe festgestellt, dass das Akzeptieren aller Singularwerte und des damit verbundenen Beitrags zu meinem Lösungsvektor zu schlechten Ergebnissen führt. Ich skaliere sie …
Durch die Nutzung unserer Website bestätigen Sie, dass Sie unsere Cookie-Richtlinie und Datenschutzrichtlinie gelesen und verstanden haben.
Licensed under cc by-sa 3.0 with attribution required.