The Born series methods, namely, traditional Born series (TBS) and convergent Born series (CBS), have been recently implemented to numerically solve the time-independent photoacoustic (PA) wave equation for an acoustically inhomogeneous source. The TBS algorithm diverges when the sound-speed mismatch is >20%, however, the CBS technique provides faithful result even beyond this limit. These protocols are iterative in nature and computationally expensive. Currently, MATLAB-based implementations are available to study PA emission from a single source mimicking a cell. However, efficient numerical implementation strategy is further needed particularly to calculate PA field from a tissue. Therefore, to develop insights, uniprocessorbased C codes were realized for these schemes. The PA field (in 2D) was computed at a distance 35 µm from a source (a light absorbing disc of radius 7.5 µm) over a frequency range from 7.32 to 512 MHz with a resolution of 7.32 MHz. The sound-speed within the source region was varied from v s = 1200, 1500 and 1800 m/s, but the same quantity for the ambient medium was fixed to v f = 1500 m/s. The C program was found to be at least 10 times slower than the corresponding MATLAB program. It may be because MATLAB inherently implements parallel computing while evaluating the forward and backward fast Fourier transforms (FFTs) associated with the Born series approaches. Multiprocessor-based FFTs and parallel nested loops are being incorporated into the C program for enhancement of its execution speed.