GCC 6.3.0 fails to compile fillCovarianceMatrix()
GCC apparently cannot deduce the template expressions:
error: expected primary-expression before '>' token
if ((*traits).config.get<std::string>("stochastic.anisotropy","none") == "none")
^