This work presents a GPU (Graphical Processing Unit) accelerated spatial domain oriented face resolution enhancement algorithm based on the homogeneity levels and relative-ratios of the pixels with respect to its surrounding pixels. The algorithm has been developed, implemented as well as tested in the MATLAB environment. MATLAB is slow in processing but at the same time a resourceful environment for the development in the area of image processing owing to its extremely rich set of functions and programmer-friendly integrated development environment. However, to compensate for the speed loss in testing and implementation phase, we have made use of GPU computing i.e. done parallelization of the algorithm on NVIDIA GPU using CUDA (Compute Unified Device Architecture) interface in the MATLAB environment. It is a simple but efficient algorithm in which kernel matrices are created encoding the homogeneity levels and relative-ratios of all pixels in surrounding four quadrants. Kernel matrices are subsequently applied to reconstruct the HR (High-Resolution) version from input LR (Low-Resolution) facial image.