2881 PetscErrorCode PetscDTReconstructPoly(PetscInt degree, PetscInt nsource, const PetscReal sourcex[],…  in PetscDTReconstructPoly()  argument
2891 …nsource, PETSC_COMM_SELF, PETSC_ERR_ARG_INCOMP, "Reconstruction degree %" PetscInt_FMT " must be l…  in PetscDTReconstructPoly()
2893 …for (i = 0; i < nsource; i++) PetscCheck(sourcex[i] < sourcex[i + 1], PETSC_COMM_SELF, PETSC_ERR_A…  in PetscDTReconstructPoly()
2897   xmax     = PetscMax(sourcex[nsource], targetx[ntarget]);  in PetscDTReconstructPoly()
2900   worksize = nsource;  in PetscDTReconstructPoly()
2901 …PetscCall(PetscMalloc4(degree + 1, &bdegrees, nsource + 1, &sourcey, nsource * (degree + 1), &Bsou…  in PetscDTReconstructPoly()
2902 …PetscCall(PetscMalloc4(nsource, &tau, nsource * (degree + 1), &Bsinv, ntarget + 1, &targety, ntarg…  in PetscDTReconstructPoly()
2903   for (i = 0; i <= nsource; i++) sourcey[i] = (sourcex[i] - center) / hscale;  in PetscDTReconstructPoly()
2905   PetscCall(PetscDTLegendreIntegrate(nsource, sourcey, degree + 1, bdegrees, PETSC_TRUE, Bsource));  in PetscDTReconstructPoly()
2906 …PetscCall(PetscDTPseudoInverseQR(nsource, nsource, degree + 1, Bsource, Bsinv, tau, nsource, work)…  in PetscDTReconstructPoly()
2911     for (j = 0; j < nsource; j++) {  in PetscDTReconstructPoly()
2913 …    for (k = 0; k < degree + 1; k++) sum += Btarget[i * (degree + 1) + k] * Bsinv[k * nsource + j];  in PetscDTReconstructPoly()
2914       R[i * nsource + j] = sum;  in PetscDTReconstructPoly()
2917     for (j = 0; j < nsource; j++) R[i * nsource + j] /= rowsum; /* normalize each row */  in PetscDTReconstructPoly()