Pytanie: Załóżmy, że masz dwa różne (uwzględnione) warunki wstępne dla symetrycznej dodatniej macierzy : i gdzie odwrotność czynników wynosi łatwy do zastosowania.ZAZAAA ≈ BT.bZA≈bT.bA \approx B^TBA ≈ CT.do,ZA≈doT.do,A \approx C^TC,B , BT., C., C.T.b,bT.,do,doT.B, B^T, C, C^T Gdy jest to...
W szczególności chcę przedłużyć pracę z B3LYP rozpoczętą od Gaussian 03, ale kontynuowałem w GAMESS-US. Energie dostarczane przez domyślne metody B3LYP nie są takie same. Dyskusja na ten temat znajduje się w podręczniku GAMESS-US (sekcja Dalsze informacje): Należy zauważyć, że B3LYP w GAMESS...
Mam PETSc Mati chciałbym oszacować jego numer
Istnieje wiele schematów FD dla równania porady dyskutować w sieci. Na przykład tutaj: http://farside.ph.utexas.edu/teaching/329/lectures/node89.html∂T∂t+u∂T∂x=0∂T∂t+u∂T∂x=0\frac{\partial T}{\partial t}+u\frac{\partial T}{\partial x}=0 Ale nie widać każdy zaproponowanie "ukryte" wiatr schemat...
Ludzie informatyki: I pierwotnie opublikowany na to pytanie w matematyce Stos Exchange i ktoś powiedział, że mogę dostać się „znacznie lepiej” odpowiedź tutaj: Jestem nowicjuszem w metodach numerycznych i Matlabie. Próbuję oszacować następującą sumę dwóch potrójnych całek (można to oczywiście...
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
Jakie podejścia są stosowane w praktyce do szacowania liczby warunków dużych rzadkich
Częstym problemem w statystyce jest obliczanie pierwiastka kwadratowego odwrotnego symetrycznej dodatniej macierzy określonej. Jaki byłby najbardziej efektywny sposób obliczenia tego? Natknąłem pewnym literaturze (które nie zostały jeszcze przeczytane) oraz jakiegoś przypadkowego kodu R tutaj ,...
Napisałem program / bibliotekę, której użyłem do uzyskania wyników w artykule. (Tutaj jest , ale moje pytanie jest ogólne). Mam testy, które regularnie uruchamiam ctest( uruchomienie zajmuje kilka minut). Aby odtworzyć niektóre tabele lub ryciny w tym artykule, muszę zbudować skrypt lub prosty...
Interesuje mnie rozwiązanie równania Poissona przy użyciu metody różnic skończonych. Chciałbym lepiej zrozumieć, jak zapisać równanie macierzowe z warunkami brzegowymi Neumanna. Czy ktoś sprawdziłby następujące, czy to prawda? Macierz różnic skończonych Równanie Poissona, ∂2)u ( x )∂x2)= d( x...
Muszę obliczyć następującą całkę: Gdziejest macierzą (kinetyka jednej cząstki i energia potencjalna wyrażona w podstawie),jest macierzą zależną od(funkcja Greena jednego ciała dla wielu ciał), a całka konturu jest lewym półkolem. Całkama bieguny na ujemnej osi rzeczywistej i ocena jest kosztowna....
Mam zestaw danych x1,x2,…,xkx1,x2),…,xkx_{1}, x_{2}, \ldots, x_{k} i chcę znaleźć parametr mmm taki, aby minimalizował sumę ∑i=1k∣∣m−xi∣∣.∑ja=1k|m-xja|.\sum_{i=1}^{k}\big|m-x_i\big|. to
O ile mi wiadomo, solwery wielosieciowe używają iteracyjnych wygładzaczy, takich jak Jacobi, Gauss-Seidel i SOR, aby tłumić błąd przy różnych częstotliwościach. Czy zamiast tego można zastosować metodę podprzestrzeni Kryłowa (jak gradient sprzężony, GMRES itp.)? Nie sądzę, że są one klasyfikowane...
Próbuję rozwiązać równanie Poissona 2D na podstawie różnic skończonych. W tym procesie otrzymuję rzadką macierz z tylko zmiennymi w każdym równaniu. Na przykład, jeśli zmienne byłyby , dyskretyzacja dałaby:555UUU Ui−1,j+Ui+1,j−4Ui,j+Ui,j−1+Ui,j+1=fi,jUi−1,j+Ui+1,j−4Ui,j+Ui,j−1+Ui,j+1=fi,jU_{i-1,j}...
Słyszałem anegdotycznie, że gdy ktoś próbuje liczbowo wykonać całkę formy ∫∞0f( x ) J0(x)dx∫0∞f(x)J0(x)dx\int_0^\infty f(x) J_0(x)\,\mathrm{d}x z gładkim i dobrze się zachowującym (np. nie jest silnie oscylacyjny, niesingularny itp.), wówczas poprawność przepisania go jakofa( x...
Istnieje system więzów liniowych . Chciałbym znaleźć ściśle dodatni wektor x > 0, który spełnia te ograniczenia. Oznacza to, że x i > 0 jest wymagane dla każdego składnika x i z x . Jak mogę użyć solwera LP, aby znaleźć tak ściśle dodatni wektor x (lub potwierdzić, że nie istnieje x )? Nie...
Przepraszam za długi post, ale chciałem załączyć wszystko, co uważałem za istotne za pierwszym razem. Czego chcę Wdrażam równoległą wersję Krystalicznych metod podprzestrzeni dla gęstych matryc. Głównie GMRES, QMR i CG. Zdałem sobie sprawę (po profilowaniu), że moja procedura DGEMV była żałosna....
Potrzebuję prostego wyjaśnienia metody wielosiatkowej lub literatury na ten temat. Znam metody iteracyjne, w tym BiCGStab, CG, GS, Jacobi i kondycjonowanie wstępne, ale jestem początkującym w metodzie wielosiatkowej. Czy ktoś może wyjaśnić to szczegółowo lub przynajmniej podać wyraźnie pseudokod...
W symulacji półprzewodników często równania są skalowane, aby miały znormalizowane wartości. Na przykład w skrajnych przypadkach gęstość elektronów w półprzewodnikach może zmieniać się o ponad 18 rzędów wielkości, a pole elektryczne może zmieniać się kształtnie, o ponad 6 (lub więcej) rzędów...
Moja grupa badawcza koncentruje się na dynamice molekularnej, która oczywiście może generować gigabajty danych w ramach jednej trajektorii, którą należy następnie przeanalizować. Kilka problemów, którymi się zajmujemy, wiąże się z korelacjami w zbiorze danych, co oznacza, że musimy śledzić duże...