W moim programie doktoranckim z zakresu obliczeń pracujemy prawie wyłącznie w C ++ i Fortran. Wygląda na to, że niektórzy profesorowie wolą jeden od drugiego. Zastanawiam się, który z nich jest „lepszy” lub czy w pewnych okolicznościach jest lepszy od
Czy dobrym pomysłem jest użycie vector<vector<double>>(przy użyciu std) do utworzenia klasy macierzy dla wysokowydajnego naukowego kodu komputerowego? Jeśli odpowiedź brzmi „nie”. Czemu?
Problem Obecnie pracuję nad symulacją elementów skończonych Naviera Stokesa i chciałbym zbadać wpływ różnych parametrów. Niektóre parametry są określone w pliku wejściowym lub za pomocą opcji wiersza poleceń; inne parametry są dostarczane jako flagi w pliku Makefile, więc mój kod musi być...
Chcesz poprawić ten post? Podaj szczegółowe odpowiedzi na to pytanie, w tym cytaty i wyjaśnienie, dlaczego Twoja odpowiedź jest poprawna. Odpowiedzi bez wystarczającej ilości szczegółów mogą być edytowane lub usuwane. Wersja 3.0 standardu MPI formalnie usunęła interfejs...
Załóżmy, że pracuję nad kodem naukowym w C ++. W niedawnej dyskusji z kolegą argumentowano, że szablony wyrażeń mogą być bardzo złe, potencjalnie czyniąc oprogramowanie kompilowalnym tylko w niektórych wersjach gcc. Podobno problem ten wpłynął na kilka kodów naukowych, jak wspomniano w podtytułach...
Szukam książek lub artykułów, postów na blogu lub innych opublikowanych materiałów, które odnoszą się konkretnie do wykorzystania nowoczesnych funkcji C ++ (semantyka ruchu, STL, iteratory, leniwa ocena itp.) W informatyce naukowej. Czy możesz coś zasugerować? Myślę, że te nowe funkcje ułatwią...
Nauczanie studentów na poziomie licencjackim kursu zatytułowanego Wprowadzenie do programowania komputerowego. Jestem trochę zmieszany. W dziedzinie fizyki obliczeniowej naukowcy używają C / C ++ lub Python lub Fortran, CUDA itp ... czas na zbudowanie bazy. Czego powinienem użyć? Wiem, że możesz...
Załóżmy, że mam funkcję, która przyjmuje jako dane wejściowe kilka wartości zmiennoprzecinkowych (pojedyncze lub podwójne), wykonuje pewne obliczenia i generuje wyjściowe wartości zmiennoprzecinkowe (także pojedyncze lub podwójne). Pracuję przede wszystkim z MSVC 2008, ale planuję także współpracę...
Szukam biblioteki tensorowej C ++, która obsługuje kod zależny od wymiarów. W szczególności muszę wykonywać operacje wzdłuż każdego wymiaru (do 3), np. Obliczać sumę ważoną. Wymiary są parametrem szablonu (a zatem stałą czasową kompilacji). Innym ograniczeniem jest to, że biblioteka powinna być...
Chcę opracować od podstaw równoległe oprogramowanie do obliczeń naukowych. Chcę się zastanowić, od którego języka zacząć. Program obejmuje odczytywanie / zapisywanie danych do plików txt i wykonywanie ciężkich obliczeń równolegle, z wieloma faktoryzacjami LU i wykorzystaniem rzadkich liniowych...
Język C ++ zapewnia ogólne programowanie i metaprogramowanie za pomocą szablonów. Techniki te znalazły zastosowanie w wielu dużych komputerowych pakietach naukowych (np. MPQC , LAMMPS , CGAL , Trilinos ). Ale co w rzeczywistości przyczynili się do obliczeń naukowych, które wykraczają poza...
Opracowuję kod do symulacji przepływu płynu z obecnymi w nim substancjami biologicznymi. Dotyczy to standardowych równań Naviera-Stokesa połączonych z niektórymi dodatkowymi modelami biologicznymi. Istnieje wiele parametrów / stałych. Napisałem funkcje do obsługi głównych obliczeń, ale mam problem...
Rzeczą, której najbardziej nie lubię w MPI, jest zajmowanie się typami danych (tj. Mapami / maskami danych), ponieważ nie pasują one tak dobrze do obiektowo zorientowanego C ++. boost::mpiobsługuje jednak tylko MPI 1.1 z ich strony internetowej: boost :: mpi to przyjazny interfejs C ++ do...
Wiele bardziej znanych bibliotek C ++ w informatyce, takich jak Eigen , Trilinos i deal.II, używa standardowego obiektu biblioteki szablonów C ++ std::complex<>do reprezentowania złożonych liczb zmiennoprzecinkowych. Gniazdo Poulson za odpowiedź na pytanie o domyślnych konstruktorów,...
Funkcja ma osobliwość zbliżoną do x = 0 . Osobliwość tę można jednak znieść: dla x = 1 należy mieć f ( x ) = 1 , ponieważ e x = ∑ k = 0 x kf:x↦(ex−1)/xf:x↦(ex−1)/xf \colon x \mapsto (e^x-1)/xx=0x=0x = 0x=1x=1x = 1f(x)=1f(x)=1f(x) = 1 a zatem (ex-1)/x=∑k=1x k - 1ex=∑k=0xkk!ex=∑k=0xkk! e^x =...
Rozwiązuję problem fizyczny za pomocą niejawnego schematu numerycznego. To prowadzi mnie do rozwiązania równania liniowego za pomocą macierzy tridiagonalnej. Algorytm kodowałem z Wikipedii. Zastanawiam się, czy istnieje wydajna biblioteka, która pozwala optymalnie rozwiązać tego rodzaju równanie....
Przenoszę istniejący kod z MATLAB do C ++ i mam liniowy system do rozwiązania (zamiast bardziej typowej formy A x = b )x A = bxZA=bxA=bA x = bZAx=bAx=b Matryca jest gęsta i ma ogólną formę, ale nie jest większa niż 1000 x 1000. Tak więc w MATLAB rozwiązaniem jest funkcja lub notacja ukośnika do...
Mam program Mathematica, który wykonuje niektóre całki w 3 lub 4 wymiarach przy użyciu tej QuasiMonteCarlometody. Problem polega na tym, że uruchomienie zajmuje denerwująco dużo czasu, do momentu, gdy niektóre z tych obliczeń nie mogą się zakończyć w maksymalnym czasie pracy dostępnym w naszym...
Artykuł „Przegląd szablonów ekspresji: Analiza wydajności obecnych metodologii” w SIAM Journal of Scientific Computing odwołuje się do biblioteki algebry liniowej „Blaze”. Nie słyszałem o tym wcześniej i nie mogę znaleźć referencji online. (Oczywiste wyszukiwania w Google zwracają powyższą wersję...
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 +=...