Recientemente, he encontrado un extraño problema con FORTRAN95. Inicialicé las variables X e Y de la siguiente manera: X=1.0 Y=0.1 Luego los agrego e imprimo el resultado: 1.10000000149012 Después de examinar las variables, parece que 0.1 no está representado con doble precisión con total...
θ θ ∼ Pr ( θ | datos ) i i { w i } θ ∼ N ( μ i , σ 2 )Pr ( θ | datos ) = ∑i = 1kwyonorte( μyo, σ2) .Pr(θ|data)=∑i=1kwiN(μi,σ2).\Pr(\theta| \text{data} ) = \sum_{i=1}^k w_i N(\mu_i, \sigma^2).θθ\thetaθ ∼ Pr ( θ | datos )θ∼Pr(θEl |datos)\theta\sim\Pr(\theta|\text{data})yoyoiyoyoi{ wyo}{wyo}\{w_i\}θ ∼...
Tengo una matriz cuadrada simétrica real densa. La dimensión es de aproximadamente 1000x1000. Necesito calcular el primer componente principal y preguntarme cuál sería el mejor algoritmo para hacerlo. Parece que MATLAB usa los algoritmos Arnoldi / Lanczos (para eigs). Pero al leer sobre ellos, no...
Resuelvo muchos PDE numéricamente, pero las matemáticas aplicadas no son mi campo. No he comprendido qué revistas de matemáticas aplicadas debería leer para estar al día con los desarrollos recientes en el campo. ¿Cuáles son buenas publicaciones para leer para mantenerse al día con los desarrollos...
Tengo una EDO: u ( 0 ) = - 1tu′= - 1000 u + s i n ( t )u′=−1000u+sin(t)u'=-1000u+sin(t) u ( 0 ) = - 11000001tu(0 0)=-11000001u(0)=-\frac{1}{1000001} Sé que este ODE particular es rígido, analíticamente. También sé que si usamos un método de paso de tiempo explícito (hacia adelante) (Euler,...
Estoy resolviendo un problema multiescala usando el método heterogéneo multiescala (HMM) . Esencialmente, mi procedimiento particular utiliza el siguiente proceso iterativo: Resuelve muchos sistemas matriciales locales. Calcule un valor de interés a partir de las soluciones de los sistemas...
En el método de campo autoconsistente de Hartree-Fock para resolver la ecuación electrónica de Schroedinger independiente del tiempo, buscamos minimizar la energía del estado fundamental, , de un sistema de electrones en un campo externo con respecto a la elección de los orbitales giratorios, { χ i...
Hace poco hice una pregunta en la misma línea para las matrices ermitas sesgadas. Inspirado por el éxito de esa pregunta, y después de golpearme la cabeza contra la pared durante un par de horas, estoy mirando la matriz exponencial de las matrices asimétricas reales. La ruta para encontrar los...
Ya he desarrollado una solución de trabajo del Método de Elementos Finitos para resolver problemas de transferencia de calor usando GPU y OpenCL usando el método de Gradiente Conjugado. La principal desventaja de este método es la gran demanda de memoria. Además, en el caso de las tarjetas...
Comencé a aprender OpenFOAM del tutorial de Cavity que se proporciona en el sitio web . Al experimentar con diferentes números de Reynolds, en la sección "2.1.8.2 Ejecutar el código", el tutorial dice que vuelva a ejecutar el solucionador porque "Es sensato aumentar el tiempo de solución". Pero...
¿En qué se diferencia la programación geométrica (generalizada) de la programación convexa general? Un programa geométrico puede transformarse en un programa convexo, y generalmente se resuelve mediante un método de punto interior. Pero, ¿cuál es la ventaja de formular directamente el problema...
Quiero saber cuáles de los solucionadores lineales clásicos (por ejemplo, Gauss-Seidel, Jacobi, SOR) están garantizados para converger para el problema donde es positivo semi- definido y, por supuesto,A x = bUNAX=siAx=bb ∈ i m ( A )UNAUNAAb ∈ i m ( A )si∈yometro(UNA)b \in im(A) (El aviso es semi...
Tengo un problema similar en la formulación de esta publicación, con algunas diferencias notables: ¿Qué métodos simples existen para muestrear adaptativamente una función 2D? Como en esa publicación: Tengo una y la evaluación de esta función es algo costosa de calcularF( x , y)F(X,y)f(x,y) A...
Como alguien que tiene que impartir cursos en ciencias computacionales, me enfrento a la antigua pregunta: ¿cómo evalúo la capacidad de los estudiantes para aprender una materia que depende de aplicaciones que son difíciles de evaluar con métodos de evaluación "estándar"? exámenes escritos u...
Sé que podemos usar técnicas de análisis matemático para probar si un IVP o BVP tiene una solución, es única y depende continuamente de los valores límite / iniciales. Para algunos PDE, especialmente los pde no lineales, es muy difícil, si no imposible, demostrar una buena postura. ¿Existe algún...
como sugiere el título, estoy tratando de calcular la integral de una función compacta (polinomio quíntico de Wendland) en un triángulo. Tenga en cuenta que el centro de la función está en algún lugar del espacio tridimensional. Integro esta función en un triángulo arbitrario, pero pequeño (...
¿Cuáles son las formas recomendadas de hacer mínimos cuadrados no lineales, min , con restricciones de recuadro l o j < = p j < = h i j ? Me parece (los tontos se apresuran) que uno podría hacer que las restricciones de la caja sean cuadráticas y minimizar ∑ i e r r i ( p ) 2 + C ∗ ∑ j t u b...
Supongamos que supiera que un generador congruencial lineal genera una secuencia de números aleatorios. Es decir, Xn + 1= ( una Xnorte+ c ) mod mXnorte+1=(unaXnorte+C)modificaciónmetrox_{n+1}=(aX_n+c) \bmod m Si me dan todo el período (o al menos una gran subsecuencia contigua de la misma), ¿cómo...
Sea un dominio Lipschitz convexo poligonalmente limitado en , deje .R 2 f ∈ L 2 ( Ω )ΩΩ\OmegaR2R2\mathbb R^2F∈ L2( Ω )F∈L2(Ω)f \in L^2(\Omega) Entonces, la solución del problema de Dirichlet en , on tiene una solución única en y está bien planteada, es decir, para alguna constante tenemos .Ω traza...
Me pregunto cómo encontrar los valores propios de una matriz dispersa en un intervalo dado [a, b] por método iterativo. Para mi comprensión personal, es más obvio usar el método del subespacio de Krylov para encontrar los valores propios extremos en lugar de los