Na przykład nVidia ma CUBLAS, który obiecuje 7-14-krotne przyspieszenie. Naiwnie nie jest to teoretyczna przepustowość żadnej z kart GPU nVidii. Jakie są wyzwania związane z przyspieszaniem algebry liniowej na GPU i czy są już dostępne szybsze trasy algebry
Chciałbym wykonać PCA na zestawie danych złożonym z około 40 000 próbek, z których każda zawiera około 10 000 funkcji. Używanie funkcji princomp Matlaba konsekwentnie zajmuje ponad pół godziny, w którym to momencie zabijam proces. Chciałbym znaleźć implementację / algorytm, który działa w mniej...
CVXOPT: http://abel.ee.ucla.edu/cvxopt/index.html OpenOpt: http://openopt.org/Welcome Jaki jest związek między nimi? Jakie są ich zalety / wady? BTW, czy jest jakaś inna wysokiej jakości biblioteka optymalizacji wypukłej ogólnego przeznaczenia dla Python / C
Czy istnieje opublikowana lista najlepszych praktyk zapewniających długowieczność kodu, z myślą o odtwarzalnych wynikach naukowych? (np. open source, praktyki dokumentacji, wybór zależności, wybór języka, maszyny wirtualne itp.). Wiedz o wszelkich badaniach (lub ich brak, przykładach /...
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...
Witaj społeczności scicomp, Pracowałem w obszarze algorytmów graficznych z wykorzystaniem frameworków takich jak NetworkX (Python), JUNG i YFiles (Java). Wchodzę teraz w obszar obliczeń równoległych i wysokowydajnych. W przypadku nowego projektu szukam biblioteki grafów C ++ z następującymi...
Twierdzenie o rozkładzie Cholesky'ego stwierdza zatem, że każda prawdziwa symetryczna macierz dodatnio określona ma rozkład Choleskiego gdzie jest dolną macierzą trójkątną.M = L L ⊤ L.MMMM=LL⊤M=LL⊤M= LL^\topL.LL Biorąc pod uwagę , wiemy już, że istnieją szybkie algorytmy do obliczania jego...
Rozwiązuję dla ogromnej rzadkiej dodatniej określonej macierzy za pomocą metody gradientu sprzężonego (CG). Czy można obliczyć wyznacznik podstawie informacji uzyskanych podczas rozwiązania?A AA x =
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...
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. Jestem zupełnie nowy w nauce obliczeniowej i szukam...
Powiedzmy, że mam funkcję , którą chcę zintegrować z czworościanem . Gdyby było arbitralne, kwadratura Gaussa byłaby dobrym rozwiązaniem, ale zdaję sobie sprawę, że jest harmoniczne. Jak bardzo można przyspieszyć kwadraturę Gaussa przy użyciu tych informacji?fa: R3)→ Rf:R3→Rf : \mathbf{R}^3 \to...
Teoretycznie, w jaki sposób czas na wykonanie obliczeń teorii funkcjonalnej gęstości (DFT) wraz z liczbą elektronów? Interesują mnie „typowe” implementacje DFT, takie jak VASP, ABINIT itp., A nie kody O
Mam niejednorodny układ liniowy A x = bAx=b Ax=b gdzie jest rzeczywistą macierzą n × n przy n ≤ 4 . Nullspace z A jest gwarancją wymiaru zerowym więc równanie ma unikalną odwrotną x = A - 1 b . Ponieważ wynik wchodzi w prawą stronę ODE, którą zamierzam rozwiązać za pomocą metody adaptacyjnej,...
Biorąc pod uwagę nieznaną funkcję , możemy ocenić jej wartość w dowolnym punkcie w jej dziedzinie, ale nie mamy jej wyrażenia. Innymi słowy, f jest dla nas jak czarna skrzynka.f:Rd→Rf:Rd→Rf:\mathbb R^d \to \mathbb Rfff Jak nazywa się problem znalezienia minimalizatora ? Jakie są metody?fff Jak...
Jak na przykład MATLAB oblicza SVD danej macierzy? Zakładam, że odpowiedź prawdopodobnie obejmuje obliczenie wektorów własnych i wartości własnych A*A'. Jeśli tak jest, chciałbym również wiedzieć, jak to
Mam problem fizyczny rządzony równaniem Poissona w dwóch wymiarach Mam pomiary dwóch składników gradientu ∂ u / ∂ x i ∂ u / ∂ y wzdłuż pewnej części granicy, Γ m , więc chciałbym narzucić ∂ u−∇2u=f(x,y),inΩ−∇2u=fa(x,y),janΩ -\nabla^2 u = f(x,y), \; in \; \Omega ∂u /...
W pracy zasadniczo działam jako niezależny konsultant. W przypadku kierownictwa i klientów muszę oszacować czas potrzebny na opracowanie oprogramowania w ramach moich badań w dziedzinie obliczeń. Jednak moje prognozy czasu są zwykle wyłączone. Wiem, że istnieją metody szacowania czasu potrzebnego...
Mam zestaw danych 3-wymiarowych punktów, dla których chciałbym zbudować siatkę, używając Pythona. Całe oprogramowanie, które widziałem, wymaga podania krawędzi. Czy w Pythonie jest program, który jako dane wejściowe przyjmuje zestaw punktów w 3D i generuje siatkę? Jeśli to możliwe, chciałbym, aby...
Począwszy od równania doradczego w formie konserwatorskiej. ut=(a(x)u)xut=(a(x)u)x u_t = (a(x)u)_x gdzie a(x)a(x)a(x) jest prędkością zależną od przestrzeni, uuu jest stężeniem gatunku, który jest zachowany. Dyskretyzacja strumienia (gdzie strumień f=a(x)uf=a(x)uf=a(x)u jest zdefiniowany na...
Obecnie opracowuję metodę dekompozycji domen dla rozwiązania problemu rozpraszania. Zasadniczo rozwiązuję iteracyjnie system BVP firmy Helmholtz. Dyskretyzuję równania metodą elementów skończonych na siatkach trójkątnych lub czworościennych. Rozwijam kod w kierunku mojej pracy doktorskiej. Zdaję...