interpolate_boundary_values (dof_handler,0, ZeroFunction<2>(), boundary_values);
12
ZeroFunction<2> tmp; interpolate_boundary_values (dof_handler,0, tmp, boundary_values);