Perguntas com a marcação «floating-point»

Eu tenho vários problemas desafiadores de otimização global não convexa para resolver. Atualmente, uso o Optimization Toolbox do MATLAB (especificamente, fmincon()com o algoritmo = 'sqp'), o que é bastante eficaz . No entanto, a maior parte do meu código está em Python, e eu adoraria fazer a otimização em Python …

77 python  optimization  nonlinear-programming  linear-algebra  linear-solver  sparse-matrix  ode  automatic-differentiation  efficiency  matrix  symbolic-computation  pde  multiphysics  operator-splitting  fluid-dynamics  pde  multigrid  pde  hyperbolic-pde  quadrature  precision  numerical-analysis  fluid-dynamics  interpolation  c  ode  testing  data-sets  pde  fluid-dynamics  hyperbolic-pde  pde  finite-element  fluid-dynamics  stability  incompressible  pde  hyperbolic-pde  fluid-dynamics  pde  precision  floating-point  convex-optimization  conditioning  pde  hyperbolic-pde  stability  implicit-methods  explicit-methods  fluid-dynamics  accuracy  pde  hyperbolic-pde  petsc  linear-algebra 

Eu tenho uma pergunta simples que é realmente difícil para o Google (além do artigo canônico O que todo cientista da computação deve saber sobre aritmética de ponto flutuante ). Quando funções como log1pou expm1devem ser usadas em vez de loge exp? Quando eles não devem ser usados? Como as …

30 floating-point 

Dada uma função real de variáveis ​​reais, existe um software disponível que possa gerar automaticamente código com precisão numérica para calcular a função em todas as entradas em uma máquina equipada com aritmética IEEE 754? Por exemplo, se a função real a ser avaliada fosse: O software consideraria o cancelamento …

25 software  floating-point  accuracy 

Um problema muito comum na cadeia de Markov Monte Carlo envolve probabilidades de computação que são soma de grandes termos exponenciais, ea1+ea2+...ea1+ea2+... e^{a_1} + e^{a_2} + ... onde os componentes de lata podem variar de muito pequeno a muito grande. Minha abordagem foi fatorar o maior termo exponencial para que:K …

24 monte-carlo  floating-point  statistics 

Dada é uma sequência finita crescente de números positivos . Qual dos dois algoritmos a seguir é melhor para calcular a soma dos números?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} ; …

22 algorithms  floating-point 

Suponha que eu tenha uma função que tome como entrada vários valores de ponto flutuante (simples ou duplo), faça alguma computação e produza valores de ponto flutuante de saída (também únicos ou duplos). Estou trabalhando principalmente com o MSVC 2008, mas também planejo trabalhar com o MinGW / GCC. Estou …

20 floating-point  c++ 

Estou tentando implementar a seguinte função no ponto flutuante de precisão dupla com baixo erro relativo : 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)) Isso é usado extensivamente em aplicativos estatísticos para adicionar probabilidades ou densidades de probabilidade representadas no espaço de log. Obviamente, exp(x)exp⁡(x)\exp(x) ou exp(y)exp⁡(y)\exp(y) poderiam facilmente estourar ou estourar, …

18 floating-point  stability  numerics 

Dadas duas implementações BLAS diferentes, podemos esperar que elas façam exatamente os mesmos cálculos de ponto flutuante e retornem os mesmos resultados? Ou pode acontecer, por exemplo, que alguém calcule um produto escalar como e um como ( x 1 y 1 + x 2 y 2 ) + ( …

17 floating-point  blas 

Eu executei o código de dinâmica molecular (MD) GROMACS em um cluster Ubuntu Linux que consiste em nós contendo 24 CPUs Intel Xeon. Meu ponto de interesse específico acaba sendo um pouco sensível à precisão aritmética de ponto flutuante, por isso tive que executar o GROMACS em precisão dupla ao …

17 hpc  floating-point  molecular-dynamics  gpu 

Eu gostaria de determinar o número teórico de FLOPs (operações de ponto flutuante) que meu computador pode executar. Alguém por favor pode me ajudar com isso. (Gostaria de comparar meu computador a alguns supercomputadores apenas para ter uma idéia da diferença entre eles)

15 floating-point 

Vamos x, yser dois números de ponto flutuante. Qual é o caminho certo para calcular a média deles? A maneira ingênua (x+y)/2pode resultar em estouros quando xe ysão muito grandes. Acho que 0.5 * x + 0.5 * ytalvez seja melhor, mas envolve duas multiplicações (o que talvez seja ineficiente), …

15 floating-point 

Estou tendo problemas para implementar uma função numericamente. Ele sofre com o fato de que, com valores de entrada grandes, o resultado é um número muito grande vezes um número muito pequeno. Não tenho certeza se o cancelamento catastrófico é o termo correto, por favor, corrija-me se for. Evidência de …

13 python  numerical-analysis  floating-point 

Eu tenho lido o livro Computer Simulation of Liquids de Allen e Tildesley. A partir da página 71, os autores discutem os vários algoritmos usados ​​para integrar as equações de movimento de Newton em simulações de dinâmica molecular (MD). A partir da página 78, os autores discutem o algoritmo Verlet, …

13 floating-point  computer-arithmetic 

Ao avaliar o número de FLOPs em uma função simples, geralmente é possível simplesmente descer a expressão que calcula os operadores aritméticos básicos. No entanto, no caso de declarações matemáticas envolvendo divisão uniforme, não se pode fazer isso e espera poder comparar com contagens de FLOP de funções com apenas …

13 performance  complexity  floating-point 

A maioria das bibliotecas de matemática possui várias versões de funções de logaritmo. Na maioria das vezes, assumimos que eles são perfeitos, mas na verdade muitos deles oferecem apenas um certo número de dígitos de precisão. Para algumas funções, existem variantes numericamente mais estáveis. Por exemplo, Fortran, R, Java e …

11 floating-point  accuracy