A Waring decomposition of a (homogeneous) polynomial f is a minimal sum of powers of linear forms expressing f. Under certain conditions, such a decomposition is unique. We discuss some algorithms to compute the Waring decomposition, which are linked to the equation of certain secant varieties and to eigenvectors of tensors. In particular we explicitly decompose a general cubic polynomial in three variables as the sum of five cubes (Sylvester Pentahedral Theorem).
- पेपर ID: 1103.0203
- शीर्षक: टेंसर के आइजेनवेक्टर और वारिंग अपघटन के लिए एल्गोरिदम
- लेखक: Luke Oeding, Giorgio Ottaviani
- वर्गीकरण: math.AG (बीजगणितीय ज्यामिति)
- प्रकाशन समय: 6 नवंबर 2012 (arXiv v2)
- पेपर लिंक: https://arxiv.org/abs/1103.0203
यह पेपर सजातीय बहुपद f के वारिंग अपघटन का अध्ययन करता है, अर्थात् f को रैखिक रूपों की घातों के न्यूनतम योग के रूप में प्रस्तुत करना। विशिष्ट शर्तों के तहत, यह अपघटन अद्वितीय है। पेपर वारिंग अपघटन की गणना के लिए एल्गोरिदम पर चर्चा करता है, जो विशिष्ट सेकेंट वैरायटी के समीकरणों और टेंसर के आइजेनवेक्टर से संबंधित हैं। विशेष रूप से, पेपर तीन चर के सामान्य घन बहुपद को पाँच घनों के योग में स्पष्ट रूप से अपघटित करता है (सिल्वेस्टर पंचमुख प्रमेय)।
इस पेपर द्वारा समाधान की जाने वाली मूल समस्या यह है: दिए गए बहुपद को देखते हुए, रैखिक रूपों की घातों के योग के रूप में इसका न्यूनतम अपघटन कैसे खोजें? इसे गणितीय रूप से वारिंग अपघटन समस्या कहा जाता है।
- सैद्धांतिक महत्व: वारिंग अपघटन बीजगणितीय ज्यामिति में एक शास्त्रीय समस्या है, जो सममित टेंसर अपघटन से घनिष्ठ रूप से संबंधित है
- अनुप्रयोग मूल्य: टेंसर अपघटन बीजगणितीय सांख्यिकी, रसायन विज्ञान, कंप्यूटर विज्ञान, विद्युत इंजीनियरिंग, तंत्रिका विज्ञान, भौतिकी और मनोमिति में व्यापक रूप से लागू होता है
- कम्प्यूटेशनल आवश्यकता: 2010 में इटली के मोनोपोली टेंसर अपघटन और अनुप्रयोग सम्मेलन की सामान्य थीम कुशल और विश्वसनीय टेंसर अपघटन एल्गोरिदम की आवश्यकता थी
- उत्प्रेरक मैट्रिक्स विधि: द्विपद रूपों के मामले में पूरी तरह सफल, लेकिन n≥2 के लिए केवल आंशिक रूप से सफल
- बर्बर विधि: समय और मेमोरी खपत विशाल है, अक्सर विफल होती है
- संख्यात्मक विधियाँ: अधिकांश टेंसर समस्याएँ अत्यंत कठिन हैं और आमतौर पर अघुलनशील हैं
पेपर का उद्देश्य बीजगणितीय ज्यामिति को एल्गोरिदम के आधार के रूप में उपयोग करना, वेक्टर बंडल तकनीकों और टेंसर आइजेनवेक्टर की अवधारणा को जोड़ते हुए, नए कुशल और मजबूत टेंसर अपघटन एल्गोरिदम विकसित करना है।
- नई एल्गोरिदम फ्रेमवर्क: कोज़ुल समतलीकरण और टेंसर आइजेनवेक्टर पर आधारित नई एल्गोरिदम (एल्गोरिदम 4) प्रस्तावित करता है, जो पेपर का मुख्य परिणाम है
- सैद्धांतिक सुधार: उत्प्रेरक मैट्रिक्स विधि की प्रयोज्यता सीमा पर Iarrobino-Kanev का सुधार (प्रमेय 2.4)
- शास्त्रीय समस्या का समाधान: सिल्वेस्टर पंचमुख प्रमेय का रचनात्मक प्रमाण और एल्गोरिदम कार्यान्वयन प्रदान करता है
- सफलता की शर्तें: एल्गोरिदम की सफलता के लिए पर्याप्त शर्तें देता है (प्रमेय 3.5 और 5.4)
- ज्यामितीय व्याख्या: सामान्यीकृत आइजेनवेक्टर संख्या पर Cartwright-Sturmfels के परिणाम के लिए नया ज्यामितीय प्रमाण प्रदान करता है
- कार्यान्वयन कोड: Macaulay2 कार्यान्वयन प्रदान करता है, जो आगे के अनुसंधान के लिए शुरुआती बिंदु है
सममित टेंसर f ∈ S^d V (d-डिग्री सजातीय बहुपद के समतुल्य) को देखते हुए, इसका वारिंग अपघटन खोजें:
f=∑i=1rci(vi)d
जहाँ v_i ∈ V 1-डिग्री रैखिक रूप हैं, c_i ∈ ℂ गुणांक हैं, r न्यूनतम संख्या है जिसके लिए यह अपघटन मौजूद है (सममित रैंक कहा जाता है)।
f ∈ S^d V के लिए, 0 ≤ a ≤ n, 1 ≤ m ≤ d-1 को ठीक करते हुए, रैखिक मानचित्र बनाएँ:
Pf:Hom(SmV,∧aV)→Hom(∧n−aV,Sd−m−1V)
जब f = v^d हो, तो परिभाषित करें:
Pvd(M)(w)=(M(vm)∧v∧w)(vd−m−1)
लेम्मा 3.3: वेक्टर v ∈ V टेंसर M का आइजेनवेक्टर है यदि और केवल यदि M ∈ ker(P_{v^d})।
परिभाषा 3.2: M ∈ Hom(S^m V, ∧^a V) को देखते हुए, वेक्टर v ∈ V को टेंसर M का आइजेनवेक्टर कहा जाता है, यदि:
M(vm)∧v=0
पेपर बीजगणितीय वैरायटी X पर वेक्टर बंडल E के प्रतिनिधित्व का उपयोग करता है, f ∈ W पर निर्भर रैखिक मानचित्र बनाता है:
Af:H0(E)→H0(E∗⊗L)∗
प्रस्ताव 4.1: यदि f = ∑v_i, Z = {v_1,...,v_r}, तो:
- H^0(I_Z ⊗ E) ⊆ ker A_f
- जब H^0(E^* ⊗ L) → H^0(E^* ⊗ L|_Z) आच्छादक हो तो समानता成立
एल्गोरिदम 4 (सामान्य टेंसर अपघटन एल्गोरिदम):
- मानचित्र A_f बनाएँ
- ker A_f की गणना करें
- ker A_f का आधार स्थान Z' खोजें
- रैखिक प्रणाली f = ∑c_i v_i^d को हल करें
f ∈ S^5 ℂ^3 के लिए:
- 18×18 ब्लॉक मैट्रिक्स P_f बनाएँ
- ker P_f की गणना करें, सामान्य तत्व M चुनें
- 2×2 उप-सारणिकों के शून्य समुच्चय के माध्यम से 7 आइजेनवेक्टर खोजें
- रैखिक प्रणाली को हल करके अद्वितीय अपघटन प्राप्त करें
f ∈ S^3 ℂ^4 के लिए:
- a=2, m=1 सेट करके P_f बनाएँ
- 9-आयामी कर्नल स्पेस की गणना करें
- 5 आइजेनवेक्टर खोजें (पंचमुख के 5 समतलों के अनुरूप)
- अद्वितीय अपघटन प्राप्त करें
प्रमेय 2.4: उत्प्रेरक मैट्रिक्स विधि की सुधारी गई सीमा
- सम डिग्री: r ≤ (n+m choose n) - n - 1
- विषम डिग्री: r ≤ (n+m-1 choose n)
प्रमेय 3.5: n=2 मामले में कोज़ुल विधि की सीमा
- यदि 2r ≤ m² + 3m + 4, तो एल्गोरिदम सफल है
- यदि 2r ≤ m² + 4m + 2, तो एल्गोरिदम अद्वितीय न्यूनतम अपघटन देता है
प्रमेय 3.4: सामान्य टेंसर M ∈ Hom(S^m V, ∧^a V) के आइजेनवेक्टर की संख्या:
- a = 1: (m^{n+1} - 1)/(m - 1)
- a = n-1: ((m+1)^{n+1} + (-1)^n)/(m + 2)
पेपर Macaulay2 कार्यान्वयन प्रदान करता है, जिसमें शामिल हैं:
- उत्प्रेरक मैट्रिक्स एल्गोरिदम: फ़ाइल "cat method.m2"
- कोज़ुल समतलीकरण एल्गोरिदम: फ़ाइल "General Kappa Method.m2"
उत्प्रेरक मैट्रिक्स विधि सफलता सीमा:
- n=2: (d=6, s≤8)
- n=3: (d=6, s≤16)
- n=4: (d=6, s≤16)
कोज़ुल विधि सफलता सीमा:
- n=2: (d=6, s≤7)
- n=3: (d=5, s≤11)
- n=4: (d=5, s≤14)
- एल्गोरिदम प्रभावकारिता: सैद्धांतिक सीमा के भीतर, एल्गोरिदम अद्वितीय वारिंग अपघटन सफलतापूर्वक खोज सकता है
- कम्प्यूटेशनल दक्षता: बर्बर विधि की तुलना में, नई एल्गोरिदम पंचमुख उदाहरण में 1 सेकंड से कम में पूरी होती है, जबकि बर्बर विधि मेमोरी और समय सीमा के कारण विफल होती है
- सीमा सटीकता: प्रयोग सैद्धांतिक सीमा की सटीकता को सत्यापित करते हैं
- पंचम वक्र: 7 पंचम घातों के योग में सफलतापूर्वक अपघटित
- पंचमुख: तीन-चर घन बहुपद को 5 घनों के योग में सफलतापूर्वक अपघटित
- परिमेय चतुर्घात वक्र: उप-उत्पाद के रूप में, 7 सामान्य बिंदुओं से गुजरने वाले अद्वितीय परिमेय चतुर्घात वक्र का नया प्रतिनिधित्व खोजा
- सिल्वेस्टर उत्प्रेरक मैट्रिक्स विधि: 19वीं सदी में विकसित, द्विपद मामले में पूरी तरह सफल
- Alexander-Hirschowitz प्रमेय: सामान्य सममित टेंसर की रैंक निर्धारित करता है
- अद्वितीयता परिणाम: Chiantini-Ciliberto आदि का कार्य
- Cartwright-Sturmfels: टेंसर आइजेनवेक्टर गणना सूत्र
- Brachat आदि: अर्ध-Hankel ऑपरेटर का उपयोग करके संख्यात्मक विधि
- Kolda-Mayo: टेंसर आइजेनवेक्टर की पुनरावृत्तिमूलक गणना विधि
- एकीकृत फ्रेमवर्क: शास्त्रीय अद्वितीयता मामलों को संभालने के लिए एकीकृत विधि प्रदान करता है
- ज्यामितीय अंतर्दृष्टि: टेंसर अपघटन को बीजगणितीय ज्यामिति में सेकेंट वैरायटी सिद्धांत से जोड़ता है
- व्यावहारिक एल्गोरिदम: कार्यान्वयन योग्य एल्गोरिदम विकसित करता है, जो विशिष्ट सीमा में सफलता की गारंटी देता है
- प्रयोज्यता सीमा: एल्गोरिदम केवल तब सफल होता है जब रैंक काफी छोटा हो और टेंसर सामान्य हो
- कम्प्यूटेशनल जटिलता: बड़े टेंसर के लिए, समस्या अभी भी कठिन है
- संख्यात्मक स्थिरता: संख्यात्मक विधियों की स्थिरता पर आगे के अनुसंधान की आवश्यकता है
- संख्यात्मक विधियाँ: GPU त्वरण के साथ आइजेनवेक्टर गणना को जोड़ना
- निम्न-रैंक सन्निकटन: मैट्रिक्स मामले की नकल करने वाली छोटी आइजेनमान उन्मूलन विधि
- अधिक सामान्य मामले: आंशिक सममित टेंसर तक विस्तार
- सैद्धांतिक नवाचार: बीजगणितीय ज्यामिति की वेक्टर बंडल तकनीकों को टेंसर अपघटन में सफलतापूर्वक लागू करता है
- विधि एकीकरण: कई शास्त्रीय समस्याओं के लिए एकीकृत उपचार फ्रेमवर्क प्रदान करता है
- पूर्ण कार्यान्वयन: पूर्ण एल्गोरिदम कार्यान्वयन और परीक्षण प्रदान करता है
- सीमा सटीकता: एल्गोरिदम सफलता के लिए सटीक सैद्धांतिक सीमा देता है
- प्रयोज्यता प्रतिबंध: एल्गोरिदम सफलता सीमा अपेक्षाकृत सीमित है
- कम्प्यूटेशनल जटिलता: उच्च-आयामी मामलों के लिए गणना अभी भी कठिन है
- संख्यात्मक समस्याएँ: परिमेय समस्याओं पर अधिक कार्य की आवश्यकता है
- सैद्धांतिक योगदान: टेंसर अपघटन सिद्धांत के लिए नया ज्यामितीय दृष्टिकोण प्रदान करता है
- व्यावहारिक मूल्य: विशिष्ट सीमा में विश्वसनीय एल्गोरिदम प्रदान करता है
- अनुवर्ती अनुसंधान: आगे की संख्यात्मक विधियों और GPU कार्यान्वयन के लिए आधार प्रदान करता है
यह विधि विशेष रूप से उपयुक्त है:
- छोटे से मध्यम आकार के सममित टेंसर अपघटन के लिए
- सैद्धांतिक अनुसंधान में जहाँ सटीक अपघटन की आवश्यकता हो
- एल्गोरिदम विकास के लिए शुरुआती बिंदु और बेंचमार्क के रूप में
यह पेपर टेंसर अपघटन क्षेत्र में महत्वपूर्ण सैद्धांतिक और एल्गोरिदमिक योगदान देता है, विशेष रूप से बीजगणितीय ज्यामिति विधियों को व्यावहारिक कम्प्यूटेशनल समस्याओं में लागू करने के संदर्भ में नई दिशा खोलता है।