matrix discretization work