ugust 2012. Crudification operator added. Switched to cudaMalloc-based memory
allocation. Added several minor functions to Piece and Poly interfaces.
January 2013. A bug fixed. Some Host<dataType> classes missed protection
against access via negative index. In some cases, summation of PiecewisePoly
classes with overlapping supports was working incorrectly.
April 2013. The code is reconstructed to reduce flow control statements in the
June 2013. Crudification operator corrected. LPoly added.