Arhn - архитектура программирования

Остановить scipy odeint при первой ошибке

Редактировать: эта проблема связана с ошибкой, исправленной в Scipy 0.15

Когда я разрабатываю и тестирую код, я могу сделать простую ошибку, например NameError. Когда я использую scipy.integrate.odeint, odeint напечатает сообщение об ошибке, но продолжит интеграцию столько раз, сколько я запрошу, поэтому я получаю много одинаковых сообщений об ошибках. Я предполагаю, что у него такое поведение, поэтому он может продолжать работу при возникновении арифметических ошибок (например, при делении на ноль), но это бесполезное поведение для ошибок программирования.

Есть ли способ заставить scipy остановиться после первого сообщения об ошибке? Было бы лучше, если бы я мог остановить его для ошибок, но не для арифметических исключений.

python scipy odeint
07.05.2013

Ответы:


1

odeint — это оболочка Python для кода C, которая вызывает LSODA (подпрограмму Fortran), которая вызывает код C для вызова моего обратного вызова Python, представляющего dy/dt. LSODA не передает исключение Python, и переход от одного бита кода C к другому сложно реализовать.

Я нашел удовлетворительное решение, просто используя ode вместо odeint. Я считаю, что ode сложнее начать использовать, но он ведет себя правильно, когда возникает исключение в Python. Приведенная ниже функция fake_odeint() — это начало создания функции, которая работает как odeint достаточно хорошо для моих целей, чтобы я мог заменить ее в своем существующем коде. Недостатком использования ode вместо odeint является то, что LSODA вызывается один раз за временной шаг; этот вызов происходит в C с odeint и медленнее в Python с ode.

import numpy as np
from scipy.integrate import ode

def fake_odeint(func, y0, t, Dfun=None):
    ig = ode(func, Dfun)
    ig.set_integrator('lsoda',
                       method='adams')
    ig.set_initial_value(y0, t=0.)
    y = []
    for tt in t:
        y.append(ig.integrate(tt))
    return np.array(y)

Я изучил механизм использования кода Fortran, C и Python, чтобы позволить odeint взаимодействовать с LSODA в этот ТАК вопрос.

25.06.2013
  • Я думаю, что использую вашу функцию fake_odeint, пока ошибка не будет решена! 02.07.2013
    • Новые материалы
    Коллекции публикаций по глубокому обучению
    Последние пару месяцев я создавал коллекции последних академических публикаций по различным подполям глубокого обучения в моем блоге https://amundtveit.com - эта публикация дает обзор 25..
    Представляем: Pepita
    Фреймворк JavaScript с открытым исходным кодом Я знаю, что недостатка в фреймворках JavaScript нет. Но я просто не мог остановиться. Я хотел написать что-то сам, со своими собственными..
    Советы по коду Laravel #2
    1-) Найти // You can specify the columns you need // in when you use the find method on a model User::find(‘id’, [‘email’,’name’]); // You can increment or decrement // a field in..
    Работа с временными рядами спутниковых изображений, часть 3 (аналитика данных)
    Анализ временных рядов спутниковых изображений для данных наблюдений за большой Землей (arXiv) Автор: Рольф Симоэс , Жильберто Камара , Жильберто Кейрос , Фелипе Соуза , Педро Р. Андраде ,..
    3 способа решить квадратное уравнение (3-й мой любимый) -
    1. Методом факторизации — 2. Используя квадратичную формулу — 3. Заполнив квадрат — Давайте поймем это, решив это простое уравнение: Мы пытаемся сделать LHS,..
    Создание VR-миров с A-Frame
    Виртуальная реальность (и дополненная реальность) стали главными модными терминами в образовательных технологиях. С недорогими VR-гарнитурами, такими как Google Cardboard , и использованием..
    Демистификация рекурсии
    КОДЕКС Демистификация рекурсии Упрощенная концепция ошеломляющей О чем весь этот шум? Рекурсия, кажется, единственная тема, от которой у каждого начинающего студента-информатика..