def trimult(U, x, uplo='U', transa='N', alpha=1., inplace=False):
    """
    b = trimult(U,x, uplo='U')


    Multiplies U x, where U is upper triangular if uplo='U'
    or lower triangular if uplo = 'L'.
    """
    if inplace:
        b = x
    else:
        b = x.copy('F')
    dtrmm_wrap(a=U, b=b, uplo=uplo, transa=transa, alpha=alpha)
    return b
b = trimult(U,x, uplo='U')


    Multiplies U x, where U is upper triangular if uplo='U'
    or lower triangular if uplo = 'L'.