J'ai une question simple qui est vraiment difficile pour Google (à part le canonique que tout informaticien devrait savoir sur le papier arithmétique à virgule flottante ). Quand faut log1p- expm1il utiliser des fonctions telles que ou au lieu de loget exp? Quand ne devraient-ils pas être utilisés?...
Quelqu'un peut-il me recommander un bon système de gestion de flux de travail (WMS), de préférence en Python? Jusqu'à présent, j'utilise GNU Make, mais cela introduit une couche de complexité que je veux éviter. Un bon WMS doit avoir les caractéristiques suivantes: s'intègre facilement avec les...
Les ingénieurs insistent souvent sur l'utilisation de méthodes localement conservatrices telles que le volume fini, la différence finie conservatrice ou les méthodes discontinues de Galerkin pour résoudre les PDE. Qu'est-ce qui peut mal tourner lorsque vous utilisez une méthode qui n'est pas...
En tant que titulaire d'un BA en physique, j'ai été quelque peu scandalisé lorsque j'ai commencé à travailler avec des simulations moléculaires. Ce fut un peu un choc de découvrir que même les simulations les plus détaillées et les plus coûteuses en calcul ne peuvent pas reproduire quantitativement...
Le problème Je travaille actuellement sur une simulation d'élément fini Navier Stokes et j'aimerais étudier les effets d'une variété de paramètres. Certains paramètres sont spécifiés dans un fichier d'entrée ou via une option de ligne de commande; d'autres paramètres sont fournis sous forme...
Contexte: Je pense que je pourrais vouloir porter du code qui calcule les produits matriciels à vecteur exponentiel en utilisant une méthode de sous-espace Krylov de MATLAB à Python. (Plus précisément, la fonction expmvp de Jitse Niesen , qui utilise un algorithme décrit dans cet article .)...
J'ai récemment rencontré une formulation du méta-phénomène : " deux c'est facile, trois c'est dur " (formulé de cette façon par Federico Poloni), qui peut être décrit comme suit: Lorsqu'un certain problème est formulé pour deux entités, il est relativement facile à résoudre; cependant, un...
Le Journal of Computational Physics a été un débouché important pour la science informatique dans le passé, et j'y ai déjà publié. Pour le bénéfice de ceux (comme moi) qui ont signé le boycott d'Elsevier , quelles revues non Elsevier seraient des endroits appropriés pour publier des articles qui...
Si j'ai une matrice inversible carrée et que je prends son déterminant, et que je trouve que , cela implique-t-il que la matrice est mal conditionnée?det(A)≈0det(A)≈0\det(A) \approx 0 L'inverse est-il également vrai? Une matrice mal conditionnée a-t-elle un déterminant presque nul? Voici quelque...
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. La version 3.0 de la norme MPI a...
En statistique et ses diverses applications, nous calculons souvent la matrice de covariance , qui est définie positive (dans les cas considérés) et symétrique, pour différentes utilisations. Parfois, nous avons besoin de l'inverse de cette matrice pour divers calculs (formes quadratiques avec cet...
Je veux modifier une matrice de transition carrée dense en place en changeant l'ordre de plusieurs de ses lignes et colonnes, en utilisant la bibliothèque numpy de python. Mathématiquement, cela correspond à la pré-multiplication de la matrice par la matrice de permutation P et à la...
Cette question a été migrée à partir de Stack Overflow car il est possible d'y répondre sur Computational Science Stack Exchange. Migré il y a 7 ans . Existe-t-il des problèmes / algorithmes célèbres dans le calcul scientifique qui ne peuvent pas être accélérés par la parallélisation? Il me semble...
Je connais et j'ai utilisé f2py2e pour envelopper un ancien code Fortran 77, mais je crois comprendre que cela ne fonctionne pas avec le code Fortran 95 plus récent. J'ai recherché ce que je devrais utiliser, et j'ai rencontré fwrap et G3 f2py, qui ne semblent pas donner d'explication sur leur état...
EDIT: Je teste si des valeurs propres ont une magnitude de un ou plus. J'ai besoin de trouver la plus grande valeur propre absolue d'une grande matrice clairsemée et non symétrique. J'ai utilisé la eigen()fonction de R , qui utilise l'algo QR d'EISPACK ou de LAPACK pour trouver toutes les valeurs...
Les développeurs de logiciels ont le choix de choisir une licence appropriée en fonction des objectifs du travail. Quelqu'un peut-il donner des recommandations / expériences sur quelle licence choisir pour le logiciel? Quels sont les avantages / inconvénients de "donner" tout le travail codé en...
Qu'est-ce qui peut mal tourner lors de l'utilisation de méthodes Krylov précondonifiées du KSP ( package de solveur linéaire de PETSc ) pour résoudre un système linéaire clairsemé tel que ceux obtenus en discrétisant et en linéarisant des équations différentielles partielles? Quelles mesures...
L'un des principaux problèmes auxquels nous devons faire face dans les simulations moléculaires est le calcul des forces dépendantes de la distance. Si nous pouvons restreindre les fonctions de force et de distance pour avoir des puissances égales de la distance de séparation , alors nous pouvons...
Je ne connais pas très bien les schémas de discrétisation courants pour les PDE. Je sais que Crank-Nicolson est un schéma populaire pour discrétiser l'équation de diffusion. Est également un bon choix pour le terme d'advection? Je suis intéressé par la résolution de l'équation...
Existe-t-il des outils spécialement conçus pour compresser les données scientifiques en virgule flottante? Si une fonction est lisse, il y a évidemment beaucoup de corrélation entre les nombres représentant cette fonction, donc les données doivent bien se compresser. Zipper / gzipper des données...