problem = Problem.from_conf(conf, init_equations=False)
problem.setup_default_output()
problem.set_equations({'eq' : conf.equations['balance']})
problem.time_update(ebcs=conf.ebcs)
state_el = problem.solve()
time_solver = problem.get_time_solver()
time_solver(save_results=True,post_process_hook=stress_strain)
ts = problem.get_timestepper()
Ans these were my solvers:
# Solvers etc.
solver_0 = {
'name' : 'ls',
'kind' : 'ls.scipy_direct',
}
solver_1 = {
'name' : 'newton',
'kind' : 'nls.newton',
'i_max' : 7,
'eps_a' : 1e-10,
'eps_r' : 1.0,
'macheps' : 1e-16,
'lin_red' : 1e-2, # Linear system error < (eps_a * lin_red).
'ls_red' : 0.1,
'ls_red_warp': 0.001,
'ls_on' : 1.1,
'ls_min' : 1e-5,
'check' : 0,
'delta' : 1e-6,
'problem' : 'nonlinear', # 'nonlinear' or 'linear' (ignore i_max)
}
solver_2 = {
'name' : 'ts',
'kind' : 'ts.simple',
't0' : 0,
't1' : 1,
'dt' : None,
'n_step' : 101, # has precedence over dt!
}