Ich frage mich, ob es eine sehr gute Mailingliste oder Google-Gruppen für Computational Science gibt, in denen wir Fragen diskutieren können, anstatt nur Fragen zu stellen und zu beantworten. Tatsächlich interessiere ich mich mehr für paralleles Rechnen und numerische Lösungen von PDEs. Aber ich...
Ich habe mich gefragt, wie ich am besten lernen kann, wie man in Fortran programmiert. Ich habe nach Büchern gesucht. Die Beispiele, die mir in Ordnung erschienen, waren "Fortran 95/2003 Explained" und "Fortran For Scientists And Engineers" . Ich weiß nicht, was ihr darüber denkt? Eine andere...
Ich schreibe ein kleines MATLAB-Paket, das bestimmte numerische Klassenprobleme löst. Es gibt 3 Stufen des Algorithmus und der Benutzer hat 5 Möglichkeiten für jede Stufe. Ich habe das gesamte Problem mit Funktionen und 3 Schalterfällen implementiert (einer für jede Algorithmusstufe). Es...
Ich habe in den letzten Tagen versucht, diesen Fehler zu beheben. Ich habe mich gefragt, ob jemand Ratschläge zum weiteren Vorgehen hat. Ich löse die Poisson-Gleichung für eine Stufenladungsverteilung (ein häufiges Problem in der Elektrostatik / Halbleiterphysik) auf einem ungleichmäßigen...
Betrachten Sie das streng konvexe uneingeschränkte OptimierungsproblemLassen Sie x_ \ text {opt} seine einzigartige Minima bezeichnen und X_0 eine gegebene erste Annäherung sein x_ \ text {opt}. Wir werden einen Vektor x eine \ epsilon- nahe Lösung von \ mathcal {O} nennen, wenn \ begin...
Diese Fragen sind in comp-sci möglicherweise etwas unangebracht. Wenn es benötigt wird, schlagen Sie bitte vor, wo es passt. Die Frage ist, wie alle Simulationsläufe effizient verwaltet werden können. Nehmen wir zum Beispiel an, eine Simulation erfordert die Festlegung von 2 Parametern, die in...
Ich habe ein großes kubisches Eigenwertproblem: ( A.0+ λ A.1+ λ2EIN2+ λ3EIN3) x =0.(EIN0+λEIN1+λ2EIN2+λ3EIN3)x=0.\left(\mathbf{A}_0 + \lambda\mathbf{A}_1 + \lambda^2\mathbf{A}_2 + \lambda^3\mathbf{A}_3\right)\mathbf{x} = 0. Ich könnte dies lösen, indem ich in ein lineares Eigenwertproblem...
Was sind die Hauptunterschiede zwischen FEM und XFEM? Wann sollten wir XFEM anstelle von FEM (nicht) verwenden und umgekehrt? Mit anderen Worten, wenn ich auf ein neues Problem stoße, wie kann ich dann feststellen, welches davon verwendet
Gibt es eine einfache Faustregel, um zu sagen, ob es sich lohnt, SOR anstelle von Gauß-Seidel durchzuführen? (und mögliche Methode zur Schätzung des Realxationsparameters )ωω\omega Ich meine, nur indem ich auf die Matrix schaue oder das Wissen über ein bestimmtes Problem, das die Matrix darstellt?...
Viele adaptive FEM-Bibliotheken verwenden erweiterte Netzdatenstrukturen, um das Hinzufügen / Entfernen von Knoten, Kanten, Dreiecken, Tetraedern usw. zu handhaben. Beispielsweise verwendet die p4est- Bibliothek Octree-Datenstrukturen zur adaptiven Netzverfeinerung . Sie würden nicht oft Oktrees...
Ich habe immer gehört, dass eine einfache Parallelisierung einer der Vorteile von DG-Methoden war, aber ich verstehe nicht wirklich, warum einer dieser Gründe nicht auch für kontinuierliches Galerkin
double trap(double func(double), double b, double a, double N) { double j; double s; double h = (b-a)/(N-1.0); //Width of trapezia double func1 = func(a); double func2; for (s=0,j=a;j<b;j+=h){ func2 = func(j+h); s = s + 0.5*(func1+func2)*h; func1 = func2; } return s; } Das Obige...
Was ist der Unterschied zwischen diesen beiden Methoden? Kann ein Problem, das mit einer Methode gelöst werden kann, mit der anderen gelöst werden? Können beide / oder einer von ihnen mit OpenMP und / oder MPI parallelisiert
Ich habe ein Projekt, in dem ich quadratische Felder verwenden muss. Speziell Zahlen der Form mita,b∈Q.a+b−3−−−√a+b−3a + b \sqrt{-3}a,b∈Qa,b∈Qa,b \in \mathbb{Q} Zum Beispiel sind hier die Primzahlen in Eisenstein-Ganzzahlen : Ich möchte keinen Salbei verwenden. Ich möchte meinen eigenen Datentyp...
In dem technischen Bericht über Galahad [1] geben die Autoren im Zusammenhang mit allgemeinen nichtlinearen Programmierproblemen an: Unserer Meinung nach gab es nie wirklich Zweifel daran, dass SQP-Methoden (Sequential Quadratic Programming) auf lange Sicht erfolgreicher sein würden [als...
Ich führe eine Zeilensuche als Teil eines Quasi-Newton-BFGS-Algorithmus durch. In einem Schritt der Zeilensuche verwende ich eine kubische Interpolation, um näher an den lokalen Minimierer heranzukommen. Sei die interessierende Funktion. Ich möchte ein so dass .x ≤ f ' ( x ≤ ) ≤ 0f: R → R , f∈...
Ich möchte eine Angelrute (oder ein Seil) modellieren, indem ich kurze Segmente verbinde. (Die Segmente können gleich lang sein, aber jedem Segment sollte eine eigene Masse zugewiesen werden.) Ein Segment beeinflusst das nächste durch das Drehmoment zwischen den Segmenten. Die Gelenke können...
Ich bin daran interessiert, die Lösung eines Lage-Systems von ODEs mit einer Krylov-Methode wie in [1] zu berechnen. Ein solches Verfahren beinhaltet Funktionen, die sich auf das Exponential beziehen (die sogenannten ). Es besteht im Wesentlichen darin, die Wirkung der Matrixfunktion zu berechnen,...
Einige ältere Bücher, die ich gesehen habe, sagen, dass die Mindestanzahl von Stufen einer expliziten Runge-Kutta-Methode einer bestimmten Reihenfolge für Bestellungen unbekannt ist . Ist das noch wahr?≥9≥9\geq 9 Welche Bibliotheken gibt es, um automatisch mit Runge-Kutta-Methoden hoher Ordnung zu...
Ich habe mich gefragt, ob jemand Erfahrung im Umgang mit Grenzen bei der Implementierung der Chebyshev-Differenzierung hat. Ich versuche derzeit, eine No-Slip-Randbedingung zu implementieren, um die inkompressiblen Navier-Stokes-Gleichungen in 3D zu lösen, um sicherzustellen, dass der Fluss an den...