A C++ implementation of the G-Scheme stiff ODE solver with multi-resolution sparse hash-table kernel lookup

कल्पना कीजिए कि आप यह भविष्यवाणी करने की कोशिश कर रहे हैं कि एक जटिल मशीन, जैसे कि कार का इंजन या रॉकेट, गर्म होने पर कैसा व्यवहार करती है। इन मशीनों के भीतर, रसायन बिजली की गति से प्रतिक्रिया कर रहे हैं, जो परमाणुओं के एक अराजक नृत्य को जन्म दे रहे हैं। इस प्रक्रिया को कंप्यूटर पर सिम्युलेट (simulate) करने के लिए, वैज्ञानिक हर एक परमाणु की गति को ट्रैक करने के लिए गणित का उपयोग करते हैं। लेकिन यहाँ एक पेंच है: कुछ परमाणु अविश्वसनीय रूप से तेज़ चलते हैं (जैसे एक हमिंगबर्ड), जबकि अन्य बहुत धीरे चलते हैं (जैसे एक स्लॉथ)। जब आप धीमे वाले परमाणुओं के पथ की गणना करने की कोशिश करते हैं, तो कंप्यूटर तेज़ वाले परमाणुओं के साथ तालमेल बिठाने की कोशिश में फंस जाता है, जिससे पूरी प्रक्रिया धीमी हो जाती है। यह "स्टिफ" (stiff) समीकरणों की समस्या है, और दशकों से, उन्हें हल करने का सबसे अच्छा तरीका बहुत भारी, धीमा, लेकिन बहुत सावधानीपूर्वक काम करने वाले गणितीय उपकरणों का उपयोग करना रहा है, जो हर छोटे कदम पर अपने काम की जाँच करते हैं।

हाल ही में, एक स्मार्ट दृष्टिकोण जिसे "जी-स्कीम" (G-Scheme) कहा गया, का आविष्कार किया गया। हर एक परमाणु की जाँच करने के बजाय, यह एक क्लब के बाउंसर की तरह काम करता है: यह उन तेज़ गति वाले परमाणुओं की पहचान करता है जो बस अपनी जगह पर घूम रहे हैं और उन्हें बाहर इंतज़ार करने के लिए कहता है, और केवल उन "सक्रिय" परमाणुओं पर ध्यान केंद्रित करता है जो वास्तव में कुछ महत्वपूर्ण कर रहे हैं। यह सिमुलेशन को बहुत तेज़ बना देता है। हालाँकि, यह जानने के लिए कि किन परमाणुओं को अनदेखा करना है, कंप्यूटर को हर एक क्षण में एक विशाल, जटिल गणना (एक विशाल पहेली को हल करने की तरह) करनी पड़ती है। यह एक ऐसे बाउंसर की तरह है जो बहुत बुद्धिमान तो है, लेकिन हर आईडी कार्ड की जाँच करने में दस मिनट लगा देता है। परिणाम यह है कि सिमुलेशन पुराने तरीके से तेज़ तो है, लेकिन फिर भी वास्तविक समय (real-time) के उपयोग के लिए बहुत धीमा है।

यहीं पर एक नया पेपर आता है, जो बाउंसर को तुरंत काम करने के लिए एक चतुर ट्रिक पेश करता है। शोधकर्ता, रिकार्डो मालपिका गैलसी और माउरो वलोरानी ने एक नई प्रोग्रामिंग भाषा (C++) में जी-स्कीम का एक सुपर-फास्ट संस्करण बनाया है और इसमें एक "लुकअप टेबल" (lookup table) सिस्टम जोड़ा है। कल्पना कीजिए कि हर बार आईडी कार्ड की जाँच करने के बजाय, बाउंसर के पास एक विशाल, व्यवस्थित फोटो एल्बम है। जब कोई व्यक्ति सामने आता है, तो बाउंसर उसकी आईडी नहीं जाँचता; वह बस उस व्यक्ति का चेहरा देखता है, एल्बम में एक मिलान वाली फोटो ढूँढता है, और तुरंत जान जाता है कि वह कौन है और कौन से नियम लागू होते हैं। यह पेपर दिखाता है कि इस "फोटो एल्बम" (एक मल्टी-रेज़ोल्यूशन स्पार्स हैश टेबल) का उपयोग करके हर बार पहेली को हल करने के बजाय उत्तर का अनुमान लगाने से, कंप्यूटर मौजूदा सर्वोत्तम तरीकों की तुलना में इन रासायनिक प्रतिक्रियाओं को 9.3 गुना तेज़ सिम्युलेट कर सकता है, यहाँ तक कि सैकड़ों अलग-अलग सामग्रियों वाले सबसे जटिल रासायनिक नुस्खों के लिए भी।

समस्या: "स्टिफ" गणितीय ट्रैफिक जाम

