Entonces entiendo (o al menos creo que lo hago) cómo se ejecuta un ciclo V. He escrito en Matlab la versión recursiva 1-D de un ciclo V. Sin embargo, cuando ejecuté mi código para FMG, mi solución no convergía. Creo que mi problema radica en mi comprensión de la parte FMG real. Lo que sé...
Necesito resolver el mismo sistema lineal disperso (300x300 a 1000x1000) con muchos lados derechos (300 a 1000). Además de este primer problema, también me gustaría resolver diferentes sistemas, pero con los mismos elementos distintos de cero (solo valores diferentes), es decir, muchos sistemas...
Estoy usando MATLAB para resolver un problema que implica resolver en cada paso de tiempo, donde cambia con el tiempo. En este momento, estoy logrando esto usando MATLAB :bA x = bAx=b\mathbf{A} \mathbf{x}=\mathbf{b}sib\mathbf{b}mldivide x = A\b Tengo la flexibilidad de hacer tantas calculaciones...
Tengo una colección de modelos computacionales que podrían describirse como autómatas celulares asíncronos. Estos modelos se parecen al modelo de Ising, pero son un poco más complicados. Parece que tales modelos se beneficiarían de ejecutarse en una GPU en lugar de una CPU. Desafortunadamente, no...
Estoy transfiriendo un código existente de MATLAB a C ++ y tengo un sistema lineal para resolver (en lugar de la forma más típica )A x = bxA=bxA=bxA=bAx=bAx=bAx=b La matriz es densa y de forma general, pero no es mayor que 1000x1000. Entonces, en MATLAB, la solución o la notación de barra diagonal...
Estoy simulando un flujo incompresible sobre un cilindro en el número de Reynold de 500. Estoy resolviendo la ecuación de Stoke usando el método de corrección de presión. Mi solución se vuelve inestable después de cierto tiempo (aproximadamente 5 segundos). He intentado refinar mi malla, stepsize...
Existen varias familias de métodos para la cuadratura numérica. Si tengo una clase específica de integrandos, ¿cómo selecciono el método ideal? ¿Cuáles son las preguntas relevantes para formular tanto sobre el integrando (por ejemplo, ¿es suave? ¿Tiene singularidades?) Como sobre el problema...
Supongamos que tengo un problema de valor límite: dud2udx2+dvdx=f in Ωd2udx2+dvdx=f in Ω\frac{d^2u}{dx^2} + \frac{dv}{dx}=f \text{ in } \Omega u=h en ∂Ωdudx+d2vdx2=g in Ωdudx+d2vdx2=g in Ω\frac{du}{dx} +\frac{d^2v}{dx^2} =g \text{ in } \Omega u=h in ∂Ωu=h in ∂Ωu=h \text{ in } \partial\Omega Mi...
Al discretizar con FEM y resolver un problema de reacción-difusión, por ejemplo, con 0 < ε ≪ 1 (perturbación singular), la solución del problema discreto exhibirá típicamente capas oscilatorias cerca del límite. Con Ω = ( 0 , 1 ) , ε = 10 - 5 y elementos finitos lineales, la solución u h...
Mi sistema es un problema de FE simétrico con multiplicadores lagrange (por ejemplo, flujo incompresible de Stokes): (ABBTC)(ABTBC)\begin{pmatrix}A & B^T \\ B & C\end{pmatrix} donde es el caso típico (incluso me he asegurado de que las ecuaciones estén numeradas para que los multiplicadores de...
Hola y gracias por darme tiempo para mirar mi pregunta. Esta es una versión actualizada de mi pregunta que publiqué anteriormente en physics.stackexchange.com Actualmente estoy estudiando un condensador Bose-Einstein Condensador de excitón 2D y tengo curiosidad sobre el estado fundamental de este...
Estoy tratando de averiguar si hay una forma más rápida de calcular todos los valores propios y vectores propios de una matriz de adyacencia muy grande y dispersa que usando scipy.sparse.linalg.eigsh Hasta donde sé, este método solo usa la dispersión y atributos de simetría de la matriz. Una matriz...
En las ecuaciones incompresibles de Navier-Stokes, ρ(ut+(u⋅∇)u)∇⋅u=−∇p+μΔu+f=0ρ(ut+(u⋅∇)u)=−∇p+μΔu+f∇⋅u=0\begin{align*} \rho\left(\mathbf{u}_t + (\mathbf{u} \cdot \nabla)\mathbf{u}\right) &= - \nabla p + \mu\Delta\mathbf{u} + \mathbf{f}\\ \nabla\cdot\mathbf{u} &= 0 \end{align*} el término de...
Estoy interesado en maximizar globalmente una función de muchos ( ) parámetros reales (resultado de una simulación compleja). Sin embargo, la función en cuestión es relativamente costosa de evaluar y requiere aproximadamente 2 días para cada conjunto de parámetros. Estoy comparando diferentes...
Estamos trabajando en un modelo bayesiano para un proceso de espacio-tiempo, y estamos utilizando un muestreador sin giro en U (NUTS) que requiere un modelo para la probabilidad logarítmica y su gradiente con respecto a los parámetros del modelo. Más sucintamente, tenemos una función de...
Sea y f ( → x ) : [ 0 , 1 ] n → C sea una función en estas variables.X⃗ = ( x1, x2, ... , xnorte) ∈ [ 0 , 1 ]nortex→=(x1,x2,…,xn)∈[0,1]n\vec{x} = (x_1, x_2, \dots, x_n) \in [0,1]^nF( x⃗ ) : [ 0 , 1 ]norte→ Cf(x→):[0,1]n→Cf(\vec{x}): [0,1]^n \to \mathbb{C} ¿Existe un esquema recursivo para esta...
Digamos que tengo versiones lentas y rápidas de algún código, y quiero informar un número de aceleración comparando los dos. Ejecuto la versión lenta veces y la versión rápida m veces, produciendo tiempos ( s 1 , ... , s n ) y ( f 1 , ... , f m ) . La forma más simple de producir una aceleración es...
Esta podría ser una pregunta más adecuada para el lado de las Recomendaciones de software de SE, sin embargo, creo que las personas que frecuentan esta parte de SE tienen más probabilidades de poder responder esta pregunta. Estoy buscando una alternativa gratuita (no solo en libertad) a Comsol...
Si tuviera que resolver un problema relativamente pequeño, es decir, un problema que puede manejarse mediante un método directo como LU, ¿el número de condición del operador lineal afecta la precisión de la solución? Uno de los problemas de investigación en los que estoy trabajando se centra en...
Acabo de leer el Capítulo 3 en "Un tutorial de múltiples cuadrículas" de Briggs / Henson / McCormick, enlace . El texto trata sobre ciclos de cuadrícula múltiple como V-cycle, mu-cycle, FMG. Lo que me llamó la atención: en la mayoría de los procedimientos iterativos, uno verifica si ha convergido...