Hi all, I am following example 13 and trying to solve a set of nonlinear PDEs using the solver of libmesh. I have noted that: >From previous Newton iterate, the variable values of "u" and "v" can be assigned by u += phi[l][qp]*navier_stokes_system.current_solution (dof_indices_u[l]); v += phi[l][qp]*navier_stokes_system.current_solution (dof_indices_v[l]); For the first nonlinear loop, "current_solution (dof_indices_u[l])" is zero, so all the initial values of "u" and "v" are assigned 0. But I would like to give the variables non-zero initial values, could you guys please tell me how should I do that? what function should be used? Thanks in advance.
Best wishes, yunfei ------------------------------------------------------------------------- This SF.Net email is sponsored by the Moblin Your Move Developer's challenge Build the coolest Linux based applications with Moblin SDK & win great prizes Grand prize is a trip for two to an Open Source event anywhere in the world http://moblin-contest.org/redirect.php?banner_id=100&url=/ _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
