function [u, v] = pm(a,b)

% [u v] = pm(a,b)
% u = (a + b)/sqrt(2);
% v = (a - b)/sqrt(2);

u = (a + b)/sqrt(2);
v = (a - b)/sqrt(2);