* grain() should compile() by default.

* compileXXX -> gatherXXX (because compile gets too many meanings)

* use the fit.grain function

* describe what extractPOT / extractMARG do

