linear-solver पर टैग किए गए जवाब

मेरे पास कई चुनौतीपूर्ण गैर-उत्तल वैश्विक अनुकूलन समस्याओं को हल करने के लिए है। वर्तमान में मैं MATLAB के ऑप्टिमाइज़ेशन टूलबॉक्स काfmincon()'sqp' उपयोग करता हूं (विशेष रूप से, एल्गोरिथ्म = ) के साथ, जो काफी प्रभावी है । हालाँकि, मेरा अधिकांश कोड पायथन में है, और मैं पायथन में भी …

77 python  optimization  nonlinear-programming  linear-algebra  linear-solver  sparse-matrix  ode  automatic-differentiation  efficiency  matrix  symbolic-computation  pde  multiphysics  operator-splitting  fluid-dynamics  pde  multigrid  pde  hyperbolic-pde  quadrature  precision  numerical-analysis  fluid-dynamics  interpolation  c  ode  testing  data-sets  pde  fluid-dynamics  hyperbolic-pde  pde  finite-element  fluid-dynamics  stability  incompressible  pde  hyperbolic-pde  fluid-dynamics  pde  precision  floating-point  convex-optimization  conditioning  pde  hyperbolic-pde  stability  implicit-methods  explicit-methods  fluid-dynamics  accuracy  pde  hyperbolic-pde  petsc  linear-algebra 

