Efficient Computation of Arbitrary-Order Directional Derivatives in Multiple Directions via Generalized Dual Numbers

तकनीकी सारांश: सामान्यीकृत द्वैत संख्याओं (Generalized Dual Numbers) के माध्यम से स्वेच्छिक-क्रम दिशात्मक अवकलजों (Arbitrary-Order Directional Derivatives) की कुशल गणना

समस्या विवरण
विभिन्न विशिष्ट दिशाओं में स्वेच्छिक-क्रम दिशात्मक अवकलजों की गणना करना बहुरेखीय विश्लेषण (multilinear analysis), उच्च-क्रम टेलर विधियों (high-order Taylor methods), संवेदनशीलता विश्लेषण (sensitivity analysis), और यांत्रिक प्रणालियों के गतिज मॉडलिंग (kinematic modeling) में एक मौलिक आवश्यकता है। जबकि प्रथम- और द्वितीय-क्रम के अवकलज मानक हैं, उभरते अनुप्रयोगों में तेजी से उच्च-क्रम विभेदन (higher-order differentiation) की मांग बढ़ रही है। पारंपरिक दृष्टिकोण महत्वपूर्ण सीमाओं का सामना करते हैं: प्रतीकात्मक विभेदन (symbolic differentiation) बड़े व्यंजकों के लिए गणनात्मक रूप से महंगा और प्रबंधित करने में कठिन हो जाता है; परिमित-अंतर सन्निकटन (finite-difference approximations) ट्रंकेशन और कैंसिलेशन त्रुटियों से ग्रस्त होते हैं; और उच्च-क्रम डेरिवेटिव टेंसर (जैसे पूर्ण हेसियन या उच्च-क्रम टेंसर) का स्पष्ट निर्माण क्रम बढ़ने के साथ तेजी से गणनात्मक रूप से निषेध योग्य (prohibitive) हो जाता है। यद्यपि मौजूदा सामान्यीकृत द्वैत-संख्या प्रारूप एक पुनरावृत्त दिशा $d^n f(v, \dots, v)$ के साथ उच्च-क्रम डेरिवेटिव को कुशलतापूर्वक कंप्यूट करते हैं, वे बिना स्पष्ट टेंसर असेंबली के आवश्यक सामान्य सममित बहुरेखीय रूप (general symmetric multilinear form) $d^n f(v_1, \dots, v_n)$ को सीधे प्रदान नहीं करते हैं।

कार्यप्रणाली
लेखक एक ढांचा प्रस्तावित करते हैं जो सामान्यीकृत द्वैत-संख्या मूल्यांकन को सममित बहुरेखीय रूपों के समावेश-अपवर्जन पुनर्निर्माण (inclusion–exclusion reconstruction) के साथ जोड़ता है।

1. सामान्यीकृत द्वैत संख्याएँ (Generalized Dual Numbers): बीजगणितीय आधार $r = \sum_{k=0}^n r_k \epsilon^k$ के रूप में सामान्यीकृत द्वैत संख्याओं का उपयोग करता है, जहाँ आधार तत्व विशिष्ट गुणन नियमों ($\epsilon^i \epsilon^j = 0$ यदि $i+j > n$, अन्यथा स्केल किया गया $\epsilon^{i+j}$) का पालन करते हैं। एक बिंदु $q$ पर $\epsilon_1 v$ द्वारा विक्षेपित द्वैत दिशा के साथ एक विश्लेतिक फलन $f$ का मूल्यांकन करने पर, $\epsilon^k$ के गुणांक में पुनरावृत्त दिशा $v$ के साथ $k$-वें क्रम का दिशात्मक अवकलज प्राप्त होता है।
2. समावेश-अपवर्जन पुनर्निर्माण (Inclusion–Exclusion Reconstruction): विविध दिशाओं $x_i$ के लिए सामान्य अवकलज $d^n f(x_1, \dots, x_n)$ की गणना करने के लिए, लेखक सममित बहुरेखीय रूप की समरूपता का लाभ उठाते हैं। दिशाओं के उपसमुच्चयों के योगों पर फलन का मूल्यांकन करके, वांछित मिश्रित अवकलज को एक समावेश-अपवर्जन सिद्धांत का उपयोग करके अलग किया जाता है:
   $$d^n f(x_1, \dots, x_n) = \frac{1}{n!} \sum_{k=1}^n (-1)^{n-k} \sum_{S_j \in \mathcal{S}_k} d^n f\left(\sum_{x_p \in S_j} x_p\right)$$
   जहाँ $\mathcal{S}_k$ आकार $k$ के उपसमुच्चय को दर्शाता है। यह समस्या को एक सामान्य $n$-वें क्रम के अवकलज की गणना से पुनरावृत्त-दिशा अवकलजों की एक परिमित संख्या के मूल्यांकन में बदल देता है, जो एकल सामान्यीकृत द्वैत मूल्यांकन के माध्यम से प्राप्त किया जा सकता है।
