High Performance Direct Gravitational N-body Simulations on Graphics Processing Unit I: An implementation in Cg