Metoda przedstawiania liczb za pomocą ustalonej liczby cyfr znaczących i wykładnika pewnej liczby podstawowej. Są scharakteryzowane w postaci
. Zazwyczaj liczby są przedstawiane w odniesieniu do podstawy = 2 (binarne).
( s i gn i fja c a n t dja gi t s ) ∗basmie x p o n e n t
Mam kilka trudnych, niewypukłych problemów globalnej optymalizacji do rozwiązania. Obecnie używam MATLAB's Optimization Toolbox (konkretnie fmincon()z algorytmem = 'sqp'), co jest dość skuteczne . Jednak większość mojego kodu znajduje się w języku Python i chciałbym również przeprowadzić optymalizację w języku Python. Czy istnieje solver NLP z powiązaniami Pythona, z którym …
Mam proste pytanie, które jest naprawdę trudne dla Google (oprócz kanonicznego Co każdy informatyk powinien wiedzieć o arytmetyki zmiennoprzecinkowej ). Kiedy należy używać funkcji takich jak log1plub expm1należy używać zamiast logi exp? Kiedy nie należy ich używać? Czym różnią się różne implementacje tych funkcji pod względem ich wykorzystania?
Biorąc pod uwagę rzeczywistą funkcję rzeczywistych zmiennych, czy jest dostępne oprogramowanie, które może automatycznie generować dokładny numerycznie kod do obliczenia funkcji dla wszystkich danych wejściowych na maszynie wyposażonej w arytmetykę IEEE 754? Na przykład, jeśli rzeczywistą funkcją do oceny były: Oprogramowanie rozważyłoby katastrofalne anulowanie i możliwe wyszukiwanie tablic wyników dla …
Bardzo częstym problemem w Markov Chain Monte Carlo jest obliczanie prawdopodobieństw, które są sumą dużych terminów wykładniczych, ea1+ea2+...ea1+ea2+... e^{a_1} + e^{a_2} + ... gdzie składniki może wahać się od bardzo małych do bardzo dużych. Moje podejście na uwzględnieniu największego wykładniczego terminu , aby:aaaK:=maxi(ai)K:=maxi(ai)K := \max_{i}(a_{i}) e ' ≡ e się …
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ę z MinGW / GCC. Programuję w C ++. Jaki jest typowy sposób programowego …
Próbuję zaimplementować następującą funkcję w zmiennoprzecinkowym podwójnej precyzji z niskim błędem względnym : logsum(x,y)=log(exp(x)+exp(y))logsum(x,y)=log(exp(x)+exp(y))\mathrm{logsum}(x,y) = \log(\exp(x) + \exp(y)) Jest to szeroko stosowane w aplikacjach statystycznych w celu dodania prawdopodobieństw lub gęstości prawdopodobieństwa, które są reprezentowane w przestrzeni dziennika. Oczywiście albo exp(x)exp(x)\exp(x) albo exp(y)exp(y)\exp(y) mogą łatwo przepełnić lub niedopełnić, co byłoby …
Biorąc pod uwagę dwie różne implementacje BLAS, czy możemy oczekiwać, że wykonają dokładnie takie same obliczenia zmiennoprzecinkowe i zwrócą te same wyniki? Lub może się zdarzyć, na przykład, że oblicza się produkt skalarny jako i jeden jako co może dać inny wynik w zmiennoprzecinkowym IEEE arytmetyka?( ( x1y1+ x2)y2)) + …
Uruchomiłem kod dynamiki molekularnej (MD) GROMACS w klastrze Ubuntu Linux składającym się z węzłów zawierających 24 procesory Intel Xeon. Moje szczególne zainteresowanie okazuje się nieco wrażliwe na zmiennoprzecinkową precyzję arytmetyczną, więc musiałem uruchomić GROMACS z podwójną precyzją, a nie z pojedynczą precyzją - pomimo wyższych kosztów obliczeniowych podwójnej precyzji. Tak …
Chciałbym określić teoretyczną liczbę operacji FLOP (operacji zmiennoprzecinkowych), jaką może wykonać mój komputer. Czy ktoś mógłby mi z tym pomóc. (Chciałbym porównać mój komputer z niektórymi superkomputerami, aby zorientować się, jaka jest między nimi różnica)
Pozwól x, ybyć dwóch liczb zmiennoprzecinkowych. Jaki jest właściwy sposób na obliczenie ich średniej? Naiwny sposób (x+y)/2może doprowadzić do przepełnienia kiedy xi ysą zbyt duże. Myślę, że 0.5 * x + 0.5 * ymoże lepiej, ale wiąże się to z dwoma multiplikacjami (co może być nieefektywne) i nie jestem pewien, …
Mam problem z implementacją funkcji numerycznie. Cierpi na tym fakt, że przy dużych wartościach wejściowych wynik jest bardzo dużą liczbą razy bardzo małą liczbą. Nie jestem pewien, czy katastrofalne anulowanie jest właściwym terminem, więc proszę mnie poprawić, jeśli tak jest. Dowody na coś nie tak: Jak mogę uniknąć oscylacji i …
Czytałem książkę Computer Simulation of Liquids autorstwa Allena i Tildesleya. Począwszy od strony 71, autorzy omawiają różne algorytmy stosowane do integracji równań ruchu Newtona w symulacjach dynamiki molekularnej (MD). Począwszy od strony 78, autorzy omawiają algorytm Verleta, który jest być może algorytmem integracji kanonicznej w MD. Stanowią one: Być może …
Oceniając liczbę FLOP w prostej funkcji, często można po prostu zejść w dół wyrażenia zestawiając podstawowe operatory arytmetyczne. Jednak w przypadku wyrażeń matematycznych obejmujących parzysty podział nie można tego zrobić i można oczekiwać, że będzie można porównać z liczbą FLOP z funkcji z tylko dodatkami i mnożeniami. Sytuacja jest jeszcze …
Większość bibliotek matematycznych ma wiele wersji funkcji logarytmicznych. Przez większość czasu zakładamy, że są idealne, ale w rzeczywistości całkiem sporo z nich oferuje określoną liczbę cyfr precyzji. W przypadku niektórych funkcji istnieją stabilniejsze numerycznie warianty. Na przykład Fortran, R, Java i C mają Math.log1p, zarówno do obliczeń log(1.0+x)(która oferuje wyższą …
Używamy plików cookie i innych technologii śledzenia w celu poprawy komfortu przeglądania naszej witryny, aby wyświetlać spersonalizowane treści i ukierunkowane reklamy, analizować ruch w naszej witrynie, i zrozumieć, skąd pochodzą nasi goście.
Kontynuując, wyrażasz zgodę na korzystanie z plików cookie i innych technologii śledzenia oraz potwierdzasz, że masz co najmniej 16 lat lub zgodę rodzica lub opiekuna.