def METHOD_NAME(pressure: np.ndarray) -> np.ndarray:
    """compute the velocity corresponding to a guessed pressure"""
    velocity = basis['u'].zeros()
    velocity[I] = solve(Aint,
                        condense(csr_matrix(A.shape),
                                 f - B.T @ pressure, I=I)[1],
                        solver=solver)
    return velocity