Same here:
```c++
const double sigma_q = 1./n_dest_mesh_elements * std::accumulate(dest_properties, dest_properties + n_dest_mesh_elements, 0,
    [const &mu](const double& sum, const double& v) { return sum + MathLib::fastpow(v - mu, 2); });
```
but I'm not sure here.