Procedure is 90% done

Missing: convolution of v, stopping criterion, etc...