Snes solver : bounds

[Hasupama]

Hi,
I want to give bounds for the solver and changed snes type to 'vinewtonrsls'. But still it does not get the snes type as vinewtonrsls. The code is posted in below.

from firedrake import *
import time as comp_time
from firedrake.petsc import PETSc

solver_phi_parameters ={"snes_type": 'vinewtonrsls',"snes_linesearch_type": 'basic',
			   "pc_factor_mat_solver_type": "mumps",
                          "snes_rtol": 1.e-4,
                          "snes_atol": 1.e-3,  
                     	 "maximum_iterations": 100,
                    		"report": True,
                     	 	"error_on_nonconvergence": True,
                     	 									 }


mesh = Mesh('domain.msh')
x = SpatialCoordinate(mesh)
Space = FunctionSpace(mesh, 'CG', 2)

dphi = TrialFunction(Space)
del_phi = TestFunction(Space)

phi = Function(Space)
F_phi = x[0]*(  8.854187E-12*phi.dx(i)*del_phi.dx(i)  )*dx
bc_list_phi=[]
bc_list_phi.append( DirichletBC(Space, 0.0, 12) )
bc_list_phi.append( DirichletBC(Space, 20e3, 11) )

min_bound= Function(Space)
max_bound = Function(Space)
min_bound.interpolate(Constant(-1e10))
max_bound.interpolate(Constant(1e100))

problem_phi = NonlinearVariationalProblem(F_phi, phi, bcs=bc_list_phi,  form_compiler_parameters={"quadrature_degree":2})
solver_phi = NonlinearVariationalSolver(problem_phi)
solver_phi.parameters.update(solver_phi_parameters)
solver_phi.solve(bounds=(min_bound,max_bound))

The error says

[ 0] Object is in wrong state
[ 0] SNES solver newtonls does not support bounds

It seems like solver parameters are not updated and has the default snes type: newtonls
Appreciate comments on this issue.

[danshapero]

Does it work if you instead do
solver_phi = NonlinearVariationalSolver(problem_phi, solver_parameters=solver_phi_parameters)
instead of using the call to parameters.update later?

[Hasupama]

Yes. It worked. Thank you very much.