I have a nonlinear eigenvalue problem for a system of equations, of the form,
-Delta u + f(u) = E * u, where E is my nonlinear eigenvalue parameter, and u and f are vector valued. I thus have the following two things to contend with: 1. E is a scalar which needs to be distributed across all the processes when the right hand side is formed 2. I would like to be able to use a da to manage the spatial and multicomponent nature of u Obviously the Vec that my nonlinear solver is going to search for has to store both bits of data. Is there a clever petsc way to handle this, or will I need to do all the indexing and broadcasting by hand? -gideon
