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,...
Estoy buscando la literatura más importante y las referencias de diapositivas para modelar los costos de energía actuales y futuros de las operaciones de punto flotante y las transferencias de datos a través de la CPU, la memoria, la red y el almacenamiento. He marcado esta pregunta como un wiki...
Soy nuevo en el paquete PETSc. Tengo un ~ 4000x4000 matriz A en formato de mercado de matriz y quiero obtener PETSc para resolver esto usando múltiples procesadores. Sé cómo resolver el sistema en un único procesador, pero no sé cómo distribuir la matriz y los vectores entre diferentes...
El problema es maxf(x) subject to Ax=bmaxf(x) subject to Ax=b\max f(\mathbf{x}) \text{ subject to } \mathbf{Ax} = \mathbf{b} donde f(x)=∑Ni=11+x4i(∑Ni=1x2i)2−−−−−−−−−−√f(x)=∑i=1N1+xi4(∑i=1Nxi2)2f(\mathbf{x}) = \sum_{i=1}^N\sqrt{1+\frac{x_i^4}{(\sum_{i=1}^{N}x_i^2)^2}} ,...
Al leer hoy sobre la teoría de las formas diferenciales, me dejó impresionado cuánto me recordaba al Método de Volumen Finito (FVM) de segundo orden. Me cuesta entender que pensar de esta manera es trivial o hay alguna conexión más profunda. Bueno, las formas diferenciales sirven para...
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...
Estoy trabajando en un simulador de fluidos 2D usando partículas de vórtice / "vortones" como se describe en Fluid Simulation for Video Games . Lo que creo que es lo mismo que el "método de vórtice discreto". Básicamente, representa el fluido con una colección de partículas con vorticidad definida,...
Actualmente estoy trabajando para resolver sistemas simétricos muy grandes (pero no positivos definidos), generados por algunos algoritmos determinados. Estas matrices tienen una buena dispersión de bloques que se puede usar para la resolución paralela. Pero no puedo decidir si debería usar un...
en matlab, tanto linsolve como mldivide se usan para resolver un sistema de ecuaciones lineales, en todos los casos determinados, sobredeterminados y subdeterminados. Al leer sus documentos, me preguntaba qué diferencias hay entre ellos. ¿Están usando casi los mismos algoritmos de factorización...
Cerrado. Esta pregunta está fuera de tema . Actualmente no está aceptando respuestas. ¿Quieres mejorar esta pregunta? Actualice la pregunta para que sea sobre el tema de Computational Science Stack Exchange. Cerrado hace 2 años . Como dice el título,...
C ++ 11 introduce la semántica de movimiento que puede, por ejemplo, mejorar el rendimiento del código en situaciones donde C ++ 03 necesitaría realizar una construcción de copia o asignación de copia. Este artículo informa que el siguiente código experimenta una velocidad 5 veces mayor cuando se...
Todos estamos familiarizados con los muchos métodos computacionales para resolver el sistema lineal estándar. A x = b .Ax=b. Ax=b. Sin embargo, tengo curiosidad si hay algún método computacional "estándar" para resolver un sistema lineal más general (dimensión finita) de la forma L A = B...
Como todos saben, el flujo incompresible no existe en realidad, es una suposición introducida para simplificar las ecuaciones de gobierno. No podemos aplicar este supuesto directamente. Generalmente el número de Mach (M <0.3 para flujo incompresible), la variación de densidad (variación de...
En el método de soluciones fabricadas (MMS) se postula una solución exacta, se sustituye en las ecuaciones y se calcula el término fuente correspondiente. La solución se usa para la verificación del código. Para las ecuaciones incomprensibles de Navier-Stokes, MMS conduce fácilmente a un término...
Estoy usando el esquema de diferencia finita de Crank-Nicolson para resolver una ecuación de calor 1D. Me pregunto si el principio máximo / mínimo de la ecuación de calor (es decir, que el máximo / mínimo ocurre en la condición inicial o en los límites) también es válido para la solución...
Estoy buscando portar algún código que resuelva un conjunto de ecuaciones diferenciales parciales (PDE) por el método de volumen finito en forma IMPLICIT (para la discretización de tiempo). Como resultado, existe un sistema tridiagonal de ecuaciones en direcciones x, y, z que se maneja mediante el...
Soy nuevo en computación científica. Estoy buscando una implementación Fortran (preferiblemente f90) de un Octree. Mi problema requiere un Octree que divide mi dominio hasta que no haya más de algunas partículas de N (o fuentes donde conozco el valor de densidad que se puede conectar en un método...
Estoy tratando de resolver la ecuación de Poisson en 2D usando el método Discontinuous Galerkin (DG) y la siguiente discretización (tengo un archivo png pero no puedo cargarlo, lo siento): Ecuación: ∇ ⋅ ( κ ∇ T) + f= 0∇⋅(κ∇T)+F=0 0\nabla \cdot( \kappa \nabla T) + f = 0 Nuevas ecuaciones: q= κ ∇...
Uso la función de kernel RBF para implementar un algoritmo de aprendizaje automático basado en kernel (KLPP), la matriz de kernel resultante está extremadamente mal condicionado. El número de condición de la norma L2 vieneKKK 1017-1064K( i , j ) = exp( - ( xyo-
Para la reducción del modelo, quiero calcular los vectores singulares izquierdos asociados a los - digamos 20 - valores singulares más grandes de una matriz , donde y . Desafortunadamente, mi matriz será densa sin ninguna estructura. N ≈ 10 6 k ≈ 10 3 AA∈RN,kA∈RN,kA \in \mathbb...