पैरारियल: Difference between revisions
(Created page with "{{Short description|Parallel algorithm from numerical analysis}} पैरारियल संख्यात्मक विश्लेषण से एक समा...") |
No edit summary |
||
| Line 19: | Line 19: | ||
इसे 2001 में [[जैक्स-लुई लायंस]], मैडे और टुरिनिसी द्वारा पेश किया गया था। तब से, यह सबसे व्यापक रूप से अध्ययन किए गए समानांतर-इन-टाइम एकीकरण तरीकों में से एक बन गया है।{{Citation needed|date=November 2018}} | इसे 2001 में [[जैक्स-लुई लायंस]], मैडे और टुरिनिसी द्वारा पेश किया गया था। तब से, यह सबसे व्यापक रूप से अध्ययन किए गए समानांतर-इन-टाइम एकीकरण तरीकों में से एक बन गया है।{{Citation needed|date=November 2018}} | ||
[[File:Parareal.svg|thumb| | [[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>).]] | ||
== समानांतर-इन-टाइम एकीकरण विधियां == | == समानांतर-इन-टाइम एकीकरण विधियां == | ||
| Line 45: | Line 45: | ||
}}</ref> | }}</ref> | ||
यह समानांतर रनगे-कुट्टा या एक्सट्रपलेशन विधियों जैसी विधि में समानता का उपयोग करने वाले दृष्टिकोणों के विपरीत है, जहां स्वतंत्र चरणों की गणना तरंग रूप विश्राम जैसी सिस्टम विधियों में समानांतर या समानांतर में की जा सकती है।<ref>{{Cite journal|last1=Iserles|first1=A.|last2=NøRSETT|first2=S. P.|date=1990-10-01|title=On the Theory of Parallel Runge—Kutta Methods|journal=IMA Journal of Numerical Analysis|language=en|volume=10|issue=4|pages=463–488|doi=10.1093/imanum/10.4.463|issn=0272-4979}}</ref><ref>{{Cite journal|last1=Ketcheson|first1=David|last2=Waheed|first2=Umair bin|date=2014-06-13|title=A comparison of high-order explicit Runge–Kutta, extrapolation, and deferred correction methods in serial and parallel|journal=Communications in Applied Mathematics and Computational Science|volume=9|issue=2|pages=175–200|doi=10.2140/camcos.2014.9.175|issn=2157-5452|arxiv=1305.6165|s2cid=15242644 }}</ref> | यह समानांतर रनगे-कुट्टा या एक्सट्रपलेशन विधियों जैसी विधि में समानता का उपयोग करने वाले दृष्टिकोणों के विपरीत है, जहां स्वतंत्र चरणों की गणना तरंग रूप विश्राम जैसी सिस्टम विधियों में समानांतर या समानांतर में की जा सकती है।<ref>{{Cite journal|last1=Iserles|first1=A.|last2=NøRSETT|first2=S. P.|date=1990-10-01|title=On the Theory of Parallel Runge—Kutta Methods|journal=IMA Journal of Numerical Analysis|language=en|volume=10|issue=4|pages=463–488|doi=10.1093/imanum/10.4.463|issn=0272-4979}}</ref><ref>{{Cite journal|last1=Ketcheson|first1=David|last2=Waheed|first2=Umair bin|date=2014-06-13|title=A comparison of high-order explicit Runge–Kutta, extrapolation, and deferred correction methods in serial and parallel|journal=Communications in Applied Mathematics and Computational Science|volume=9|issue=2|pages=175–200|doi=10.2140/camcos.2014.9.175|issn=2157-5452|arxiv=1305.6165|s2cid=15242644 }}</ref> | ||
== इतिहास == | == इतिहास == | ||
| Line 111: | Line 109: | ||
| s2cid = 6361754 | | s2cid = 6361754 | ||
}}</ref> | }}</ref> | ||
== एल्गोरिथम == | == एल्गोरिथम == | ||
| Line 162: | Line 158: | ||
यदि प्रारंभिक मूल्य समस्या पीडीई के विवेकाधिकार से उत्पन्न होती है, <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> | यदि प्रारंभिक मूल्य समस्या पीडीई के विवेकाधिकार से उत्पन्न होती है, <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>.]] | |||
[[File:Parareal Animation.ogg|thumb| | |||
===स्पीडअप=== | ===स्पीडअप=== | ||
| Line 214: | Line 209: | ||
जैसे-जैसे पुनरावृत्तियों की संख्या बढ़ती है, space <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" /> | जैसे-जैसे पुनरावृत्तियों की संख्या बढ़ती है, space <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>, संकर विधि में समानांतर दक्षता पर बाध्य हो जाता है | वर्णक्रमीय विलंबित सुधार (एसडीसी) के साथ पैरारियल के संयोजन पर आधारित बेहतर समानांतर दक्षता वाली एक विधि <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>, संकर विधि में समानांतर दक्षता पर बाध्य हो जाता है | ||
| Line 221: | Line 214: | ||
<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> पीएफएएसएसटी के प्रदर्शन का अध्ययन पीईपीसी के लिए किया गया है, जो जूलिच सुपरकंप्यूटिंग सेंटर में विकसित एक बार्न्स-हट सिमुलेशन|बार्न्स-हट ट्री कोड आधारित कण सॉल्वर है। IBM ब्लू जीन/पी सिस्टम JUGENE पर सभी 262,144 कोर का उपयोग करके सिमुलेशन से पता चला कि PFASST स्थानिक वृक्ष समानांतरीकरण की संतृप्ति से परे अतिरिक्त गति उत्पन्न कर सकता है।<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> | साथ <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> पीएफएएसएसटी के प्रदर्शन का अध्ययन पीईपीसी के लिए किया गया है, जो जूलिच सुपरकंप्यूटिंग सेंटर में विकसित एक बार्न्स-हट सिमुलेशन|बार्न्स-हट ट्री कोड आधारित कण सॉल्वर है। IBM ब्लू जीन/पी सिस्टम JUGENE पर सभी 262,144 कोर का उपयोग करके सिमुलेशन से पता चला कि PFASST स्थानिक वृक्ष समानांतरीकरण की संतृप्ति से परे अतिरिक्त गति उत्पन्न कर सकता है।<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 XBraid] लाइब्रेरी [[ लॉरेंस लिवरमोर राष्ट्रीय प्रयोगशाला ]] द्वारा विकसित की जा रही है। | मल्टीग्रिड रिडक्शन इन टाइम मेथड (एमजीआरआईटी) विभिन्न स्मूथर्स का उपयोग करके कई स्तरों पर मल्टीग्रिड-इन-टाइम एल्गोरिदम के रूप में पैरारियल की व्याख्या को सामान्यीकृत करता है।<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 XBraid] लाइब्रेरी [[ लॉरेंस लिवरमोर राष्ट्रीय प्रयोगशाला ]] द्वारा विकसित की जा रही है। | ||
| Line 231: | Line 222: | ||
==संदर्भ== | ==संदर्भ== | ||
{{Reflist}} | {{Reflist}} | ||
== बाहरी संबंध == | == बाहरी संबंध == | ||
*[https://www.parallel-in-time.org/ parallel-in-time.org] | *[https://www.parallel-in-time.org/ parallel-in-time.org] | ||
Revision as of 15:21, 25 July 2023
पैरारियल संख्यात्मक विश्लेषण से एक समानांतर एल्गोरिदम है और प्रारंभिक मूल्य समस्याओं के समाधान के लिए उपयोग किया जाता है।[1] इसे 2001 में जैक्स-लुई लायंस, मैडे और टुरिनिसी द्वारा पेश किया गया था। तब से, यह सबसे व्यापक रूप से अध्ययन किए गए समानांतर-इन-टाइम एकीकरण तरीकों में से एक बन गया है।[citation needed]
समानांतर-इन-टाइम एकीकरण विधियां
उदाहरण के विपरीत रनगे-कुट्टा विधियां|रंज-कुट्टा या रैखिक मल्टीस्टेप विधि|मल्टी-स्टेप विधियां, पैरारियल में कुछ गणनाएं समानांतर कंप्यूटिंग की जा सकती हैं और पैरारियल इसलिए समानांतर-इन-टाइम एकीकरण विधि का एक उदाहरण है। जबकि ऐतिहासिक रूप से आंशिक अंतर समीकरणों के संख्यात्मक आंशिक अंतर समीकरणों को समानांतर करने के अधिकांश प्रयास स्थानिक विवेकीकरण पर केंद्रित हैं, एक्सास्केल कंप्यूटिंग की चुनौतियों को देखते हुए, अस्थायी विवेकीकरण के समानांतर तरीकों को संख्यात्मक विश्लेषण सॉफ्टवेयर की सूची में समरूपता बढ़ाने के संभावित तरीके के रूप में पहचाना गया है।[3] क्योंकि पैरारियल समानांतर में कई समय चरणों के लिए संख्यात्मक समाधान की गणना करता है, इसे चरणों में समानांतर विधि के रूप में वर्गीकृत किया गया है।[4] यह समानांतर रनगे-कुट्टा या एक्सट्रपलेशन विधियों जैसी विधि में समानता का उपयोग करने वाले दृष्टिकोणों के विपरीत है, जहां स्वतंत्र चरणों की गणना तरंग रूप विश्राम जैसी सिस्टम विधियों में समानांतर या समानांतर में की जा सकती है।[5][6]
इतिहास
पैरारियल को समय विधि में मल्टीग्रिड विधि या समय अक्ष के साथ प्रत्यक्ष एकाधिक शूटिंग विधि दोनों के रूप में प्राप्त किया जा सकता है।[7] दोनों विचार, समय में मल्टीग्रिड और साथ ही समय एकीकरण के लिए मल्टीपल शूटिंग को अपनाना, 1980 और 1990 के दशक में वापस चले गए।[8][9] पैरारियल एक व्यापक रूप से अध्ययन की जाने वाली विधि है और विभिन्न अनुप्रयोगों के लिए इसका उपयोग और संशोधन किया गया है।[10] प्रारंभिक मूल्य समस्याओं के समाधान को समानांतर करने के विचार और भी पुराने हैं: समानांतर-इन-टाइम एकीकरण विधि का प्रस्ताव करने वाला पहला पेपर 1964 में सामने आया था।[11]
एल्गोरिथम
समस्या
लक्ष्य प्रपत्र की प्रारंभिक मूल्य समस्या को हल करना है
दाहिने हाथ की ओर इसे एक सुचारू (संभवतः अरैखिक) फ़ंक्शन माना जाता है। यह रेखा दृष्टिकोण की विधि में आंशिक अंतर समीकरण के स्थानिक विवेक के अनुरूप भी हो सकता है। हम इस समस्या को अस्थायी आधार पर हल करना चाहते हैं समान दूरी वाले बिंदु , कहाँ और . इस विवेक का पालन करते हुए हमें समय के टुकड़ों से युक्त एक विभाजित समय अंतराल प्राप्त होता है के लिए .
इसका उद्देश्य संख्यात्मक सन्निकटन की गणना करना है सटीक समाधान के लिए एक सीरियल टाइम-स्टेपिंग विधि (जैसे रनगे-कुट्टा) का उपयोग करना जिसमें उच्च संख्यात्मक सटीकता (और इसलिए उच्च कम्प्यूटेशनल लागत) है। हम इस पद्धति को फाइन सॉल्वर कहते हैं , जो प्रारंभिक मान को प्रसारित करता है समय पर एक टर्मिनल मान के लिए समय पर . लक्ष्य का उपयोग करके समाधान की गणना (उच्च संख्यात्मक सटीकता के साथ) करना है जैसे कि हम प्राप्त करते हैं
इस (और सबसे पहले समानांतर में हल करने का प्रयास करने का कारण) समाधान के साथ समस्या यह है कि वास्तविक समय में गणना करना कम्प्यूटेशनल रूप से संभव नहीं है।
यह कैसे काम करता है
प्रारंभिक मूल्य समस्या को हल करने के लिए एकल प्रोसेसर का उपयोग करने के बजाय (जैसा कि शास्त्रीय समय-चरण विधियों के साथ किया जाता है), पैरारियल इसका उपयोग करता है प्रोसेसर. का उद्देश्य उपयोग करना है प्रोसेसर को हल करना होगा समानांतर में छोटी प्रारंभिक मूल्य समस्याएं (प्रत्येक समय स्लाइस पर एक)। उदाहरण के लिए, संदेश पासिंग इंटरफ़ेस आधारित कोड में, प्रक्रियाओं की संख्या होगी, जबकि ओपनएमपी आधारित कोड में, थ्रेड (कंप्यूटिंग) की संख्या के बराबर होगी।
पैरारियल इस प्रारंभिक मूल्य समस्या को समानांतर में हल करने के लिए दूसरी टाइम-स्टेपिंग विधि का उपयोग करता है, जिसे मोटे सॉल्वर के रूप में जाना जाता है . मोटे सॉल्वर ठीक सॉल्वर की तरह ही काम करता है, लंबाई के समय अंतराल पर प्रारंभिक मूल्य का प्रसार करता है हालाँकि, यह इससे कहीं कम संख्यात्मक सटीकता पर ऐसा करता है (और इसलिए बहुत कम कम्प्यूटेशनल लागत पर)। एक मोटे सॉल्वर का होना जो कि फाइन सॉल्वर की तुलना में बहुत कम कम्प्यूटेशनल रूप से महंगा है, पैरारियल के साथ समानांतर गति प्राप्त करने की कुंजी है।
अब से, हम समय पर पैरारियल समाधान को निरूपित करेंगे और पुनरावृत्ति द्वारा .
शून्य पुनरावृत्ति
सबसे पहले, मोटे सॉल्वर को पूरे समय अंतराल पर क्रमिक रूप से चलाएं समाधान के लिए अनुमानित प्रारंभिक अनुमान की गणना करने के लिए:
अनुवर्ती पुनरावृत्तियाँ
इसके बाद, सबसे अद्यतित समाधान मानों से, समानांतर में, प्रत्येक समय स्लाइस पर बढ़िया सॉल्वर चलाएं:
अब प्रेडिक्टर-करेक्टर का उपयोग करके क्रमिक रूप से पैरारियल समाधान मानों को अपडेट करें:
इस स्तर पर, कोई यह निर्धारित करने के लिए एक स्टॉपिंग मानदंड का उपयोग कर सकता है कि क्या समाधान मान अब प्रत्येक पुनरावृत्ति में नहीं बदल रहे हैं। उदाहरण के लिए, कोई इसकी जांच करके यह जांच सकता है कि क्या
और कुछ सहनशीलता . यदि यह मानदंड संतुष्ट नहीं है, तो बाद के पुनरावृत्तियों को समानांतर में फाइन सॉल्वर और फिर भविष्यवक्ता-सुधारक को लागू करके चलाया जा सकता है। हालाँकि, एक बार जब मानदंड संतुष्ट हो जाता है, तो कहा जाता है कि एल्गोरिदम अभिसरण हो गया है पुनरावृत्तियाँ ध्यान दें कि अन्य रोक मानदंड मौजूद हैं और पैरारियल में सफलतापूर्वक परीक्षण किया गया है।
टिप्पणियाँ
पैरारियल को उस समाधान को पुन: पेश करना चाहिए जो फाइन सॉल्वर के क्रमिक अनुप्रयोग द्वारा प्राप्त किया जाता है और अधिकतम में परिवर्तित हो जाएगा पुनरावृत्तियाँ[7]हालाँकि, पैरारियल को स्पीडअप प्रदान करने के लिए, इसे समय स्लाइस की संख्या की तुलना में काफी कम संख्या में पुनरावृत्तियों में परिवर्तित करना होगा, अर्थात। .
पैरारियल पुनरावृत्ति में, कम्प्यूटेशनल रूप से महंगा मूल्यांकन पर समानांतर में निष्पादित किया जा सकता है प्रसंस्करण इकाइयाँ। इसके विपरीत, की निर्भरता पर इसका मतलब है कि मोटे सुधार की गणना क्रमिक क्रम में की जानी है।
आमतौर पर, रंज-कुट्टा विधि का कुछ रूप मोटे और बारीक इंटीग्रेटर दोनों के लिए चुना जाता है, जहां निम्न क्रम का हो सकता है और इससे बड़े समय के चरण का उपयोग किया जा सकता है . यदि प्रारंभिक मूल्य समस्या पीडीई के विवेकाधिकार से उत्पन्न होती है, मोटे स्थानिक विवेक का भी उपयोग कर सकते हैं, लेकिन यह अभिसरण पर नकारात्मक प्रभाव डाल सकता है जब तक कि उच्च क्रम के प्रक्षेप का उपयोग नहीं किया जाता है।[12]
स्पीडअप
कुछ मान्यताओं के तहत, अमदहल के पैरारियल के नियम के लिए एक सरल सैद्धांतिक मॉडल प्राप्त किया जा सकता है।[13] हालाँकि अनुप्रयोगों में ये धारणाएँ बहुत अधिक प्रतिबंधात्मक हो सकती हैं, फिर भी मॉडल पैरारियल के साथ स्पीडअप प्राप्त करने में शामिल ट्रेडऑफ़ को चित्रित करने के लिए उपयोगी है।
सबसे पहले, मान लें कि हर बार स्लाइस करें बिल्कुल शामिल है फाइन इंटीग्रेटर और के चरण मोटे इंटीग्रेटर के चरण. इसमें विशेष रूप से यह धारणा शामिल है कि सभी समय स्लाइस समान लंबाई के होते हैं और मोटे और बारीक इंटीग्रेटर दोनों पूर्ण सिमुलेशन पर एक स्थिर चरण आकार का उपयोग करते हैं। दूसरा, द्वारा निरूपित करें और क्रमशः बारीक और मोटे तरीकों के एक चरण के लिए आवश्यक कंप्यूटिंग समय, और मान लें कि दोनों स्थिर हैं। यह आम तौर पर बिल्कुल सच नहीं है जब टेम्पोरल डिस्क्रेटाइजेशन#इम्प्लिसिट टाइम इंटीग्रेशन विधि का उपयोग किया जाता है, क्योंकि तब रनटाइम इटरेटिव विधि द्वारा आवश्यक पुनरावृत्तियों की संख्या के आधार पर भिन्न होता है।
इन दो मान्यताओं के तहत, ठीक विधि को एकीकृत करने का रनटाइम खत्म हो गया है समय के टुकड़ों को इस प्रकार प्रतिरूपित किया जा सकता है
पैरारियल का उपयोग करने का रनटाइम प्रसंस्करण इकाइयाँ और प्रदर्शन पुनरावृत्तियाँ है
पैरारियल का स्पीडअप तो है
ये दो सीमाएँ मोटे तरीके को चुनने में किए जाने वाले व्यापार को दर्शाती हैं: एक ओर, इसे सस्ता होना होगा और/या पहली सीमा को जितना संभव हो उतना बड़ा बनाने के लिए बहुत बड़े समय के कदम का उपयोग करना होगा, दूसरी ओर पुनरावृत्तियों की संख्या दूसरे बाउंड को बड़ा रखने के लिए इसे नीचे रखना होगा। विशेष रूप से, स्पीडअप#अतिरिक्त विवरण|पैरारियल की समानांतर दक्षता सीमित है
यह आवश्यक पुनरावृत्तियों की संख्या के व्युत्क्रम से है।
काल्पनिक eigenvalues के लिए अस्थिरता
पैरारियल के वेनिला संस्करण में काल्पनिक आइगेनवैल्यूज़ एवं आइगेनवेक्टर्स के साथ समस्याएं हैं।[7]यह आम तौर पर केवल अंतिम पुनरावृत्तियों की ओर ही परिवर्तित होता है, अर्थात दृष्टिकोण , और स्पीडअप हमेशा एक से छोटा होगा. तो या तो पुनरावृत्तियों की संख्या छोटी है और पैरारियल अस्थिर है या, यदि पैरारियल को स्थिर बनाने के लिए पर्याप्त बड़ा है, कोई गति संभव नहीं है। इसका यह भी अर्थ है कि हाइपरबोलिक आंशिक अंतर समीकरण समीकरणों के लिए पैरारियल आमतौर पर अस्थिर है।[14] भले ही गैंडर और वांडेवेले द्वारा औपचारिक विश्लेषण केवल निरंतर गुणांक के साथ रैखिक समस्याओं को कवर करता है, समस्या तब भी उत्पन्न होती है जब पैरारियल को गैर-रेखीय नेवियर-स्टोक्स समीकरणों पर लागू किया जाता है जब चिपचिपाहट गुणांक बहुत छोटा हो जाता है और रेनॉल्ड्स संख्या बहुत बड़ी हो जाती है।[15] पैरारियल को स्थिर करने के लिए विभिन्न दृष्टिकोण मौजूद हैं,[16][17][18] इनमें से एक क्रायलोव-सबस्पेस संवर्धित पैरारियल है।
वेरिएंट
ऐसे कई एल्गोरिदम हैं जो सीधे तौर पर आधारित हैं या कम से कम मूल पैरारियल एल्गोरिदम से प्रेरित हैं।
क्रायलोव-सबस्पेस ने पैरारियल को बढ़ाया
प्रारंभ में ही यह माना गया कि रैखिक समस्याओं के लिए सूचना सूक्ष्म विधि द्वारा उत्पन्न की जाती है मोटे तरीके की सटीकता में सुधार के लिए इसका उपयोग किया जा सकता है .[17]मूल रूप से, यह विचार समानांतर अंतर्निहित समय-एकीकरणकर्ता PITA के लिए तैयार किया गया था,[19] यह एक विधि है जो पैरारियल से निकटता से संबंधित है लेकिन सुधार कैसे किया जाता है इसमें थोड़ा अंतर है। प्रत्येक पुनरावृत्ति में परिणाम मूल्यों के लिए गणना की जाती है के लिए . इस जानकारी के आधार पर, सदिश स्थल
प्रत्येक पैरारियल पुनरावृत्ति के बाद परिभाषित और अद्यतन किया जाता है।[20] के रूप में निरूपित करें ओर्थोगोनल प्रक्षेपण से को . फिर, मोटे तरीके को बेहतर इंटीग्रेटर से बदलें .
जैसे-जैसे पुनरावृत्तियों की संख्या बढ़ती है, space बढ़ेगा और संशोधित प्रचारक अधिक सटीक हो जाएगा. इससे तेजी से अभिसरण हो सकेगा. पैरारियल का यह संस्करण रैखिक अतिशयोक्तिपूर्ण आंशिक अंतर समीकरणों को भी स्थिर रूप से एकीकृत कर सकता है।[21] कम आधार पद्धति पर आधारित अरेखीय समस्याओं का विस्तार भी मौजूद है।[18]
हाइब्रिड पैरारियल वर्णक्रमीय आस्थगित सुधार
वर्णक्रमीय विलंबित सुधार (एसडीसी) के साथ पैरारियल के संयोजन पर आधारित बेहतर समानांतर दक्षता वाली एक विधि [22] एम. मिनियन द्वारा प्रस्तावित किया गया है।[23] यह बेहतर समानांतर दक्षता के लिए लचीलेपन का त्याग करते हुए मोटे और बारीक इंटीग्रेटर के विकल्प को एसडीसी तक सीमित कर देता है। की सीमा के बजाय , संकर विधि में समानांतर दक्षता पर बाध्य हो जाता है
साथ धारावाहिक एसडीसी आधार पद्धति के पुनरावृत्तियों की संख्या और आमतौर पर समानांतर संकर विधि के पुनरावृत्तियों की अधिक संख्या। नॉनलाइनियर मल्टीग्रिड विधि में उपयोग की जाने वाली पूर्ण सन्निकटन योजना को जोड़कर पैरारियल-एसडीसी हाइब्रिड को और बेहतर बनाया गया है। इससे अंतरिक्ष और समय में समानांतर पूर्ण सन्निकटन योजना (पीएफएएसएसटी) का विकास हुआ।[24] पीएफएएसएसटी के प्रदर्शन का अध्ययन पीईपीसी के लिए किया गया है, जो जूलिच सुपरकंप्यूटिंग सेंटर में विकसित एक बार्न्स-हट सिमुलेशन|बार्न्स-हट ट्री कोड आधारित कण सॉल्वर है। IBM ब्लू जीन/पी सिस्टम JUGENE पर सभी 262,144 कोर का उपयोग करके सिमुलेशन से पता चला कि PFASST स्थानिक वृक्ष समानांतरीकरण की संतृप्ति से परे अतिरिक्त गति उत्पन्न कर सकता है।[25]
समय में मल्टीग्रिड कटौती (एमजीआरआईटी)
मल्टीग्रिड रिडक्शन इन टाइम मेथड (एमजीआरआईटी) विभिन्न स्मूथर्स का उपयोग करके कई स्तरों पर मल्टीग्रिड-इन-टाइम एल्गोरिदम के रूप में पैरारियल की व्याख्या को सामान्यीकृत करता है।[26] यह एक अधिक सामान्य दृष्टिकोण है लेकिन मापदंडों की एक विशिष्ट पसंद के लिए यह पैरारियल के बराबर है। एमजीआरआईटी को लागू करने वाली XBraid लाइब्रेरी लॉरेंस लिवरमोर राष्ट्रीय प्रयोगशाला द्वारा विकसित की जा रही है।
पैराएक्सप
ParaExp, Parareal के भीतर घातीय इंटीग्रेटर्स का उपयोग करता है।[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.