2977   const PetscReal *gllnodes = nodes;  in PetscGaussLobattoLegendreElementLaplacianCreate()  local
2988     x  = gllnodes[j];  in PetscGaussLobattoLegendreElementLaplacianCreate()
2999         A[j][j] = 2. / (3. * (1. - gllnodes[j] * gllnodes[j]) * Lpj * Lpj);  in PetscGaussLobattoLegendreElementLaplacianCreate()
3001         x  = gllnodes[r];  in PetscGaussLobattoLegendreElementLaplacianCreate()
3010 …tscReal)p) * (((PetscReal)p) + 1.) * Lpj * Lpr * (gllnodes[j] - gllnodes[r]) * (gllnodes[j] - glln…  in PetscGaussLobattoLegendreElementLaplacianCreate()
3015     x  = gllnodes[j];  in PetscGaussLobattoLegendreElementLaplacianCreate()
3024 …-1., p) / (((PetscReal)p) * (((PetscReal)p) + 1.) * Lpj * (1. + gllnodes[j]) * (1. + gllnodes[j]));  in PetscGaussLobattoLegendreElementLaplacianCreate()
3028     x  = gllnodes[j];  in PetscGaussLobattoLegendreElementLaplacianCreate()
3038 …j] = 4. / (((PetscReal)p) * (((PetscReal)p) + 1.) * Lpj * (1. - gllnodes[j]) * (1. - gllnodes[j]));  in PetscGaussLobattoLegendreElementLaplacianCreate()
3097   const PetscReal *gllnodes = nodes;  in PetscGaussLobattoLegendreElementGradientCreate()  local
3118       PetscCall(PetscDTComputeJacobi(0., 0., p, gllnodes[i], &Li));  in PetscGaussLobattoLegendreElementGradientCreate()
3119       PetscCall(PetscDTComputeJacobi(0., 0., p, gllnodes[j], &Lj));  in PetscGaussLobattoLegendreElementGradientCreate()
3120       if (i != j) A[i][j] = Li / (Lj * (gllnodes[i] - gllnodes[j]));  in PetscGaussLobattoLegendreElementGradientCreate()