¿Qué puede salir mal cuando se utilizan los métodos de Krylov preconditonado de KSP ( paquete de solucionador lineal de PETSc ) para resolver un sistema lineal disperso como los obtenidos al discretizar y linealizar ecuaciones diferenciales parciales? ¿Qué pasos puedo tomar para determinar qué...
Esta pregunta se migró de Mathematics Stack Exchange porque se puede responder en Computational Science Stack Exchange. Migrado hace 6 años . Estoy a punto de comenzar a trabajar en una biblioteca de software de solucionadores numéricos de EDO , y estoy luchando sobre...
¿Cuál es la biblioteca más rápida para realizar la triangulación delaunay de conjuntos con millones si tiene puntos 3D? ¿También hay versiones de GPU disponibles? Desde el otro lado, tener la teselación de voronoi del mismo conjunto de puntos, ¿ayudaría (en términos de rendimiento) a obtener la...
Uno de los principales problemas con los que tenemos que lidiar en las simulaciones moleculares es el cálculo de las fuerzas dependientes de la distancia. Si podemos restringir las funciones de fuerza y distancia para tener potencias iguales de la distancia de separación , entonces podemos...
Pregunta relacionada: Estado del sistema operativo Mac en informática científica y HPC Un número significativo de paquetes de software en ciencia computacional están escritos en Fortran, y Fortran no va a desaparecer. También se requiere un compilador Fortran para construir otros paquetes de...
Estoy buscando una herramienta para visualizar gráficos de enlaces direccionales muy grandes. Actualmente tengo ~ 2 millones de nodos con ~ 10 millones de bordes. He intentado algunas cosas diferentes, pero la mayoría toma horas incluso para hacer gráficos de 100k nodos Lo que he intentado: pasé...
En este comentario escribí: ... integrador predeterminado de SciPy, que supongo que solo usa métodos simplécticos. en el que me refiero a SciPy's odeint, que utiliza un "método no rígido (Adams)" o un "método rígido (BDF)". Según la fuente : def odeint(func, y0, t, args=(), Dfun=None,...
Necesito evaluar numéricamente la integral a continuación: ∫∞0sinc′(xr)rE(r)−−−−√dr∫0∞sinc′(xr)rE(r)dr\int_0^\infty \mathrm{sinc}'(xr) r \sqrt{E(r)} dr donde E(r)=r4(λκ2+r2−−−−−−√)−ν−5/2K−ν−5/2(λκ2+r2−−−−−−√)E(r)=r4(λκ2+r2)−ν−5/2K−ν−5/2(λκ2+r2)E(r) = r^4 (\lambda\sqrt{\kappa^2+r^2})^{-\nu-5/2}...
¿Qué consideraciones debo hacer al elegir entre BFGS y el gradiente conjugado para la optimización? La función que estoy tratando de ajustar con estas variables son funciones exponenciales; sin embargo, la función objetivo real implica la integración, entre otras cosas, y es muy costosa si eso...
¿Cómo se puede resolver numéricamente el problema del cuerpo n gravitatorio en paralelo? ¿Es posible el equilibrio precisión-complejidad? ¿Cómo influye la precisión en la calidad del modelo?
No entiendo el comportamiento diferente de la ecuación de advección-difusión cuando aplico diferentes condiciones de contorno. Mi motivación es la simulación de una cantidad física real (densidad de partículas) bajo difusión y advección. La densidad de partículas debe conservarse en el interior a...
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...
No estoy ejecutando ningún código paralelo en este momento, pero preveo ejecutar código paralelo en el futuro usando un híbrido de OpenMP y MPI. Los depuradores han sido herramientas invaluables para mí al ejecutar proyectos en serie. ¿Alguien puede recomendar un depurador paralelo (o múltiples...
En términos generales, he escuchado a analistas numéricos pronunciar la opinión de que "Por supuesto, matemáticamente hablando, el tiempo es solo otra dimensión, pero aún así, el tiempo es especial" ¿Cómo justificar esto? ¿En qué sentido es el tiempo especial para la ciencia...
Según tengo entendido, hay dos categorías principales de métodos iterativos para resolver sistemas lineales de ecuaciones: Métodos estacionarios (Jacobi, Gauss-Seidel, SOR, Multigrid) Métodos del subespacio de Krylov (gradiente conjugado, GMRES, etc.) Entiendo que la mayoría de los métodos...
Por lo que puedo decir, los dos grandes marcos genéricos de software de ciencia computacional del Departamento de Energía de EE . UU . Son PETSc y Trilinos . Parecen similares a primera vista, más allá de las diferencias en el lenguaje (C versus C ++). ¿Cuáles son las principales diferencias entre...
Supongamos que estoy trabajando en un código científico en C ++. En una discusión reciente con un colega, se argumentó que las plantillas de expresión podrían ser algo realmente malo, lo que podría hacer que el software sea compilable solo en ciertas versiones de gcc. Supuestamente, este problema...
Descripción del experimento: En la interpolación de Lagrange, la ecuación exacta se muestrea en puntos (orden polinomial ) y se interpola en 101 puntos. Aquí es variado de 2 a 64. Cada vez , y parcelas de error se preparan. Se ve que, cuando la función se muestrea en puntos equidistantes, gotas...
Cuando se pasa de la forma fuerte de una PDE a la forma FEM, parece que siempre se debe hacer esto al indicar primero la forma variacional. Para hacer esto, multiplique la forma fuerte por un elemento en algún espacio (Sobolev) e integre sobre su región. Esto lo puedo aceptar. Lo que no entiendo es...
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...