Je calcule toujours les choses à la main, mais maintenant mes camarades deviennent méchants et font beaucoup d'exercices répétitifs impliquant simplement de brancher des choses comme l'expression ci-dessus. Je suis particulièrement intéressé par les logiciels open source tels que Python ou R pour...
Je développe un code pour simuler l'écoulement des fluides avec des substances biologiques présentes dans l'écoulement. Cela implique les équations de Navier-Stokes standard couplées à certains modèles biologiques supplémentaires. Il existe de nombreux paramètres / constantes. J'ai écrit des...
Le langage C ++ fournit une programmation générique et une métaprogrammation via des modèles. Ces techniques ont trouvé leur chemin dans de nombreux progiciels de calcul scientifique à grande échelle (par exemple, MPQC , LAMMPS , CGAL , Trilinos ). Mais qu'est-ce qu'ils ont réellement contribué au...
Je sais que la méthode de Nitsche est une méthode très intéressante car elle permet de prendre en compte les conditions aux limites de type Dirichlet ou le contact avec les conditions aux limites de frottement de façon faible sans utiliser de multiplicateurs de Lagrange. Et son avantage, qui est de...
J'ai besoin d'écrire quelques données à partir d'un calcul, qui seront lues plus tard par Paraview (fichier .vtu ou vtk). En ce qui concerne la taille du fichier, dois-je opter pour le format ASCII ou le format
Supposons que j'ai une fonction FFf et que je veuille trouver XXx tel que F( x ) ≈ 0F(X)≈0f(x)\approx 0 . Je pourrais utiliser la méthode Newton-Raphson. Mais cela nécessite que je connaisse la fonction dérivée F′( x )F′(X)f'(x) . Une expression analytique pour FFf peut ne pas être disponible. Par...
J'ai une fonction 2D non convexe bornée dont j'aimerais trouver le minimum. La fonction est assez fluide. L'évaluer coûte cher. Une erreur acceptable est d'environ 3% du domaine de la fonction dans chaque axe. J'ai essayé d'exécuter l'implémentation de l'algorithme DIRECT dans la bibliothèque...
Récemment, j'ai comparé différents solveurs non linéaires de scipy et j'ai été particulièrement impressionné par l' exemple de Newton-Krylov dans le livre de recettes Scipy dans lequel ils résolvent une équation d'équation différentielle de second ordre avec un terme de réaction non linéaire dans...
De nombreuses méthodes numériques pour les EDP hyperboliques sont basées sur l'utilisation de solveurs de Riemann. De tels solveurs sont essentiels pour capturer avec précision les ondes de choc. Il existe une gamme de ces solveurs disponibles pour les systèmes les mieux étudiés (par exemple, les...
Les méthodes numériques de résolution des EDP (ou ODE) se divisent en deux grandes catégories: les méthodes explicites et implicites. Les méthodes implicites permettent des pas de temps stables plus grands mais nécessitent plus de travail par étape. Pour les EDP hyperboliques, la sagesse courante...
Supposons que j'ai un maillage 2D constitué de triangles non chevauchants , et un ensemble de points { p i } M i = 1 ⊂ ∪ N k = 1 T K . Quelle est la meilleure façon de déterminer dans quel triangle chacun des points se
Étant donné une évolution PDE ut= A u + B uut=UNEu+Buu_t = Au + Bu où sont des opérateurs différentiels (éventuellement non linéaires) qui ne commutent pas, une approche numérique courante consiste à alterner entre la résolutionA , BUNE,BA,B ut= A uut=UNEuu_t = Au et ut= B u .ut=Bu.u_t = Bu. La...
La chose que je n'aime pas le plus à propos de MPI concerne les types de données (c'est-à-dire les cartes / masques de données) car ils ne correspondent pas très bien au C ++ orienté objet. boost::mpine prend en charge que MPI 1.1, cependant, à partir de leur site Web: boost :: mpi est une...
Je suis dans un domaine qui ne fait pas nécessairement beaucoup de travail HPC, et quand il le rencontre, c'est souvent le résultat de chercheurs d'autres domaines qui explorent de nouvelles applications à leurs méthodes et autres. Ce que cela signifie principalement, c'est qu'il n'est jamais...
Je suis curieux de savoir quels sont les bons algorithmes numériques pour évaluer la fonction hypergéométrique généralisée (ou série), définie comme pFq( un1, … , Unp; b1, … , Bq; z) = ∑k = 0∞( un1)k⋯ ( ap)k( b1)k⋯ ( bq)kzkk
Considérons avec presque singulier, ce qui signifie qu'il existe une valeur propre de qui est très petite. Le critère d'arrêt habituel d'une méthode itérative est basé sur le résiduel et concerne les itérations peuvent s'arrêter lorsque avec n le numéro d'itération. Mais dans le cas que nous...
Je fais une tâche de classification de texte avec R, et j'obtiens une matrice de termes de document avec la taille 22490 par 120 000 (seulement 4 millions d'entrées non nulles, moins de 1% d'entrées). Maintenant, je veux réduire la dimensionnalité en utilisant PCA (analyse en composantes...
Il s'agit probablement d'une question de niveau étudiant, mais je ne peux pas vraiment m'arranger. Pourquoi est-il plus précis d'utiliser des grilles non uniformes dans les méthodes numériques? Je pense dans le contexte d'une méthode de différence finie pour l'EDP de la forme . Et supposons que je...
J'utilise Valgrind + Callgrind pour profiler un solveur que j'ai écrit. Comme l'indique le manuel de l'utilisateur de Valgrind, j'ai compilé mon code avec les options de débogage du compilateur: "Sans informations de débogage, les meilleurs outils Valgrind seront capables de deviner à quelle...
J'aimerais être indépendant des logiciels commerciaux pour mon travail scientifique. Je trouve une dépendance à des packages commerciaux tels que Matlab et ses boîtes à outils insatisfaisante, car je ne sais pas si j'aurai accès à Matlab à l'avenir, et parce que je n'aime pas la langue. Par...