function [pioDist] = pioDistance(pioA, pioB) pioDist = sum(sum( (pioA-pioB).^2 ));