3. मिश्रित आंशिक अवकलज (Mixed Partial Derivatives): इन्हें मानक आधार वेक्टरों को दिशाओं के रूप में चुनकर पुनः प्राप्त किया जाता है।
4. गतिज मात्राएँ (Kinematic Quantities): यह ढांचा सदिश-मान वाले स्थिति फलनों $r(t) = f(q(t))$ के स्वेच्छिक-क्रम समय अवकलजों की गणना के लिए विस्तारित होता है। वेक्टर दिशात्मक अवकलजों के रूप में कुल अवकलजों को व्यक्त करके और आंशिक बेल पॉलिनोमिअल्स (partial Bell polynomials) द्वारा व्यवस्थित मल्टीवेरिएट फाई डी ब्रूनुइल (Faà di Bruno) विस्तारों का उपयोग करके, यह विधि व्यवस्थित रूप से वेग, त्वरण, जर्क (jerk), स्नैप (snap), और उच्च-क्रम गतिज मात्राओं को बिना प्रतीकात्मक विभेदन के उत्पन्न करती है।
5. टेलर सीरीज़ विधियाँ (TSM): स्वचालित रूप से उत्पन्न उच्च-क्रम समय अवकलजों की उपलब्धता साधारण अंतर समीकरणों (ODEs) के सिस्टम को हल करने के लिए टेलर सीरीज़ विधियों के कार्यान्वयन को सक्षम बनाती है।

प्रमुख योगदान और कार्यान्वयन

• एल्गोरिद्मिक ढांचा: कई विशिष्ट दिशाओं के साथ स्वेच्छей-क्रम दिशात्मक अवकलजों को स्पष्ट रूप से उच्च-क्रम डेरिवेटिव टेंसरों का निर्माण किए बिना कंप्यूट करने की एक विधि।
• एकीकृत गतिज गणना: द्वैत संख्याओं और बेल पॉलिनोमिअल्स कॉम्बिनेटोरिक्स का उपयोग करके स्थिति मैपिंग से सीधे उच्च-क्रम गतिज मात्राओं (मानक वेग और त्वरण से परे) की गणना करने के लिए एक व्यवस्थित प्रक्रिया।
• सॉफ्टवेयर कार्यान्वयन: इस कार्यप्रणाली को आधुनिक फोर्ट्रान (Fortran) में DNAOAD इंफ्रास्ट्रक्चर के भीतर लागू किया गया है, जिसे फोर्ट्रान पैकेज मैनेजर पारिस्थितिकी तंत्र के साथ संगतता के लिए संशोधित किया गया है। कोड ओपन-सोर्स है और Zenodo पर संग्रहीत है।
• दक्षता: यह दृष्टिकोण सभी आंशिक अवकलजों को सूचीबद्ध करने से जुड़े कॉम्बिनेटोरियल विस्फोट (combinatorial explosion) से बचता है, और इसके बजाय चयनित दिशात्मक अवकलजों पर ध्यान केंद्रित करता है।

परिणाम
लेखक कई संख्यात्मक उदाहरणों के माध्यम से कार्यप्रणाली को मान्य करते हैं:

