Science computationnelle

Questions-réponses pour les scientifiques utilisant des ordinateurs pour résoudre des problèmes scientifiques

6
Quelle est la meilleure façon de faire des recherches reproductibles si vous avez besoin de bibliothèques propriétaires?
La recherche reproductible en calcul vise à mettre le code nécessaire pour générer les résultats dans un document de calcul à la disposition d'autres chercheurs afin qu'ils puissent exécuter ce code pour reproduire les résultats dans ce document. Je voudrais que toutes mes recherches soient reproductibles, mais je me heurte …
6
Comment écrire du code dimensionnellement agnostique?
Je me retrouve souvent à écrire du code très similaire pour les versions à une, deux et trois dimensions d'une opération / algorithme donné. La maintenance de toutes ces versions peut devenir fastidieuse. La génération de code simple fonctionne assez bien, mais il semble que l'on pense qu'il doit y …
3
Quelles directives dois-je utiliser lors de la recherche de bonnes méthodes de préconditionnement pour un problème spécifique?
Pour la solution de grands systèmes linéaires utilisant des méthodes itératives, il est souvent intéressant d'introduire un préconditionnement, par exemple résoudre à la place M - 1 ( A x = b ) , où M est ici utilisé pour le préconditionnement à gauche du système. Typiquement, nous devrions avoir …
1
Pourquoi Octrees est-il utilisé pour la décomposition d'espace multipolaire?
Dans la plupart (toutes?) Des implémentations de la méthode multipolaire rapide (FMM), les octrees sont utilisés pour décomposer le domaine concerné. Théoriquement, les octrees fournissent une borne volumétrique simple, qui est utile pour prouver le temps d'exécution O (n) d'un FMM. Au-delà de cette justification théorique, y a-t-il des avantages …
18 algorithms 
2
Qu'est-ce qu'un pseudo pas de temps?
En lisant de la littérature sur les solveurs PDE, je suis tombé sur le terme pseudo pas de temps aujourd'hui. Il semble que ce soit un terme courant, mais je n'ai pas réussi à trouver une bonne définition ou un article d'introduction. Par conséquent: Qu'est-ce qu'un pseudo pas de temps …
2
F2Py avec des tableaux de formes attribuables et supposés
Je voudrais utiliser f2pyavec le Fortran moderne. En particulier, j'essaie de faire fonctionner l'exemple de base suivant. Ceci est le plus petit exemple utile que j'ai pu générer. ! alloc_test.f90 subroutine f(x, z) implicit none ! Argument Declarations ! real*8, intent(in) :: x(:) real*8, intent(out) :: z(:) ! Variable Declarations …
18 python  fortran 
1
Annulation catastrophique dans LogSum
J'essaie d'implémenter la fonction suivante en virgule flottante double précision avec une faible erreur relative : logsum(x,y)=log(exp(x)+exp(y))logsum(x,y)=log⁡(exp⁡(x)+exp⁡(y))\mathrm{logsum}(x,y) = \log(\exp(x) + \exp(y)) Ceci est largement utilisé dans les applications statistiques pour ajouter des probabilités ou des densités de probabilité qui sont représentées dans l'espace logarithmique. Bien sûr, ou pourrait facilement déborder …
1
Comment les ondelettes peuvent-elles être appliquées à la PDE?
Je voudrais savoir comment les méthodes d'ondelettes peuvent être appliquées à la PDE, mais malheureusement je ne connais pas de bonne ressource pour en savoir plus sur ce sujet. Il semble que de nombreuses introductions aux ondelettes se concentrent sur la théorie de l'interpolation, par exemple, l'assemblage d'un signal par …
18 pde  wavelet 
4
Existe-t-il une bibliothèque à usage général pour le raffinement de maillage adaptatif à grille structurée?
Voulez-vous améliorer ce post? Fournissez des réponses détaillées à cette question, y compris des citations et une explication de la raison pour laquelle votre réponse est correcte. Les réponses sans suffisamment de détails peuvent être modifiées ou supprimées. Le raffinement adaptatif du maillage (AMR) est une technique courante pour traiter …
4
Quelles méthodes peuvent garantir que les quantités physiques restent positives tout au long d'une simulation PDE?
Les quantités physiques comme la pression, la densité, l'énergie, la température et la concentration doivent toujours être positives, mais les méthodes numériques calculent parfois des valeurs négatives au cours du processus de solution. Ce n'est pas correct car les équations calculent des valeurs complexes ou infinies (écrasant généralement le code). …
En utilisant notre site, vous reconnaissez avoir lu et compris notre politique liée aux cookies et notre politique de confidentialité.
Licensed under cc by-sa 3.0 with attribution required.