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.