कम्प्यूटेशनल विज्ञान

पृष्ठभूमि मैं तरल सिद्धांत से ऑर्नस्टीन-ज़र्निक समीकरण के एक संस्करण को हल कर रहा हूं । संक्षेप में, समस्या को निर्धारित बिंदु समस्या को हल करने के रूप में दर्शाया जा सकता है , जहां एक पूर्णांक-बीजीय ऑपरेटर है और समाधान फ़ंक्शन (OZ प्रत्यक्ष सहसंबंध समारोह) है। मैं पिकार्ड पुनरावृत्ति …

13 iterative-method  convergence 

मेरे पास लगभग n ~ 100000 नोड्स के साथ दो ग्राफ हैं। दोनों ग्राफ़ में, प्रत्येक नोड ठीक 3 अन्य नोड्स से जुड़ा होता है, ताकि आसन्न मैट्रिक्स सममित और बहुत विरल हो। मुश्किल हिस्सा यह है कि मुझे आसन्न मैट्रिक्स के सभी eigenvalues ​​की आवश्यकता है, लेकिन eigenvectors की …

13 linear-algebra  sparse-matrix  eigenvalues  matrix 

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

13 spectral-method  precision  condition-number 

मैं Hesthaven / Warburton पुस्तक का उपयोग करके DG-FEM विधियों के पीछे सिद्धांत सीख रहा हूं और 'संख्यात्मक प्रवाह' की भूमिका के बारे में थोड़ा भ्रमित हूं। मैं माफी माँगता हूँ अगर यह एक बुनियादी सवाल है, लेकिन मैंने देखा है और इसका संतोषजनक उत्तर नहीं मिला है। पर विचार …

13 finite-element  discontinuous-galerkin 

मैं एक FEM सिमुलेशन विकसित कर रहा हूं। शुरुआती परीक्षण के लिए, मैं सरल स्व-लिखित मेशर और मेष ग्राफ के विज़ुअलाइज़ेशन का उपयोग करूंगा। लेकिन मैं एक मौजूदा मेज़र द्वारा उत्पन्न डेटा का उपयोग करने के लिए अपना कार्यक्रम तैयार करना चाहता हूं और इसे मौजूदा विज़ुअलाइज़ेशन टूल में आउटपुट …

13 finite-element  mesh  graph-theory 

लहर समीकरण में: c2∇⋅∇u(x,t)−∂2u(x,t)∂t2=f(x,t)c2∇⋅∇u(x,t)−∂2u(x,t)∂t2=f(x,t)c^2 \nabla \cdot \nabla u(x,t) - \frac{\partial^2 u(x,t)}{\partial t^2} = f(x,t) एकीकृत करने से पहले हम पहले परीक्षण फ़ंक्शन v(x,t)v(x,t)v(x,t) से गुणा क्यों करते हैं ?

13 finite-element 

हमें फॉर्म की समस्या के साथ शुरू करते हैं (L+k2)u=0(L+k2)u=0(\mathcal{L} + k^2) u=0 दी गई सीमा स्थितियों ( ड्यूरिचलेट , न्यूमैन , रॉबिन , आवधिक , बलोच-आवधिक ) के एक सेट के साथ । यह कुछ ऑपरेटर एलL\mathcal{L} के लिए कुछ जियोमेट्री और सीमा स्थितियों के तहत आइजेनवल और ईजीनवेक्टर्स …

13 pde  finite-element  eigensystem  verification 

मल्टीग्रेड (एमजी) का उपयोग एक रेखीय प्रणाली को हल करने के लिए किया जा सकता है एक प्रारंभिक अनुमान निर्माण और लिए निम्नलिखित दोहराते हुए : अभिसरण तक:एक्स 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 

मैंने अजगर 3 (सुन्न का उपयोग करके) में एक पिछड़े-यूलर सॉल्वर को लागू किया है। अपनी सुविधा के लिए और एक अभ्यास के रूप में, मैंने एक छोटा सा कार्य भी लिखा, जो कि ढाल के परिमित अंतर की गणना करता है ताकि मुझे हमेशा जेकबिएन को विश्लेषणात्मक रूप से …

13 pde  stability  numpy  scipy  newton-method 

मुझे एक फ़ंक्शन को संख्यात्मक रूप से लागू करने में समस्या हो रही है। यह इस तथ्य से ग्रस्त है कि बड़े इनपुट मूल्यों पर परिणाम बहुत बड़ी संख्या में बहुत छोटी संख्या है। मुझे यकीन नहीं है कि अगर विपत्तिपूर्ण रद्दीकरण सही शब्द है तो कृपया मुझे सही करें …

13 python  numerical-analysis  floating-point 

हमें एक आयाम में एक चिकनी प्रारंभिक हालत और गर्मी समीकरण पर विचार करें: ∂tu=∂xxu∂tu=∂xxu \partial_t u = \partial_{xx} u खुला अंतराल में ]0,1[]0,1[]0,1[ , और हमें लगता है कि हम परिमित अंतर के साथ संख्यानुसार इसे हल करना चाहते हैं। मुझे पता है कि मेरी समस्या को अच्छी तरह …

13 pde  boundary-conditions  parabolic-pde 

मैं जिस तरह से matlab संख्यात्मक एकीकरण बनाम Scipy संभालता है पर कुछ निराशा का अनुभव कर रहा हूँ। मैं अपने परीक्षण कोड में निम्नलिखित अंतरों का पालन करता हूं: मतलाब का संस्करण मेरे अजगर के मुकाबले औसतन 24 गुना तेज चलता है ! मैथलैब का संस्करण बिना किसी चेतावनी …

13 matlab  python  quadrature  numba 

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

13 linear-solver  sparse-matrix  complexity  multigrid 

मान लें कि मैं http://www.kc.gov/users/computational-systems/edison/configuration पर 4k घंटों के लिए 100k कोर पर एक सुपर कंप्यूटर संगणना चलाता हूं , नेटवर्क पर लगभग 4 PB डेटा का आदान-प्रदान करता हूं और I / 4 TB का प्रदर्शन करता हूं। ओ गणना सभी पूर्णांक है, इसलिए परिणाम या तो सही या …

13 error-estimation 

मैं अपने द्वारा संख्यात्मक रूप से पीडीई को हल करने के बारे में जानने की कोशिश कर रहा हूं। मैं कुछ समय के लिए परिमित अंतर विधि (एफडीएम) के साथ शुरुआत कर रहा हूं क्योंकि मैंने सुना है कि एफडीएम पीडीई के लिए कई संख्यात्मक तरीकों का मूल है। अब …

13 pde  finite-difference  boundary-conditions