Matrix-free is explained in step-67 https://dealii.org/current/doxygen/deal.II/step_67.html Modify `system_legendre_mpi` to use above approach.