Recientemente, vi esta pregunta: cómo medir el error de un método de diferencia finita Soy estudiante de ciencias de la simulación y desafortunadamente, para mí, no está del todo claro qué norma usar en qué contexto. Muy a menudo, usamos la norma euclidiana o la norma L2, pero ¿por qué uno elige...
Tengo parte de un problema que se describe en la ecuación de conservación del momento: ∂ρ∂t+ 1pecadoθ∂∂θ( ρ u sinθ ) = 0∂ρ∂t+1sinθ∂∂θ(ρusinθ)=0\frac{\partial \rho}{\partial t} + \frac{1}{\sin\theta} \frac{\partial}{\partial \theta}(\rho u \sin \theta) =0 Donde y ρ = f ( θ , t ) (velocidad...
He visto muchas publicaciones en revistas de física computacional que utilizan diferentes métricas para el desempeño de su código. Especialmente para el código GPGPU, parece que hay una gran variedad de resultados de tiempo que la gente publica. En particular, he visto: Comparaciones de...
Tengo un montón de procesos MPI que intercambian mensajes de solicitud de ida y vuelta. Los procesos no saben qué otros procesos les enviarán mensajes, o cuántos. Dada esta situación, quiero una manera eficiente de saber si todos los demás procesos se consideran terminados enviando mensajes. Esto...
La ecuación de calor se discretiza en el espacio con FV (o FEM), y se obtiene una ecuación semi-discreta (sistema de EDO). Este enfoque, conocido como el método de líneas , permite cambiar fácilmente de una discretización temporal a otra, sin duplicación de código. En particular, puede reutilizar...
Me disculpo de antemano si esta pregunta es tonta. Necesito calcular la raíz de u−f(u)=0u−f(u)=0\begin{equation} u -f(u) =0 \end{equation} Donde es un vector real yf ( u ) es una función con valor de vector real. Comencé con el método de Newton (que funcionó), pero luego me di cuenta de que un...
He estado leyendo sobre la implementación de la multiplicación de matriz densa cuando la matriz no cabe en el caché. Uno de los gráficos que he visto (diapositiva 9 de estas diapositivas ) muestra caídas repentinas en el rendimiento utilizando el algoritmo ingenuo. Estas caídas son de alrededor del...
Recientemente heredé un gran cuerpo de código heredado que resuelve un problema muy rígido y transitorio. Me gustaría demostrar que los tamaños de paso espacial y temporal son lo suficientemente pequeños como para que la naturaleza cualitativa de la solución calculada no cambie si se redujeran. En...
Quiero hacer que mi programa de Python sea rápido usando cython, ¡pero mi ciclo interno todavía está haciendo llamadas lentas de Python al generador de números aleatorios! Hace varios años, alguien con el apoyo de Sage planteó este mismo problema y parecía que no había una buena solución en ese...
Suponiendo que se realiza una simulación numérica directa, ¿cuál es un buen método para estimar el impacto de la pequeña escala a gran escala en la dinámica de fluidos? Por ejemplo, ¿es pertinente comparar dos corridas con diferentes tamaños de cuadrícula o dos corridas con diferentes viscosidades?...
Necesito resolver el sistema de dos ecuaciones diferenciales parciales acopladas numéricamente. ∂x1∂t∂x2∂t=c1∇2x1+f1(x1,x2)=c2∇2x2+K∂x1∂t∂x1∂t=c1∇2x1+f1(x1,x2)∂x2∂t=c2∇2x2+K∂x1∂t\begin{align} \frac{\partial x_1}{\partial t} &= c_1\nabla ^2 x_1 + f_1(x_1,x_2)\\ \frac{\partial x_2}{\partial t} &=...
Estoy tratando de aprender más sobre el mallado de volumen y he decidido implementar un simple mallador de volumen. La estrategia que he elegido es subdividir mi espacio usando un octree, refinado basado en algunos criterios. Me he asegurado de que mi octree esté equilibrado, de modo que la...
¿Existe una implementación de código abierto y segura para subprocesos de solucionadores de optimización convexa en C / C ++? Algunas bibliotecas como NLopt , Ipopt , OPT ++ no cumplen mis requisitos. OPT ++ e Ipopt no son seguros para subprocesos, y NLopt no parece tener una forma específica /...
Tengo dos variables ky tcomo funciones de otras dos variables p1y p2. También sé sus valores máximos. No tengo ninguna expresión analítica para esto. Quiero encontrar los valores de ky tcuáles son los más cercanos a sus valores máximos. ¿Hay alguna manera de optimizar el k = f1(p1, p2)y t = f2(p1,...
Estoy buscando ejemplos de las ecuaciones de Helmholtz y Biharmonic en coordenadas cartesianas con soluciones exactas, para poder comparar mis soluciones numéricas con ellas. Pude encontrar bastantes ejemplos en Internet, donde el problema con las condiciones de contorno se definió con precisión....
Estoy tratando de volver a implementar redes neuronales en Python. Implementé la función de costo y el algoritmo de retropropagación correctamente. Los he comprobado ejecutando su código equivalente de octava. Pero cuando intento usar la scipy.optimize.fmin_cgfunción, las iteraciones tardan mucho...
¿Cuál es la forma más eficiente de obtener una matriz dispersa compleja de mi código Fortran a PETSc? Entiendo que esto depende del problema, por lo que traté de dar tantos detalles relevantes como sea posible a continuación. He estado jugando con el solucionador de valores propios FEAST [1] para...
Implementé un Navier Stokes basado en la discretización FEM y PETSc para resolver el sistema lineal de ecuaciones. Para crear un procedimiento de solución eficiente, sigo el documento "Precondicionamiento eficiente de las ecuaciones de Navier-Stokes linealizadas para flujo incompresible" (Silvester...
Si quiero usar una matriz grande, digamos mpz_t A [100000], recibí "Error de segmentación (núcleo volcado)" durante mi compilación. ¿Hay alguna manera más fácil de resolver
Estoy tratando de entender el algoritmo de renacimiento de la integración de Monte Carlo de VEGAS ( publicación original ( preimpresión de LKlevin) y notas de implementación ). Intentaré explicar primero lo que creo que entendí y luego plantear mis preguntas. Para simplificar, supongamos que...