fixed diffusivity
fixed compile -- as-cran works