Niejawne metody to metody z krokiem czasu, które używają inwersji w każdym kroku. Pozwala to na znacznie lepsze właściwości stabilności niż metody jawne, chociaż w niektórych przypadkach wiąże się to z poważną karą prędkości. Przykłady ukrytych metod obejmują Backward Euler i Crank-Nicholson.
Czy ktoś ma rekomendacje dotyczące użytecznej, szybkiej biblioteki macierzy C ++? Mam na myśli to, co użyteczne: Obiekty macierzy mają intuicyjny interfejs (np .: mogę używać wierszy i kolumn podczas indeksowania) Mogę zrobić wszystko z klasą macierzy, co mogę zrobić z LAPACK i BLAS Łatwy w...
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ć...
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
Metody numeryczne rozwiązywania PDE (lub ODE) dzielą się na dwie szerokie kategorie: metody jawne i niejawne. Metody niejawne pozwalają na większe stabilne kroki czasowe, ale wymagają więcej pracy na krok. W przypadku hiperbolicznych PDE powszechnie uważa się, że metody niejawne zwykle się nie...
Symulacja wielofizyczna obejmuje sprzężenie wielu „fizyków”, często z różnymi skalami przestrzeni i / lub czasu. Ponadto kody pojedynczej fizyki są często pisane przez różne zespoły. Najczęściej stosowaną techniką sprzęgania jest podział operatora pierwszego rzędu, ale ma on słabą dokładność i...
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...
Próbuję zrozumieć niektóre wyniki i doceniłbym kilka ogólnych uwag na temat rozwiązywania problemów nieliniowych. Równanie Fishera (nieliniowe PDE z dyfuzją reakcyjną), ut= dux x+ βu ( 1 - u ) = F.( u )ut=reuxx+βu(1-u)=fa(u) u_t = du_{xx} + \beta u (1 - u) = F(u) w formie dyskretnej, u′jot=...
Chcę zaimportować kod, który rozwiązuje zbiór równań różniczkowych cząstkowych (PDE) metodą objętości skończonej w formie IMPLICIT (dla dyskretyzacji czasu). W wyniku tego powstaje tridiagonalny układ równań w kierunkach x, y, z, który jest obsługiwany przez schemat ADI / TDMA. Nie mogę znaleźć...
Jaka jest dokładnie różnica między jawnym MES a niejawnym MES? Zgodnie z postem tutaj wydaje się, że jedyną różnicą jest to, czy stosowana jest integracja czasowa niejawna czy jawna. Jak pamiętam z jednej książki, którą przeczytałem, niejawna MES jest miejscem, w którym masa nie jest zbita w...