1d_diffusion_solver