Lapack: DSYEVR गैर-ऑर्थोगोनल वैक्टर लौटाता है

हाय मिरो, बग 126 को गिटहब मुद्दे के रूप में जोड़ने के लिए धन्यवाद। यह भी खूब रही। Osni ( @oamarques ) #100 पर काम कर रहे हैं। अपना परीक्षण कोड हमारे साथ साझा करने के लिए धन्यवाद। हां, लैपैक के लिए एक नए परीक्षण सूट पर काम करने की योजना है और किसी भी खराब मैट्रिक्स का स्वागत है। ओस्नी के पास ऐसे कई मैट्रिसेस हैं और परीक्षण के लिए बुनियादी ढांचे के साथ-साथ इनका एक संग्रह विकसित करना लक्ष्य होगा। किसी भी स्थिति में, बग 126 को GitHub मुद्दे के रूप में जोड़ने के लिए धन्यवाद। यह भी खूब रही। चीयर्स, जूलियन।

यह बग डीएसटीईएमआर में प्रयुक्त विभाजन मानदंड से संबंधित है, जो संबंधित स्रोत फ़ाइल में लाइन 606 से शुरू होता है (जैसा कि LAPACK 3.7.0 में है)। इनपुट वैरिएबल TRYRAC के आधार पर, सबरूटीन DLARRR INFO=1 लौटाता है (ऐसा लगता है कि यह DLARRR के अपेक्षित व्यवहार के विपरीत है), जो बाद में इनपुट के रूप में पारित (विशेष) ट्राइडिगोनल मैट्रिक्स के हानिकारक विभाजन को ट्रिगर करता है। _अगर हम DLARRR की उपेक्षा करते हैं और "विभाजन मानदंड सेट करें" से ठीक पहले IINFO=0 सेट करते हैं, तो लौटाए गए वैक्टर ऑर्थोगोनल होते हैं, DLARRF के नए संस्करण का उपयोग करते समय, समस्या 100 देखें। हालांकि, यह इस बग का उचित समाधान नहीं है: हम यह पता लगाने की जरूरत है कि DLARRR का सही व्यवहार क्या होना चाहिए।_)

नमस्ते, इस बग को हल करने से संबंधित कोई खबर?

@oamarques क्या अंतर्निहित DSTEMR के बजाय DSTEGR का उपयोग इस समस्या को कम करता है?

@uihsnv : अच्छा सवाल। STEGR का अब कोई रखरखाव नहीं किया जाता है। इसके साथ मुद्दे हैं। हम केवल एसटीईएमआर पर काम कर रहे हैं और उसका रखरखाव कर रहे हैं। लेकिन, हाँ, यह देखना एक अच्छा विचार होगा कि इस समस्या पर STEGR क्या करता है।

@langou ओह! मुझे नहीं पता था!! मैं अपने कोड में STEGR के कुछ उदाहरणों का उपयोग कर रहा हूं। क्या मुझे एसटीईएमआर का उपयोग करने के लिए स्विच करना चाहिए?

साथ ही, मैंने @miroi द्वारा साझा किए गए परीक्षण में कुछ संशोधन किए हैं, जिसे यदि वह स्वीकार करता है, तो आप अंतर देख सकते हैं।

मैं अपने कोड में STEGR के कुछ उदाहरणों का उपयोग कर रहा हूं। क्या मुझे एसटीईएमआर का उपयोग करने के लिए स्विच करना चाहिए?

मुझे भी ऐसा ही लगता है। मुझे उम्मीद है कि @oamarques जल्द ही हमारे साथ इस पर कार्रवाई करेगा।

साथ ही, मैंने @miroi द्वारा साझा किए गए परीक्षण में कुछ संशोधन किए हैं, जिसे यदि वह स्वीकार करता है, तो आप अंतर देख सकते हैं।

संशोधित परीक्षण को @oamarques . के साथ साझा करना बहुत अच्छा होगा

