Описание тега odeint

Odeint is a modern C++ library for numerically solving Ordinary Differential Equations.
1 ответ

Есть ли инструмент для анализа влияния переменных

Когда вы пишете программу с большим количеством кода, трудно определить, какие значения оказывают большое влияние на ваш конечный результат. В моем случае у меня есть несколько дифференциальных уравнений, которые я решаю с odeint. Потребовалось бы м…
2 ответа

Несколько экземпляров scipy.integrate.ode

Я хотел бы использовать экземпляры scipy.integrate.ode (или scipy.integrate.odeint) в нескольких потоках (по одному на каждое ядро ​​ЦП) для решения нескольких IVP одновременно. Однако в документации сказано: "Этот интегратор не является повторно вх…
15 дек '15 в 14:23
1 ответ

Scipy.ode "Ошибка проверки неоднократно"

Я пытаюсь решить систему связанных сложных ODE в Python, используя scipy.ode с интегратором zvode. Но это сообщение об ошибке появляется, когда я запускаю код. ZVODE-- At T(=R1) and step size H(=R2), the error test failed repeatedly or with abs(H) =…
19 июн '18 в 15:30
1 ответ

Обновление старого кода Python

Довольно нубский вопрос, поэтому, пожалуйста, потерпите меня, потому что я новичок в Python. Ниже приведен код об обыкновенных дифференциальных уравнениях, генерирующих автоматические комбинации ОДУ. Я пытался выполнить этот код на Python 3.6.3 и Sp…
20 ноя '17 в 12:39
0 ответов

Переменная временного шага в ODEINT

Сначала позвольте мне извиниться, если это простой вопрос, или на него уже был дан ответ. Правда в том, что я даже не знаю, как провести поиск по моему вопросу. Допустим, у меня есть следующий набор связанных ODE: dN / dt = (hN (t) + i Q (t)) * G (t…
19 июн '15 в 17:21
1 ответ

TypeError: model() принимает ровно 3 аргумента (дано 5)

odeint с последующей настройкой работает нормально; import numpy as np import matplotlib.pyplot as plt from scipy.integrate import odeint v0 = 10.0 k1 = 0.5 k2 = 0.35 def model(x,t): dx0 = v0 - k1*x[0] dx1 = k1*x[0] - k2*x[1] return [dx0, dx1] time …
29 мар '16 в 04:46
1 ответ

Тяга и жесткие ODE решатели в одном

Я изучаю параметризованную систему жестких ОДУ, и мне нужно решить ее для многих значений ее параметров. Я работаю в C/C++/CUDA и хочу решить эти системы параллельно. Мне было интересно, могу ли я использовать Odeint для этого. Я нашел этот пример, …
28 июн '15 в 09:02
0 ответов

Сложность использования scipy.integrate.odeint для решения радиально-волновой функции водорода

Я пытаюсь использовать scipy.integrate.odeint для определения волновой функции атома водорода, но, похоже, она не дает правильной формы решения. Вот функция, которая должна реализовывать уравнение SE. def eq(y,r,l,E_nl): [u_nl, v_nl] = y #constants …
30 окт '16 в 12:15
1 ответ

Управление выходным пределом в ODEINT Сципи

У меня есть следующая упрощенная настройка псевдокода, которая включает систему ODE, которую я пытаюсь решить с помощью Сципи. from scipy.integrate import odeint def diff_func(y, time, parms): # Do stuff with parms that depends upon y and t. new_par…
1 ответ

Возьмите переменные из функции define в python

У меня есть вопрос по извлечению переменных в def Например, def func(z,W): A = z[1] B = z[2] C = z[3] Var = k*B/C dAdW = A*3 dBdW = B*2/Var dCdW = C**2*Var**2 return dAdW, dBdW, dCdw init = [0,0,1] W = np.linspace(0,100,1000) sol = odeint(func, init…
30 мар '18 в 01:26
0 ответов

Ошибка в смежных взглядах броненосца?

Я столкнулся с очень странным поведением броненосца с использованием метода subvec. Вот короткий фрагмент, иллюстрирующий проблему. Возникает в рамках обычая operator() метод, используемый boost::numeric::odeint, X является вектором столбца длиной 1…
22 дек '17 в 14:37
1 ответ

Кажется, что функция Python odeint не работает

Я хочу изучить движение заряженной частицы во время путешествия через магнитное поле, моделируя его с помощью Python. Я попытался использовать функцию odeint из scipy.integrate, и она не работает так, как я ожидал. Вот что я ожидал, учитывая начальн…
03 мар '18 в 10:54
1 ответ

Параметры шаблона `boost::numeric::odeint::runge_kutta-X` совместимы с CUDA/OpenMP

Подпись типа для класса stepper, который я использую, суммирована здесь: http://www.boost.org/doc/libs/1_56_0/libs/numeric/odeint/doc/html/boost/numeric/odeint/runge_kutta_dopri5.html Это может быть реализовано следующим образом: boost::numeric::ode…
12 апр '16 в 05:27
0 ответов

Использование Eigen::MatrixXd в качестве состояния в odeint error stepper

Я новичок здесь, поэтому я надеюсь правильно задать свой вопрос!:) Я пытаюсь реализовать решатель для диссипативного уравнения Линдблада с odeint в C++, используя собственную библиотеку для обработки матриц. С постоянным размером шагов все работает …
17 июн '16 в 08:25
2 ответа

Скудный одеинт со сложными начальными значениями

Мне нужно решить систему ODE, определяемую комплексным доменом, со сложными начальными значениями. scipy.integrate.odeint не работает на сложных системах. Я рассказываю о том, как разрезать мою систему в реальной и мнимой части и решать отдельно, но…
11 ноя '13 в 15:40
1 ответ

Интеграция Runge kutta для отклика подвеса

Я должен использовать интеграцию Runge Kutta, чтобы решить и построить реакцию ускорения подвески транспортного средства для определенного интервала времени. Когда на графике изображен отклик, он должен иметь несколько колебаний, но из-за демпфирова…
01 ноя '16 в 07:53
1 ответ

Как решить это обыкновенное дифференциальное уравнение с помощью библиотеки C++ Boost odeint

Используя библиотеку Odeint Boost C++, как решить следующие обыкновенные дифференциальные уравнения движения, z '' = -n² · z. Аналитическое решение вышеупомянутого ОДУ z(t) = (z0'/n)·sin(n·t) + z0·cos(n·t). Попытка решения согласно комментарию (17 я…
17 янв '16 в 19:52
1 ответ

Не могу использовать неявный решатель из OdeInt правильно в Xcode 7.3.1

Мне нужно решить некоторые жесткие обыкновенные дифференциальные уравнения, поэтому я узнал о неявных методах, которые предоставляет Odeint. Однако, когда я копирую и соблюдаю в Xcode пример в документации Odeint, http://headmyshoulder.github.io/ode…
07 окт '16 в 22:15
1 ответ

Сценарий Python занимает более 20 часов

У меня есть система из 16 связанных дифференциальных уравнений, и я использую пакет Scipy.integrate.odeint. Код должен работать от 0 до 5 гигов, а затем наносить на график данные, природа дифференциальных уравнений заключается в том, что они чрезвыч…
1 ответ

Интегральная система управления не ведет себя должным образом

Вчера я разместил здесь вопрос: ValueError и odepack.error с использованием integrate.odeint(), на который, как мне показалось, был успешно получен ответ. Однако с тех пор я заметил несколько вещей. При запуске этой программы она не стремится к жела…
08 дек '13 в 00:41