Preguntas etiquetadas con numerics

En mi campo de investigación, la especificación de los errores experimentales es comúnmente aceptada y las publicaciones que no los proporcionan son muy criticadas. Al mismo tiempo, a menudo encuentro que los resultados de los cálculos numéricos se proporcionan sin tener en cuenta los errores numéricos, a pesar de que …

40 error-estimation  numerics 

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 

Estoy paralelizando el código para resolver numéricamente un modelo de equilibrio de población de 5 dimensiones. Actualmente tengo un muy buen código paralelo MPICH2 en FORTRAN, pero a medida que aumentamos los valores de los parámetros, las matrices se vuelven demasiado grandes para ejecutarse en modo de memoria distribuida. Tengo …

14 parallel-computing  fortran  numerics  mpi  compiling 

Cada matriz real se puede reducir a un formato Schur real de usando un ortogonal similiary transformada . Aquí la matriz es una forma cuasi triangular con 1 por 1 o 2 por 2 bloques en la diagonal principal. Cada 1 por 1 bloque corresponde a un valor propio real …

10 eigenvalues  eigensystem  lapack  numerics  scalapack 

Los algoritmos existentes para resolver ODE manejan funciones , dondey∈Rn. Pero en muchos sistemas físicos, la ecuación diferencial es autónoma, entoncesdyreyret= f( y, t )reyret=F(y,t)\frac{dy}{dt} = f(y, t)y∈ Rnortey∈Rnortey \in \mathbb R^n,y∈Rn, con latfuera. Con este supuesto simplificador, ¿qué mejoras se pueden ver en los métodos numéricos existentes? Por ejemplo, …

10 ode  numerics 

Existen muchos métodos numéricos bien conocidos para resolver ecuaciones del tipo por ejemplo, método de bisección, método de Newton, etc.f( x ) =0,x ∈ Rnorte,f(X)=0 0,X∈Rnorte, f(x) = 0, \quad x \in \mathbb{R}^n, En mi aplicación, se calcula con un método estocástico (el resultado es un promedio).F( x )F(X)f(x) ¿Hay …

10 monte-carlo  nonlinear-equations  stochastic  numerics 

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 fijo. ¿Hay …

10 floating-point  numerics 

Cuando tenemos un modelo numérico que representa un sistema físico real, y que exhibe caos (por ejemplo, modelos de dinámica de fluidos, modelos climáticos), ¿cómo podemos saber que el modelo está funcionando como debería? No podemos comparar dos conjuntos de resultados del modelo directamente, porque incluso pequeños cambios en las …

10 testing  numerics 

Estoy trabajando con funciones que, en general, son mucho más suaves y se comportan mejor en el espacio log-log --- así que ahí es donde realizo la interpolación / extrapolación, etc., y eso funciona muy bien. ¿Hay alguna manera de integrar estas funciones numéricas en el espacio log-log? es decir, …

10 numerics  integration 

Soy doble licenciado en informática y matemáticas. Amo los dos temas. Estoy pensando en tomar una carrera de posgrado, tal vez en informática científica. ¿Cuál es la verdadera diferencia entre la computación científica y el análisis numérico? ¿Se estudian como carreras?

9 numerics  education  terminology  career-development 

Antecedentes: solo he creado una solución numérica de trabajo para 2d Navier-Stokes, para un curso. Era una solución para el flujo de la cavidad impulsada por la tapa. Sin embargo, el curso discutió una serie de esquemas para discretizaciones espaciales y discretizaciones de tiempo. También he tomado más cursos de …

9 finite-difference  fluid-dynamics  machine-learning  numerics 

¿Cuáles son las principales diferencias entre FEM y XFEM? ¿Cuándo deberíamos (no) usar XFEM intead de FEM y viceversa? En otras palabras, cuando encuentro un nuevo problema, ¿cómo puedo saber cuál de ellos debo usar?

9 finite-element  numerical-analysis  adaptive-mesh-refinement  numerics 

Necesito calcular la integral yo= ∫R0 0F( r ) Jnorte( zn mrR) rdrI=∫0Rf(r)Jn(znmrR)rdrI = \int_0^R f(r)J_n\left(\frac{z_{nm}r}{R}\right)rdr donde es el n t h para funciones de Bessel del primer tipo, z n m es su m t h cero y f ( r ) es una función real que es algo …

8 quadrature  integration  special-functions  numerics 

Me gustaría calcular, como parte de la solución de la ecuación de Laplace utilizando el Método rápido multipolar, la segunda derivada de las funciones de leyenda asociadas del primer tipo . Específicamente, estoy buscando implementaciones de C o simplemente la relación de recurrencia correcta para usar para escribir la función …

8 special-functions  numerics 

Según Wolfram Alpha y el sistema de álgebra computacional Sage , la siguiente identidad se mantiene: cos(arctan(l1−l2d))=11+(l1−l2)2d2−−−−−−−−√cos⁡(arctan⁡(l1−l2d))=11+(l1−l2)2d2 \cos\left(\arctan\left(\frac{l_1-l_2}{d}\right)\right) = \frac{1}{\sqrt{1 + \frac{(l_1-l_2)^2}{d^2}}} Sin embargo, cuando traté de verificarlo con un ejemplo arbitrario en NumPy, noté una diferencia bastante grande en los valores reales calculados por ambos lados de la identidad. …

8 python  numerics