Calculating epistatic interactions between genomic variants in studies incorporating complex endophenotypes is a computationally challenging problem that requires emphasis on accelerating and parallelizing the code and on workload distribution efficiency. Accelerating and scaling this process will enable the detection of epistasis in many existing GWAS datasets, from both the biomedical and agricultural areas.