विरल रैखिक सिस्टम अनुप्रयोगों में बढ़ती आवृत्ति के साथ बदल जाते हैं। इन प्रणालियों को हल करने के लिए किसी को चुनने के लिए बहुत सी दिनचर्याएँ हैं। उच्चतम स्तर पर, प्रत्यक्ष (जैसे विरल गॉसियन उन्मूलन या चोल्स्की अपघटन, विशेष आदेश एल्गोरिदम, और मल्टीफ्रंटल तरीकों के साथ) और पुनरावृत्त (जैसे …

49 linear-algebra  linear-solver  sparse-matrix 

मेरी जानकारी के लिए, रैखिक समीकरणों की एक प्रणाली को हल करने के 4 तरीके हैं (अधिक होने पर मुझे सही करें): यदि सिस्टम मैट्रिक्स एक पूर्ण-रैंक वर्ग मैट्रिक्स है, तो आप Cramer के नियम का उपयोग कर सकते हैं; सिस्टम मैट्रिक्स के व्युत्क्रम या छद्म बिंदु की गणना करें; …

31 linear-algebra  linear-solver  iterative-method 

मेरे पास और । विरल है और साथ बहुत बड़ा है (कई मिलियन के आदेश पर हो सकता है।) G एक n \ n मीटर लंबा मैट्रिक्स है जिसमें m छोटा ( 1 \ lt m \ lt 1000 ) और प्रत्येक कॉलम हो सकता है केवल एक ही है …

22 linear-algebra  linear-solver  sparse-matrix  krylov-method 

वहाँ कई अलग-अलग पुस्तकालय हैं जो समीकरणों की एक विरल रैखिक प्रणाली को हल करते हैं, हालांकि मुझे यह पता लगाना मुश्किल है कि अंतर क्या हैं। : जहां तक मेरा यह बता सकें कि तीन प्रमुख संकुल हैं Trilinos , PETSc , और इंटेल MKL । वे सभी मैट्रिक्स …

21 linear-algebra  linear-solver  libraries 

मुझे उन मामलों में दिलचस्पी है जहाँ Conjugate ग्रेडिएंट GMRES विधि से बहुत बेहतर काम करता है। सामान्य तौर पर, सीजी एसपीडी (सममित-सकारात्मक-निश्चित) के कई मामलों में बेहतर विकल्प है क्योंकि इसमें सीजी के लिए अभिसरण दर पर कम भंडारण और सैद्धांतिक बाध्यता की आवश्यकता होती है जो जीएमआरईएस से …

17 linear-solver  conjugate-gradient  gmres 

मैं विशेष रूप से विषम हेल्महोल्त्ज़ के लिए बहुस्तरीय उलटा-आधारित आईएलयू पूर्व- क्रम के धारावाहिक प्रदर्शन से बहुत प्रभावित हूं, लेकिन मुझे कोई भी खुला स्रोत कार्यान्वयन नहीं मिल पा रहा है। विशेष रूप से, ILUPACK बायनेरिज़ को स्वतंत्र रूप से शिक्षाविदों के लिए उपलब्ध कराता है, लेकिन ऐसा नहीं …

15 pde  algorithms  linear-solver  libraries  preconditioning 

मान लीजिए कि हमारे पास एक रैखिक प्रणाली है और हम इसके कंडीशनिंग के बारे में कुछ नहीं जानते हैं और समाधान के बारे में कोई प्रारंभिक जानकारी नहीं है। हम नेत्रहीन रूप से गाऊसी उन्मूलन लागू करते हैं और कुछ समाधान प्राप्त करते हैं । क्या यह निर्धारित करना …

14 linear-solver  condition-number  conditioning 

मल्टीग्रेड (एमजी) का उपयोग एक रेखीय प्रणाली को हल करने के लिए किया जा सकता है एक प्रारंभिक अनुमान निर्माण और लिए निम्नलिखित दोहराते हुए : अभिसरण तक:एक्स 0 आई = 0 , 1 ।।अ x = खAx=bAx=bएक्स0x0x_0मैं = 0 , 1 ।।i=0,1..i=0,1.. अवशिष्टआरमैं= बी - ए एक्समैंri=b−Axir_i = b-Ax_i …

13 linear-algebra  linear-solver  multigrid  krylov-method 

मुझे आश्चर्य है कि अगर थॉमस एल्गोरिथ्म सबसे तेज़ तरीका है (निश्चित रूप से?) एक सममित तिरछे वर्चस्व को हल करने के लिए एल्गोरिदमिक जटिलता (LAPACK आदि जैसे कार्यान्वयन पैकेज की तलाश में नहीं) के संदर्भ में। मुझे पता है कि थॉमस एल्गोरिथ्म और मल्टीग्रिड दोनों जटिलता हैं, लेकिन शायद …

13 linear-solver  sparse-matrix  complexity  multigrid 

एकल परिशुद्धता फ़्लोटिंग पॉइंट संख्या आधी मेमोरी लेती हैं और आधुनिक मशीनों (यहां तक ​​कि GPU पर भी ऐसा लगता है) का संचालन उनके साथ दोहरी सटीकता की तुलना में लगभग दोगुनी गति से किया जा सकता है। कई FDTD कोड जो मैंने विशेष रूप से एकल परिशुद्धता अंकगणित और …

13 matrix  linear-solver  precision 

Gmres या bicgstab की तरह विधि में यह एक पूर्ववर्ती के रूप में एक और krylov विधि का उपयोग करने के लिए आकर्षक हो सकता है। आखिरकार वे मैट्रिक्स-मुक्त तरीके से और समानांतर वातावरण में लागू करना आसान है। उदाहरण के लिए, एक कौल कुछ का उपयोग करते हैं (मान …

13 linear-algebra  linear-solver  preconditioning  krylov-method  gmres 

मैं MATLAB से C ++ में एक मौजूदा कोड पोर्ट कर रहा हूं और (बजाय अधिक विशिष्ट रूप से ) हल करने के लिए एक रैखिक प्रणाली हैxA=bएक्सए=खxA=bAx=bएएक्स=खAx=b मैट्रिक्स , घने और सामान्य रूप का है, लेकिन 1000x1000 से बड़ा नहीं है। तो MATLAB में, समाधान फ़ंक्शन या फ़ॉरवर्ड-स्लैश नोटेशन …

12 linear-solver  c++  lapack 

मुझे एक ही विरल रेखीय प्रणाली (300x300 से 1000x1000) को कई दाहिने हाथों (300 से 1000) के साथ हल करने की आवश्यकता है। इस पहली समस्या के अलावा, मैं विभिन्न प्रणालियों को भी हल करना चाहूंगा, लेकिन एक ही गैर-शून्य तत्वों (बस अलग-अलग मूल्यों) के साथ, यह निरंतर स्पार्सिटी पैटर्न …

12 linear-algebra  petsc  linear-solver  sparse-matrix 

मैं एक समस्या को हल करने के लिए MATLAB का उपयोग कर रहा हूं जिसमें प्रत्येक टाइमस्टेप पर को हल करना शामिल है , जहां समय के साथ बदल जाता है। अभी, मैं MATLAB का उपयोग करके इसे पूरा कर रहा हूं :बीअ x = खAx=b\mathbf{A} \mathbf{x}=\mathbf{b}खb\mathbf{b}mldivide x = A\b …

12 linear-algebra  matlab  linear-solver  matrix  linear-system