क्या पेट्स ने कभी स्पार्स मैट्रिक्स गणित के लिए लैप पुस्तकालयों का उपयोग किया है?


10

क्या PETSc को बाहरी BLAS / LAPACK लाइब्रेरी के साथ संकलित करने से विरल मैट्रिस पर प्रदर्शन प्रभावित होता है, या क्या यह केवल घना मैट्रिक्स गणित के लिए उन पुस्तकालयों का उपयोग करता है?


आप पूर्वव्यापी के लिए एलयू का उपयोग कर सकते हैं।
stali

जवाबों:


10

पेट्सक कुछ वेक्टर प्राइमरी के लिए BLAS का उपयोग करता है, लेकिन ये आम तौर पर मेमोरी बैंडविड्थ द्वारा सीमित होते हैं और "ऑप्टिमाइज़ेशन" में अधिक भिन्नता नहीं होती है, इसलिए यह बहुत अधिक प्रदर्शन अंतर नहीं करता है।

यह कुछ विश्लेषण के लिए लैपैक का भी उपयोग करता है, जैसे लैंक्ज़ोस या अर्नोल्डी के आइजनवेल्स और विलक्षण मूल्यों का अनुमान, लेकिन ये आम तौर पर प्रदर्शन-संवेदनशील नहीं होते हैं।

घने "स्तर 3" ऑपरेशन आम तौर पर केवल एक प्रदर्शन-संवेदनशील संदर्भ में दिखाई देते हैं, जब तीसरे पक्ष के पुस्तकालयों (जैसे MUMPS, SuperLU, UMFPACK) से विरल प्रत्यक्ष सॉल्वर का उपयोग करते हैं, जिस स्थिति में अंततः घने समस्याओं का कारण बनता है जो कि लाभ के लिए पर्याप्त रूप से बड़ी होती है बीएलएएस बुला रहा है।

यदि आप बड़ी समस्याओं के लिए लागू इन विरल प्रत्यक्ष सॉल्वरों पर भरोसा करते हैं, तो यह एक ट्यून किए गए BLAS कार्यान्वयन के साथ निर्माण के लायक है, अन्यथा यह बहुत कम अंतर बनाता है।


जेड, वितरित मैट्रिसिस के मामले में, पेट्सन ईजन / विलक्षण मूल्यों के आकलन के लिए लैपैक का उपयोग कैसे करता है? क्या यह तब ब्लॉक-वार दिनचर्या या स्कैल्प है?
अलेक्जेंडर

अर्नोल्डी पुनरावृति पुनरावृत्तियों की संख्या के बराबर आयाम का एक हेसेनबर्ग मैट्रिक्स उत्पन्न करता है (30 या 100 कहते हैं)। हेसेनबर्ग मैट्रिक्स के eigen- या एकवचन मूल्यों को LAPACK का उपयोग करके प्रत्येक प्रक्रिया पर अनावश्यक रूप से गणना की जाती है। यह आकार में छोटा होने के कारण इसे समानांतर में करने की कोशिश में बहुत तेज है और डेटा पहले से ही वितरित किया गया है। लैंक्ज़ोस के लिए भी ऐसा ही है।
जेड ब्राउन
हमारी साइट का प्रयोग करके, आप स्वीकार करते हैं कि आपने हमारी Cookie Policy और निजता नीति को पढ़ और समझा लिया है।
Licensed under cc by-sa 3.0 with attribution required.