A series of bug fixes including: 1) gradient calculation error;
2) addressing boundary edges interpolation in certain rare cases;
3) setting scalar values (if requested) in one fill_n operation;
4) correct treatment of normals.

Also the documentation was improved, and a new .cxx test case was
added.