308 …ymin, ymax, *bins, *values, *weights, binSize, binLeft, binRight, maxHeight, mean, var, totwt = 0.;  in PetscDrawHGDraw()  local
336   var       = 0.0;  in PetscDrawHGDraw()
349       var += values[p] * values[p] * weights[p];  in PetscDrawHGDraw()
360         var = (var - totwt * mean * mean) / totwt;  in PetscDrawHGDraw()
365         if (numValues > 1) var = (var - numValues * mean * mean) / (PetscReal)(numValues - 1);  in PetscDrawHGDraw()
366         else var = 0.0;  in PetscDrawHGDraw()
370       PetscCall(PetscSNPrintf(title, 256, "Mean: %g  Var: %g", (double)mean, (double)var));  in PetscDrawHGDraw()
410           var += values[p] * values[p] * weights[p];  in PetscDrawHGDraw()
423         var = (var - totwt * mean * mean) / totwt;  in PetscDrawHGDraw()
428         if (numValues > 1) var = (var - numValues * mean * mean) / (PetscReal)(numValues - 1);  in PetscDrawHGDraw()
429         else var = 0.0;  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
513   var       = 0.0;  in PetscDrawHGView()
521       var += values[p] * values[p] * weights[p];  in PetscDrawHGView()
553           var += values[p] * values[p] * weights[p];  in PetscDrawHGView()
572       var = (var - totwt * mean * mean) / totwt;  in PetscDrawHGView()
577       if (numValues > 1) var = (var - numValues * mean * mean) / (PetscReal)(numValues - 1);  in PetscDrawHGView()
578       else var = 0.0;  in PetscDrawHGView()
582     PetscCall(PetscViewerASCIIPrintf(viewer, "Mean: %g  Var: %g\n", (double)mean, (double)var));  in PetscDrawHGView()