Jak mogę zastąpić metodę Eulera czwartym rzędem Runge-Kutty, aby określić ruch swobodnego spadku o nie stałej wielkości grawitacyjnej (np. Swobodny spadek z 10 000 km nad ziemią)? Do tej pory napisałem prostą integrację metodą Eulera: while() { v += getMagnitude(x) * dt; x += v * dt; time +=...
C ++ 11 wprowadza semantykę przenoszenia, która może na przykład poprawić wydajność kodu w sytuacjach, w których C ++ 03 musiałby wykonać konstrukcję kopii lub przypisanie kopii. W tym artykule opisano, że następujący kod przyspiesza 5-krotnie po skompilowaniu z C + 11: vector<vector<int>...
Chciałbym użyć metody Runge-Kutta 8. rzędu (89) w niebiańskiej aplikacji do mechaniki / astrodynamiki, napisanej w C ++, za pomocą komputera z systemem Windows. Dlatego zastanawiam się, czy ktoś zna dobrą bibliotekę / implementację, która jest udokumentowana i może być używana bezpłatnie? Jest w...
Mam własną małą procedurę integracji numerycznej (kwadraturę), która jest adaptacją C ++ programu ALGOL opublikowanego przez Bulirsch & Stoer w 1967 r. (Numerische Mathematik, 9, 271-278). Chciałbym uaktualnić do bardziej nowoczesnego (adaptacyjnego) algorytmu i zastanawiać się, czy istnieją...
Obecnie pracuję nad metodą rozwiązywania równań różniczkowych zwaną kolokacją podstawa-splajn. Mam problem z budowaniem metody budowania splajnu dowolnego porządku z relacją przy czym warunek początkowy i mam problem nawet z początkiem tego problemu, ponieważ jest on rekurencyjny, może...
Zastanawiam się więc, czy powinienem zawracać sobie głowę nauką języka Python. Po rozmowie z moimi profesorami Matlab wydaje się być wspólnym językiem stosowanym w matematyce stosowanej / naukach komputerowych w odniesieniu do środowiska akademickiego ; podczas gdy w przemyśle moi profesorowie...
Mam duży problem wartości własnej sześciennej: (A0+λA1+λ2A2+λ3A3)x=0.(A0+λA1+λ2A2+λ3A3)x=0.\left(\mathbf{A}_0 + \lambda\mathbf{A}_1 + \lambda^2\mathbf{A}_2 + \lambda^3\mathbf{A}_3\right)\mathbf{x} = 0. Mógłbym to rozwiązać, przechodząc na liniowy problem wartości własnych, ale spowodowałoby to,...
Obecnie próbuję rozwiązać problem minimalizacji ograniczonej nieliniowo, jak zaimplementowano w funkcji matlab „fmincon”. Moje oczekiwania to minimalizacja (fun1, x0, uB, lB, fun2), gdzie x0 to stan początkowy, fun1 to funkcja, którą należy zminimalizować, uB to górne granice, lB to dolne granice,...
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; } Powyżej...
Do moich prac badawczych opracowałem małą bibliotekę C ++ mającą na celu ułatwienie komunikacji między kodami obliczeniowymi C ++ a Octave / Matlab (gdy ten ostatni jest wykorzystywany do celów przetwarzania końcowego). Chciałbym wydać taką bibliotekę na licencji GPL na darmowej stronie...
Większość mojego programowania to jednorazowe kody badawcze w C na własny użytek. Nigdy nie rozpowszechniałem żadnego kodu poza bliskimi współpracownikami. Opracowałem algorytm, który publikuję w czasopiśmie naukowym. Chcę podać kod źródłowy i być może kod wykonywalny w internetowym suplemencie do...
Zasadniczo FEM wydaje się problemem, który jest w zasadzie „rozwiązany”. Istnieje wiele potężnych frameworków, takich jak Trilinos, PETSc, FEniCS, Libmesh lub MOOSE. Łączy je jedno: są wyjątkowo „ciężkie”. Po pierwsze, instalacja zwykle jest bardzo bolesna. Po drugie, ich interfejs / API jest...