J'ai eu un aperçu de l'analyse numérique (principalement, des méthodes numériques comme la recherche de racines, les équations quadratiques et d'autres choses préliminaires) dans ma classe de calcul, mais maintenant, je me retrouve à vouloir plus de sophistication dans mon travail. Existe-t-il un...
Il existe de nombreux schémas FD pour l'équation d'advection discuter sur le web. Par exemple ici: http://farside.ph.utexas.edu/teaching/329/lectures/node89.html∂T∂t+ u ∂T∂X= 0∂T∂t+u∂T∂X=0\frac{\partial T}{\partial t}+u\frac{\partial T}{\partial x}=0 Mais je n'ai vu personne proposer un schéma de...
quelles sont les méthodes de pointe pour la solution numérique des ODE à côté droit discontinu? Je suis surtout intéressé par les fonctions du côté droit par morceaux, par exemple le signe. J'essaie de résoudre l'équation d'un type suivant: X˙v˙= v= { ( | Fexterne| - | Ffriction| )signe(...
Je suis intéressé à résoudre l'équation de Poisson en utilisant l'approche des différences finies. Je voudrais mieux comprendre comment écrire l'équation matricielle avec les conditions aux limites de Neumann. Quelqu'un pourrait-il revoir ce qui suit, est-ce exact? La matrice des différences finies...
J'exécute des simulations de dynamique moléculaire (MD) en utilisant plusieurs progiciels, comme Gromacs et DL_POLY. Gromacs prend désormais en charge les algorithmes de décomposition des particules et de décomposition de domaine. Par défaut, les simulations Gromacs utilisent la décomposition de...
Existe-t-il un algorithme (efficace) pour sélectionner un sous-ensemble de points dans un ensemble de points ( ) de telle sorte qu'ils "couvrent" la plus grande partie de la zone (sur tous les sous-ensembles possibles de taille )?MMMNNNM< NM<NM < NMMM Je suppose que les points sont dans un...
Laissez x, ydeux nombres à virgule flottante. Quelle est la bonne façon de calculer leur moyenne? La façon naïve (x+y)/2peut entraîner des refoulements quand xet ysont trop grandes. Je pense que c'est 0.5 * x + 0.5 * ypeut-être mieux, mais cela implique deux multiplications (ce qui est peut-être...
Y a-t-il des cas où Hartree-Fock n'est pas une bonne approximation pour calculer la géométrie d'équilibre lorsque la molécule n'est pas en rupture de liaison?
Je suis nouveau dans les simulations de dynamique moléculaire (MD). Quelle est la complexité d'une simulation de dynamique moléculaire en termes de temps de simulation? En d'autres termes, si je veux augmenter le temps simulé de 10 nanosecondes à 20 nanosecondes, à quoi puis-je m'attendre en termes...
Les méthodes multigrilles résolvent généralement les problèmes de Dirichlet aux niveaux (par exemple, le point Jacobi ou Gauss-Seidel). Lors de l'utilisation de méthodes d'éléments finis continues, il est beaucoup moins coûteux d'assembler de petits problèmes de Neumann que d'assembler de petits...
Je travaille assez fortement en biologie mathématique / épidémiologie, où la plupart des travaux de modélisation / science informatique sont encore dominés par des ensembles d'OD, certes parfois assez élaborés. L'un des avantages de ces modèles est qu'ils sont assez faciles à décrire et à...
Nos simulations de dynamique du plasma produisent souvent trop d'informations. Au cours des simulations, nous enregistrons diverses propriétés physiques sur une grille (x, y, z, t) qui est aussi grande que (8192x1024x1024x1500), pour au moins 10 propriétés. Ces informations sont traitées une fois...
Je compare deux méthodes itératives pour inverser des matrices carrées aléatoires. Étant donné que les matrices sont aléatoires, chaque cas de test prend à la fois différentes quantités d'itérations et différents temps écoulés. Ma question est, en plus du temps CPU moyen, est la valeur moyenne des...
Je veux calculer le spectre ( toutes les valeurs propres) d'une grande matrice clairsemée (des centaines de milliers de lignes). C'est dur. Je suis prêt à me contenter d'une approximation. Existe-t-il des méthodes d'approximation pour ce faire? Bien que j'espère une réponse générale à cette...
Si je comprends bien, étant donné qu'une solution à un programme linéaire se produit toujours au sommet de son ensemble réalisable polyédrique (si une solution existe et que la valeur optimale de la fonction objective est limitée par le bas, en supposant un problème de minimisation), comment une...
Je viens de commencer à jouer avec FEniCS. Je résout Poisson avec des éléments de troisième ordre et je souhaite visualiser les résultats. Cependant, lorsque j'utilise plot (u), la visualisation n'est qu'une interpolation linéaire des résultats. J'obtiens la même chose lorsque je produis en VTK....
En pratique, le temps d'exécution de la résolution numérique d'un IVP est souvent dominé par la durée de l'évaluation du côté droit (RHS) . Supposons donc que toutes les autres opérations sont instantanées (c'est-à-dire sans coût de calcul). Si le temps d'exécution global pour résoudre l'IVP est...
Je suis en train de résoudre un système de deux PDE couplés dans deux dimensions spatiales et en temps de calcul. Étant donné que les évaluations de fonctions sont coûteuses, je voudrais utiliser une méthode en plusieurs étapes (initialisée à l'aide de Runge-Kutta 4-5). La méthode Adams-Bashforth...
La plupart des bibliothèques C ++ les plus connues en science informatique telles que Eigen , Trilinos et deal.II utilisent l'objet de bibliothèque d'en-tête de modèle C ++ standard std::complex<>, pour représenter des nombres à virgule flottante complexes. Dans la réponse de Jack Poulson à...
Quelqu'un ici a-t-il utilisé le calcul scientifique en double précision avec des GPU de nouvelle génération (par exemple K20) via Python? Je sais que cette technologie évolue rapidement, mais quelle est la meilleure façon de le faire actuellement? Le GPU est hors de portée des bibliothèques...