c++ पर टैग किए गए जवाब

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

12 linear-solver  c++  lapack 

SIAM जर्नल ऑफ साइंटिफिक कंप्यूटिंग में पेपर "एक्सप्रेशन टेम्प्लेट्स रिविजिटेड: ए परफॉर्मेंस एनालिसिस ऑफ करंट मेथडोलोजी" "ब्लेज़" लीनियर अलजेब्रा लाइब्रेरी का संदर्भ देता है। मैंने इसके बारे में पहले नहीं सुना है, और ऑनलाइन संदर्भ खोजने के लिए प्रतीत नहीं हो सकता है। (स्पष्ट Google खोजें उपरोक्त पेपर को वापस …

12 linear-algebra  reference-request  c++ 

मेरे पास एक गणितज्ञ कार्यक्रम है जो QuasiMonteCarloविधि का उपयोग करके 3 या 4 आयामों में कुछ अभिन्न प्रदर्शन करता है । समस्या यह है कि इसे चलाने के लिए एक कष्टप्रद समय लगता है, इस बिंदु पर जहां इनमें से कुछ गणना हमारे एचपीसी क्लस्टर पर उपलब्ध अधिकतम नौकरी …

12 libraries  monte-carlo  c++  mathematica 

मैं अंतर्निहित संख्यात्मक योजना का उपयोग करके एक शारीरिक समस्या को हल कर रहा हूं। यह मुझे त्रिदलीय मैट्रिक्स के साथ एक रैखिक समीकरण को हल करने की ओर ले जाता है। मैंने इस एल्गोरिथम को विकिपीडिया से कोडित किया है । मुझे आश्चर्य है कि अगर एक कुशल पुस्तकालय …

12 linear-algebra  libraries  c++ 

मैं रनवे-कुट्टा 4 के क्रम से यूलर विधि को कैसे बदल सकता हूं ताकि निरंतर गुरुत्वाकर्षण परिमाण में मुक्त गिरावट गति का निर्धारण किया जा सके (उदाहरण के लिए जमीन से 10 000 किमी ऊपर से गिरना)। अब तक मैंने यूलर विधि द्वारा सरल एकीकरण लिखा है: while() { v …

11 numerical-analysis  c++  runge-kutta  integral-equations 

मुझे एक मैट्रिक्स जो कि सममित, उल्टा, सकारात्मक निश्चित और घना है। मुझे परीक्षण करने की आवश्यकता है अगर जहां जे सभी मैट्रिक्स हैं।१२ × १२12×1212 \times 12क्यूQQdet ( Q ) = det ( 12 मैं- क्यू - जे)( 1 )det(Q)=det(12I−Q−J)(1)\det(Q) = \det(12I-Q-J) \; \; (1)जेJJ मैं वर्तमान में आर्मडिलो …

11 linear-algebra  matrix  c++  graph-theory  c 

न्यूमेरिकल इंटीग्रेशन (क्वाडरेचर) के लिए मेरे पास अपना छोटा सबरूटीन है, जो कि 1967 में बुलिरस्च एंड स्टोअर द्वारा प्रकाशित ALGOL प्रोग्राम का C ++ रूपांतरण है (न्यूमेरिस मैथिक, 9, 27-2-278)। मैं एक अधिक आधुनिक (अनुकूली) एल्गोरिथ्म में अपग्रेड करना चाहूंगा और आश्चर्यचकित करूंगा कि क्या कोई (फ्री) सी ++ …

10 c++  quadrature 

मैं एक celestial यांत्रिकी / astrodynamics एप्लिकेशन में Runge-Kutta 8th ऑर्डर विधि (89) का उपयोग करना चाहता हूं, जो C ++ में लिखा गया है, विंडोज मशीन का उपयोग कर रहा है। इसलिए मुझे आश्चर्य है कि अगर किसी को एक अच्छी लाइब्रेरी / कार्यान्वयन पता है जो प्रलेखित और …

10 c++  reference-request  libraries  c  runge-kutta 

C ++ 11 चालित शब्दार्थ का परिचय देता है, उदाहरण के लिए, उन परिस्थितियों में कोड प्रदर्शन में सुधार कर सकता है जहाँ C ++ 03 को कॉपी निर्माण या कॉपी असाइनमेंट करने की आवश्यकता होगी। यह लेख बताता है कि निम्नलिखित कोड 5x गति का अनुभव करता है जब …

10 finite-element  performance  c++ 

फिलहाल मैं एक अंतर समीकरण हल करने की विधि पर काम कर रहा हूं, जिसे आधार-स्पलाइन कोलोलेशन कहा जाता है। मुझे जो परेशानी हो रही है, वह साथ, एक मनमाना क्रम बनाने के लिए एक विधि का निर्माण कर रहा है। के साथ प्रारंभिक स्थिति और मुझे इस समस्या से …

10 c++  b-spline 

मूल रूप से, FEM एक समस्या है जो बहुत अधिक "हल" है। कई शक्तिशाली चौखटे मौजूद हैं, जैसे त्रिलीनोस, पेट्सक, फेनीसीएस, लिबमेश या एमओओएसई। एक चीज जो उनके पास आम है: वे हैं बेहद "भारी-वजन"। सबसे पहले, स्थापना सामान्य रूप से सुपर दर्दनाक है। दूसरा, उनका इंटरफ़ेस / एपीआई मोटा …

9 finite-element  python  c++  petsc  fenics 

मुझे एक बड़ी क्यूबिक स्वदेशी समस्या है: (ए0+ λए1+λ2ए2+λ3ए3) x =0।(ए0+λए1+λ2ए2+λ3ए3)एक्स=0।\left(\mathbf{A}_0 + \lambda\mathbf{A}_1 + \lambda^2\mathbf{A}_2 + \lambda^3\mathbf{A}_3\right)\mathbf{x} = 0. मैं एक रेखीय eigenvalue समस्या में परिवर्तित करके इसे हल कर सकता था, लेकिन इसका परिणाम रूप में एक प्रणाली के रूप में होगा:32323^2 ⎡⎣⎢-ए0000मैं000मैं⎤⎦⎥⎡⎣⎢एक्सyz⎤⎦⎥= λ⎡⎣⎢ए1मैं0ए20मैंए300⎤⎦⎥⎡⎣⎢एक्सyz⎤⎦⎥,[-ए0000मैं000मैं][एक्सyz]=λ[ए1ए2ए3मैं000मैं0][एक्सyz],\begin{bmatrix} -\mathbf{A}_0 & 0 & 0 …

9 c++  eigensystem  eigenvalues 

double trap(double func(double), double b, double a, double N) { double j; double s; double h = (b-a)/(N-1.0); //Width of trapezia double func1 = func(a); double func2; for (s=0,j=a;j<b;j+=h){ func2 = func(j+h); s = s + 0.5*(func1+func2)*h; func1 = func2; } return s; } ऊपर की func()सीमाओं के बीच 1 …

9 c++  performance 

मैं वर्तमान में matlab "fmincon" फ़ंक्शन में कार्यान्वित किए गए गैर-समसामयिक कम से कम समस्या को हल करने का प्रयास कर रहा हूं। मेरी अपेक्षाएँ हैं, कम से कम (fun1, x0, uB, lB, fun2) जहाँ x0 प्रारंभिक अवस्था है, fun1 वह कार्य है जिसे कम से कम करने की आवश्यकता …

9 c++  nonlinear-programming 

इसलिए मैं इस बात पर बहस कर रहा हूं कि मुझे पायथन सीखने में परेशान होना चाहिए या नहीं। मेरे प्राध्यापकों के साथ बोलने से लेकर, मतलाब सामान्य भाषा में प्रयुक्त गणित / कम्प्यूटेशनल विज्ञान में प्रयुक्त होने वाली भाषा प्रतीत होती है जहाँ तक अकादमिक का संबंध है; उद्योग …

9 matlab  python  c++