Preguntas etiquetadas con floating-point

Tengo que resolver varios problemas desafiantes de optimización global no convexo. Actualmente uso la Caja de herramientas de optimización de MATLAB (específicamente, fmincon()con algoritmo = 'sqp'), que es bastante eficaz . Sin embargo, la mayor parte de mi código está en Python, y me encantaría...

Tengo una pregunta simple que es realmente difícil para Google (además del canónico Lo que todo informático debe saber sobre el papel de aritmética de punto flotante ). ¿Cuándo deberían usarse funciones como log1po en expm1lugar de logy exp? ¿Cuándo no deberían ser utilizados? ¿Cómo difieren las...

Dada la función real de las variables reales, ¿existe un software disponible que pueda generar automáticamente un código numérico exacto para calcular la función sobre todas las entradas en una máquina equipada con aritmética IEEE 754? Por ejemplo, si la función real a evaluar fuera: El...

Un problema muy común en Markov Chain Monte Carlo implica calcular probabilidades que son la suma de términos exponenciales grandes, ea1+ea2+...ea1+ea2+... e^{a_1} + e^{a_2} + ... donde los componentes de lata pueden variar de muy pequeños a muy grandes. Mi enfoque ha sido factorizar el término...

Se da una secuencia finita creciente de números positivos . ¿Cuál de los siguientes dos algoritmos es mejor para calcular la suma de los números?z1, z2, . . . . . znortez1,z2,.....znortez_{1} ,z_{2},.....z_{n} s=0; for \ i=1:n s=s + z_{i} ; end O: s=0; for \ i=1:n s=s + z_{n-i+1} ;...

Supongamos que tengo una función que toma como entrada varios valores de punto flotante (simple o doble), hace algunos cálculos y produce valores de punto flotante de salida (también simple o doble). Estoy trabajando principalmente con MSVC 2008, pero también planeo trabajar con MinGW / GCC. Estoy...

Estoy tratando de implementar la siguiente función en coma flotante de doble precisión con un error relativo bajo : 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)) Esto se usa ampliamente en aplicaciones estadísticas para agregar...

Dadas dos implementaciones de BLAS diferentes, ¿podemos esperar que hagan exactamente los mismos cálculos de coma flotante y devuelvan los mismos resultados? ¿O puede suceder, por ejemplo, que uno calcule un producto escalar como y uno como por lo que posiblemente da un resultado diferente en...

He ejecutar el código de dinámica molecular (MD) GROMACS en un clúster Linux Ubuntu que consta de nodos que contienen 24 CPUs Intel Xeon. Mi punto de interés particular resulta ser algo sensible a la precisión aritmética de coma flotante, por lo que he tenido que ejecutar GROMACS con doble...

Me gustaría determinar el número teórico de FLOP (operaciones de punto flotante) que puede hacer mi computadora. ¿Puede alguien ayudarme con esto? (Me gustaría comparar mi computadora con algunas supercomputadoras solo para tener una idea de la diferencia entre

Dejado x, yser dos números de punto flotante. ¿Cuál es la forma correcta de calcular su media? La forma ingenua (x+y)/2puede dar lugar a desbordamientos cuando xy yson demasiado grandes. Creo que 0.5 * x + 0.5 * ytal vez sea mejor, pero implica dos multiplicaciones (que tal vez sea ineficiente), y...

He estado leyendo el libro Computer Simulation of Liquids de Allen y Tildesley. A partir de la página 71, los autores analizan los diversos algoritmos que se utilizan para integrar las ecuaciones de movimiento de Newton en simulaciones de dinámica molecular (DM). A partir de la página 78, los...

Tengo problemas para implementar una función numéricamente. Sufre el hecho de que a valores de entrada grandes el resultado es un número muy grande multiplicado por un número muy pequeño. No estoy seguro de si la cancelación catastrófica es el término correcto, así que corríjame si es así....

Al evaluar el número de FLOP en una función simple, a menudo se puede bajar la expresión que cuenta los operadores aritméticos básicos. Sin embargo, en el caso de enunciados matemáticos que involucran una división pareja, uno no puede hacer esto y esperar poder comparar con los conteos FLOP de...

La mayoría de las bibliotecas matemáticas tienen varias versiones de funciones de logaritmo. La mayoría de las veces suponemos que son perfectos, pero en realidad muchos de ellos solo ofrecen un cierto número de dígitos de precisión. Para algunas funciones, hay variantes numéricamente más...

Tengo una función numérica que f(x, y)devuelve un número de coma flotante doble que implementa alguna fórmula y quiero verificar que sea correcta contra las expresiones analíticas para todas las combinaciones de los parámetros xy yque me interese. ¿Cuál es la forma correcta de comparar lo calculado...

Veo muy pocas bibliotecas / paquetes de computación de punto no flotante. Dadas las diversas imprecisiones de la representación de punto flotante, surge la pregunta de por qué no hay al menos algunos campos en los que esta mayor precisión pueda valer las complejidades de trabajar con punto...

He leido eso (1) Las operaciones mal condicionadas deben realizarse antes que las bien condicionadas. Como ejemplo, uno debe calcular x z- yzXz-yzxz-yz como ( x - y) z(X-y)z(x-y)z ya que la resta está mal condicionada mientras que la multiplicación no. Sin embargo, un análisis de error de...

Tengo un modelo considerable (~ 5000 líneas) escrito en C. Es un programa en serie, sin generación de números aleatorios en ninguna parte. Utiliza la biblioteca FFTW para funciones que usan FFT: no conozco los detalles de la implementación de FFTW, pero supongo que las funciones allí también son...

Recientemente, he encontrado un extraño problema con FORTRAN95. Inicialicé las variables X e Y de la siguiente manera: X=1.0 Y=0.1 Luego los agrego e imprimo el resultado: 1.10000000149012 Después de examinar las variables, parece que 0.1 no está representado con doble precisión con total...