Mój system to symetryczny problem FE z mnożnikami lagrange (np. Nieściśliwy przepływ Stokesa): (ABBTC)(ABTBC)\begin{pmatrix}A & B^T \\ B & C\end{pmatrix} gdzie jest typowym przypadkiem (upewniłem się nawet, że równania są ponumerowane, aby mnożniki Lagrange'a pojawiały się na końcu). System jest...
Oto typowe użycie wskaźników funkcji w C. Chciałbym zrobić coś podobnego w Fortranie. Mam kilka pomysłów, ale chciałbym wiedzieć, czy istnieje jakiś kanoniczny sposób. Wskaźniki funkcji i konteksty przekazane przez użytkownika są przechowywane, a następnie wywoływane później. typedef...
Obliczenia równoległe są często modelowane przy użyciu deterministycznej lokalnej szybkości obliczeń, obciążenia związanego z opóźnieniami i przepustowości sieci. W rzeczywistości są one zmienne przestrzennie i niedeterministyczne. Badania takie jak Skinner i Kramer (2005) obserwują rozkłady...
O ile wolniej jest Petsc4py vs C / C ++ / Fortran? Zdaję sobie sprawę, że będzie to w znacznym stopniu zależeć od wykonywanego kodu, ale co z czymś prostym, takim jak produkt
Dano mi macierz , która jest symetryczna, odwracalna, dodatnia określona i gęsta. Muszę sprawdzić, czy gdzie J jest macierzą wszystkich jedynek.Q det ( Q ) = det ( 12 I - Q - J )12×1212×1212 \times 12QQQJdet ( Q ) = det ( 12 I- Q - J)( 1 )det(Q)=det(12I−Q−J)(1)\det(Q) = \det(12I-Q-J) \; \;...
Czy kompilacja PETSc z zewnętrzną biblioteką BLAS / LAPACK znacząco wpływa na wydajność na rzadkich macierzach, czy używa tych bibliotek tylko do matematyki gęstej
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...
Który jest obliczany szybciej, lub lub ? , i są liczbami dodatnimi z .zaloguj a c b √ababa^blogaclogac\log_a c abcb>1c√bcb\sqrt[b]{c}aaabbbcccb>1b>1b>1 Jakiego rodzaju algorytmów użyjesz w porównaniu? Jakie są ich zawiłości? Na przykład, gdy lub c ≈ a bc≡abc≡abc \equiv a^bc≈abc≈abc...
Istnieje wiele dobrze znanych metod numerycznych rozwiązywania równań typu np. metoda bisekcji, metoda Newtona itp.f(x)=0,x∈Rn,f(x)=0,x∈Rn, f(x) = 0, \quad x \in \mathbb{R}^n, W mojej aplikacji oblicza się metodą stochastyczną (wynik jest średnią).f(x)f(x)f(x) Czy są jakieś metody rozwiązywania...
Szukam wiarygodnych referencji określających, ile zasobów superkomputery wydają na koordynację, a nie na faktyczną pracę związaną z zadaniami. Zasoby mogą być dostępną mocą obliczeniową, ale nawet waty wydają się być prawidłową jednostką. Wierzę, że jeden z moich profesorów lub podręczników...
Jestem nowy w pakiecie PETSc. Mam macierz A ~ 4000x4000 w formacie matrycowym i chcę, aby PETSc rozwiązał ten problem za pomocą wielu procesorów. Wiem, jak rozwiązać system na jednym procesorze, ale nie wiem, jak rozdzielić macierz i wektory między różne procesory. Czy istnieje prosty zestaw...
Mam spory model (~ 5000 linii) napisany w C. Jest to program szeregowy, bez generowania liczb losowych w dowolnym miejscu. Wykorzystuje bibliotekę FFTW dla funkcji korzystających z FFT - nie znam szczegółów implementacji FFTW, ale zakładam, że funkcje w niej również są deterministyczne (poprawcie...
W ostatnich latach pojawiło się kilka projektów bibliotek / oprogramowania, które oferują jakąś formę równoległego zarządzania pamięcią współużytkowaną opartą na danych ogólnych. Główną ideą jest to, że zamiast pisać kod z jawnym wątkiem, programiści implementują swoje algorytmy jako zadania...
Jestem bardzo nowy w programowaniu GPGPU, więc proszę wybacz mi, jeśli pytanie nie jest szczególnie odpowiednie. Z tego, co rozumiem, programowanie na GPU to bardzo skomplikowana praca inżynierska w porównaniu do zwykłego programowania procesora. Należy bardzo uważać na problemy z rozbieżnościami,...
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...
Zanurzam się w fascynujący świat analizy elementów skończonych i chciałbym rozwiązać duży problem termomechaniczny (tylko mechaniczny efekt termiczny , brak sprzężenia zwrotnego).→→\rightarrow W przypadku problemu mechanicznego zrozumiałem już odpowiedź Geoffa , że muszę użyć iteracyjnego solvera...
Muszę rozwiązać Ax = b, ale zdaję sobie sprawę, że nawet jeśli jest on rzadki, przechowywanie współczynników macierzy mojego problemu zajmie zbyt dużo pamięci. Zastanawiam się teraz nad użyciem metody bez macierzy, ponieważ te same współczynniki pojawiają się w macierzy dużo czasu, więc mogłem użyć...
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...
Jak zilustrowano w odpowiedzi Jeda Browna na Koszty wyszukiwania w porównaniu z obliczeniami , użycie wektoryzowanych vs niewektoryzowanych operacji zmiennoprzecinkowych skutkuje znacznie szybszym kodem. Wiele współczesnych kompilatorów twierdzi, że może wykonywać automatyczną wektoryzację. Jak...