रसायन विज्ञान की दुनिया में, विशेष रूप से जब चीजें जल रही होती हैं या विस्फोट हो रहे होते हैं, तो प्रतिक्रियाएं बहुत अलग समय के पैमाने पर होती हैं। कुछ रासायनिक बंधन एक ट्रिलियनवें सेकंड में टूट जाते हैं, जबकि अन्य सेकंडों तक चलते हैं। जब आप इसे सिम्युलेट करने के लिए कंप्यूटर प्रोग्राम लिखने की कोशिश करते हैं, तो कंप्यूटर को बहुत छोटे कदम उठाने पड़ते हैं ताकि वह तेज़ धमाकों को मिस न कर दे। इसे एक "स्टिफ" (stiff) सिस्टम कहा जाता है।

वर्षों तक, इसे संभालने के लिए मानक उपकरण CVODE नामक एक सॉल्वर था। CVODE को एक बहुत ही सावधान, सूक्ष्म अकाउंटेंट (लेखाकार) के रूप में सोचें। यह हर एक संख्या की जाँच करता है, पूरे सिस्टम की पुनर्गणना करता है, और हर कदम पर अपने काम की दोबारा जाँच करता है। यह अविश्वसनीय रूप से सटीक है, लेकिन धीमा है। यदि आप सैकड़ों रासायनिक प्रजातियों (विभिन्न प्रकार के अणुओं) वाले जटिल इंजन का सिमुलेशन करना चाहते हैं, तो CVODE बहुत समय लेता है क्योंकि यह बार-बार वही भारी काम कर रहा होता है।

जी-स्कीम: स्मार्ट बाउंसर

कुछ साल पहले, वैज्ञानिकों ने जी-स्कीम विकसित किया। यह एक अलग प्रकार का सॉल्वर है। एक सावधान अकाउंटेंट होने के बजाय, यह एक स्मार्ट बाउंसर है। यह रासायनिक प्रणाली को देखता है और महसूस करता है, "हे, इनमें से अधिकांश तेज़ गति वाले अणु वास्तव में अपनी जगह पर कंपन कर रहे हैं; वे वास्तव में प्रतिक्रिया के परिणाम को नहीं बदल रहे हैं।"

जी-स्कीम इन "तेज़ लेकिन थके हुए" मोड की पहचान करने के लिए कंप्यूटेशनल सिंगुलर पर्टरबेशन (CSP) नामक एक गणितीय तकनीक का उपयोग करता है। यह प्रभावी रूप से कहता है, "हमें इन तेज़ चालों को स्टेप-बाय-स्टेप ट्रैक करने की ज़रूरत नहीं है। हम बस उनके प्रभाव का अनुमान लगा सकते हैं और उन कुछ 'सक्रिय' अणुओं पर ध्यान केंद्रित कर सकते हैं जो प्रतिक्रिया को आगे बढ़ा रहे हैं।" यह कंप्यूटर को बहुत बड़े कदम उठाने की अनुमति देता है, जिससे वह तेज़ परमाणुओं के ट्रैफिक जाम को छोड़ देता है।

हालाँकि, एक समस्या थी। यह जानने के लिए कि किन अणुओं को अनदेखा करना है, जी-स्कीम को सिमुलेशन के हर एक चरण में एक विशाल, जटिल गणना (एक जैकोबियन मैट्रिक्स का आइजनडिकम्पोजिशन) करनी पड़ती थी। यह एक ऐसे बाउंसर की तरह था जो बहुत बुद्धिमान है लेकिन हर आईडी कार्ड चेक करने में 10 मिनट लेता है। छोटे समस्याओं के लिए, यह ठीक था। लेकिन सैकड़ों प्रजातियों वाले बड़े, वास्तविक रासायनिक मॉडलों के लिए, वह "चेक" इतना लंबा था कि जी-स्कीम पुराने, धीमे CVODE से बहुत अधिक तेज़ नहीं था।

नया समाधान: "फोटो एल्बम" लुकअप

लेखकों ने एक सरल प्रश्न पूछा: क्या हमें वास्तव में हर बार उस विशाल पहेली को हल करने की आवश्यकता है?

उन्होंने महसूस किया कि एक रासायनिक प्रतिक्रिया में, सिस्टम अक्सर एक ही स्थिति से बार-बार गुजरता है। यदि तापमान और दबाव समान हैं, तो "बाउंसर" संभवतः यह निर्णय लेने के बारे में वही निर्णय लेगा कि किन अणुओं को अनदेखा करना है। तो, उन निर्णयों को सहेज कर क्यों नहीं रखा जाता?

उन्होंने एक मल्टी-रेज़ोल्यूशन स्पार्स हैश टेबल पेश किया। इसे एक विशाल, सुपर-व्यवस्थित फोटो एल्बम के रूप में सोचें।

