308 …ymin, ymax, *bins, *values, *weights, binSize, binLeft, binRight, maxHeight, mean, var, totwt = 0.;  in PetscDrawHGDraw()  local
335   mean      = 0.0;  in PetscDrawHGDraw()
348       mean += values[p] * weights[p];  in PetscDrawHGDraw()
359         mean /= totwt;  in PetscDrawHGDraw()
360         var = (var - totwt * mean * mean) / totwt;  in PetscDrawHGDraw()
364         mean /= numValues;  in PetscDrawHGDraw()
365         if (numValues > 1) var = (var - numValues * mean * mean) / (PetscReal)(numValues - 1);  in PetscDrawHGDraw()
370       PetscCall(PetscSNPrintf(title, 256, "Mean: %g  Var: %g", (double)mean, (double)var));  in PetscDrawHGDraw()
409           mean += values[p] * weights[p];  in PetscDrawHGDraw()
422         mean /= totwt;  in PetscDrawHGDraw()
423         var = (var - totwt * mean * mean) / totwt;  in PetscDrawHGDraw()
427         mean /= numValues;  in PetscDrawHGDraw()
428         if (numValues > 1) var = (var - numValues * mean * mean) / (PetscReal)(numValues - 1);  in PetscDrawHGDraw()
433       PetscCall(PetscSNPrintf(title, 256, "Mean: %g  Var: %g", (double)mean, (double)var));  in PetscDrawHGDraw()
494   PetscReal xmax, xmin, *bins, *values, *weights, binSize, binLeft, binRight, mean, var, totwt = 0.;  in PetscDrawHGView()  local
512   mean      = 0.0;  in PetscDrawHGView()
520       mean += values[p] * weights[p];  in PetscDrawHGView()
552           mean += values[p] * weights[p];  in PetscDrawHGView()
571       mean /= totwt;  in PetscDrawHGView()
572       var = (var - totwt * mean * mean) / totwt;  in PetscDrawHGView()
576       mean /= numValues;  in PetscDrawHGView()
577       if (numValues > 1) var = (var - numValues * mean * mean) / (PetscReal)(numValues - 1);  in PetscDrawHGView()
582     PetscCall(PetscViewerASCIIPrintf(viewer, "Mean: %g  Var: %g\n", (double)mean, (double)var));  in PetscDrawHGView()