Observations at sub-mm/mm wavelengths allow us to probe the solids in the interior of protoplanetary disks, where the bulk of the dust is located and planet formation is expected to occur. However, the actual size of dust grains is still largely unknown due to the limited angular resolution and sensitivity of past observations. The upgraded VLA and, especially, the ALMA observatories provide now powerful tools to resolve grain growth in disks, making the time ripe for developing a multi-wavelength analysis of sub-mm/mm observations of disks. In my contribution I will present a novel analysis method for multi-wavelength ALMA/VLA observations which, based on the self-consistent modelling of the sub-mm/mm disk continuum emission, allows us to constrain simultaneously the size distribution of dust grains and the disk's physical structure (Tazzari et al. 2016, A&A 588 A53). I will also present the recent analysis of spatially resolved ALMA Band 7 observations of a large sample of disks in the Lupus star forming region, from which we obtained a tentative evidence of a disk size-disk mass correlation (Tazzari et al. 2017, arXiv:1707.01499). Finally, I will introduce galario, a GPU Accelerated Library for the Analysis of Radio Interferometry Observations. Fitting the observed visibilities in the uv-plane is computationally demanding: with galario we solve this problem for the current as well as for the full-science ALMA capabilities by leveraging on the computing power of GPUs, providing the computational breakthrough needed to fully exploit the new wealth of information delivered by ALMA.