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.