options_.gradient method = {12,22} two points formula.
options_.gradient method = {13,23} three points formula.
options_.gradient method = {15,25} four points formula.

Still need to decide what to do with (local) option zgrad.