Should VecSqrt() for complex numbers take the square root of the absolute value (as it does for real numbers) or should be actually be the square root? Currently it uses absolute values for all cases. Barry