2010-12-01 18 views
8

मैं वर्तमान में एक सी ++ पर काम कर रहा हूँ - बड़े, विरल रेखीय बीजगणित समस्याओं के लिए आधारित पुस्तकालय (हाँ, मैं इस तरह के कई पुस्तकालयों मौजूद पता है, लेकिन मैं यह जानने के लिए अपने खुद के रोलिंग कर रहा हूँ ज्यादातर पुनरावृत्त हलकों, स्पैस भंडारण कंटेनर, आदि के बारे में ..)।परीक्षण मैट्रिक्स/पुनरावृत्ति रैखिक के लिए सिस्टम के लिए खोज रहे solver

मैं बिंदु जहां मैंने अपने अन्य प्रोग्रामिंग परियोजनाओं के भीतर मेरी समाधानकर्ताओं उपयोग कर रहा हूँ, और समस्याओं है कि मेरे अपने नहीं हैं के खिलाफ समाधानकर्ताओं परीक्षण करना चाहते हैं के लिए कर रहा हूँ। मुख्य रूप से, मैं सममित स्पैस सिस्टम के खिलाफ परीक्षण करना चाहता हूं जो सकारात्मक निश्चित हैं।

Matrix Market UF Sparse Matrix Collection

कहा जा रहा है, मैं अभी तक अच्छा परीक्षण मैट्रिक्स है कि पूरे प्रणाली प्रणाली मैट्रिक्स और आरएचएस शामिल की किसी भी स्रोत नहीं मिला है: मैं इस तरह के रूप में इस तरह के सिस्टम मैट्रिक्स के लिए कई स्रोतों मिल गया है। परिणाम जांचने के लिए यह बहुत अच्छा होगा। जहाँ मैं इस तरह के पूर्ण सिस्टम पा सकते हैं, या वैकल्पिक रूप पर कोई सुझाव, मैं क्या प्रणाली मैट्रिक्स मैं ऑनलाइन प्राप्त कर सकते हैं के लिए एक "अच्छा" आरएचएस उत्पन्न करने के लिए कर सकता है? मैं वर्तमान में यादृच्छिक मूल्यों, या सभी के साथ एक मैट्रिक्स भर रहा हूं, लेकिन संदेह है कि यह जरूरी नहीं है कि यह सबसे अच्छा तरीका है।

+0

"ऐसे कई पुस्तकालयों मौजूद": नहीं वास्तव में (कम से कम मूल रूप से सी के लिए लिखा ++) यूएफ विरल मैट्रिक्स संग्रह के लिए मैं इस लिंक के बजाय सुझाव देना चाहेंगे। बड़े स्पैस मैट्रिस से निपटने वाले फोर्टन पुस्तकालयों के लिए स्वच्छ रैपर लिखना ईमानदार होने के लिए पहले से ही एक चुनौती है। –

+0

हालांकि, मुझे याद है कि कुछ शोध पत्रों में बीमार स्थिति वाले परीक्षण मामलों के कुछ संदर्भ हैं, लेकिन आईआईआरसी वे एसपीडी स्पैस मैट्रिस के लिए नहीं थे। आपकी स्थिति में परीक्षण के मामलों का निर्माण करने का एक आसान तरीका एक यादृच्छिक एन एक्स पी मैट्रिक्स एम लेना है, इसे अपने स्वयं के हस्तांतरण से गुणा करें, और कुछ लैम्ब्डा * पहचान जोड़ें ताकि यह अचूक हो। लेकिन यह स्पैर मैट्रिस का उत्पादन नहीं करेगा। –

+0

इसके अलावा, यूएफ स्पैर मैट्रिक्स संग्रह में क्या गलत है? कुछ यादृच्छिक आरएचएस लेना मेरे लिए बिल्कुल ठीक लगता है। –

उत्तर

0

