cmat_reg and sqr roots now put right
ft_mats are the same - modulo a transpose ??