How to define laplacian in Flux or Zygote?

I want to solve a PDE using a neural network and I need to calculate the laplacian of my network.