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 hacer la optimización también en …

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 

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 diferentes implementaciones de …

30 floating-point 

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 software consideraría la cancelación …

25 software  floating-point  accuracy 

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 exponencial más grande para …

24 monte-carlo  floating-point  statistics 

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} …

22 algorithms  floating-point 

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 programando en C …

20 floating-point  c++ 

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 probabilidades o densidades de probabilidad que se representan en el espacio logarítmico. Por supuesto, exp(x)exp⁡(x)\exp(x) o exp(y)exp⁡(y)\exp(y) …

18 floating-point  stability  numerics 

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 el punto flotante IEEE …

17 floating-point  blas 

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 …

17 hpc  floating-point  molecular-dynamics  gpu 

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 …

15 floating-point 

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 ellas)

15 floating-point 

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í. Evidencia de …

13 python  numerical-analysis  floating-point 

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 autores …

13 floating-point  computer-arithmetic 

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 …

13 performance  complexity  floating-point 

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 estables. Por ejemplo, Fortran, R, Java …

11 floating-point  accuracy