Comment puis-je remplacer la méthode Euler par Runge-Kutta 4e ordre pour déterminer le mouvement de chute libre dans une amplitude gravitationnelle non constante (par exemple, chute libre à partir de 10 000 km au-dessus du sol)? Jusqu'à présent, j'ai écrit une intégration simple par la méthode...
C ++ 11 introduit la sémantique de déplacement qui peut, par exemple, améliorer les performances du code dans les situations où C ++ 03 devrait effectuer une construction de copie ou une affectation de copie. Cet article signale que le code suivant connaît une accélération de 5 fois lorsqu'il est...
Je voudrais utiliser la méthode Runge-Kutta de 8e ordre (89) dans une application de mécanique / astrodynamique céleste, écrite en C ++, à l'aide d'une machine Windows. Par conséquent, je me demande si quelqu'un connaît une bonne bibliothèque / implémentation documentée et gratuite à utiliser?...
J'ai mon propre petit sous-programme d'intégration numérique (quadrature), qui est une adaptation en C ++ d'un programme ALGOL publié par Bulirsch & Stoer en 1967 (Numerische Mathematik, 9, 271-278). Je voudrais passer à un algorithme plus moderne (adaptatif) et me demander s'il existe des...
En ce moment, je travaille sur une méthode de résolution d'équations différentielles appelée collocation base-spline. Ce avec quoi j'ai du mal à construire une méthode pour construire une spline d'ordre arbitraire, avec la relation avec la condition initiale B 1 i (x)={ 1Bk + 1je( x ) = x - xjeXk +...
J'ai donc débattu pour savoir si je devrais ou non prendre la peine d'apprendre Python. De parler avec mes professeurs, Matlab semble être le langage commun utilisé dans les mathématiques appliquées / science informatique en ce qui concerne le monde universitaire ; tandis que dans l' industrie ,...
J'ai un gros problème de valeur propre cubique: ( A0+ λ A1+λ2UNE2+λ3UNE3) x =0.(UNE0+λUNE1+λ2UNE2+λ3UNE3)X=0.\left(\mathbf{A}_0 + \lambda\mathbf{A}_1 + \lambda^2\mathbf{A}_2 + \lambda^3\mathbf{A}_3\right)\mathbf{x} = 0. Je pourrais résoudre cela en convertissant en un problème de valeur propre...
J'essaie actuellement de résoudre le problème de minimisation contraint non linéaire tel qu'implémenté dans la fonction matlab "fmincon". Mes attentes sont, minimiser (fun1, x0, uB, lB, fun2) où x0 est l'état initial, fun1 est une fonction qui doit être minimisée, uB sont des limites supérieures,...
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; } Ce qui précède est mon...
Pour mes travaux de recherche, j'ai développé une petite bibliothèque C ++ visant à faciliter la communication entre les codes de calcul C ++ et Octave / Matlab (lorsque ce dernier est utilisé à des fins de post-traitement). Je voudrais publier une telle bibliothèque sous GPL sur un site...
La plupart de mes programmes sont des codes de recherche uniques en C pour mon propre usage. Je n'ai jamais distribué de code à des collaborateurs proches. J'ai développé un algorithme que je publie dans une revue scientifique. Je veux fournir le code source et peut-être le code exécutable dans le...
Fondamentalement, FEM semble être un problème à peu près «résolu». Il existe de nombreux frameworks puissants, comme Trilinos, PETSc, FEniCS, Libmesh ou MOOSE. Une chose qu'ils ont en commun: ils sont extrêmement "lourds". Tout d'abord, l'installation est normalement super douloureuse....