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

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

14 finite-element  fenics 

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

14 ode  numerical-analysis 

मेरे पास 3D में लाखों डेटा बिंदुओं पर चलने वाला डेटासेट है। मैं जो गणना कर रहा हूं, उसके लिए मुझे एक त्रिज्या में प्रत्येक डेटा बिंदु पर पड़ोसी (रेंज खोज) की गणना करने की आवश्यकता है, एक फ़ंक्शन को फिट करने का प्रयास करें, फिट के लिए त्रुटि की …

14 algorithms  data-structures  nearest-neighbors  spatial-data 

मैं पायथन में निम्नलिखित अभिव्यक्ति को लागू करना चाहता हूं: जहां और आकार , और संख्यात्मक सारणी हैं। आकार का एक अरै सरणी । आकार लगभग 10000 तक हो सकता है, और फ़ंक्शन एक आंतरिक लूप का हिस्सा है जिसका मूल्यांकन कई बार किया जाएगा, इसलिए गति महत्वपूर्ण है।x y …

14 python  numpy 

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

14 linear-solver  condition-number  conditioning 

मान लें कि मेरे पास एक रैखिक प्रणाली Ax=bAx=bA x = b , जो सभी लिए एक उपयुक्त क्रायलोव विधि (जैसे सीजी या जीएमआरईएस) का उपयोग करके जल्दी से परिवर्तित होती है bbb। यदि BBB कम रैंक साथ एक मैट्रिक्स है rrr, तो सिस्टम पर एक ही क्रायलोव विधि (A+B)x=b(A+B)x=b(A …

14 krylov-method 

मैं एक प्रतिलिपि प्रस्तुत करने योग्य कागज लिख रहा हूं, और कागज में कम्प्यूटेशनल परिणाम हैं जो एक पायथन स्क्रिप्ट (एक समान MATLAB स्क्रिप्ट लगभग समान परिणाम उत्पन्न करता है) द्वारा उत्पन्न होते हैं। मुझे लगता है कि पेपर को पाठकों के लिए समझना आसान होगा यदि वे कोड में …

14 python  software  reproducibility 

प्रारंभिक सीमा मूल्य PDEs के संख्यात्मक समाधान में, अंतरिक्ष में समानता को नियोजित करना बहुत आम है । यह समय के विवेकाधिकार में समानता के कुछ प्रकार को नियोजित करने के लिए बहुत कम आम है , और यह समानता आमतौर पर बहुत अधिक सीमित है। मुझे कोड और प्रकाशित …

14 pde  parallel-computing  time-integration 

मैं पीडीई को हल करने के लिए परिमित अंतर विधियों का उपयोग करते समय सीमा की स्थिति का चयन करने में मदद करने के लिए कुछ संसाधनों को खोजने की कोशिश कर रहा हूं। वर्तमान में मेरे पास जितनी भी किताबें और नोट्स हैं, वे सभी समान हैं: सीमाओं की …

14 pde  finite-difference  boundary-conditions  advection 

मैं एक बड़े विरल मैट्रिक्स (सैकड़ों हजारों पंक्तियों) के स्पेक्ट्रम ( सभी eigenvalues) की गणना करना चाहता हूं । यह कठिन है। मैं एक अनुमान के लिए समझौता करने को तैयार हूं। क्या ऐसा करने के लिए अनुमानित तरीके हैं? जबकि मैं इस प्रश्न के सामान्य उत्तर की आशा करता …

14 linear-algebra  eigensystem  graph-theory  approximation 

मैंने एक विशेष PDE सिस्टम को हल करने के लिए फोरट्रान प्रोग्राम को कोड करने पर (तरल प्रवाह / दहन का वर्णन करते हुए) पिछले कुछ महीने बिताए हैं। मैंने नवीनतम-मानक फोरट्रान का उपयोग करने का प्रयास किया और नए ओआरओपी कैपबिलिटीज़ आधुनिक फोरट्रान के पास हैं। मैं अपने दम …

14 fortran 

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

14 pde  finite-difference 

रिमेज़ एल्गोरिथ्म एक बहुप्रतिक्षित पुनरावृत्त दिनचर्या है जो न्यूनतम पॉलिसियों में एक बहुपद द्वारा एक फ़ंक्शन को अनुमानित करता है। लेकिन, निक ट्रेफेथेन के रूप में [1] इसके बारे में कहते हैं: इनमें से अधिकांश [कार्यान्वयन] कई साल पीछे चले जाते हैं और वास्तव में, उनमें से अधिकांश सामान्य सर्वश्रेष्ठ …

14 optimization  approximation 

जब मुझे उच्च आदेश केंद्र अंतर सन्निकटन का उपयोग करना हो तो मुझे एक समस्या है: (−ui+2,j+16ui+1,j−30ui,j+16ui−1,j−ui−2,j12)(−ui+2,j+16ui+1,j−30ui,j+16ui−1,j−ui−2,j12)\left(\frac{-u_{i+2,j}+16u_{i+1,j}-30u_{i,j}+16u_{i-1,j}-u_{i-2,j}}{12}\right) पोइसन समीकरण के लिए (uxx+uyy=0)(uxx+uyy=0)(u_{xx}+u_{yy}=0) एक वर्गाकार डोमेन में जिसमें सीमा की स्थितियाँ हैं: Δ एक्स = Δ y = 0.1u(0,y)=u(x,0)=u(x,1)=0,u(1,y)=sinπyu(0,y)=u(x,0)=u(x,1)=0,u(1,y)=sin⁡πyu(0,y)=u(x,0)=u(x,1)=0,u(1,y)=\sin \pi y Δx=Δy=0.1Δx=Δy=0.1\Delta{x}=\Delta{y}=0.1 जब मैं डोमेन के अंदरूनी बिंदुओं का मूल्य प्राप्त …

14 pde  finite-difference 

जहां कहीं भी मैंने देखा है, पेट्सक ट्यूटोरियल / दस्तावेज आदि कहते हैं कि यह रैखिक बीजगणित के लिए उपयोगी है और आमतौर पर निर्दिष्ट करता है कि विरल प्रणालियों को लाभ होगा। घने मैट्रिस के बारे में क्या? मैं घने ए के लिए को हल करने के बारे में …

14 linear-algebra  petsc  blas