क्या scipy.integrate.odeint के लिए आंतरिक गणनाओं को आउटपुट करना संभव है

यह संभव है, आप अपने eom फ़ंक्शन के रिटर्न वैल्यू का उपयोग नहीं कर सकते। इसलिए आपको eom से डेटा की तस्करी के किसी अन्य तरीके की आवश्यकता है। ऐसा करने के कई, कई अलग-अलग तरीके हैं। शायद वैश्विक चर का उपयोग करना सबसे आसान है:

import scipy.integrate as spi

count = 0

def pend(t, y):
    global count

    theta, omega = y
    dydt = [omega, -.25*omega - 5*np.sin(theta)]

    count += 1
    return dydt

sol = spi.solve_ivp(pend, [0, 10], [np.pi - 0.1, 0.0])
print(count)

आउटपुट:

182

साथ ही, ध्यान दें कि मैंने उपरोक्त कोड में odeint के बजाय solve_ivp का उपयोग किया है। odeint डॉक्स कहते हैं कि नया कोड लिखते समय, अब आपको पुराने odeint के बजाय solve_ivp का उपयोग करना चाहिए।

यदि यह मेरा अपना कोड होता, तो शायद मैं अपने फ़ंक्शन के आंशिक संस्करण में एक संचायक ऑब्जेक्ट पास करके कार्य पूरा करता:

class Acc:
    def __init__(self):
        self.x = 0

    def __str__(self):
        return str(self.x)

def pend_partial(acc):
    def pend(t, y):
        theta, omega = y
        dydt = [omega, -.25*omega - 5*np.sin(theta)]

        acc.x += 1
        return dydt
    return pend

count = Acc()
sol = spi.solve_ivp(pend_partial(count), [0, 10], [np.pi - 0.1, 0.0])
print(count)

आउटपुट:

182

हालाँकि, यदि आप केवल एक छोटी स्क्रिप्ट या कुछ और लिख रहे हैं, तो आपको शायद सरल global दृष्टिकोण का उपयोग करना चाहिए। यह इसके लिए एक बहुत अच्छा उपयोग मामला है।