This is a port of the proximal gradient method from `https://github.com/pmelchior/proxmin`. Original code credits to the authors of that package, errors are my own.