Pertanyaan yang diberi tag «c++»

Saya porting kode yang ada dari MATLAB ke C ++ dan memiliki sistem linier untuk menyelesaikan (daripada bentuk yang lebih khas A x = b )xA=bxA=bxA=bAx=bAx=bAx=b Matriks padat, dan berbentuk umum, tetapi tidak lebih besar dari 1000x1000. Jadi dalam MATLAB, solusinya ditemukan oleh fungsi atau notasi forward-slashAAAmrdivide(b,A)x = b/A; Bagaimana …

12 linear-solver  c++  lapack 

Makalah "Expression Templates Revisited: A Performance Analysis of Current Metodologi" di SIAM Journal of Scientific Computing merujuk pada perpustakaan aljabar linear "Blaze". Saya belum pernah mendengarnya sebelumnya, dan sepertinya tidak dapat menemukan referensi online. (Pencarian google yang jelas memberikan kertas di atas kembali.) Jadi apa perpustakaan ini dan di mana …

12 linear-algebra  reference-request  c++ 

Saya memiliki program Mathematica yang melakukan beberapa integral dalam 3 atau 4 dimensi menggunakan QuasiMonteCarlometode ini. Masalahnya adalah, butuh waktu lama yang mengganggu untuk dijalankan, ke titik di mana beberapa perhitungan ini tidak dapat diselesaikan dalam waktu kerja maksimum yang tersedia di kluster HPC kami. Jadi saya sedang mempertimbangkan untuk …

12 libraries  monte-carlo  c++  mathematica 

Saya memecahkan masalah fisik menggunakan skema numerik implisit. Ini menuntun saya untuk memecahkan persamaan linear dengan matriks tridiagonal. Saya telah memberi kode algoritma ini dari Wikipedia. Saya ingin tahu apakah ada perpustakaan yang efisien yang memungkinkan untuk menyelesaikan persamaan jenis ini dengan cara yang dioptimalkan. Catatan penting adalah bahwa matriks …

12 linear-algebra  libraries  c++ 

Bagaimana saya bisa mengganti metode Euler dengan urutan Runge-Kutta ke-4 untuk menentukan gerakan jatuh bebas dalam besaran gravitasi tidak konstan (mis. Jatuh bebas dari 10.000 km di atas tanah)? Sejauh ini saya menulis integrasi sederhana dengan metode Euler: while() { v += getMagnitude(x) * dt; x += v * dt; …

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

12×1212×1212 \times 12QQQJdet(Q)=det(12I−Q−J)(1)det(Q)=det(12I−Q−J)(1)\det(Q) = \det(12I-Q-J) \; \; (1)JJJ Saat ini saya melakukan ini dengan perpustakaan armadillo tetapi ternyata terlalu lambat. Masalahnya adalah bahwa saya perlu melakukan ini untuk satu triliun matriks dan ternyata menghitung dua penentu adalah hambatan dari program saya. Karena itu saya punya dua pertanyaan Apakah ada trik …

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

Saya memiliki subrutin kecil saya sendiri untuk integrasi numerik (quadrature), yang merupakan adaptasi C ++ dari program ALGOL yang diterbitkan oleh Bulirsch & Stoer pada tahun 1967 (Numerische Mathematik, 9, 271-278). Saya ingin meningkatkan ke algoritma yang lebih modern (adaptif) dan bertanya-tanya apakah ada perpustakaan C ++ (gratis) yang menyediakannya. …

10 c++  quadrature 

Saya ingin menggunakan metode urutan ke-8 Runge-Kutta (89) dalam aplikasi mekanika / astrodinamika surgawi, yang ditulis dalam C ++, menggunakan mesin Windows. Karena itu saya bertanya-tanya apakah ada yang tahu perpustakaan / implementasi yang baik yang didokumentasikan dan bebas untuk digunakan? Tidak apa-apa jika ditulis dalam C, selama tidak ada …

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

C ++ 11 memperkenalkan langkah semantik yang dapat, misalnya, meningkatkan kinerja kode dalam situasi di mana C ++ 03 perlu melakukan konstruksi salin atau tugas penyalinan. Artikel ini melaporkan bahwa kode berikut mengalami kecepatan 5x saat dikompilasi dengan C + 11: vector<vector<int> > V; for(int k = 0; k < …

10 finite-element  performance  c++ 

Saat ini saya sedang mengerjakan metode penyelesaian persamaan diferensial yang disebut kolokasi basis-spline. Apa yang saya mengalami masalah adalah membangun metode untuk membangun spline pesanan sewenang-wenang, dengan hubungan dengan kondisi awal B 1 i (x)={ 1Bk+1i(x)=x−xixk+i−xiBki+xk+i+1−xxk+i+1−xi+1Bki+1(x)Bik+1(x)=x−xixk+i−xiBik+xk+i+1−xxk+i+1−xi+1Bi+1k(x) B^{k+1}_{i}(x)= \frac{x-x_i}{x_{k+i}-x_i}B^k_i + \frac{x_{k+i+1}-x}{x_{k+i+1}-x_{i+1}}B^k_{i+1}(x) dan saya mengalami kesulitan bahkan mulai dengan masalah ini, karena …

10 c++  b-spline 

Pada dasarnya, FEM tampaknya menjadi masalah yang cukup banyak "diselesaikan". Ada banyak kerangka kerja yang kuat yang ada, seperti Trilinos, PETSc, FEniCS, Libmesh atau MOOSE. Satu hal yang mereka miliki bersama: Mereka sangat "berat". Pertama, instalasi biasanya sangat menyakitkan. Kedua, antarmuka / API mereka tebal dan berat - Anda harus …

9 finite-element  python  c++  petsc  fenics 

Saya memiliki masalah nilai eigen kubik besar: ( A0+ λ A1+ λ2SEBUAH2+ λ3SEBUAH3) x =0.(SEBUAH0+λSEBUAH1+λ2SEBUAH2+λ3SEBUAH3)x=0.\left(\mathbf{A}_0 + \lambda\mathbf{A}_1 + \lambda^2\mathbf{A}_2 + \lambda^3\mathbf{A}_3\right)\mathbf{x} = 0. Aku bisa mengatasi ini dengan mengkonversi ke masalah nilai eigen linear tetapi akan menghasilkan sistem sebagai besar:32323^2 ⎡⎣⎢- A0000saya000saya⎤⎦⎥⎡⎣⎢xyz⎤⎦⎥= λ ⎡⎣⎢SEBUAH1saya0SEBUAH20sayaSEBUAH300⎤⎦⎥⎡⎣⎢xyz⎤⎦⎥,[-SEBUAH0000saya000saya][xyz]=λ[SEBUAH1SEBUAH2SEBUAH3saya000saya0][xyz],\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; } Di atas adalah kode C ++ …

9 c++  performance 

Saat ini saya mencoba untuk menyelesaikan masalah minimisasi terbatas nonlinear seperti yang diterapkan dalam fungsi "fmincon" matlab. Harapan saya adalah, meminimalkan (fun1, x0, uB, lB, fun2) di mana x0 adalah keadaan awal, fun1 adalah fungsi yang perlu diminimalkan, uB adalah batas atas, lB adalah batas bawah dan fungsi fun2 adalah …

9 c++  nonlinear-programming 

Jadi saya sudah berdebat apakah saya harus repot-repot belajar Python. Dari berbicara dengan profesor saya, Matlab tampaknya menjadi bahasa umum yang digunakan dalam matematika terapan / ilmu komputasi sejauh menyangkut akademisi ; sementara di industri , profesor saya (terutama mereka yang telah bekerja di industri) mengatakan bahwa belajar c ++ …

9 matlab  python  c++