Perguntas com a marcação «ode»

Eu tenho vários problemas desafiadores de otimização global não convexa para resolver. Atualmente, uso o Optimization Toolbox do MATLAB (especificamente, fmincon()com o algoritmo = 'sqp'), o que é bastante eficaz . No entanto, a maior parte do meu código está em Python, e eu adoraria fazer a otimização em Python …

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 

Atualmente, estou procurando métodos paralelos para integração com ODE. Há muita literatura nova e antiga por aí descrevendo uma ampla variedade de abordagens, mas não encontrei nenhuma pesquisa recente ou artigo de visão geral descrevendo o tópico em geral. Existe o livro de Burrage [1], mas tem quase 20 anos …

39 parallel-computing  ode  numerical-analysis 

Esta pergunta foi migrada do Mathematics Stack Exchange porque pode ser respondida no Computational Science Stack Exchange. Migrou há 6 anos . Estou prestes a começar a trabalhar em uma biblioteca de software de solucionadores numéricos de ODE e estou tendo dificuldades em formular testes para as implementações do solucionador. …

26 ode 

Como o problema gravitacional do corpo-n pode ser resolvido numericamente em paralelo? É possível compensar a complexidade da precisão? Como a precisão influencia a qualidade do modelo?

25 algorithms  ode  parallel-computing  complexity  precision 

Neste comentário eu escrevi: ... integrador SciPy padrão, que suponho usar apenas métodos simpléticos. no qual estou me referindo ao SciPy's odeint, que usa um "método não rígido (Adams)" ou um "método rígido (BDF)". Segundo a fonte : def odeint(func, y0, t, args=(), Dfun=None, col_deriv=0, full_output=0, ml=None, mu=None, rtol=None, atol=None, …

25 ode  time-integration  integration  differential-equations 

Eu estava curioso para saber por que os métodos Runge-Kutta de alta ordem (ou seja, maiores que 4) quase nunca são discutidos / empregados (pelo menos que eu saiba). Entendo que isso requer maior tempo computacional por etapa (por exemplo, RK14 com etapa incorporada de 12ª ordem ), mas existem …

17 ode  runge-kutta 

Considere um IVP para o sistema ODE , . Geralmente esse problema é considerado rígido quando a matriz de Jacobi \ frac {\ parcial f} {\ parcial y} (x_0, y_0) tem ambos os valores próprios com parte real negativa muito grande e valores próprios com parte real negativa muito pequena …

17 ode  stiffness 

Desejo simular o comportamento de um sistema semelhante ao pêndulo duplo. O sistema é um manipulador de robôs de 2 graus de liberdade que não é acionado e, portanto, se comportará principalmente como um pêndulo duplo afetado pela gravidade. A única diferença principal com um pêndulo duplo é que ele …

15 matlab  ode  numerical-modelling 

Quais são os métodos de última geração para solução numérica de EDOs com o lado direito descontínuo? Estou interessado principalmente nas funções do lado direito, como por exemplo, sinal. Estou tentando resolver a equação do seguinte tipo: x˙v˙=v={(|Fexternal|−|Ffriction|)sign(Fexternal)0:|Fexternal|<|Ffriction|:otherwisex˙=vv˙={(|Fexternal|−|Ffriction|)sign⁡(Fexternal):|Fexternal|<|Ffriction|0:otherwise\begin{align*} \dot x &= v\\ \dot v &= \begin{cases} (|F_\text{external}| - |F_\text{friction}|) \mathop{\rm …

15 ode 

Na prática, o tempo de execução de resolver numericamente um IVP é frequentemente dominada pela duração da avaliação do lado direito (RHS) . Vamos, portanto, assumir que todas as outras operações são instantâneas (ou seja, sem custo computacional). Se o tempo de execução geral para resolver o IVP for limitado, …

14 ode  numerical-analysis 

Trabalho bastante pesadamente em biologia / epidemiologia matemática, onde a maior parte do trabalho de modelagem / ciência computacional ainda é dominada por conjuntos de EDOs, reconhecidamente, às vezes, conjuntos bastante elaborados deles. Um dos pontos positivos desses modelos é que eles são fáceis de descrever e replicar. Uma tabela …

14 ode  computational-biology 

Eu gostaria de descobrir como resolver sistemas de ODEs em GPUs, em uma configuração 'trivialmente paralelizável'. Por exemplo, fazendo uma análise de sensibilidade com 512 conjuntos de parâmetros diferentes. Idealmente, eu quero resolver o ODE com um solucionador de tempo inteligente adaptável como CVODE, em vez de um tempo fixo …

14 ode  gpu 

Olá, obrigado por dar tempo para analisar minha pergunta. Esta é uma versão atualizada da minha pergunta, que eu postei anteriormente em physics.stackexchange.com Atualmente, estou estudando um spinor de exciton 2D Bose-Einstein Condensate e estou curioso sobre o estado fundamental deste sistema. O método matemático de chegar ao estado fundamental …

12 matlab  ode  quantum-mechanics 

Gostaria de saber: qual é o melhor algoritmo para resolver dvocêdt= A ududt=Au\begin{equation} \frac{du}{dt} = Au \end{equation} OndeUMAAAé umamatrizn×n × nn×nn\times nreal. A não é explicitamente dependente do tempo, geralmente escasso, mas não necessariamente em faixas. Seus autovalores possuem partes reais não positivas. A também é diagonalizável, mas pode ser …

12 linear-algebra  ode 

Eu tenho usado o GSL como base de muitas de minhas simulações, mas é um pouco exagerado para meus propósitos e define seu próprio tipo complexo por motivos legados. Em vez de codificar meu próprio solucionador Runge-Kutta ODE, o que provavelmente não seria muito eficiente, há algum solucionador ODE de …

12 c  ode