Fast distributed phononic band-structure calculations through a GPU accelerated mixed-variational formulation