Mam układ równań nieliniowych, które chcę rozwiązać numerycznie:nnn f = ( f 1 , … , f n )f(x)=af(x)=a\mathbf{f}(\mathbf{x})=\mathbf{a} f=(f1,…,fn)x=(x1,…,xn)f=(f1,…,fn)x=(x1,…,xn)\mathbf{f}=(f_1,\dots,f_n)\quad\mathbf{x}=(x_1,\dots,x_n) Ten system ma wiele cech, które sprawiają, że jest...
Korzystałem z SVD Intela MKL ( dgesvdprzez SciPy) i zauważyłem, że wyniki są znacząco różne, kiedy zmieniam precyzję pomiędzy float32i float64kiedy moja matryca jest źle uwarunkowana / nie w pełni ustawiona. Czy istnieje przewodnik dotyczący minimalnej ilości regularyzacji, którą powinienem dodać,...
Napisałem kodu jednocześnie temu, która polegała obliczyć bez korzystania z funkcji biblioteki. Wczoraj sprawdzałem stary kod i starałem się, aby był jak najszybciej (i poprawiony). Oto moja dotychczasowa próba:log(x)log(x)log(x) const double ee = exp(1); double series_ln_taylor(double n){ /* n =...
Nie mogę sobie wyobrazić, że pierwszy pomyślę o następującym problemie, więc będę zadowolony z referencji (ale zawsze doceniamy pełną, szczegółową odpowiedź): Załóżmy, że masz symetryczny dodatni określony . jest uważane za bardzo duże, więc trzymanie w pamięci jest niemożliwe. Możesz jednak...
Czytam artykuł [1], w którym rozwiązują następujące równanie nieliniowe przy użyciu metod różnic skończonych. Analizują także stabilność schematów za pomocą analizy stabilności von Neumanna. Jednak, jak zdają sobie sprawę autorzy, dotyczy to tylko liniowych PDE. Autorzy więc to poprzez...
Interesuje mnie obliczenie rozwiązania systemu lage ODE przy użyciu metody krylova jak w [1]. Taka metoda obejmuje funkcje związane z wykładniczym (tzwφφ\varphi-Funkcje). Zasadniczo polega na obliczeniu działania funkcji macierzowej przez zbudowanie podprzestrzeni Kryłowa za pomocą iteracji...
Przeprowadzam wyszukiwanie linii w ramach quasi-Newtona algorytmu BFGS. W jednym kroku wyszukiwania linii używam interpolacji sześciennej, aby zbliżyć się do lokalnego minimalizatora. Niech będzie funkcją będącą przedmiotem zainteresowania. Chcę znaleźć takie, że .f:R→R,f∈C1f:R→R,f∈C1f : R...
Jakie są główne różnice między FEM a XFEM? Kiedy powinniśmy (nie) korzystać z XFEM zamiast FEM i odwrotnie? Innymi słowy, kiedy napotykam nowy problem, skąd mogę wiedzieć, z którego z nich