How to use gulp to calculate defect formation energy

Dear Sir/Madam,
I calculated that Zn was doped with ZrO2 of cell expansion 222, and I am now randomly doped with two Zn, I want to compare the defect formation energy of different doping methods, please ask how to calculate the defect can be compared?
Best regards,