Python equivalent of MATLAB’s “ismember” function

After many attempts trying optimize code, it seems that one last resource would be to attempt to run the code below using multiple cores. I don’t know exactly how to convert/re-structure my code so that it can run much faster using multiple cores. I will appreciate if I could get guidance to achieve the end goal. The end goal is to be able to run this code as fast as possible for arrays A and B where each array holds about 700,000 elements. Here is the code using small arrays. The 700k element arrays are commented out.