  Geometric phases have been used in NMR, to implement controlled phase shift
gates for quantum information processing, only in weakly coupled systems in
which the individual spins can be identified as qubits. In this work, we
implement controlled phase shift gates in strongly coupled systems, by using
non-adiabatic geometric phases, obtained by evolving the magnetization of
fictitious spin-1/2 subspaces, over a closed loop on the Bloch sphere. The
dynamical phase accumulated during the evolution of the subspaces, is refocused
by a spin echo pulse sequence and by setting the delay of transition selective
pulses such that the evolution under the homonuclear coupling makes a complete
$2\pi$ rotation. A detailed theoretical explanation of non-adiabatic geometric
phases in NMR is given, by using single transition operators. Controlled phase
shift gates, two qubit Deutsch-Jozsa algorithm and parity algorithm in a
qubit-qutrit system have been implemented in various strongly dipolar coupled
systems obtained by orienting the molecules in liquid crystal media.
