THESE/codes/wave/allcode/pm4inv.m

   1 function [u, v, q, r] = pm4inv(a, b, c, d)
   2
   3 % [u, v, q, r] = pm4(a, b, c, d);
   4 % u = (a + b - c - d)/2;
   5 % v = (a - b + c + d)/2;
   6 % q = (a + b - c + d)/2;
   7 % r = (a + b + c - d)/2;
   8
   9 u = ( a + b + c + d)/2;
  10 v = (-a - b + c + d)/2;
  11 q = (-a + b - c + d)/2;
  12 r = (-a + b + c - d)/2;