I have a case where the heatmap of PureCN optimal purity/ploidy solutions for a sample shows the maximum likelihood solution (marked with a 1 on the heatmap) is located at what appears to NOT be a heatmap peak. I examined the heatmap numbers (used PureCN:::.matrixTotalPloidyToTumorPloidy() to get the heatmap matrix), and I find that the matrix number located at the #1 solution peak doesn't match either the log.likelihood value or the total.log.likelihood value of the #1 solution. Does the matrix factor in other numbers?