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) *...
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>...
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?...
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...
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)={
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...
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...
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...
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...
Untuk pekerjaan penelitian saya, saya telah mengembangkan perpustakaan C ++ kecil yang bertujuan memfasilitasi komunikasi antara kode komputasi C ++ dan Octave / Matlab (ketika yang terakhir digunakan untuk tujuan pasca-pemrosesan). Saya ingin merilis perpustakaan seperti itu di bawah GPL di...
Sebagian besar pemrograman saya adalah kode penelitian satu kali dalam C untuk saya gunakan sendiri. Saya belum pernah mendistribusikan kode apa pun selain dari kolaborator dekat. Saya telah mengembangkan suatu algoritma yang saya terbitkan dalam jurnal ilmiah. Saya ingin memberikan kode sumber dan...
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,...