I have to choose between the PETSc and Trilinos libraries for a project that involves the solution of non-linear systems of equations.  From their web sites they both mention support for GPUs, without giving many details (e.g. here and here).  I would like to know which of these libraries is better for GPGPU. In particular, I am looking for a library that, on the basis of the same calling source code, can use either the CPU or the GPU on the basis of a pre-compiler switch or, better, dynamically, on the basis of the time taken to run the computations and to move the data in and out the graphics memory.
PETSc or Trilinos for GPU?
petsc;nonlinear equations;gpu;trilinos