code: def _symbolic_product_helper():
      
      from sympy import symbols, Matrix
  
      D11, D12, D13, D21, D22, D23, D31, D32, D33 =\
          symbols('D11, D12, D13, D21, D22, D23, D31, D32, D33')
  
      D = Matrix([[D11, D12, D13], [D21, D22, D23], [D31, D32, D33]])
      grad = Matrix([['dx', 'dy', 'dz']]).T
      div = grad.T
  
      a = div * D * grad
  
      print(a[0])
docs: Use SymPy to generate the 3D products for diffusion_stencil_3d.