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) …
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 nazywa się problem znalezienia rozwiązania …
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ć …
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 \mathbf{R}T.⊂ R3)T⊂R3T \subset \mathbf{R}^3fafaffafaf Na przykład, jeśli byłaby zamiast …
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 (N).
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 multimodalne, ale analiza wydajności wydaje się zawsze wykorzystywać albo rozkłady deterministyczne, albo gaussowskie (nie tylko …
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 współczynnika Cholesky'ego .LM.MMLLL Załóżmy teraz, że otrzymałem prostokątną macierz , …
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 dobrego punktu wyjścia. Rozumiem, że nie ma obiektywnie najlepszego języka, ale chciałbym …
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 / ∂x∂u/∂x\partial{u}/\partial{x}∂u / ∂y∂u/∂y\partial{u}/\partial{y}ΓmΓm\Gamma_m …
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ą …
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 funkcjami: ma intuicyjny interfejs, który umożliwia tworzenie algorytmów obsługuje operacje dynamiczne: np. …
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 = bZAx=bAx=bZAZAAZAZAA
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 ++?
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 niż 10 minut. Jaki byłby najszybszy …
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 liniowej?
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.