हैलो, इस बग से संबंधित कोई खबर? और इसे "दस्तावेज़ीकरण" क्यों लेबल किया गया था?

फिर भी, ताजा लैपैक sthis त्रुटि देता है (gfortran 4.9.2 के साथ परीक्षण किया गया, https://github.com/miroi/lapack-dsyevr-test):

~~~
DIRAC के विकर्णीकरण दिनचर्या के परीक्षण कार्यक्रम में आपका स्वागत है!
पढ़ें एन = 9 एनजेड = 1
विकर्ण के लिए मैट्रिक्स पढ़ा गया था
आरएस सबरूटीन के लिए सरणियाँ आवंटित
ईसपैक का आरएस स्वदेशी मूल्य:
1 -1.50000000000000062
2 -1.5000000000000002
3 -1.4999999999999982
4 0.49999999999999772
5 0.499999999999999856
6 0.49999999999999994
7 0.5000000000000033
8 0.5000000000000122
9 2.499999999999969
*ईआईएसपैक रुपये*
यू^{+} ए यू - ईपीएस?= 0> मानदंड/डायग:0.1912डी-15 मानदंड/ऑफडिएग:0.1008डी-15
यू^{+} यू-आई?= 0> मानदंड/डायग:0.2344डी-15 मानदंड/ऑफडिएग:0.7239डी-16यू यू^{+} - मैं ?= 0> मानदंड/डायग:0.9869डी-16 मानदंड/ऑफडिएग:0.9320डी-16
लैपैक डीएसवाईईवीआर eigenvalues:
1 -1.50000000000000056
2 -1.5000000000000009
3 -1.4999999999999973
4 0.49999999999999523
5 0.499999999999999867
6 0.499999999999999939
7 0.499999999999999950
8 0.50000000000000100
9 2.4999999999999862
*लैपैक डीएसवाईईवीआर*
यू^{+} ए यू-ईपीएस?= 0> मानदंड/डायग:0.1807डी-14 मानदंड/ऑफडिएग:0.9012डी-07
यू^{+} यू-आई?= 0> मानदंड/डायग:0.2591डी-15 मानदंड/ऑफडिएग:0.1802डी-06यू यू^{+} - मैं ?= 0> मानदंड/डायग:0.7543डी-06 मानदंड/ऑफडिएग:0.4951डी-06
LAPACK DSTEGR eigenvalues:
1 -1.50000000000000056
2 -1.5000000000000009
3 -1.4999999999999973
4 0.49999999999999523
5 0.499999999999999867
6 0.499999999999999939
7 0.499999999999999950
8 0.50000000000000100
9 2.4999999999999862
*लैपैक डीएसवाईटीआरडी+डीएसटीईजीआर*
यू^{+} ए यू - ईपीएस?= 0> मानदंड/डायग:0.1617डी+01 मानदंड/ऑफडिएग:0.1644डी+00
यू^{+} यू-आई?= 0> मानदंड/डायग:0.1234डी-15 मानदंड/ऑफडिएग:0.1802डी-06यू यू^{+} - मैं ?= 0> मानदंड/ डायग:0.5411डी -06 मानदंड/ ऑफडिएग:0.2301डी-06

~~~

इसे फिर से पढ़ने पर, ऐसा लगता है कि इसे ठीक करना # 100 को ठीक करने पर टिका है जो एक अनंत लूप को ठीक करने के प्रयास के कारण हुआ था (महसूस कर रहा हूं कि मैं हाल ही में उन लोगों के लिए लंगड़ा कम-तकनीकी समाधान प्रस्तावित करने में विशेषज्ञता प्राप्त कर रहा हूं) और पूरी बात एक अनंत लूप में चली गई एक आदर्श समाधान क्या हो सकता है, इसके बारे में संदेह। क्या किसी ने 3.5.0 पर वापस जाकर यह जांचा कि #100 में बताए गए फ़्लैग वेरिएबल को बदलने से पहले यह मामला कैसा व्यवहार करता है?