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...
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...
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...
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,...
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...
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...
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...
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 de y cada 2 por 2 bloque corresponde...
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,...
¿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
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...
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
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,...
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...
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 similar a J n (pero no el mismo,...