Bilinear forms with trace functions over arbitrary sets, and applications to Sato-Tate

We prove non-trivial upper bounds for general bilinear forms with trace functions of bountiful sheaves, where the supports of two variables can be arbitrary subsets in $\mathbf{F}_p$ of suitable sizes. This essentially recovers the P\'olya-Vinogradov range, and also applies to symmetric powers of Kloosterman sums and Frobenius traces of elliptic curves. In the case of hyper-Kloosterman sums, we can beat the P\'olya-Vinogradov barrier by combining additive combinatorics with a deep result of Kowalski, Michel and Sawin on sum-products of Kloosterman sheaves. Two Sato-Tate distributions of Kloosterman sums and Frobenius traces of elliptic curves in sparse families are also concluded.