Pytania oznaczone «petsc»

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ć...

Co może pójść nie tak, gdy zastosuje się wstępne metody Kryłowa z KSP ( pakiet solvera liniowego PETSc ) do rozwiązania rzadkiego układu liniowego, takiego jak te uzyskane przez dyskretyzację i linearyzację równań różniczkowych cząstkowych? Jakie kroki mogę podjąć, aby ustalić, co jest nie tak z...

O ile mogę stwierdzić, dwie duże ogólne ramy oprogramowania do obliczeń informatycznych Departamentu Energii w USA to PETSc i Trilinos . Na pierwszy rzut oka wydają się podobne, poza różnicami językowymi (C kontra C ++). Jakie są główne różnice między tymi dwoma ramami i jakie czynniki powinny...

Używam nieliniowego pakietu solvera PETSc SNES do rozwiązania układu równań nieliniowych uzyskanych przez dyskretyzację równania różniczkowego cząstkowego. Jak mogę ustalić, dlaczego solver nie jest zbieżny i co mogę zrobić, aby pomyślnie rozwiązać moje

Chciałbym wiedzieć, czy istnieje szybki sposób na obliczenie odległości euklidesowej dwóch wektorów w oktawie. Wydaje się, że nie ma do tego żadnej specjalnej funkcji, więc czy powinienem po prostu użyć formuły z

Planuję użyć multigrid do obliczenia wartości własnych i wektorów, i zauważyłem, że PETSc ma wsparcie wysokiego poziomu dla multigrid. Dokumentacja PETSc mówi, że ta część PETSc nie powinna być używana, ponieważ wkrótce zostanie wymieniona. Które inne biblioteki mają obsługę wysokiego poziomu dla...

Mam PETSc Mati chciałbym oszacować jego numer

Bawiłem się PETSc i zauważyłem, że kiedy uruchamiam mój program z więcej niż jednym procesem przez MPI , wydaje się, że działa jeszcze wolniej ! Jak mogę sprawdzić, co się

Gdziekolwiek widziałem, tutorial / dokumenty PETSc itp. Mówią, że jest to przydatne w algebrze liniowej i zwykle określa, że ​​korzyści będą miały systemy rzadkie. Co z gęstymi matrycami? Jestem zaniepokojony rozwiązywania dla gęstej A .A x = bAx=bAx=bZAZAA W Fortranie napisałem własny kod dla CG...

Zwykle piszę kod seryjny, a kiedy to robię, piszę testy jednostkowe za pomocą szkieletu testowego w stylu xUnit (MATLAB xUnit, PyUnit / nose lub szkielet testowy Google C ++). Na podstawie pobieżnej wyszukiwarki Google nie widziałem wiele na temat tego, jak praktykujący testują kod jednostkowy,...

Pytanie: Jakie metody są dostępne w celu dokładnego i wydajnego obliczenia struktury rzadkości matrycy elementów skończonych? Informacje: Pracuję nad solwerem Poissona Równania Ciśnienia, stosując metodę Galerkina z kwadratową podstawą Lagrange'a, napisaną w C, i używam PETSc do rzadkiego...

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...

Muszę rozwiązać ten sam rzadki układ liniowy (300 x 300 do 1000 x 1000) z wieloma prawymi bokami (300 do 1000). Oprócz tego pierwszego problemu chciałbym również rozwiązać różne systemy, ale z tymi samymi niezerowymi elementami (tylko różne wartości), to znaczy wiele rzadkich systemów o stałym...

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...

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

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

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...

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,...

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ć...

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...