मैं अभी तक उपयोग नहीं किया है, मैं करने के बारे में हूँ, लेकिन GiNAC सबसे अच्छी बात मैं सी ++ के लिए मिल गया है की तरह लगता है। यह सीएएस के लिए मेपल के पीछे उपयोग की जाने वाली लाइब्रेरी है, मुझे इसके प्रदर्शन के बारे में पता नहीं है।

http://www.ginac.de/

0

यह निर्दिष्ट करने के लिए समस्याओं आप किस प्रकार की सुलझाने कर रहे हैं ... विभिन्न समस्याओं वैधता की जांच के लिए किसी भी काम का हो करने के लिए विभिन्न आरएचएस की आवश्यकता होगी अच्छी तरह से करना होगा ..... मैं क्या सुझाव देंगे DUNE Numerics की तरह कुछ परियोजनाओं (मैं इस अधिकार पर काम कर रहा हूँ), FENICS, deal.ii जो पहले से ही मैट्रिक हल करने के लिए उपयोग कर रहे हैं समाधानकर्ताओं से कुछ उदाहरण कोड प्राप्त है ... आम तौर पर वे कुछ कार्यक्षमता उत्पादन के लिए किसी तरह का में अपने मैट्रिक्स होगा फ़ाइल (ड्यून न्यूमेरिक्स में matlab-compliant फ़ाइलों में आउटपुट मैट्रिस और आरएचएस आउटपुट करने की कार्यक्षमता है)।

यह आप तो अपने समाधानकर्ताओं को खिलाने सकते हैं .. और उसके बाद फिर उनके पुस्तकालयों कार्यक्षमता का उपयोग उत्पादन डेटा बनाने के लिए (जैसे DUNE अंकीय एक VTK स्वरूप का उपयोग करता) ... यही था, तो आप डेटा का विश्लेषण करने का प्रयास करेंगे शक्तिशाली उपकरण का उपयोग .....

आपको उन पुस्तकालयों को संकलित करने और उपयोग करने के बारे में कुछ सीखना पड़ सकता है ... लेकिन यह बहुत कुछ नहीं है ... और मेरा मानना ​​है कि आपको जो कार्यक्षमता मिलेगी वह लायक होगी समय निवेश किया गया ......

मुझे लगता है कि आपके पुस्तकालयों का परीक्षण करने के लिए भी एक अच्छी तरह से परिभाषित और उचित जटिल समस्या पर्याप्त होनी चाहिए .... वास्तव में दो Ax = b समस्याओं के लिए एक और कुल्हाड़ी = cbx (eigenvalue समस्याओं) के लिए एक और ....

1

मैं एक दाएँ हाथ की ओर वेक्टर एक पूर्वनिर्धारित 'लक्ष्य' से प्राप्त का उपयोग कर सुझाव है कि समाधान एक्स:

b = A*x 

तो फिर तुम solver से एक लक्ष्य समाधान, एक्स, और एक परिणामस्वरूप समाधान, एक्स है,। इसका मतलब है कि आप त्रुटि (लक्ष्य के अंतर और परिणामी समाधान) के साथ-साथ अवशिष्ट (ए * एक्स - बी) की तुलना कर सकते हैं।

ध्यान दें कि एक पुनरावर्तक सॉल्वर के सावधानीपूर्वक मूल्यांकन के लिए आपको यह भी विचार करना होगा कि प्रारंभिक एक्स के लिए क्या उपयोग करना है।

मैट्रिक्स की ऑनलाइन संग्रह मुख्य रूप से बाएं हाथ की ओर मैट्रिक्स होते हैं, लेकिन कुछ दाएँ हाथ-पक्षों को शामिल करते हैं और भी कुछ समाधान वैक्टर भी .:

http://www.cise.ufl.edu/research/sparse/matrices/rhs.txt

वैसे है,

http://www.cise.ufl.edu/research/sparse/matrices/

 संबंधित मुद्दे

  • कोई संबंधित समस्या नहीं^_^