1915   PetscReal     *x, *w, *xw, *ww;  in PetscDTGaussTensorQuadrature()  local
1934     PetscCall(PetscMalloc1(npoints, &ww));  in PetscDTGaussTensorQuadrature()
1935     PetscCall(PetscDTGaussQuadrature(npoints, a, b, x, ww));  in PetscDTGaussTensorQuadrature()
1937       for (PetscInt c = 0; c < Nc; ++c) w[i * Nc + c] = ww[i];  in PetscDTGaussTensorQuadrature()
1938     PetscCall(PetscFree(ww));  in PetscDTGaussTensorQuadrature()
1942     PetscCall(PetscMalloc2(npoints, &xw, npoints, &ww));  in PetscDTGaussTensorQuadrature()
1943     PetscCall(PetscDTGaussQuadrature(npoints, a, b, xw, ww));  in PetscDTGaussTensorQuadrature()
1948         for (PetscInt c = 0; c < Nc; ++c) w[(i * npoints + j) * Nc + c] = ww[i] * ww[j];  in PetscDTGaussTensorQuadrature()
1951     PetscCall(PetscFree2(xw, ww));  in PetscDTGaussTensorQuadrature()
1955     PetscCall(PetscMalloc2(npoints, &xw, npoints, &ww));  in PetscDTGaussTensorQuadrature()
1956     PetscCall(PetscDTGaussQuadrature(npoints, a, b, xw, ww));  in PetscDTGaussTensorQuadrature()
1963 …PetscInt c = 0; c < Nc; ++c) w[((i * npoints + j) * npoints + k) * Nc + c] = ww[i] * ww[j] * ww[k];  in PetscDTGaussTensorQuadrature()
1967     PetscCall(PetscFree2(xw, ww));  in PetscDTGaussTensorQuadrature()