def test_root_diagbroyden(self): res = root(F, F.xin, method='diagbroyden', options={ 'nit': 11, 'jac_options': { 'alpha': 1 } }) assert_(nonlin.norm(res.x) < 1e-8) assert_(nonlin.norm(res.fun) < 1e-8)
def test_root_linearmixing(self): res = root(F, F.xin, method='linearmixing', options={ 'nit': 60, 'jac_options': { 'alpha': 0.5 } }) assert_(nonlin.norm(res.x) < 1e-7) assert_(nonlin.norm(res.fun) < 1e-7)
def test_root_excitingmixing(self): res = root(F, F.xin, method='excitingmixing', options={ 'nit': 20, 'jac_options': { 'alpha': 0.5 } }) assert_(nonlin.norm(res.x) < 1e-5) assert_(nonlin.norm(res.fun) < 1e-5)
def test_root_broyden2(self): res = root(F, F.xin, method='broyden2', options={ 'nit': 12, 'jac_options': { 'alpha': 1 } }) assert_(nonlin.norm(res.x) < 1e-9) assert_(nonlin.norm(res.fun) < 1e-9)
def test_root_anderson(self): res = root(F, F.xin, method='anderson', options={ 'nit': 12, 'jac_options': { 'alpha': 0.03, 'M': 5 } }) assert_(nonlin.norm(res.x) < 0.33)
def test_diagbroyden(self): x = nonlin.diagbroyden(F, F.xin, iter=11, alpha=1) assert_(nonlin.norm(x) < 1e-8) assert_(nonlin.norm(F(x)) < 1e-8)
def test_exciting(self): x = nonlin.excitingmixing(F, F.xin, iter=20, alpha=0.5) assert_(nonlin.norm(x) < 1e-5) assert_(nonlin.norm(F(x)) < 1e-5)
def test_linearmixing(self): x = nonlin.linearmixing(F, F.xin, iter=60, alpha=0.5) assert_(nonlin.norm(x) < 1e-7) assert_(nonlin.norm(F(x)) < 1e-7)
def test_anderson(self): x = nonlin.anderson(F, F.xin, iter=12, alpha=0.03, M=5) assert_(nonlin.norm(x) < 0.33)
def test_broyden2(self): x = nonlin.broyden2(F, F.xin, iter=12, alpha=1) assert_(nonlin.norm(x) < 1e-9) assert_(nonlin.norm(F(x)) < 1e-9)