Preguntas etiquetadas con ode

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 hacer la optimización también en …

77 python  optimization  nonlinear-programming  linear-algebra  linear-solver  sparse-matrix  ode  automatic-differentiation  efficiency  matrix  symbolic-computation  pde  multiphysics  operator-splitting  fluid-dynamics  pde  multigrid  pde  hyperbolic-pde  quadrature  precision  numerical-analysis  fluid-dynamics  interpolation  c  ode  testing  data-sets  pde  fluid-dynamics  hyperbolic-pde  pde  finite-element  fluid-dynamics  stability  incompressible  pde  hyperbolic-pde  fluid-dynamics  pde  precision  floating-point  convex-optimization  conditioning  pde  hyperbolic-pde  stability  implicit-methods  explicit-methods  fluid-dynamics  accuracy  pde  hyperbolic-pde  petsc  linear-algebra 

Actualmente estoy buscando métodos paralelos para la integración de ODE. Existe una gran cantidad de literatura nueva y antigua que describe una amplia gama de enfoques, pero no he encontrado encuestas recientes ni artículos de descripción general que describan el tema en general. Existe el libro de Burrage [1], pero …

39 parallel-computing  ode  numerical-analysis 

Esta pregunta se migró de Mathematics Stack Exchange porque se puede responder en Computational Science Stack Exchange. Migrado hace 6 años . Estoy a punto de comenzar a trabajar en una biblioteca de software de solucionadores numéricos de EDO , y estoy luchando sobre cómo formular pruebas para las implementaciones …

26 ode 

¿Cómo se puede resolver numéricamente el problema del cuerpo n gravitatorio en paralelo? ¿Es posible el equilibrio precisión-complejidad? ¿Cómo influye la precisión en la calidad del modelo?

25 algorithms  ode  parallel-computing  complexity  precision 

En este comentario escribí: ... integrador predeterminado de SciPy, que supongo que solo usa métodos simplécticos. en el que me refiero a SciPy's odeint, que utiliza un "método no rígido (Adams)" o un "método rígido (BDF)". Según la fuente : def odeint(func, y0, t, args=(), Dfun=None, col_deriv=0, full_output=0, ml=None, mu=None, …

25 ode  time-integration  integration  differential-equations 

Solo tenía curiosidad por saber por qué los métodos Runge – Kutta de alto orden (es decir, mayores de 4) casi nunca se discuten / emplean (al menos que yo sepa). Entiendo que requiere un mayor tiempo de cómputo por paso (por ejemplo, RK14 con paso incorporado de 12 ° …

17 ode  runge-kutta 

Considere un IVP para el sistema ODE , y ( x 0 ) = y 0 . Más comúnmente, este problema se considera rígido cuando la matriz de Jacobi ∂ fy′=f(x,y)y′=f(x,y)y'=f(x,y)y(x0)=y0y(x0)=y0y(x_0)=y_0tienedosvalores propios con gran parte real negativa y valores propios con parte real negativa muy pequeña (que considero sólo el …

17 ode  stiffness 

Deseo simular el comportamiento de un sistema de doble péndulo. El sistema es un robot manipulador de 2 grados de libertad que no se acciona y, por lo tanto, se comportará principalmente como un doble péndulo afectado por la gravedad. La única diferencia principal con un péndulo doble es que …

15 matlab  ode  numerical-modelling 

¿Cuáles son los métodos más modernos para la solución numérica de EDO con lado derecho discontinuo? Me interesan principalmente las funciones del lado derecho liso por partes, por ejemplo, signo. Estoy tratando de resolver la ecuación de un tipo siguiente: X˙v˙= v= { ( | FexternoEl | - | FfricciónEl …

15 ode 

Me gustaría desarrollar sistemas de resolución de EDO en GPU, en un entorno 'trivialmente paralelo'. Por ejemplo, haciendo un análisis de sensibilidad con 512 conjuntos de parámetros diferentes. Idealmente, quiero resolver ODE con un solucionador de paso de tiempo inteligente y adaptable como CVODE, en lugar de un paso de …

15 ode  gpu 

En la práctica, el tiempo de ejecución de resolver numéricamente un IVP menudo está dominado por la duración de la evaluación del lado derecho (RHS). Supongamos, por lo tanto, que todas las demás operaciones son instantáneas (es decir, sin costo computacional). Si el tiempo de ejecución general para resolver el …

14 ode  numerical-analysis 

Trabajo bastante en biología matemática / epidemiología, donde la mayor parte del trabajo de modelado / ciencia computacional todavía está dominado por conjuntos de EDO, es cierto que a veces son conjuntos bastante elaborados de ellos. Una de las ventajas de estos modelos es que son bastante fáciles de describir …

14 ode  computational-biology 

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 sistema. El método matemático para llegar al estado …

12 matlab  ode  quantum-mechanics 

Me pregunto: ¿cuál es el mejor algoritmo para resolver donde es una matriz real . A no es explícitamente dependiente del tiempo, generalmente escaso pero no necesariamente en bandas. Sus valores propios tienen partes reales no positivas. A también es diagonalizable, pero puede ser demasiado grande para que una diagonalización …

12 linear-algebra  ode 

He estado usando GSL como base de muchas de mis simulaciones, pero es un poco exagerado para mis propósitos y define su propio tipo complejo por razones heredadas. En lugar de codificar mi propio solucionador ODE Runge-Kutta, que probablemente no sería muy eficiente, ¿hay algún solucionador ODE de código abierto …

12 c  ode