We present a fully non-perturbative determination of a relativistic heavy quark action’s parameters
on the CLS $n_f=2+1$ Wilson-clover ensembles using neural networks. We then
further illustrate the applicability of such an approach for lattice NRQCD bottom quarks, and finally investigate some physics quantities under our tuning. In particular, we look at the excited spectrum of bottomonia, a popular $ud\bar{b}\bar{b}$ tetraquark candidate, and the not-yet observed bottom-strange cousins of the exotic $J^P=0^+$ $D_{s0}^*(2317)$ and $J^P=1^+$ $D_{s1}(2460)$ mesons.