Questions marquées «floating-point»

J'ai plusieurs problèmes d'optimisation globale non convexe difficiles à résoudre. Actuellement, j'utilise la boîte à outils Optimization de MATLAB (en particulier, fmincon()avec algorithm = 'sqp'), ce qui est assez efficace . Cependant, la majeure partie de mon code est en Python et j'aimerais...

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?...

Étant donné une fonction réelle de variables réelles, existe-t-il un logiciel capable de générer automatiquement un code numérique précis pour calculer la fonction sur toutes les entrées d'une machine équipée de l'arithmétique IEEE 754? Par exemple, si la fonction réelle à évaluer était: Le...

Un problème très courant dans Markov Chain Monte Carlo implique le calcul de probabilités qui sont la somme de grands termes exponentiels, eune1+ eune2+ . . .eune1+eune2+... e^{a_1} + e^{a_2} + ... où les composants d' peut varier de très petit à très grand. Mon approche a été de factoriser le plus...

Etant donné une suite finie croissante de nombres positifs . Lequel des deux algorithmes suivants est le meilleur pour calculer la somme des nombres?z1, z2, . . . . . znz1,z2,.....znz_{1} ,z_{2},.....z_{n} s=0; for \ i=1:n s=s + z_{i} ; end Ou: s=0; for \ i=1:n s=s + z_{n-i+1} ; end À mon avis, il...

Supposons que j'ai une fonction qui prend en entrée plusieurs valeurs à virgule flottante (simple ou double), effectue certains calculs et produit des valeurs à virgule flottante en sortie (également simples ou doubles). Je travaille principalement avec MSVC 2008, mais je prévois également de...

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...

Étant donné deux implémentations BLAS différentes, pouvons-nous nous attendre à ce qu'elles effectuent exactement les mêmes calculs en virgule flottante et renvoient les mêmes résultats? Ou peut-il arriver, par exemple, que l'on calcule un produit scalaire comme et un comme donnant ainsi...

J'ai exécuté le code de dynamique moléculaire (MD) GROMACS sur un cluster Linux Ubuntu composé de nœuds contenant 24 processeurs Intel Xeon. Mon point d'intérêt particulier se révèle être quelque peu sensible à la précision arithmétique à virgule flottante, j'ai donc dû exécuter GROMACS en double...

Je voudrais déterminer le nombre théorique de FLOP (opérations en virgule flottante) que mon ordinateur peut effectuer. Quelqu'un pourrait m'aider avec ça. (Je voudrais comparer mon ordinateur à certains superordinateurs juste pour avoir une idée de la différence entre

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...

J'ai lu le livre Computer Simulation of Liquids d'Allen et Tildesley. À partir de la page 71, les auteurs discutent des divers algorithmes utilisés pour intégrer les équations de mouvement de Newton dans les simulations de dynamique moléculaire (MD). À partir de la page 78, les auteurs discutent de...

J'ai du mal à implémenter une fonction numériquement. Il souffre du fait qu'à de grandes valeurs d'entrée, le résultat est un très grand nombre fois un très petit nombre. Je ne sais pas si l'annulation catastrophique est le terme correct, alors corrigez-moi si c'est le cas. Preuve que quelque chose...

Lors de l'évaluation du nombre de FLOP dans une fonction simple, on peut souvent simplement descendre l'expression en comptant les opérateurs arithmétiques de base. Cependant, dans le cas d'instructions mathématiques impliquant une division paire, on ne peut pas le faire et s'attendre à pouvoir...

La plupart des bibliothèques mathématiques ont un certain nombre de versions de fonctions logarithmes. La plupart du temps, nous supposons qu'ils sont parfaits, mais en fait, beaucoup d'entre eux n'offrent qu'un certain nombre de chiffres de précision. Pour certaines fonctions, il existe des...

Je vois très peu de bibliothèques / packages informatiques à virgule non flottante. Étant donné les diverses inexactitudes de la représentation en virgule flottante, la question se pose de savoir pourquoi il n'y a pas au moins certains domaines où cette précision accrue pourrait valoir la...

J'ai une fonction numérique f(x, y)renvoyant un double nombre à virgule flottante qui implémente une formule et je veux vérifier qu'elle est correcte par rapport aux expressions analytiques pour toutes les combinaisons de paramètres xet yque je suis intéressé. Quelle est la bonne façon de comparer...

J'ai lu ça (1) Les opérations mal conditionnées doivent être effectuées avant celles bien conditionnées. Par exemple, on devrait calculer comme puisque la soustraction est mal conditionnée tandis que la multiplication ne l'est pas.xz−yzxz−yzxz-yz(x−y)z(x−y)z(x-y)z Cependant, une analyse d'erreur de...

J'ai un modèle non négligeable (~ 5000 lignes) écrit en C. C'est un programme en série, sans génération de nombre aléatoire nulle part. Il utilise la bibliothèque FFTW pour les fonctions utilisant FFT - Je ne connais pas les détails de l'implémentation FFTW, mais je suppose que les fonctions y sont...

Récemment, j'ai rencontré un problème bizarre avec FORTRAN95. J'ai initialisé les variables X et Y comme suit: X=1.0 Y=0.1 Plus tard, je les additionne et j'imprime le résultat: 1.10000000149012 Après examen des variables, il semble que 0,1 ne soit pas représenté en double précision avec une...