This allows int-vectors to be multiplied by float-scalars, producing
accurate int values (e.g. 4.2 * 10 == 42).

- separate data types for input arrays of binary operations
  (necessary for the equality comparison of two arrays of different
  type)
- added separate template parameter for ResultType, which defaults to
  the RealType
  (the goal is to allow sum of bool vectors into int instead of bool)

Refactoring the Dirichlet and the Neumann boundary conditions.

RowLengths was renamed to CompressedRowsLengths.
OpenMP support wads added to matrix-vector multiplication and to scalar
products.
CSR format can use Cusparse for SpMV.

Min. interations feature was fixed in the GMRES method.
Timers for explicit and semi-implicit timestepping were added.

Adding --fast-build optin for the build script to speed-up compilation.