public static function Trxp(array $r, array $p, array &$trp) {
    $tr = [];

    /* Transpose of matrix r. */
    IAU::Tr($r, $tr);

    /* Matrix tr * vector p -> vector trp. */
    IAU::Rxp($tr, $p, $trp);

    return;
  }