Wie kann ich die Euler-Methode durch Runge-Kutta 4. Ordnung ersetzen, um die Bewegung des freien Falls in nicht konstanter Gravitationsgröße (z. B. freier Fall aus 10 000 km Höhe) zu bestimmen? Bisher habe ich eine einfache Integration nach der Euler-Methode geschrieben: while() { v +=...
Ich habe meine eigene kleine Subroutine für die numerische Integration (Quadratur), eine C ++ - Adaption eines ALGOL-Programms, das 1967 von Bulirsch & Stoer veröffentlicht wurde (Numerische Mathematik, 9, 271-278). Ich möchte auf einen moderneren (adaptiven) Algorithmus upgraden und mich...
C ++ 11 führt eine Verschiebungssemantik ein, die beispielsweise die Codeleistung in Situationen verbessern kann, in denen C ++ 03 eine Kopierkonstruktion oder eine Kopierzuweisung durchführen müsste. Dieser Artikel berichtet, dass der folgende Code beim Kompilieren mit C + 11 eine 5-fache...
Im Moment arbeite ich an einer Differentialgleichungslösungsmethode namens Basis-Spline-Kollokation. Ich habe Probleme damit, eine Methode zum Erstellen eines Splines beliebiger Ordnung mit der Beziehung zu erstellen. mit dem Anfangsbedingung und ich habe Probleme, selbst mit diesem Problem zu...
Ich möchte die Runge-Kutta-Methode 8. Ordnung (89) in einer in C ++ geschriebenen Anwendung für Himmelsmechanik / Astrodynamik unter Verwendung eines Windows-Computers verwenden. Daher frage ich mich, ob jemand eine gute Bibliothek / Implementierung kennt, die dokumentiert und kostenlos zu...
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...
Grundsätzlich scheint FEM ein Problem zu sein, das so ziemlich "gelöst" ist. Es gibt zahlreiche leistungsstarke Frameworks wie Trilinos, PETSc, FEniCS, Libmesh oder MOOSE. Eines haben sie gemeinsam: Sie sind extrem "schwer". Erstens ist die Installation normalerweise sehr schmerzhaft. Zweitens ist...
Ich verwende Armadillo, um sehr intensive Matrixmultiplikationen mit Seitenlängen von , wobei bis zu 20 oder sogar mehr kann. Ich verwende Armadillo mit OpenBLAS für die Matrixmultiplikation, was in parallelen Kernen sehr gute Arbeit zu leisten scheint, außer dass ich ein Problem mit dem...
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...
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...
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...
Gibt es Präprozessoranweisungen, die zur Verwendung der Polylog-Funktion verwendet werden könnten? Oder ist es in cmath enthalten? Wenn ja, nennen Sie es Li oder Polylog? EDIT: Was ich wirklich versuche, ist einen analytischen Wert für das unbestimmte Integral der Funktion anzugeben x3ex-...
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...
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...
Ich suche nach einer praktischen, kostenlosen C ++ - Bibliothek für Diagramme, die einen Löser für das MCS-Problem (Maximum Common Subgraph) enthält. Ich kenne die Boost Graph Library und LEMON , aber keiner enthält einen
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...
Mein aktuelles Projekt ist die Neuprogrammierung eines Proteinfaltungsmodells, bei dem Tausende von ODEs in C ++ gelöst werden. Ich habe einige Stopp- und Startfortschritte gemacht, während ich den Solver schreibe, der vollständig auf der GPU ausgeführt werden soll. Ich habe es endlich integriert,...
Ich bin ein Doktorand in Scientific Computing und habe in den letzten Monaten viel Zeit damit verbracht, Python und C ++ richtig zu lernen. Ich habe das Gefühl, dass ich C ++ gut gelernt habe und mit Python tun kann, was ich will, wenn ich ein gutes Nachschlagewerk habe. Ich kenne MATLAB auch...