Preguntas etiquetadas con stability

Tengo que resolver varios problemas desafiantes de optimización global no convexo. Actualmente uso la Caja de herramientas de optimización de MATLAB (específicamente, fmincon()con algoritmo = 'sqp'), que es bastante eficaz . Sin embargo, la mayor parte de mi código está en Python, y me encantaría...

El teorema de equivalencia de Lax establece que la consistencia y la estabilidad de un esquema numérico para un problema de valor inicial lineal es una condición necesaria y suficiente para la convergencia. Pero para problemas no lineales, los métodos numéricos pueden converger de manera muy...

Estoy tratando de implementar la siguiente función en coma flotante de doble precisión con un error relativo bajo : logsum(x,y)=log(exp(x)+exp(y))logsum(x,y)=log⁡(exp⁡(x)+exp⁡(y))\mathrm{logsum}(x,y) = \log(\exp(x) + \exp(y)) Esto se usa ampliamente en aplicaciones estadísticas para agregar...

Los métodos numéricos para resolver PDEs (u ODE) se dividen en dos grandes categorías: métodos explícitos e implícitos. Los métodos implícitos permiten pasos de tiempo estables más grandes pero requieren más trabajo por paso. Para las PDE hiperbólicas, la sabiduría común es que los métodos...

Me gustaría saber si hay una manera rápida de calcular la distancia euclidiana de dos vectores en Octave. Parece que no hay una función especial para eso, así que ¿debería usar la fórmula con

Me encontré con un comentario desconcertante en el periódico. PJ van der Houwen, El desarrollo de los métodos de Runge-Kutta para ecuaciones diferenciales parciales, Appl. Num. Matemáticas. 20: 261, 1996 En las líneas 8ff en la página 264, van der Houwen escribe: "Para los polinomios de...

Después de hacer algunas matemáticas relacionadas con la estabilidad de los elementos en el problema de 3D Stokes, me sorprendió un poco darme cuenta de que PAG2- P1PAG2-PAG1P_2-P_1 no es estable para una malla tetraédrica arbitraria. Más precisamente, en caso de que tenga un elemento donde todos...

Estoy trabajando en resolver las ecuaciones de poroelasticidad unidimensionales acopladas (modelo de biot), dadas como: ∂−(λ+2μ)∂2u∂x2+∂p∂x=0−(λ+2μ)∂2u∂x2+∂p∂x=0-(\lambda+ 2\mu) \frac{\partial^2 u}{\partial x^2} + \frac{\partial p}{\partial x} =

He implementado un solucionador de Euler hacia atrás en Python 3 (usando numpy). Para mi propia conveniencia y como ejercicio, también escribí una pequeña función que calcula una aproximación de diferencia finita del gradiente para que no siempre tenga que determinar el jacobiano analíticamente...

Los altos flujos de números de Reynolds producen capas límite muy delgadas. Si se usa la resolución de pared en la simulación de Eddy grande, la relación de aspecto puede ser del orden de 10610610^6 . Muchos métodos se vuelven inestables en este régimen porque la constante inf-sup se degrada como...

Suponga que tengo una función de valor real de algunas variables que quiero evaluar numéricamente. En general, la fórmula para puede contener productos, racionales, funciones transitorias, etc., y será demasiado tiempo para investigar analíticamente su estabilidad numérica. O al menos tomará mucho...

En este momento tengo un código que usa el algoritmo de Crank-Nicholson, pero creo que me gustaría pasar a un algoritmo de orden superior para pasar el tiempo. Sé que el algoritmo de Crank-Nicholson es estable en el dominio en el que quiero trabajar, pero me preocupa que otros algoritmos no lo...

He leido eso (1) Las operaciones mal condicionadas deben realizarse antes que las bien condicionadas. Como ejemplo, uno debe calcular x z- yzXz-yzxz-yz como ( x - y) z(X-y)z(x-y)z ya que la resta está mal condicionada mientras que la multiplicación no. Sin embargo, un análisis de error de...

He estado usando SVD de Intel MKL (a dgesvdtravés de SciPy) y noté que los resultados son significativamente diferentes cuando cambio la precisión entre float32y float64cuando mi matriz está mal condicionada / no tiene rango completo. ¿Existe una guía sobre la cantidad mínima de regularización que...

Estoy leyendo un artículo [1] donde resuelven la siguiente ecuación no lineal utilizando métodos de diferencia finita. También analizan la estabilidad de los esquemas utilizando el análisis de estabilidad de Von Neumann. Sin embargo, como los autores se dan cuenta, esto solo es aplicable a PDE...

Antecedentes: El método multiplicador de Lagrange se ha empleado en numerosos campos, como problemas de contacto, interfaces de materiales, transformación de fase, restricciones rígidas o deslizamiento a lo largo de las interfaces. Es bien sabido que una mala elección o diseño del espacio...