1. सेटअप (ऑफलाइन): सिमुलेशन शुरू करने से पहले, वे एक ट्रेनिंग सेशन चलाते हैं जहाँ वे प्रतिक्रिया का सिमुलेशन करते हैं और हजारों अलग-अलग स्थितियों के लिए "बाउंसर के निर्णयों" (कर्नेल सेट) को इस एल्बम में सहेजते हैं।
2. लुकअप (ऑनलाइन): वास्तविक सिमुलेशन के दौरान, उस विशाल पहेली को हल करने के बजाय, कंप्यूटर वर्तमान स्थिति (तापमान, दबाव, आदि) लेता है, उसे एक सरल कोड में छोटा करता है, और एल्बम में उसे ढूँढता है।
3. मल्टी-रेज़ोल्यूशन ट्रिक: एल्बम केवल एक बड़ी सूची नहीं है। इसमें विवरण के विभिन्न स्तर हैं। यदि कंप्यूटर को "हाई-डेफिनिशन" सेक्शन में सटीक मिलान नहीं मिलता है, तो वह "मीडियम-डेफिनिशन" सेक्शन को चेक करता है, फिर "लो-डेफिनिशन" सेक्शन को। यह सुनिश्चित करता है कि भले ही कंप्यूटर ऐसी स्थिति का सामना करे जो उसने पहले नहीं देखी है, फिर भी वह एक थोड़े अलग राज्य से एक अच्छा-खासा अनुमान पा सके, बजाय इसके कि वह फंस जाए।

परिणाम: सिमुलेशन की गति बढ़ाना

शोधकर्ताओं ने इस प्रणाली का एक नया, हाई-परफॉर्मेंस संस्करण C++ (एक ऐसी भाषा जो अपनी कच्ची गति के लिए जानी जाती है) में बनाया और इसकी तुलना मानक CVODE सॉल्वर से की। उन्होंने n-हेप्टेन (एक ईंधन जिसका उपयोग जेट इंजनों में किया जाता है) के लिए 33 अलग-अलग रासायनिक तंत्रों का उपयोग किया, जो 56 प्रजातियों वाले सरल मॉडल से लेकर 654 प्रजातियों वाले एक विशाल, विस्तृत मॉडल तक विस्तृत हैं।

यहाँ उन्होंने क्या पाया:

• 100% हिट रेट: हर एक टेस्ट में, "फोटो एल्बम" के पास एक उत्तर तैयार था। कंप्यूटर को कभी भी ऑनलाइन उस विशाल पहेली को हल करने के लिए रुकना नहीं पड़ा। उसने 100% बार टेबल से उत्तर प्राप्त किया।
• भारी गति लाभ: नया सिस्टम प्रतिस्पर्धा से काफी तेज़ था।
  • सबसे छोटे मॉडल (56 प्रजातियां) के लिए, यह CVODE से 2.7 गुना तेज़ था।
  • सबसे बड़े, सबसे जटिल मॉडल (654 प्रजातियां) के लिए, यह 9.3 गुना तेज़ था।
  • एक विशिष्ट टेस्ट केस में (80-प्रजाति मॉडल), नए सॉल्वर ने काम 0.028 सेकंड में पूरा कर दिया, जबकि मानक CVODE ने 0.093 सेकंड लिए। पुराने जी-स्कीम (जो अभी भी पहेली हल कर रहा था) की तुलना में, नया वाला 157 गुना तेज़ था।
• सटीकता: लेखकों ने यह जाँचने के लिए सावधानी बरती कि क्या शॉर्टकट लेने से परिणाम गलत हो गए। उन्होंने "इग्निशन डिले" (यह कितनी देर में जलना शुरू होता है) और रसायनों की अंतिम अवस्था को मापा। उन्होंने पाया कि नया तरीका पुराने तरीके जितना ही सटीक है। त्रुटियां इतनी कम थीं कि वे स्वयं रासायनिक मॉडल को सरल बनाने से होने वाली त्रुटियों की तुलना में नगण्य थीं।
• स्केलिंग: जैसे-जैसे रासायनिक मॉडल बड़े होते गए, नया सॉल्वर अपेक्षाकृत अधिक तेज़ होता गया। इसे चलाने में लगने वाला समय CVODE के मुकाबले बहुत धीमी गति से बढ़ा।

यह क्यों मायने रखता है

यह पेपर साबित करता है कि आपको रासायनिक सिमुलेशन में गति और सटीकता के बीच चुनाव करने की आवश्यकता नहीं है। एक स्मार्ट गणितीय रणनीति (जी-स्कीम) को एक चतुर डेटा लुकअप सिस्टम (हैश टेबल) के साथ जोड़कर, उन्होंने एक ऐसा टूल बनाया है जो जटिल दहन प्रक्रियाओं (combustion processes) को लगभग तुरंत सिम्युलेट कर सकता है।

लेखकों ने यह भी दिखाया कि यह नया C++ संस्करण एक "प्रोडक्शन-ग्रेड" टूल है, जिसका अर्थ है कि यह वास्तविक इंजीनियरिंग के लिए पर्याप्त मजबूत है। उन्होंने इस सॉफ़्टवेयर को ओपन-सोर्स भी बनाया है, ताकि अन्य वैज्ञानिक बेहतर इंजन, स्वच्छ ईंधन या सुरक्षित अग्नि शमन प्रणालियों को डिजाइन करने के लिए इसका उपयोग कर सकें। मुख्य बात यह है कि पिछली गणनाओं को दोबारा करने के बजाय उन्हें याद रखकर, हम उन सिमुलेशन को अनलॉक कर सकते हैं जो पहले व्यावहारिक होने के लिए बहुत धीमे थे।