• उच्च-आयामी दिशात्मक अवकलज: आयाम $D \in \{100, 1000, 2000, 3000\}$ के लिए एक साइनुसोइडल समस्या फलन का सातवाँ-क्रम का दिशात्मक अवकलज कंप्यूट किया गया। सभी मामले मानक हार्डवेयर पर 0.1 सेकंड से कम समय में पूरे हुए, जो उस मापनीयता (scalability) को प्रदर्शित करते हैं जहाँ स्पष्ट टेंसर निर्माण अव्यवहार्य होता।
• मिश्रित आंशिक अवकलज: ढांचे ने सफलतापूर्वक विशिष्ट उच्च-क्रम मिश्रित आंशिक अवकलजों (जैसे $\partial^5 f / \partial x^2 \partial y \partial z^2$) को कंप्यूट किया, जो इंडेक्स मैट्रिसेस को निर्दिष्ट करने के लिए एक समर्पित फंक्शन (SPD) का उपयोग करता है।
• गतिज मात्राएँ: एक RCR रोबोट मैनिपुलेटर के लिए पाँचवें-क्रम की गतिज मात्रा की गणना की गई, जो जटिल यांत्रिक प्रणालियों पर इसके अनुप्रयोग को दर्शाती है।
• ODE एकीकरण: साधारण अंतर समीकरणों (ODEs) के सिस्टम को हल करने के लिए एक चौथा-क्रम का टेलर सीरीज़ मेथड लागू किया गया, जिसने अपेक्षित व्यवहारों से मेल खाने वाले संख्यात्मक समाधान प्रस्तुत किए, जो समय-एकीकरण कार्यों के लिए विधि की उपयोगिता को प्रदर्शित करता है।

महत्व और दावे
लेखक इस कार्य को पूर्ण टेंसर-आधारित ऑटोमैटिक डिफरेंशिएशन के प्रतिस्थापन के रूप में नहीं, बल्कि निर्धारित दिशाओं के साथ चयनित उच्च-क्रम डेरिवेटिव को कुशलतापूर्वक कंप्यूट करने के लिए एक लक्षित उपकरण के रूप में प्रस्तुत करते हैं। इसका महत्व निम्नलिखित में निहित है:

1. एक अंतराल को भरना: पुनरावृत्त-दिशा मूल्यांकन से सामान्य सममित बहुरेखीय रूपों को पुनः प्राप्त करने के लिए एक व्यावहारिक तंत्र प्रदान करना, एक ऐसी क्षमता जो पहले कुशल द्वैत-संख्या ढांचों में लुप्त थी।
2. मापनीयता (Scalability): उच्च-आयामी स्थानों में उच्च-क्रम डेरिवेटिव की गणना को सक्षम बनाना जहाँ स्पष्ट टेंसर गणना कॉम्बिनेटोरियल रूप से निषेध योग्य है।
3. व्यवस्थित गतिज विज्ञान (Systematic Kinematics): मानक वेग और त्वरण पदों से आगे बढ़कर, स्वेच्छिक-क्रम गतिज विश्लेषण के लिए एक एकीकृत दृष्टिकोण प्रदान करना।
4. सुलभता: वैज्ञानिक कंप्यूटिंग के लिए उच्च-क्रम विभेदन हेतु एक मजबूत, पारदर्शी और पुनरुत्पादक फोर्ट्रान कार्यान्वयन प्रदान करना, जो उच्च-क्रम AD उपकरणों के वर्तमान परिदृश्य में एक सीमा को संबोधित करता है।

लेखक निष्कर्ष निकालते हैं कि यह फॉर्मूलेशन संख्यात्मक एकीकरण और संवेदनशीलता विश्लेषण के लिए एक कठोर और आसानी से सामान्यीकरण योग्य विकल्प प्रदान करता है, विशेष रूप से उन परिदृश्यों में जिनमें विशिष्ट दिशाओं के साथ उच्च-क्रम विस्तार की आवश्यकता होती है।