We develop a novel approach to constrain the Hubble parameter H ₀ and the primordial power spectrum amplitude A _s using type Ia supernovae (SNIa) data. By considering SNIa as tracers of the peculiar velocity field, we can model their distance and their covariance as a function of cosmological parameters without the need of calibrators like Cepheids; this yields a new independent probe of the large-scale structure based on SNIa data without distance anchors. Crucially, we implement a differentiable pipeline in JAX , including efficient emulators and affine sampling, reducing inference time from years to hours on a single GPU. We first validate our method on mock datasets, demonstrating that we can constrain H ₀ and log 10 ¹⁰ A _s within 10% and 15%, respectively, using 𝒪(10 ³ ) SNIa. We then test our pipeline with SNIa from an N -body simulation, obtaining 6%-level unbiased constraints on H ₀ with a moderate noise level. We finally apply our method to Pantheon+ data, constraining H ₀ at the 15% level without Cepheids when fixing A _s to its Planck value. On the other hand, we obtain 20%-level constraints on log 10 ¹⁰ A _s in agreement with Planck when including Cepheids in the analysis. In light of upcoming observations of low redshift SNIa from the Zwicky Transient Facility and the Vera Rubin Legacy Survey of Space and Time, surveys for which our method will develop its full potential, we make our veloce ( https://github.com/dpiras/veloce ) code publicly available.