third order gradient optimization w/ option for sub-selecting the stats. still working on memory consumption issues for scaling this up further