+ plogp(T)) / float(T)</span>
        <span class="c1"># hyz = (-(plogp(hist.sum(axis=0))).sum() + plogp(T)) / float(T)</span>
        <span class="c1"># hz = (-(plogp(hist.sum(axis=0).sum(axis=0))).sum() + plogp(T)) / float(T)</span>

        <span class="c1"># Multivariate X, Y version</span>
        <span class="n">plogp</span> <span class="o">=</span> <span class="n">_plogp_vector</span><span class="p">(</span><span class="n">T</span><span class="p">)</span>
        <span class="n">hxyz</span> <span cla