test[x_]:=Boole[PrimeQ[x,GaussianIntegers->True]];
DistributeDefinitions[test];
data=ParallelTable[test[x+ y I],{x,-50,50},{y,-50,50}];
ArrayPlot[data]