मैं कुछ घने, बीमार हालत वाले मैट्रिसेस को विकर्ण करने की कोशिश कर रहा हूं। मशीन परिशुद्धता में, परिणाम गलत हैं (नकारात्मक आईजेनवेल्स लौटाते हैं, eigenvectors के पास अपेक्षित समरूपता नहीं है)। मैंने मैथमेटिका के ईगेंससिस्टम [] कार्य पर स्विच किया ताकि मनमानी परिशुद्धता का लाभ उठाया जा सके, लेकिन गणना बेहद धीमी है। मैं किसी भी समाधान के लिए खुला हूं। क्या ऐसे पैकेज / एल्गोरिदम हैं जो अच्छी तरह से वातानुकूलित समस्याओं के अनुकूल हैं? मैं पूर्वशर्त पर विशेषज्ञ नहीं हूं, इसलिए मुझे यकीन नहीं है कि यह कितना मदद कर सकता है। अन्यथा, मैं सोच सकता हूं कि सभी समानांतर रूप से मनमाने ढंग से सटीक eigenvalue सॉल्वर हैं, लेकिन मैं गणितज्ञ, MATLAB और C ++ से परे किसी भी चीज से परिचित नहीं हूं।
समस्या पर कुछ पृष्ठभूमि देने के लिए, मैट्रिसेस बड़े हैं, लेकिन विशाल नहीं (4096x4096 से 32768x32768 सबसे अधिक)। वे वास्तविक, सममित हैं, और eigenvalues 0 और 1 (अनन्य) के बीच बंधे हुए हैं, कई eigenvalues 0 के बहुत करीब हैं और 1 के करीब कोई भी नहीं है। मैट्रिक्स अनिवार्य रूप से एक कन्वेंशन ऑपरेटर है। मुझे अपने सभी मैट्रिसेस को तिरस्कृत करने की आवश्यकता नहीं है, लेकिन जितना बड़ा मैं जा सकता हूं, उतना बेहतर होगा। मेरे पास कई प्रोसेसर के साथ कंप्यूटिंग क्लस्टर तक पहुंच है और कंप्यूटिंग क्षमताओं को वितरित किया है।
धन्यवाद