पैरारियल: Difference between revisions
No edit summary |
No edit summary |
||
(7 intermediate revisions by 3 users not shown) | |||
Line 18: | Line 18: | ||
}}</ref> | }}</ref> | ||
इसे 2001 में [[जैक्स-लुई लायंस]], मैडे और टुरिनिसी द्वारा प्रस्तुत किया गया था। | इसे 2001 में [[जैक्स-लुई लायंस]], मैडे और टुरिनिसी द्वारा प्रस्तुत किया गया था। तब से यह समय एकीकरण में सबसे व्यापक रूप से अध्ययन किए जाने वाले समानांतर प्रणालियों में से एक बन गया है।{{Citation needed|date=November 2018}} | ||
[[File:Parareal.svg|thumb|342x342px|पैरारियल में पहले पुनरावृत्ति का चित्रण (मूल संस्करण से अनुकूलित)<ref>{{Cite journal|last1=Pentland|first1=Kamran|last2=Tamborrino|first2=Massimiliano|last3=Samaddar|first3=Debasmita|last4=Appel|first4=Lynton|date=2022|title=Stochastic parareal: an application of probabilistic methods to time-parallelization|journal=SIAM Journal on Scientific Computing|language=en|pages=S82–S102|doi=10.1137/21M1414231|s2cid=235485544 }}</ref>).]] | [[File:Parareal.svg|thumb|342x342px|पैरारियल में पहले पुनरावृत्ति का चित्रण (मूल संस्करण से अनुकूलित)<ref>{{Cite journal|last1=Pentland|first1=Kamran|last2=Tamborrino|first2=Massimiliano|last3=Samaddar|first3=Debasmita|last4=Appel|first4=Lynton|date=2022|title=Stochastic parareal: an application of probabilistic methods to time-parallelization|journal=SIAM Journal on Scientific Computing|language=en|pages=S82–S102|doi=10.1137/21M1414231|s2cid=235485544 }}</ref>).]] | ||
== समानांतर | == समय में समानांतर एकीकरण विधियां == | ||
उदाहरण के विपरीत रंज-कुट्टा या [[रैखिक मल्टीस्टेप विधि|मल्टी-स्टेप विधियां]], पैरारियल में कुछ गणनाएं [[समानांतर कंप्यूटिंग]] में की जा सकती हैं और इसलिए पैरारियल समय एकीकरण विधि में समानांतर का एक उदाहरण है। जबकि ऐतिहासिक रूप से आंशिक अंतर समीकरणों के संख्यात्मक आंशिक अंतर समीकरणों को समानांतर करने के अधिकांश प्रयास स्थानिक विवेकीकरण पर केंद्रित हैं, [[एक्सास्केल कंप्यूटिंग]] की चुनौतियों को देखते हुए, [[अस्थायी विवेक|अस्थायी विवेकीकरण]] के समानांतर प्रणालियों को [[संख्यात्मक विश्लेषण सॉफ्टवेयर की सूची]] में समरूपता बढ़ाने के संभावित विधि के रूप में पहचाना गया है।<ref>{{Cite report | उदाहरण के विपरीत रंज-कुट्टा या [[रैखिक मल्टीस्टेप विधि|मल्टी-स्टेप विधियां]], पैरारियल में कुछ गणनाएं [[समानांतर कंप्यूटिंग]] में की जा सकती हैं और इसलिए पैरारियल समय एकीकरण विधि में समानांतर का एक उदाहरण है। जबकि ऐतिहासिक रूप से आंशिक अंतर समीकरणों के संख्यात्मक आंशिक अंतर समीकरणों को समानांतर करने के अधिकांश प्रयास स्थानिक विवेकीकरण पर केंद्रित हैं, [[एक्सास्केल कंप्यूटिंग]] की चुनौतियों को देखते हुए, [[अस्थायी विवेक|अस्थायी विवेकीकरण]] के समानांतर प्रणालियों को [[संख्यात्मक विश्लेषण सॉफ्टवेयर की सूची]] में समरूपता बढ़ाने के संभावित विधि के रूप में पहचाना गया है।<ref>{{Cite report | ||
Line 64: | Line 64: | ||
| citeseerx = 10.1.1.154.6042 | | citeseerx = 10.1.1.154.6042 | ||
}}</ref> | }}</ref> | ||
समय में मल्टीग्रिड के साथ-साथ समय एकीकरण के लिए मल्टीपल शूटिंग को अपनाने के दोनों विचार 1980 और 1990 के दशक से चले आ रहे हैं।<ref>{{cite book | |||
| last = Hackbusch | | last = Hackbusch | ||
| first = Wolfgang | | first = Wolfgang | ||
Line 86: | Line 87: | ||
| doi = 10.1016/S0167-8191(06)80013-X | | doi = 10.1016/S0167-8191(06)80013-X | ||
}}</ref> | }}</ref> | ||
पैरारियल व्यापक रूप से अध्ययन की जाने वाली विधि है और विभिन्न अनुप्रयोगों के लिए इसका उपयोग और संशोधन किया गया है।<ref>{{cite conference | पैरारियल व्यापक रूप से अध्ययन की जाने वाली विधि है और विभिन्न अनुप्रयोगों के लिए इसका उपयोग और संशोधन किया गया है।<ref>{{cite conference | ||
| last = Gander | | last = Gander | ||
Line 99: | Line 101: | ||
| doi-access = free | | doi-access = free | ||
}}</ref> | }}</ref> | ||
प्रारंभिक मूल्य समस्याओं के समाधान को समानांतर करने के विचार और भी | |||
प्रारंभिक मूल्य समस्याओं के समाधान को समानांतर करने के विचार और भी प्राचीन हैं: समय में समानांतर एकीकरण विधि का प्रस्ताव करने वाला पहला पेपर 1964 में सामने आया था।<ref> | |||
{{cite journal | {{cite journal | ||
| last = Nievergelt | | last = Nievergelt | ||
Line 115: | Line 118: | ||
===समस्या=== | ===समस्या=== | ||
लक्ष्य प्रपत्र की प्रारंभिक मूल्य समस्या | लक्ष्य प्रपत्र की प्रारंभिक मूल्य समस्या का समाधान करना है | ||
<math> \frac{\mathrm{d} u}{\mathrm{d}t} = f(t,u) \quad \text{over} \quad t \in [t_0,T] \quad \text{with} \quad u(t_0) = u^0.</math> | <math> \frac{\mathrm{d} u}{\mathrm{d}t} = f(t,u) \quad \text{over} \quad t \in [t_0,T] \quad \text{with} \quad u(t_0) = u^0.</math> | ||
दाहिने हाथ की ओर <math>f</math> को एक सुचारू (संभवतः अरेखीय) फलन माना जाता है। यह रेखा दृष्टिकोण की विधि में आंशिक अंतर समीकरण के स्थानिक विवेक के अनुरूप भी हो सकता है। हम इस समस्या को समान दूरी वाले बिंदु <math>(t_0,t_1,\ldots,t_N)</math>, जहां <math>t_{j+1} = t_j + \Delta T </math> और <math>\Delta T = (T-t_0)/N</math>, <math>N+1</math> के अस्थायी आधार पर समाधान करना चाहते हैं। इस विवेक को आगे बढ़ाते हुए हमें एक विभाजित समय अंतराल प्राप्त होता है जिसमें <math>j = 0,\ldots,N-1 </math> के लिए समय स्लाइस <math>[t_j,t_{j+1}] </math> सम्मिलित होता है। | |||
इसका उद्देश्य सीरियल टाइम-स्टेपिंग विधि (जैसे रनगे-कुट्टा) का उपयोग करके त्रुटिहीन समाधान <math>u(t_j)</math> के लिए संख्यात्मक अनुमान <math>U_j</math> की गणना करना है जिसमें उच्च संख्यात्मक शुद्धता (और इसलिए उच्च कम्प्यूटेशनल लागत) है। हम इस विधि को फाइन सॉल्वर <math>\mathcal{F}</math> के रूप में संदर्भित करते हैं, जो समय <math>t_j</math> पर प्रारंभिक मान <math>U_j</math> को समय <math>t_{j+1}</math> पर टर्मिनल मान <math>U_{j+1}</math> तक प्रसारित करता है। लक्ष्य <math>\mathcal{F} </math> का उपयोग करके समाधान (उच्च संख्यात्मक शुद्धता के साथ) की गणना करना है, जो हमें प्राप्त होता है | |||
<math> U_{j+1} = \mathcal{F}(t_j,t_{j+1},U_j), \quad \text{where} \quad U_0 = u^0.</math> | <math> U_{j+1} = \mathcal{F}(t_j,t_{j+1},U_j), \quad \text{where} \quad U_0 = u^0.</math> | ||
इस (और सबसे पहले समानांतर में | इस (और सबसे पहले समानांतर में समाधान करने का प्रयास करने का कारण) समाधान के साथ समस्या यह है कि वास्तविक समय में गणना करना कम्प्यूटेशनल रूप से संभव नहीं है। | ||
===यह कैसे काम करता है=== | ===यह कैसे काम करता है=== | ||
प्रारंभिक मूल्य समस्या | प्रारंभिक मूल्य समस्या का समाधान करने के लिए एकल प्रोसेसर का उपयोग करने के अतिरिक्त (जैसा कि पारंपरिक समय-चरण विधियों के साथ किया जाता है), पैरारियल <math> N</math> प्रोसेसर का उपयोग करता है। इसका उद्देश्य समानांतर में <math>N </math> छोटी प्रारंभिक मूल्य समस्याओं (प्रत्येक समय स्लाइस पर एक) का समाधान करने के लिए <math> N</math> प्रोसेसर का उपयोग करना है। उदाहरण के लिए, [[संदेश पासिंग इंटरफ़ेस]] आधारित कोड में, <math>N </math> प्रक्रियाओं की संख्या होगी, जबकि [[ओपनएमपी]] आधारित कोड में, <math>N</math> [[थ्रेड (कंप्यूटिंग)]] की संख्या के बराबर होगी। | ||
पैरारियल इस प्रारंभिक मूल्य समस्या को समानांतर में समाधान करने के लिए दूसरी टाइम-स्टेपिंग विधि का उपयोग करता है, जिसे कोर्स सॉल्वर <math>\mathcal{G}</math> के रूप में जाना जाता है। | |||
कोर्स सॉल्वर ठीक सॉल्वर की तरह ही काम करता है, लंबाई <math>\Delta T</math> के समय अंतराल पर प्रारंभिक मूल्य का प्रचार करता है, चूंकि यह <math>\mathcal{F}</math> (और इसलिए बहुत कम कम्प्यूटेशनल लागत पर) की तुलना में बहुत कम संख्यात्मक शुद्धता पर ऐसा करता है। एक कोर्स सॉल्वर का होना जो कि फाइन सॉल्वर की तुलना में बहुत कम कम्प्यूटेशनल रूप से एक्सपेंसिव है, पैरारियल के साथ समानांतर गति प्राप्त करने की कुंजी है। | |||
अब से, हम समय | अब से, हम समय <math> t_j</math> और पुनरावृत्ति <math> k</math> पर पैरारियल समाधान को <math> U^k_j</math> द्वारा निरूपित करते हैं। | ||
====शून्य पुनरावृत्ति==== | ====शून्य पुनरावृत्ति==== | ||
सबसे पहले, | सबसे पहले, समाधान के अनुमानित प्रारंभिक अनुमान की गणना करने के लिए कोर्स सॉल्वर को पूरे समय अंतराल <math>[t_0,T] </math> पर क्रमिक रूप से चलाएं: | ||
<math>U^0_{j+1} = \mathcal{G}(t_j,t_{j+1},U^0_j), \quad j=0,\ldots,N-1. </math> | <math>U^0_{j+1} = \mathcal{G}(t_j,t_{j+1},U^0_j), \quad j=0,\ldots,N-1. </math> | ||
Line 141: | Line 145: | ||
====अनुवर्ती पुनरावृत्तियाँ==== | ====अनुवर्ती पुनरावृत्तियाँ==== | ||
इसके बाद, सबसे अद्यतित समाधान मानों से, समानांतर में, प्रत्येक समय स्लाइस पर | इसके बाद, सबसे अद्यतित समाधान मानों से, समानांतर में, प्रत्येक समय स्लाइस पर अच्छे सॉल्वर चलाएं: | ||
<math> \mathcal{F}(t_j, t_{j+1},U^{k-1}_j), \quad j=0, \ldots, N-1.</math> | <math> \mathcal{F}(t_j, t_{j+1},U^{k-1}_j), \quad j=0, \ldots, N-1.</math> | ||
अब प्रेडिक्टर-करेक्टर का उपयोग करके क्रमिक रूप से पैरारियल समाधान मानों को अपडेट करें: | अब प्रेडिक्टर-करेक्टर का उपयोग करके क्रमिक रूप से पैरारियल समाधान मानों को अपडेट करें: | ||
<math> U_{j+1}^{k} = \mathcal{G}(t_j, t_{j+1},U^{k}_j) + \mathcal{F}(t_j, t_{j+1},U^{k-1}_j) - \mathcal{G}(t_j, t_{j+1},U^{k-1}_j), \quad j=0, \ldots, N-1.</math> | <math> U_{j+1}^{k} = \mathcal{G}(t_j, t_{j+1},U^{k}_j) + \mathcal{F}(t_j, t_{j+1},U^{k-1}_j) - \mathcal{G}(t_j, t_{j+1},U^{k-1}_j), \quad j=0, \ldots, N-1.</math> | ||
इस स्तर पर, कोई यह निर्धारित करने के लिए स्टॉपिंग मानदंड का उपयोग कर सकता है कि क्या समाधान मान अब प्रत्येक पुनरावृत्ति में नहीं बदल रहे हैं। उदाहरण के लिए, कोई इसकी जांच करके यह जांच सकता है कि क्या | इस स्तर पर, कोई यह निर्धारित करने के लिए स्टॉपिंग मानदंड का उपयोग कर सकता है कि क्या समाधान मान अब प्रत्येक पुनरावृत्ति में नहीं बदल रहे हैं। उदाहरण के लिए, कोई इसकी जांच करके यह जांच सकता है कि क्या | ||
<math>| U^{k}_j - U^{k-1}_j | < \varepsilon \quad \forall \ j \leq N,</math> | <math>| U^{k}_j - U^{k-1}_j | < \varepsilon \quad \forall \ j \leq N,</math> | ||
और कुछ सहनशीलता <math>\varepsilon > 0</math> | और कुछ सहनशीलता <math>\varepsilon > 0</math> है। यदि यह मानदंड संतुष्ट नहीं है, तो बाद के पुनरावृत्तियों को समानांतर में फाइन सॉल्वर और फिर भविष्यवक्ता-सुधारक को लागू करके चलाया जा सकता है। चूँकि, एक बार जब मानदंड संतुष्ट हो जाता है, तो कहा जाता है कि एल्गोरिदम <math> k \leq N </math> पुनरावृत्तियों में परिवर्तित हो गया है। ध्यान दें कि अन्य रोक मानदंड उपस्थित हैं और पैरारियल में सफलतापूर्वक परीक्षण किया गया है। | ||
====टिप्पणियाँ==== | ====टिप्पणियाँ==== | ||
पैरारियल को उस समाधान को पुन: प्रस्तुत करना चाहिए जो फाइन सॉल्वर के क्रमिक अनुप्रयोग द्वारा प्राप्त किया | पैरारियल को उस समाधान को पुन: प्रस्तुत करना चाहिए जो फाइन सॉल्वर के क्रमिक अनुप्रयोग द्वारा प्राप्त किया गया है और अधिकतम <math>N</math> पुनरावृत्तियों में परिवर्तित होगा।<ref name="gander2007" /> चूँकि, पैरारियल को स्पीडअप प्रदान करने के लिए, इसे समय स्लाइस की संख्या अर्थात। <math>k \ll N</math> की तुलना में काफी कम संख्या में पुनरावृत्तियों में परिवर्तित करना होगा। | ||
पैरारियल पुनरावृत्ति में, | पैरारियल पुनरावृत्ति में, <math>\mathcal{F}(t_j, t_{j+1},U^{k-1}_j)</math> का कम्प्यूटेशनल रूप से एक्सपेंसिव मूल्यांकन <math>N</math> प्रसंस्करण इकाइयों पर समानांतर में किया जा सकता है। इसके विपरीत, <math>\mathcal{G}(t_j, t_{j+1},U^{k}_j)</math> पर <math>U^{k}_{j+1}</math> की निर्भरता का अर्थ है कि कोर्स सुधार की गणना क्रमिक क्रम में की जानी है। | ||
इसके विपरीत, | |||
सामान्यतः, रंज-कुट्टा विधि का कुछ रूप कोर्स और बारीक इंटीग्रेटर दोनों के लिए चुना जाता है, जहां <math>\mathcal{G}</math> कम क्रम का हो सकता है और <math>\mathcal{F}</math> की तुलना में बड़े समय चरण का उपयोग कर सकता है। | |||
[[File:Parareal Animation.ogg|thumb|216x216px|पैरारियल एल्गोरिथम का विज़ुअलाइज़ेशन। यहां | यदि प्रारंभिक मूल्य समस्या पीडीई के विवेकाधिकार से उत्पन्न होती है, तो <math>\mathcal{G}</math> एक कोर्स स्थानिक विवेकीकरण का भी उपयोग कर सकता है, किन्तु यह अभिसरण पर नकारात्मक प्रभाव डाल सकता है जब तक कि उच्च क्रम इंटरपोलेशन का उपयोग नहीं किया जाता है।<ref>{{Cite journal|last=Ruprecht|first=Daniel|date=2014-12-01|title=स्थानिक मोटेपन के साथ पैरारियल का अभिसरण|journal=Proceedings in Applied Mathematics and Mechanics|language=en|volume=14|issue=1|pages=1031–1034|doi=10.1002/pamm.201410490|s2cid=26356904 |issn=1617-7061|url=http://eprints.whiterose.ac.uk/90536/1/paper.pdf}}</ref> | ||
[[File:Parareal Animation.ogg|thumb|216x216px|पैरारियल एल्गोरिथम का विज़ुअलाइज़ेशन। यहां कोर्स प्रचारक को लेबल किया गया है <math>\bar{\varphi}</math> जबकि अच्छे प्रचारक का लेबल लगा हुआ है <math>\varphi</math>.]] | |||
===स्पीडअप=== | ===स्पीडअप=== | ||
कुछ मान्यताओं के | कुछ मान्यताओं के अनुसार, पैरारियल की गति के लिए एक सरल सैद्धांतिक मॉडल प्राप्त किया जा सकता है।<ref>{{cite journal | ||
| last = Minion | | last = Minion | ||
| first = Michael L. | | first = Michael L. | ||
Line 176: | Line 182: | ||
| doi-access= free | | doi-access= free | ||
}}</ref> | }}</ref> | ||
चूँकि अनुप्रयोगों में ये धारणाएँ बहुत अधिक प्रतिबंधात्मक हो सकती हैं, फिर भी मॉडल पैरारियल के साथ स्पीडअप प्राप्त करने में सम्मिलित ट्रेडऑफ़ को दर्शाने के लिए उपयोगी है। | |||
इन दो मान्यताओं के | सबसे पहले, मान लें कि हर बार स्लाइस <math>[t_j, t_{j+1}]</math> में फाइन इंटीग्रेटर के बिल्कुल <math>N_f</math> चरण और कोर्स इंटीग्रेटर के <math>N_c</math> चरण होते हैं। इसमें विशेष रूप से यह धारणा सम्मिलित है कि सभी समय स्लाइस समान लंबाई के होते हैं और कोर्स और फाइन इंटीग्रेटर दोनों पूर्ण सिमुलेशन पर एक स्थिर चरण आकार का उपयोग करते हैं। दूसरा, क्रमशः सूक्ष्म और कोर्स प्रणालियों के एक चरण के लिए आवश्यक कंप्यूटिंग समय को <math>\tau_f</math> और <math>\tau_c</math> द्वारा निरूपित करें, और मान लें कि दोनों स्थिर हैं। यह सामान्यतः बिल्कुल सच नहीं है जब एक अंतर्निहित विधि का उपयोग किया जाता है, क्योंकि तब रनटाइम पुनरावृत्त सॉल्वर द्वारा आवश्यक पुनरावृत्तियों की संख्या के आधार पर भिन्न होता है। | ||
इन दो मान्यताओं के अनुसार, <math>P</math> टाइम स्लाइस पर एकीकृत करने वाली फाइन विधि के रनटाइम को इस प्रकार मॉडल किया जा सकता है | |||
<math> c_{\text{fine}} = N N_{f} \tau_f. </math> | <math> c_{\text{fine}} = N N_{f} \tau_f. </math> | ||
<math>P</math> प्रसंस्करण इकाइयों का उपयोग करके और <math>k</math> पुनरावृत्तियों को निष्पादित करके पैरारियल का रनटाइम है | |||
<math> c_{\text{parareal}} = (k+1) N N_{c} \tau_c + k N_{f} \tau_f. </math> | <math> c_{\text{parareal}} = (k+1) N N_{c} \tau_c + k N_{f} \tau_f. </math> | ||
पैरारियल का स्पीडअप तो है | पैरारियल का स्पीडअप तो है | ||
<math> S_{p} = \frac{c_{\text{fine}}}{c_{\text{parareal}}} = \frac{1}{ (k+1) \frac{N_c}{N_f} \frac{\tau_c}{\tau_f} + \frac{k}{N}} \leq \min\left\{ \frac{N_f \tau_f}{N_c \tau_c}, \frac{N}{k} \right\}.</math> | <math> S_{p} = \frac{c_{\text{fine}}}{c_{\text{parareal}}} = \frac{1}{ (k+1) \frac{N_c}{N_f} \frac{\tau_c}{\tau_f} + \frac{k}{N}} \leq \min\left\{ \frac{N_f \tau_f}{N_c \tau_c}, \frac{N}{k} \right\}.</math> | ||
ये दो सीमाएँ | |||
विशेष रूप से, | ये दो सीमाएँ कोर्स विधियों को चुनने में किए जाने वाले व्यापार को दर्शाती हैं: एक ओर, इसे सस्ता होना होगा और/या पहली सीमा को जितना संभव हो उतना बड़ा बनाने के लिए बहुत बड़े समय के कदम का उपयोग करना होगा, दूसरी ओर दूसरी सीमा को बड़ा रखने के लिए पुनरावृत्तियों की संख्या <math>k</math> को कम रखना होगा। | ||
विशेष रूप से, पैरारियल की समानांतर दक्षता सीमित है | |||
<math> E_{p} = \frac{S_p}{N} \leq \frac{1}{k},</math> | <math> E_{p} = \frac{S_p}{N} \leq \frac{1}{k},</math> | ||
यह आवश्यक पुनरावृत्तियों की संख्या के व्युत्क्रम से है। | यह आवश्यक पुनरावृत्तियों की संख्या के व्युत्क्रम से है। | ||
=== काल्पनिक | === काल्पनिक आइजेनवैल्यूज़ के लिए अस्थिरता === | ||
पैरारियल के वेनिला संस्करण में काल्पनिक [[आइगेनवैल्यूज़ एवं आइगेनवेक्टर्स]] के साथ समस्याएं हैं।<ref name="gander2007" />यह | पैरारियल के वेनिला संस्करण में काल्पनिक [[आइगेनवैल्यूज़ एवं आइगेनवेक्टर्स]] के साथ समस्याएं हैं।<ref name="gander2007" /> यह सामान्यतः केवल अंतिम पुनरावृत्तियों की ओर ही परिवर्तित होता है, अर्थात जैसे-जैसे <math>k</math> <math>N</math> क्वे निकट पहुंचता है, और स्पीडअप <math> S_p</math> सदैव से छोटा होता है। तो या तो पुनरावृत्तियों की संख्या छोटी है और पैरारियल अस्थिर है या, यदि <math>k</math> पैरारियल को स्थिर बनाने के लिए पर्याप्त बड़ा है, कोई गति संभव नहीं है। इसका यह भी अर्थ है कि हाइपरबोलिक आंशिक अंतर समीकरण समीकरणों के लिए पैरारियल सामान्यतः अस्थिर है।<ref>{{Cite book|title=पैरारियल एल्गोरिथम की स्थिरता|last1=Staff|first1=Gunnar Andreas|last2=Rønquist|first2=Einar M.|date=2005-01-01|publisher=Springer Berlin Heidelberg|isbn=9783540225232|editor-last=Barth|editor-first=Timothy J.|series=Lecture Notes in Computational Science and Engineering|pages=449–456|language=en|doi=10.1007/3-540-26825-1_46|editor-last2=Griebel|editor-first2=Michael|editor-last3=Keyes|editor-first3=David E.|editor-last4=Nieminen|editor-first4=Risto M.|editor-last5=Roose|editor-first5=Dirk|editor-last6=Schlick|editor-first6=Tamar|editor-last7=Kornhuber|editor-first7=Ralf|editor-last8=Hoppe|editor-first8=Ronald|editor-last9=Périaux|editor-first9=Jacques}}</ref> भले ही गैंडर और वांडेवेले द्वारा औपचारिक विश्लेषण केवल निरंतर गुणांक के साथ रैखिक समस्याओं को कवर करता है, समस्या तब भी उत्पन्न होती है जब पैरारियल को गैर-रेखीय नेवियर-स्टोक्स समीकरणों पर लागू किया जाता है जब चिपचिपाहट गुणांक बहुत छोटा हो जाता है और [[रेनॉल्ड्स संख्या]] बहुत बड़ी हो जाती है।<ref>{{Cite book|title=रेनॉल्ड्स संख्या के आधार पर नेवियर-स्टोक्स समीकरणों के लिए पैरारियल का अभिसरण|last1=Steiner|first1=Johannes|last2=Ruprecht|first2=Daniel|last3=Speck|first3=Robert|last4=Krause|first4=Rolf|date=2015-01-01|publisher=Springer International Publishing|isbn=9783319107042|editor-last=Abdulle|editor-first=Assyr|series=Lecture Notes in Computational Science and Engineering|pages=195–202|language=en|doi=10.1007/978-3-319-10705-9_19|editor-last2=Deparis|editor-first2=Simone|editor-last3=Kressner|editor-first3=Daniel|editor-last4=Nobile|editor-first4=Fabio|editor-last5=Picasso|editor-first5=Marco|citeseerx = 10.1.1.764.6242}}</ref> पैरारियल को स्थिर करने के लिए विभिन्न दृष्टिकोण उपस्थित हैं,<ref>{{Cite journal|last1=Dai|first1=X.|last2=Maday|first2=Y.|date=2013-01-01|title=पहले और दूसरे क्रम की हाइपरबोलिक प्रणालियों के लिए समय विधि में स्थिर पैरारियल|journal=SIAM Journal on Scientific Computing|volume=35|issue=1|pages=A52–A78|doi=10.1137/110861002|arxiv=1201.1064 |bibcode=2013SJSC...35A..52D |s2cid=32336370|issn=1064-8275|url=https://semanticscholar.org/paper/beaf448f4acc6ee6a6f31a2fc8de69805ba57ff4}}</ref><ref name=":0">{{Cite journal|last1=Farhat|first1=Charbel|last2=Cortial|first2=Julien|last3=Dastillung|first3=Climène|last4=Bavestrello|first4=Henri|date=2006-07-30|title=रैखिक संरचनात्मक गतिशील प्रतिक्रियाओं की निकट-वास्तविक समय भविष्यवाणी के लिए समय-समानांतर अंतर्निहित इंटीग्रेटर्स|journal=International Journal for Numerical Methods in Engineering|language=en|volume=67|issue=5|pages=697–724|doi=10.1002/nme.1653|issn=1097-0207|bibcode=2006IJNME..67..697F|s2cid=121254646 }}</ref><ref name=":1">{{Cite book|title=पैरारियल विधि को तेज करने और स्थिर करने के लिए कम आधार विधियों के उपयोग पर|last1=Chen|first1=Feng|last2=Hesthaven|first2=Jan S.|last3=Zhu|first3=Xueyu|date=2014-01-01|publisher=Springer International Publishing|isbn=9783319020891|editor-last=Quarteroni|editor-first=Alfio|series=MS&A - Modeling, Simulation and Applications|pages=187–214|language=en|doi=10.1007/978-3-319-02090-7_7|editor-last2=Rozza|editor-first2=Gianluigi|url = http://infoscience.epfl.ch/record/190666}}</ref> क्रायलोव-सबस्पेस संवर्धित पैरारियल में से एक है। | ||
== वेरिएंट == | == वेरिएंट == | ||
ऐसे कई एल्गोरिदम हैं जो | ऐसे कई एल्गोरिदम हैं जो सामान्यतः आधारित हैं या कम से कम मूल पैरारियल एल्गोरिदम से प्रेरित हैं। | ||
=== क्रायलोव-सबस्पेस ने पैरारियल को बढ़ाया === | === क्रायलोव-सबस्पेस ने पैरारियल को बढ़ाया === | ||
प्रारंभ में ही यह माना गया कि रैखिक समस्याओं के लिए सूचना | प्रारंभ में ही यह माना गया कि रैखिक समस्याओं के लिए सूचना फाइन विधि <math>\mathcal{F}_{\delta t}</math> द्वारा उत्पन्न जानकारी का उपयोग कोर्स विधि <math>\mathcal{G}_{\Delta t}</math> की शुद्धता में सुधार के लिए किया जा सकता है।<ref name=":0" /> मूल रूप से, यह विचार समानांतर अंतर्निहित समय-एकीकरणकर्ता पीआईटीए के लिए तैयार किया गया था,<ref>{{Cite journal|last1=Farhat|first1=Charbel|last2=Chandesris|first2=Marion|date=2003-11-07|title=Time-decomposed parallel time-integrators: theory and feasibility studies for fluid, structure, and fluid–structure applications|journal=International Journal for Numerical Methods in Engineering|language=en|volume=58|issue=9|pages=1397–1434|doi=10.1002/nme.860|issn=1097-0207|bibcode=2003IJNME..58.1397F|s2cid=61667246 }}</ref> यह विधि है जो पैरारियल से निकटता से संबंधित है किन्तु सुधार कैसे किया जाता है इसमें थोड़ा अंतर है। प्रत्येक पुनरावृत्ति <math>k</math> में परिणाम <math>\mathcal{F}_{\delta t}(U^k_j)</math> की गणना <math>j=0, \ldots, N-1</math> के लिए मान <math>u^k_j \in \mathbb{R}^d</math> के लिए की जाती है। इस जानकारी के आधार पर, [[ सदिश स्थल |सदिश स्थल]] | ||
<math> | <math> | ||
S_k := \left\{ U^{k'}_j : 0 \leq k' \leq k, j=0, \ldots, N-1 \right\} </math> | S_k := \left\{ U^{k'}_j : 0 \leq k' \leq k, j=0, \ldots, N-1 \right\} </math> | ||
जैसे-जैसे पुनरावृत्तियों की संख्या बढ़ती है, | प्रत्येक पैरारियल पुनरावृत्ति के बाद परिभाषित और अद्यतन किया जाता है।<ref>{{Cite journal|last1=Gander|first1=M.|last2=Petcu|first2=M.|title=रैखिक समस्याओं के लिए क्रायलोव उप-स्थान संवर्धित पैरारियल एल्गोरिदम का विश्लेषण|journal=ESAIM: Proceedings|volume=25|pages=114–129|doi=10.1051/proc:082508|year=2008|doi-access=free}}</ref> <math> \mathbb{R}^d</math> से <math> S_k</math> तक ओर्थोगोनल प्रक्षेपण को <math> P_k</math> के रूप में निरूपित करें। फिर, मोटे तरीके को बेहतर इंटीग्रेटर केडी से बदलें। | ||
=== हाइब्रिड पैरारियल | |||
<math> \mathcal{K}_{\Delta t}(u) = \mathcal{F}_{\delta t}(P_k u) + \mathcal{G}_{\Delta t}((I-P_k)u)</math>. | |||
जैसे-जैसे पुनरावृत्तियों की संख्या बढ़ती है, स्पेस <math> S_k</math> बढ़ेगा और संशोधित प्रचारक <math> \mathcal{K}_{\Delta t}</math> अधिक त्रुटिहीन हो जाएगा। इससे तेजी से अभिसरण हो सकता है। पैरारियल का यह संस्करण रैखिक अतिशयोक्तिपूर्ण आंशिक अंतर समीकरणों को भी स्थिर रूप से एकीकृत कर सकता है।<ref>{{Cite journal|last1=Ruprecht|first1=D.|last2=Krause|first2=R.|date=2012-04-30|title=एक रैखिक ध्वनिक-संवहन प्रणाली का स्पष्ट समानांतर-समय एकीकरण|journal=Computers & Fluids|volume=59|pages=72–83|doi=10.1016/j.compfluid.2012.02.015|arxiv=1510.02237|s2cid=15703896 }}</ref> कम आधार पद्धति पर आधारित अरेखीय समस्याओं का विस्तार भी उपस्थित है।<ref name=":1" /> | |||
=== हाइब्रिड पैरारियल स्पेक्ट्रल डैफर्ड करेक्शन === | |||
स्पेक्ट्रल डैफर्ड करेक्शन (एसडीसी) के साथ पैरारियल के संयोजन पर आधारित उत्तम समानांतर दक्षता वाली विधि <ref>{{Cite journal|last1=Dutt|first1=Alok|last2=Greengard|first2=Leslie|last3=Rokhlin|first3=Vladimir|date=2000-06-01|title=साधारण विभेदक समीकरणों के लिए वर्णक्रमीय आस्थगित सुधार विधियाँ|journal=BIT Numerical Mathematics|language=en|volume=40|issue=2|pages=241–266|doi=10.1023/A:1022338906936|s2cid=43449672 |issn=0006-3835}}</ref> एम. मिनियन द्वारा प्रस्तावित किया गया है।<ref>{{Cite journal|last=Minion|first=Michael|date=2011-01-05|title=एक संकर पैरारियल वर्णक्रमीय आस्थगित सुधार विधि|journal=Communications in Applied Mathematics and Computational Science|volume=5|issue=2|pages=265–301|doi=10.2140/camcos.2010.5.265|issn=2157-5452|doi-access=free}}</ref> यह उत्तम समानांतर दक्षता के लिए लचीलेपन का त्याग करते हुए कोर्स और फाइन इंटीग्रेटर के विकल्प को एसडीसी तक सीमित कर देता है। <math>1/k</math> की सीमा के अतिरिक्त, हाइब्रिड विधि में समानांतर दक्षता पर सीमा | |||
<math>E_p \leq \frac{k_s}{k_p}</math> | <math>E_p \leq \frac{k_s}{k_p}</math> | ||
=== | बन जाती है, जिसमें <math>k_s</math> सीरियल एसडीसी बेस विधि के पुनरावृत्तियों की संख्या होती है और <math>k_p</math> सामान्यतः समानांतर हाइब्रिड विधि के पुनरावृत्तियों की अधिक संख्या होती है। नॉनलाइनियर मल्टीग्रिड विधि में उपयोग की जाने वाली पूर्ण सन्निकटन योजना को जोड़कर पैरारियल-एसडीसी हाइब्रिड को और बेहतर बनाया गया है। इससे अंतरिक्ष और समय में समानांतर पूर्ण सन्निकटन योजना (पीएफएएसएसटी) का विकास हुआ।<ref>{{Cite journal|last1=Emmett|first1=Matthew|last2=Minion|first2=Michael|date=2012-03-28|title=आंशिक अंतर समीकरणों के लिए समय में एक कुशल समानांतर विधि की ओर|journal=Communications in Applied Mathematics and Computational Science|volume=7|issue=1|pages=105–132|doi=10.2140/camcos.2012.7.105|issn=2157-5452|doi-access=free}}</ref> पीएफएएसएसटी के प्रदर्शन का अध्ययन पीईपीसी के लिए किया गया है, जो जूलिच सुपरकंप्यूटिंग सेंटर में विकसित बार्न्स-हट ट्री कोड आधारित कण सॉल्वर है। आईबीएम ब्लू जीन/पी सिस्टम जुगीन पर सभी 262,144 कोर का उपयोग करके सिमुलेशन से पता चला कि पीएफएएसएसटी स्थानिक वृक्ष समानांतरीकरण की संतृप्ति से परे अतिरिक्त गति उत्पन्न कर सकता है।<ref>{{Cite book|last1=Speck|first1=R.|last2=Ruprecht|first2=D.|last3=Krause|first3=R.|last4=Emmett|first4=M.|last5=Minion|first5=M.|last6=Winkel|first6=M.|last7=Gibbon|first7=P.|date=2012-11-01|title=एक विशाल अंतरिक्ष-समय समानांतर एन-बॉडी सॉल्वर|journal=High Performance Computing, Networking, Storage and Analysis (SC), 2012 International Conference for|pages=1–11|doi=10.1109/SC.2012.6|isbn=978-1-4673-0805-2|s2cid=1620219 }}</ref> | ||
=== मल्टीग्रिड रिडक्शन इन टाइम मेथड (एमजीआरआईटी) === | |||
मल्टीग्रिड रिडक्शन इन टाइम मेथड (एमजीआरआईटी) विभिन्न स्मूथर्स का उपयोग करके कई स्तरों पर मल्टीग्रिड-इन-टाइम एल्गोरिदम के रूप में पैरारियल की व्याख्या को सामान्यीकृत करता है।<ref>{{Cite journal|last1=Falgout|first1=R.|last2=Friedhoff|first2=S.|last3=Kolev|first3=T.|last4=MacLachlan|first4=S.|last5=Schroder|first5=J.|date=2014-01-01|title=मल्टीग्रिड के साथ समानांतर समय एकीकरण|journal=SIAM Journal on Scientific Computing|volume=36|issue=6|pages=C635–C661|doi=10.1137/130944230|bibcode=2014SJSC...36C.635F |issn=1064-8275|citeseerx=10.1.1.701.2603}}</ref> यह अधिक सामान्य दृष्टिकोण है किन्तु मापदंडों की विशिष्ट पसंद के लिए यह पैरारियल के बराबर है। एमजीआरआईटी को लागू करने वाली [http://computation.llnl.gov/projects/parallel-time-integration-multigrid एक्सब्रैड] लाइब्रेरी [[ लॉरेंस लिवरमोर राष्ट्रीय प्रयोगशाला |लॉरेंस लिवरमोर राष्ट्रीय प्रयोगशाला]] द्वारा विकसित की जा रही है। | |||
=== पैराईएक्सपी === | |||
पैराईएक्सपी, पैरारियल के अन्दर घातीय इंटीग्रेटर्स का उपयोग करता है।<ref>{{Cite journal|last1=Gander|first1=M.|last2=Güttel|first2=S.|date=2013-01-01|title=PARAEXP: A Parallel Integrator for Linear Initial-Value Problems|journal=SIAM Journal on Scientific Computing|volume=35|issue=2|pages=C123–C142|doi=10.1137/110856137|bibcode=2013SJSC...35C.123G |issn=1064-8275|citeseerx=10.1.1.800.5938}}</ref> रैखिक समस्याओं तक सीमित रहते हुए, यह लगभग इष्टतम समानांतर गति उत्पन्न कर सकता है। | |||
==संदर्भ== | ==संदर्भ== | ||
Line 227: | Line 240: | ||
== बाहरी संबंध == | == बाहरी संबंध == | ||
*[https://www.parallel-in-time.org/ parallel-in-time.org] | *[https://www.parallel-in-time.org/ parallel-in-time.org] | ||
[[Category: | [[Category:All articles with unsourced statements]] | ||
[[Category:Articles with unsourced statements from November 2018]] | |||
[[Category:CS1 English-language sources (en)]] | |||
[[Category:Created On 23/07/2023]] | [[Category:Created On 23/07/2023]] | ||
[[Category:Lua-based templates]] | |||
[[Category:Machine Translated Page]] | |||
[[Category:Pages with script errors]] | |||
[[Category:Templates Vigyan Ready]] | |||
[[Category:Templates that add a tracking category]] | |||
[[Category:Templates that generate short descriptions]] | |||
[[Category:Templates using TemplateData]] | |||
[[Category:कम्प्यूटेशनल विज्ञान]] | |||
[[Category:संख्यात्मक अंतर समीकरण]] | |||
[[Category:संख्यात्मक विश्लेषण]] | |||
[[Category:समानांतर कंप्यूटिंग]] |
Latest revision as of 12:52, 28 July 2023
पैरारियल संख्यात्मक विश्लेषण से एक समानांतर एल्गोरिदम है और प्रारंभिक मूल्य समस्याओं के समाधान के लिए उपयोग किया जाता है।[1]
इसे 2001 में जैक्स-लुई लायंस, मैडे और टुरिनिसी द्वारा प्रस्तुत किया गया था। तब से यह समय एकीकरण में सबसे व्यापक रूप से अध्ययन किए जाने वाले समानांतर प्रणालियों में से एक बन गया है।[citation needed]
समय में समानांतर एकीकरण विधियां
उदाहरण के विपरीत रंज-कुट्टा या मल्टी-स्टेप विधियां, पैरारियल में कुछ गणनाएं समानांतर कंप्यूटिंग में की जा सकती हैं और इसलिए पैरारियल समय एकीकरण विधि में समानांतर का एक उदाहरण है। जबकि ऐतिहासिक रूप से आंशिक अंतर समीकरणों के संख्यात्मक आंशिक अंतर समीकरणों को समानांतर करने के अधिकांश प्रयास स्थानिक विवेकीकरण पर केंद्रित हैं, एक्सास्केल कंप्यूटिंग की चुनौतियों को देखते हुए, अस्थायी विवेकीकरण के समानांतर प्रणालियों को संख्यात्मक विश्लेषण सॉफ्टवेयर की सूची में समरूपता बढ़ाने के संभावित विधि के रूप में पहचाना गया है।[3]
क्योंकि पैरारियल समानांतर में कई समय चरणों के लिए संख्यात्मक समाधान की गणना करता है, इसे चरणों में समानांतर विधि के रूप में वर्गीकृत किया गया है।[4]
यह समानांतर रंज-कुट्टा या एक्सट्रपलेशन विधियों जैसी विधि में समानता का उपयोग करने वाले दृष्टिकोणों के विपरीत है, जहां स्वतंत्र चरणों की गणना तरंग रूप विश्राम जैसी सिस्टम विधियों में समानांतर या समानांतर में की जा सकती है।[5][6]
इतिहास
पैरारियल को समय विधि में मल्टीग्रिड विधि या समय अक्ष के साथ प्रत्यक्ष एकाधिक शूटिंग विधि दोनों के रूप में प्राप्त किया जा सकता है।[7]
समय में मल्टीग्रिड के साथ-साथ समय एकीकरण के लिए मल्टीपल शूटिंग को अपनाने के दोनों विचार 1980 और 1990 के दशक से चले आ रहे हैं।[8][9]
पैरारियल व्यापक रूप से अध्ययन की जाने वाली विधि है और विभिन्न अनुप्रयोगों के लिए इसका उपयोग और संशोधन किया गया है।[10]
प्रारंभिक मूल्य समस्याओं के समाधान को समानांतर करने के विचार और भी प्राचीन हैं: समय में समानांतर एकीकरण विधि का प्रस्ताव करने वाला पहला पेपर 1964 में सामने आया था।[11]
एल्गोरिथम
समस्या
लक्ष्य प्रपत्र की प्रारंभिक मूल्य समस्या का समाधान करना है
दाहिने हाथ की ओर को एक सुचारू (संभवतः अरेखीय) फलन माना जाता है। यह रेखा दृष्टिकोण की विधि में आंशिक अंतर समीकरण के स्थानिक विवेक के अनुरूप भी हो सकता है। हम इस समस्या को समान दूरी वाले बिंदु , जहां और , के अस्थायी आधार पर समाधान करना चाहते हैं। इस विवेक को आगे बढ़ाते हुए हमें एक विभाजित समय अंतराल प्राप्त होता है जिसमें के लिए समय स्लाइस सम्मिलित होता है।
इसका उद्देश्य सीरियल टाइम-स्टेपिंग विधि (जैसे रनगे-कुट्टा) का उपयोग करके त्रुटिहीन समाधान के लिए संख्यात्मक अनुमान की गणना करना है जिसमें उच्च संख्यात्मक शुद्धता (और इसलिए उच्च कम्प्यूटेशनल लागत) है। हम इस विधि को फाइन सॉल्वर के रूप में संदर्भित करते हैं, जो समय पर प्रारंभिक मान को समय पर टर्मिनल मान तक प्रसारित करता है। लक्ष्य का उपयोग करके समाधान (उच्च संख्यात्मक शुद्धता के साथ) की गणना करना है, जो हमें प्राप्त होता है
इस (और सबसे पहले समानांतर में समाधान करने का प्रयास करने का कारण) समाधान के साथ समस्या यह है कि वास्तविक समय में गणना करना कम्प्यूटेशनल रूप से संभव नहीं है।
यह कैसे काम करता है
प्रारंभिक मूल्य समस्या का समाधान करने के लिए एकल प्रोसेसर का उपयोग करने के अतिरिक्त (जैसा कि पारंपरिक समय-चरण विधियों के साथ किया जाता है), पैरारियल प्रोसेसर का उपयोग करता है। इसका उद्देश्य समानांतर में छोटी प्रारंभिक मूल्य समस्याओं (प्रत्येक समय स्लाइस पर एक) का समाधान करने के लिए प्रोसेसर का उपयोग करना है। उदाहरण के लिए, संदेश पासिंग इंटरफ़ेस आधारित कोड में, प्रक्रियाओं की संख्या होगी, जबकि ओपनएमपी आधारित कोड में, थ्रेड (कंप्यूटिंग) की संख्या के बराबर होगी।
पैरारियल इस प्रारंभिक मूल्य समस्या को समानांतर में समाधान करने के लिए दूसरी टाइम-स्टेपिंग विधि का उपयोग करता है, जिसे कोर्स सॉल्वर के रूप में जाना जाता है।
कोर्स सॉल्वर ठीक सॉल्वर की तरह ही काम करता है, लंबाई के समय अंतराल पर प्रारंभिक मूल्य का प्रचार करता है, चूंकि यह (और इसलिए बहुत कम कम्प्यूटेशनल लागत पर) की तुलना में बहुत कम संख्यात्मक शुद्धता पर ऐसा करता है। एक कोर्स सॉल्वर का होना जो कि फाइन सॉल्वर की तुलना में बहुत कम कम्प्यूटेशनल रूप से एक्सपेंसिव है, पैरारियल के साथ समानांतर गति प्राप्त करने की कुंजी है।
अब से, हम समय और पुनरावृत्ति पर पैरारियल समाधान को द्वारा निरूपित करते हैं।
शून्य पुनरावृत्ति
सबसे पहले, समाधान के अनुमानित प्रारंभिक अनुमान की गणना करने के लिए कोर्स सॉल्वर को पूरे समय अंतराल पर क्रमिक रूप से चलाएं:
अनुवर्ती पुनरावृत्तियाँ
इसके बाद, सबसे अद्यतित समाधान मानों से, समानांतर में, प्रत्येक समय स्लाइस पर अच्छे सॉल्वर चलाएं:
अब प्रेडिक्टर-करेक्टर का उपयोग करके क्रमिक रूप से पैरारियल समाधान मानों को अपडेट करें:
इस स्तर पर, कोई यह निर्धारित करने के लिए स्टॉपिंग मानदंड का उपयोग कर सकता है कि क्या समाधान मान अब प्रत्येक पुनरावृत्ति में नहीं बदल रहे हैं। उदाहरण के लिए, कोई इसकी जांच करके यह जांच सकता है कि क्या
और कुछ सहनशीलता है। यदि यह मानदंड संतुष्ट नहीं है, तो बाद के पुनरावृत्तियों को समानांतर में फाइन सॉल्वर और फिर भविष्यवक्ता-सुधारक को लागू करके चलाया जा सकता है। चूँकि, एक बार जब मानदंड संतुष्ट हो जाता है, तो कहा जाता है कि एल्गोरिदम पुनरावृत्तियों में परिवर्तित हो गया है। ध्यान दें कि अन्य रोक मानदंड उपस्थित हैं और पैरारियल में सफलतापूर्वक परीक्षण किया गया है।
टिप्पणियाँ
पैरारियल को उस समाधान को पुन: प्रस्तुत करना चाहिए जो फाइन सॉल्वर के क्रमिक अनुप्रयोग द्वारा प्राप्त किया गया है और अधिकतम पुनरावृत्तियों में परिवर्तित होगा।[7] चूँकि, पैरारियल को स्पीडअप प्रदान करने के लिए, इसे समय स्लाइस की संख्या अर्थात। की तुलना में काफी कम संख्या में पुनरावृत्तियों में परिवर्तित करना होगा।
पैरारियल पुनरावृत्ति में, का कम्प्यूटेशनल रूप से एक्सपेंसिव मूल्यांकन प्रसंस्करण इकाइयों पर समानांतर में किया जा सकता है। इसके विपरीत, पर की निर्भरता का अर्थ है कि कोर्स सुधार की गणना क्रमिक क्रम में की जानी है।
सामान्यतः, रंज-कुट्टा विधि का कुछ रूप कोर्स और बारीक इंटीग्रेटर दोनों के लिए चुना जाता है, जहां कम क्रम का हो सकता है और की तुलना में बड़े समय चरण का उपयोग कर सकता है।
यदि प्रारंभिक मूल्य समस्या पीडीई के विवेकाधिकार से उत्पन्न होती है, तो एक कोर्स स्थानिक विवेकीकरण का भी उपयोग कर सकता है, किन्तु यह अभिसरण पर नकारात्मक प्रभाव डाल सकता है जब तक कि उच्च क्रम इंटरपोलेशन का उपयोग नहीं किया जाता है।[12]
स्पीडअप
कुछ मान्यताओं के अनुसार, पैरारियल की गति के लिए एक सरल सैद्धांतिक मॉडल प्राप्त किया जा सकता है।[13]
चूँकि अनुप्रयोगों में ये धारणाएँ बहुत अधिक प्रतिबंधात्मक हो सकती हैं, फिर भी मॉडल पैरारियल के साथ स्पीडअप प्राप्त करने में सम्मिलित ट्रेडऑफ़ को दर्शाने के लिए उपयोगी है।
सबसे पहले, मान लें कि हर बार स्लाइस में फाइन इंटीग्रेटर के बिल्कुल चरण और कोर्स इंटीग्रेटर के चरण होते हैं। इसमें विशेष रूप से यह धारणा सम्मिलित है कि सभी समय स्लाइस समान लंबाई के होते हैं और कोर्स और फाइन इंटीग्रेटर दोनों पूर्ण सिमुलेशन पर एक स्थिर चरण आकार का उपयोग करते हैं। दूसरा, क्रमशः सूक्ष्म और कोर्स प्रणालियों के एक चरण के लिए आवश्यक कंप्यूटिंग समय को और द्वारा निरूपित करें, और मान लें कि दोनों स्थिर हैं। यह सामान्यतः बिल्कुल सच नहीं है जब एक अंतर्निहित विधि का उपयोग किया जाता है, क्योंकि तब रनटाइम पुनरावृत्त सॉल्वर द्वारा आवश्यक पुनरावृत्तियों की संख्या के आधार पर भिन्न होता है।
इन दो मान्यताओं के अनुसार, टाइम स्लाइस पर एकीकृत करने वाली फाइन विधि के रनटाइम को इस प्रकार मॉडल किया जा सकता है
प्रसंस्करण इकाइयों का उपयोग करके और पुनरावृत्तियों को निष्पादित करके पैरारियल का रनटाइम है
पैरारियल का स्पीडअप तो है
ये दो सीमाएँ कोर्स विधियों को चुनने में किए जाने वाले व्यापार को दर्शाती हैं: एक ओर, इसे सस्ता होना होगा और/या पहली सीमा को जितना संभव हो उतना बड़ा बनाने के लिए बहुत बड़े समय के कदम का उपयोग करना होगा, दूसरी ओर दूसरी सीमा को बड़ा रखने के लिए पुनरावृत्तियों की संख्या को कम रखना होगा।
विशेष रूप से, पैरारियल की समानांतर दक्षता सीमित है
यह आवश्यक पुनरावृत्तियों की संख्या के व्युत्क्रम से है।
काल्पनिक आइजेनवैल्यूज़ के लिए अस्थिरता
पैरारियल के वेनिला संस्करण में काल्पनिक आइगेनवैल्यूज़ एवं आइगेनवेक्टर्स के साथ समस्याएं हैं।[7] यह सामान्यतः केवल अंतिम पुनरावृत्तियों की ओर ही परिवर्तित होता है, अर्थात जैसे-जैसे क्वे निकट पहुंचता है, और स्पीडअप सदैव से छोटा होता है। तो या तो पुनरावृत्तियों की संख्या छोटी है और पैरारियल अस्थिर है या, यदि पैरारियल को स्थिर बनाने के लिए पर्याप्त बड़ा है, कोई गति संभव नहीं है। इसका यह भी अर्थ है कि हाइपरबोलिक आंशिक अंतर समीकरण समीकरणों के लिए पैरारियल सामान्यतः अस्थिर है।[14] भले ही गैंडर और वांडेवेले द्वारा औपचारिक विश्लेषण केवल निरंतर गुणांक के साथ रैखिक समस्याओं को कवर करता है, समस्या तब भी उत्पन्न होती है जब पैरारियल को गैर-रेखीय नेवियर-स्टोक्स समीकरणों पर लागू किया जाता है जब चिपचिपाहट गुणांक बहुत छोटा हो जाता है और रेनॉल्ड्स संख्या बहुत बड़ी हो जाती है।[15] पैरारियल को स्थिर करने के लिए विभिन्न दृष्टिकोण उपस्थित हैं,[16][17][18] क्रायलोव-सबस्पेस संवर्धित पैरारियल में से एक है।
वेरिएंट
ऐसे कई एल्गोरिदम हैं जो सामान्यतः आधारित हैं या कम से कम मूल पैरारियल एल्गोरिदम से प्रेरित हैं।
क्रायलोव-सबस्पेस ने पैरारियल को बढ़ाया
प्रारंभ में ही यह माना गया कि रैखिक समस्याओं के लिए सूचना फाइन विधि द्वारा उत्पन्न जानकारी का उपयोग कोर्स विधि की शुद्धता में सुधार के लिए किया जा सकता है।[17] मूल रूप से, यह विचार समानांतर अंतर्निहित समय-एकीकरणकर्ता पीआईटीए के लिए तैयार किया गया था,[19] यह विधि है जो पैरारियल से निकटता से संबंधित है किन्तु सुधार कैसे किया जाता है इसमें थोड़ा अंतर है। प्रत्येक पुनरावृत्ति में परिणाम की गणना के लिए मान के लिए की जाती है। इस जानकारी के आधार पर, सदिश स्थल
प्रत्येक पैरारियल पुनरावृत्ति के बाद परिभाषित और अद्यतन किया जाता है।[20] से तक ओर्थोगोनल प्रक्षेपण को के रूप में निरूपित करें। फिर, मोटे तरीके को बेहतर इंटीग्रेटर केडी से बदलें।
.
जैसे-जैसे पुनरावृत्तियों की संख्या बढ़ती है, स्पेस बढ़ेगा और संशोधित प्रचारक अधिक त्रुटिहीन हो जाएगा। इससे तेजी से अभिसरण हो सकता है। पैरारियल का यह संस्करण रैखिक अतिशयोक्तिपूर्ण आंशिक अंतर समीकरणों को भी स्थिर रूप से एकीकृत कर सकता है।[21] कम आधार पद्धति पर आधारित अरेखीय समस्याओं का विस्तार भी उपस्थित है।[18]
हाइब्रिड पैरारियल स्पेक्ट्रल डैफर्ड करेक्शन
स्पेक्ट्रल डैफर्ड करेक्शन (एसडीसी) के साथ पैरारियल के संयोजन पर आधारित उत्तम समानांतर दक्षता वाली विधि [22] एम. मिनियन द्वारा प्रस्तावित किया गया है।[23] यह उत्तम समानांतर दक्षता के लिए लचीलेपन का त्याग करते हुए कोर्स और फाइन इंटीग्रेटर के विकल्प को एसडीसी तक सीमित कर देता है। की सीमा के अतिरिक्त, हाइब्रिड विधि में समानांतर दक्षता पर सीमा
बन जाती है, जिसमें सीरियल एसडीसी बेस विधि के पुनरावृत्तियों की संख्या होती है और सामान्यतः समानांतर हाइब्रिड विधि के पुनरावृत्तियों की अधिक संख्या होती है। नॉनलाइनियर मल्टीग्रिड विधि में उपयोग की जाने वाली पूर्ण सन्निकटन योजना को जोड़कर पैरारियल-एसडीसी हाइब्रिड को और बेहतर बनाया गया है। इससे अंतरिक्ष और समय में समानांतर पूर्ण सन्निकटन योजना (पीएफएएसएसटी) का विकास हुआ।[24] पीएफएएसएसटी के प्रदर्शन का अध्ययन पीईपीसी के लिए किया गया है, जो जूलिच सुपरकंप्यूटिंग सेंटर में विकसित बार्न्स-हट ट्री कोड आधारित कण सॉल्वर है। आईबीएम ब्लू जीन/पी सिस्टम जुगीन पर सभी 262,144 कोर का उपयोग करके सिमुलेशन से पता चला कि पीएफएएसएसटी स्थानिक वृक्ष समानांतरीकरण की संतृप्ति से परे अतिरिक्त गति उत्पन्न कर सकता है।[25]
मल्टीग्रिड रिडक्शन इन टाइम मेथड (एमजीआरआईटी)
मल्टीग्रिड रिडक्शन इन टाइम मेथड (एमजीआरआईटी) विभिन्न स्मूथर्स का उपयोग करके कई स्तरों पर मल्टीग्रिड-इन-टाइम एल्गोरिदम के रूप में पैरारियल की व्याख्या को सामान्यीकृत करता है।[26] यह अधिक सामान्य दृष्टिकोण है किन्तु मापदंडों की विशिष्ट पसंद के लिए यह पैरारियल के बराबर है। एमजीआरआईटी को लागू करने वाली एक्सब्रैड लाइब्रेरी लॉरेंस लिवरमोर राष्ट्रीय प्रयोगशाला द्वारा विकसित की जा रही है।
पैराईएक्सपी
पैराईएक्सपी, पैरारियल के अन्दर घातीय इंटीग्रेटर्स का उपयोग करता है।[27] रैखिक समस्याओं तक सीमित रहते हुए, यह लगभग इष्टतम समानांतर गति उत्पन्न कर सकता है।
संदर्भ
- ↑ Lions, Jacques-Louis; Maday, Yvon; Turinici, Gabriel (2015). "A "parareal" in time discretization of PDE's" (PDF). Comptes Rendus de l'Académie des Sciences, Série I. 332 (7): 661–668. Bibcode:2001CRASM.332..661L. doi:10.1016/S0764-4442(00)01793-6.
- ↑ Pentland, Kamran; Tamborrino, Massimiliano; Samaddar, Debasmita; Appel, Lynton (2022). "Stochastic parareal: an application of probabilistic methods to time-parallelization". SIAM Journal on Scientific Computing (in English): S82–S102. doi:10.1137/21M1414231. S2CID 235485544.
- ↑ Jack Dongarra; Jeffrey Hittinger; John Bell; Luis Chacon; Robert Falgout; Michael Heroux; Paul Hovland; Esmond Ng; Clayton Webster; Stefan Wild (March 2014). Applied Mathematics Research for Exascale Computing (PDF) (Report). US Department of Energy. Retrieved August 29, 2015.
- ↑ Burrage, Kevin (1997). "Parallel methods for ODEs". Advances in Computational Mathematics. 7 (1–2): 1–31. doi:10.1023/A:1018997130884. S2CID 15778447.
- ↑ Iserles, A.; NøRSETT, S. P. (1990-10-01). "On the Theory of Parallel Runge—Kutta Methods". IMA Journal of Numerical Analysis (in English). 10 (4): 463–488. doi:10.1093/imanum/10.4.463. ISSN 0272-4979.
- ↑ Ketcheson, David; Waheed, Umair bin (2014-06-13). "A comparison of high-order explicit Runge–Kutta, extrapolation, and deferred correction methods in serial and parallel". Communications in Applied Mathematics and Computational Science. 9 (2): 175–200. arXiv:1305.6165. doi:10.2140/camcos.2014.9.175. ISSN 2157-5452. S2CID 15242644.
- ↑ 7.0 7.1 7.2 Gander, Martin J.; Vandewalle, Stefan (2007). "Analysis of the Parareal Time‐Parallel Time‐Integration Method". SIAM Journal on Scientific Computing. 29 (2): 556–578. CiteSeerX 10.1.1.154.6042. doi:10.1137/05064607X.
- ↑ Hackbusch, Wolfgang (1985). Parabolic multi-grid methods. pp. 189–197. ISBN 9780444875976. Retrieved August 29, 2015.
{{cite book}}
:|journal=
ignored (help) - ↑ Kiehl, Martin (1994). "Parallel multiple shooting for the solution of initial value problems". Parallel Computing. 20 (3): 275–295. doi:10.1016/S0167-8191(06)80013-X.
- ↑ Gander, Martin J. (2015). 50 years of Time Parallel Time Integration. Contributions in Mathematical and Computational Sciences. Vol. 9 (1 ed.). Springer International Publishing. doi:10.1007/978-3-319-23321-5. ISBN 978-3-319-23321-5.
- ↑ Nievergelt, Jürg (1964). "Parallel methods for integrating ordinary differential equations". Communications of the ACM. 7 (12): 731–733. doi:10.1145/355588.365137. S2CID 6361754.
- ↑ Ruprecht, Daniel (2014-12-01). "स्थानिक मोटेपन के साथ पैरारियल का अभिसरण" (PDF). Proceedings in Applied Mathematics and Mechanics (in English). 14 (1): 1031–1034. doi:10.1002/pamm.201410490. ISSN 1617-7061. S2CID 26356904.
- ↑ Minion, Michael L. (2010). "A Hybrid Parareal Spectral Deferred Corrections Method". Communications in Applied Mathematics and Computational Science. 5 (2): 265–301. doi:10.2140/camcos.2010.5.265.
- ↑ Staff, Gunnar Andreas; Rønquist, Einar M. (2005-01-01). Barth, Timothy J.; Griebel, Michael; Keyes, David E.; Nieminen, Risto M.; Roose, Dirk; Schlick, Tamar; Kornhuber, Ralf; Hoppe, Ronald; Périaux, Jacques (eds.). पैरारियल एल्गोरिथम की स्थिरता. Lecture Notes in Computational Science and Engineering (in English). Springer Berlin Heidelberg. pp. 449–456. doi:10.1007/3-540-26825-1_46. ISBN 9783540225232.
- ↑ Steiner, Johannes; Ruprecht, Daniel; Speck, Robert; Krause, Rolf (2015-01-01). Abdulle, Assyr; Deparis, Simone; Kressner, Daniel; Nobile, Fabio; Picasso, Marco (eds.). रेनॉल्ड्स संख्या के आधार पर नेवियर-स्टोक्स समीकरणों के लिए पैरारियल का अभिसरण. Lecture Notes in Computational Science and Engineering (in English). Springer International Publishing. pp. 195–202. CiteSeerX 10.1.1.764.6242. doi:10.1007/978-3-319-10705-9_19. ISBN 9783319107042.
- ↑ Dai, X.; Maday, Y. (2013-01-01). "पहले और दूसरे क्रम की हाइपरबोलिक प्रणालियों के लिए समय विधि में स्थिर पैरारियल". SIAM Journal on Scientific Computing. 35 (1): A52–A78. arXiv:1201.1064. Bibcode:2013SJSC...35A..52D. doi:10.1137/110861002. ISSN 1064-8275. S2CID 32336370.
- ↑ 17.0 17.1 Farhat, Charbel; Cortial, Julien; Dastillung, Climène; Bavestrello, Henri (2006-07-30). "रैखिक संरचनात्मक गतिशील प्रतिक्रियाओं की निकट-वास्तविक समय भविष्यवाणी के लिए समय-समानांतर अंतर्निहित इंटीग्रेटर्स". International Journal for Numerical Methods in Engineering (in English). 67 (5): 697–724. Bibcode:2006IJNME..67..697F. doi:10.1002/nme.1653. ISSN 1097-0207. S2CID 121254646.
- ↑ 18.0 18.1 Chen, Feng; Hesthaven, Jan S.; Zhu, Xueyu (2014-01-01). Quarteroni, Alfio; Rozza, Gianluigi (eds.). पैरारियल विधि को तेज करने और स्थिर करने के लिए कम आधार विधियों के उपयोग पर. MS&A - Modeling, Simulation and Applications (in English). Springer International Publishing. pp. 187–214. doi:10.1007/978-3-319-02090-7_7. ISBN 9783319020891.
- ↑ Farhat, Charbel; Chandesris, Marion (2003-11-07). "Time-decomposed parallel time-integrators: theory and feasibility studies for fluid, structure, and fluid–structure applications". International Journal for Numerical Methods in Engineering (in English). 58 (9): 1397–1434. Bibcode:2003IJNME..58.1397F. doi:10.1002/nme.860. ISSN 1097-0207. S2CID 61667246.
- ↑ Gander, M.; Petcu, M. (2008). "रैखिक समस्याओं के लिए क्रायलोव उप-स्थान संवर्धित पैरारियल एल्गोरिदम का विश्लेषण". ESAIM: Proceedings. 25: 114–129. doi:10.1051/proc:082508.
- ↑ Ruprecht, D.; Krause, R. (2012-04-30). "एक रैखिक ध्वनिक-संवहन प्रणाली का स्पष्ट समानांतर-समय एकीकरण". Computers & Fluids. 59: 72–83. arXiv:1510.02237. doi:10.1016/j.compfluid.2012.02.015. S2CID 15703896.
- ↑ Dutt, Alok; Greengard, Leslie; Rokhlin, Vladimir (2000-06-01). "साधारण विभेदक समीकरणों के लिए वर्णक्रमीय आस्थगित सुधार विधियाँ". BIT Numerical Mathematics (in English). 40 (2): 241–266. doi:10.1023/A:1022338906936. ISSN 0006-3835. S2CID 43449672.
- ↑ Minion, Michael (2011-01-05). "एक संकर पैरारियल वर्णक्रमीय आस्थगित सुधार विधि". Communications in Applied Mathematics and Computational Science. 5 (2): 265–301. doi:10.2140/camcos.2010.5.265. ISSN 2157-5452.
- ↑ Emmett, Matthew; Minion, Michael (2012-03-28). "आंशिक अंतर समीकरणों के लिए समय में एक कुशल समानांतर विधि की ओर". Communications in Applied Mathematics and Computational Science. 7 (1): 105–132. doi:10.2140/camcos.2012.7.105. ISSN 2157-5452.
- ↑ Speck, R.; Ruprecht, D.; Krause, R.; Emmett, M.; Minion, M.; Winkel, M.; Gibbon, P. (2012-11-01). एक विशाल अंतरिक्ष-समय समानांतर एन-बॉडी सॉल्वर. pp. 1–11. doi:10.1109/SC.2012.6. ISBN 978-1-4673-0805-2. S2CID 1620219.
{{cite book}}
:|journal=
ignored (help) - ↑ Falgout, R.; Friedhoff, S.; Kolev, T.; MacLachlan, S.; Schroder, J. (2014-01-01). "मल्टीग्रिड के साथ समानांतर समय एकीकरण". SIAM Journal on Scientific Computing. 36 (6): C635–C661. Bibcode:2014SJSC...36C.635F. CiteSeerX 10.1.1.701.2603. doi:10.1137/130944230. ISSN 1064-8275.
- ↑ Gander, M.; Güttel, S. (2013-01-01). "PARAEXP: A Parallel Integrator for Linear Initial-Value Problems". SIAM Journal on Scientific Computing. 35 (2): C123–C142. Bibcode:2013SJSC...35C.123G. CiteSeerX 10.1.1.800.5938. doi:10.1137/110856137. ISSN 1064-8275.