Abstract. The BPS/CFT correspondence predicts that the partition function and other BPS observables in a 3d N=2 gauge theory can be mapped to states in a q-deformed CFT, i.e. a theory whose symmetry algebra is q-Virasoro. Following this logic, we derive q-Virasoro constraints for the generating functions of 1/2-BPS Wilson loops for a family of 3d N=2 gauge theories. Under certain assumptions on the Chern-Simons level and the number of flavors we show that a complete solution can be obtained via a recursion formula. Furthermore, we show that the generating function admits a semi-classical expansion around an Hermitian matrix model with polynomial potential of degree equal to the number of flavors. Finally, we compute expectation values of Macdonald polynomials and a character expansion for the generating function. The talk is based on https://arxiv.org/abs/2007.10354.