Computational Science newton-method

3

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 sqrt?

18 octave discretization nonlinear-equations newton-method visualization fluid-dynamics mesh-generation finite-element finite-volume optimization algorithms approximation fluid-dynamics navier-stokes comsol modeling optimization sparse-matrix matrix condition-number visualization matlab quadrature blas intel-mkl finite-element gpu discontinuous-galerkin mathematica optimization convex-optimization algorithms reference-request matlab statistics finite-element numerical-analysis petsc molecular-dynamics machine-learning statistics visualization open-source statistics image-processing visualization python petsc finite-element fluid-dynamics stability navier-stokes incompressible

2

Czy można rozwiązać nieliniowe PDE bez iteracji Newtona-Raphsona?

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= L …

15 numerical-analysis nonlinear-equations implicit-methods newton-method explicit-methods

1

Czy przybliżony jakobian ze skończonymi różnicami może powodować niestabilność w metodzie Newtona?

Zaimplementowałem solver z Eulerem wstecznym w Pythonie 3 (używając numpy). Dla własnej wygody i jako ćwiczenie napisałem również małą funkcję, która oblicza przybliżoną różnicę skończoną różnicy gradientu, aby nie zawsze musiałem określać analitycznie jakobian (jeśli to w ogóle możliwe!). Korzystając z opisów podanych w Ascher i Petzold 1998 , napisałem …

13 pde stability numpy scipy newton-method

2

Strategie metody Newtona, gdy jakobian przy rozwiązaniu jest osobliwy

Próbuję rozwiązać następujący układ równań dla zmiennych i (wszystkie pozostałe są stałymi):x 2P., x1P.,x1P,x_1x2)x2)x_2 A ( 1 - P)2)- k1x1= 0A P.2)- k2)x2)= 0( 1 - P) ( r1+ x1)4L.1- P( r1+ x2))4L.2)= 0ZA(1-P.)2)-k1x1=0ZAP.2)-k2)x2)=0(1-P.)(r1+x1)4L.1-P.(r1+x2))4L.2)=0\frac{A(1-P)}{2}-k_1x_1=0 \\ \frac{AP}{2}-k_2x_2=0 \\ \frac{(1-P)(r_1+x_1)^4}{L_1}-\frac{P(r_1+x_2)^4}{L_2}=0 Widzę, że mogę przekształcić ten układ równań w pojedyncze równanie jednej zmiennej …

12 optimization convergence newton-method

3

Metody rozwiązywania nieliniowych układów doradczo-dyfuzyjnych poza Newton-Raphson?

Pracuję nad projektem, w którym mam dwie domeny sprzężone za pomocą adv-diff za pomocą odpowiednich terminów źródłowych (jedna domena dodaje masę, druga odejmuje masę). Dla zwięzłości modeluję je w stanie ustalonym. Równania to standardowe równanie transportowo-dyfuzyjne z terminem źródłowym wyglądającym tak: ∂c1∂t=0=F1+Q1(c1,c2)∂c2∂t=0=F2+Q2(c1,c2)∂c1∂t=0=F1+Q1(c1,c2)∂c2∂t=0=F2+Q2(c1,c2) \frac{\partial c_1}{\partial t} = 0 = \mathcal{F}_1 + …

9 nonlinear-equations newton-method advection-diffusion coupling

Pytania otagowane jako newton-method