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

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

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

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

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} ; end Na minha...

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

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

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 ) + ( x 3 y 3 + x 4( ( x1 1y1 1+...

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

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

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), e não...

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

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

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

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

Eu tenho uma função numérica f(x, y)retornando um número de ponto flutuante duplo que implementa alguma fórmula e quero verificar se ele está correto em relação às expressões analíticas para todas as combinações de parâmetros xe no yqual estou interessado. Qual é a maneira correta de comparar os...

Eu vejo muito poucas bibliotecas / pacotes de computação de ponto não flutuante. Dadas as várias imprecisões da representação de ponto flutuante, surge a questão de por que não existem pelo menos alguns campos em que essa precisão aumentada pode valer os meandros de trabalhar com ponto...

Eu li isso (1) Operações mal condicionadas devem ser realizadas antes das bem condicionadas. Como exemplo, deve-se calcular xz−yzxz−yzxz-yz como (x−y)z(x−y)z(x-y)z pois a subtração está mal condicionada enquanto a multiplicação não. Entretanto, uma análise de erro de primeira ordem de ambos...

Eu tenho um modelo considerável (~ 5000 linhas) escrito em C. É um programa serial, sem geração de números aleatórios em nenhum lugar. Utiliza a biblioteca FFTW para funções que utilizam a FFT - não conheço os detalhes da implementação da FFTW, mas presumo que as funções nela também sejam...

Recentemente, encontrei um problema bizarro com o FORTRAN95. Inicializei as variáveis ​​X e Y da seguinte maneira: X=1.0 Y=0.1 Mais tarde eu os adiciono e imprimo o resultado: 1.10000000149012 Após examinar as variáveis, parece que 0,1 não é representado em dupla precisão com total precisão....