<svg xmlns="http://www.w3.org/2000/svg" width="24" height="24" viewBox="0 0 24 24" fill="none" stroke="currentColor" stroke-width="2" stroke-linecap="round" stroke-linejoin="round" class="feather feather-triangle"><path d="M10.29 3.86L1.82 18a2 2 0 0 0 1.71 3h16.94a2 2 0 0 0 1.71-3L13.71 3.86a2 2 0 0 0-3.42 0z"></path></svg> type='application/atom+xml'/>
<link rel='vcs-git' href='https://git.libcamera.org/libcamera/libcamera.git' title='libcamera/libcamera.git Git repository'/>
</head>
<body>
<div id='cgit'><table id='header'>
<tr>
<td class='logo' rowspan='2'><a href='/'><img src='/logo.png' alt='cgit logo'/></a></td>
<td class='main'><a href='/'>index</a> : <a title='libcamera/libcamera.git' href='/libcamera/libcamera.git/'>libcamera/libcamera.git</a></td><td class='form'><form method='get'>
<input type='hidden' name='id' value='b3befe2ef03794f5d8c2241e5809cb7dca8735e7'/><select name='h' onchange='this.form.submit();'>
<option value='master'>master</option>
</select> <input type='submit' value='switch'/></form></td></tr>
<tr><td class='sub'>libcamera official repository</td><td class='sub right'>git repository hosting on libcamera.org</td></tr></table>
<table class='tabs'><tr><td>
<a href='/libcamera/libcamera.git/?h=v0.2.0'>summary</a><a href='/libcamera/libcamera.git/refs/?h=v0.2.0&id=b3befe2ef03794f5d8c2241e5809cb7dca8735e7'>refs</a><a href='/libcamera/libcamera.git/log/src/ipa/raspberrypi/controller/histogram.cpp?h=v0.2.0'>log</a><a class='active' href='/libcamera/libcamera.git/tree/src/ipa/raspberrypi/controller/histogram.cpp?h=v0.2.0&id=b3befe2ef03794f5d8c2241e5809cb7dca8735e7'>tree</a><a href='/libcamera/libcamera.git/commit/src/ipa/raspberrypi/controller/histogram.cpp?h=v0.2.0&id=b3befe2ef03794f5d8c2241e5809cb7dca8735e7'>commit</a><a href='/libcamera/libcamera.git/diff/src/ipa/raspberrypi/controller/histogram.cpp?h=v0.2.0&id=b3befe2ef03794f5d8c2241e5809cb7dca8735e7'>diff</a></td><td class='form'><form class='right' method='get' action='/libcamera/libcamera.git/log/src/ipa/raspberrypi/controller/histogram.cpp'>
<input type='hidden' name='h' value='v0.2.0'/><input type='hidden' name='id' value='b3befe2ef03794f5d8c2241e5809cb7dca8735e7'/><select name='qt'>
<option value='grep'>log msg</option>
<option value='author'>author</option>
<option value='committer'>committer</option>
<option value='range'>range</option>
</select>
<input class='txt' type='search' size='10' name='q' value=''/>
<input type='submit' value='search'/>
</form>
</td></tr></table>
<div class='path'>path: <a href='/libcamera/libcamera.git/tree/?h=v0.2.0&id=b3befe2ef03794f5d8c2241e5809cb7dca8735e7'>root</a>/<a href='/libcamera/libcamera.git/tree/src?h=v0.2.0&id=b3befe2ef03794f5d8c2241e5809cb7dca8735e7'>src</a>/<a href='/libcamera/libcamera.git/tree/src/ipa?h=v0.2.0&id=b3befe2ef03794f5d8c2241e5809cb7dca8735e7'>ipa</a>/<a href='/libcamera/libcamera.git/tree/src/ipa/raspberrypi?h=v0.2.0&id=b3befe2ef03794f5d8c2241e5809cb7dca8735e7'>raspberrypi</a>/<a href='/libcamera/libcamera.git/tree/src/ipa/raspberrypi/controller?h=v0.2.0&id=b3befe2ef03794f5d8c2241e5809cb7dca8735e7'>controller</a>/<a href='/libcamera/libcamera.git/tree/src/ipa/raspberrypi/controller/histogram.cpp?h=v0.2.0&id=b3befe2ef03794f5d8c2241e5809cb7dca8735e7'>histogram.cpp</a></div><div class='content'>blob: 9916b3ed7f717a82ca35a3b451e15ef0d921282a (<a href='/libcamera/libcamera.git/plain/src/ipa/raspberrypi/controller/histogram.cpp?h=v0.2.0&id=b3befe2ef03794f5d8c2241e5809cb7dca8735e7'>plain</a>)
<table summary='blob content' class='blob'>
<tr><td class='linenumbers'><pre><a id='n1' href='#n1'>1</a>
<a id='n2' href='#n2'>2</a>
<a id='n3' href='#n3'>3</a>
<a id='n4' href='#n4'>4</a>
<a id='n5' href='#n5'>5</a>
<a id='n6' href='#n6'>6</a>
<a id='n7' href='#n7'>7</a>
<a id='n8' href='#n8'>8</a>
<a id='n9' href='#n9'>9</a>
<a id='n10' href='#n10'>10</a>
<a id='n11' href='#n11'>11</a>
<a id='n12' href='#n12'>12</a>
<a id='n13' href='#n13'>13</a>
<a id='n14' href='#n14'>14</a>
<a id='n15' href='#n15'>15</a>
<a id='n16' href='#n16'>16</a>
<a id='n17' href='#n17'>17</a>
<a id='n18' href='#n18'>18</a>
<a id='n19' href='#n19'>19</a>
<a id='n20' href='#n20'>20</a>
<a id='n21' href='#n21'>21</a>
<a id='n22' href='#n22'>22</a>
<a id='n23' href='#n23'>23</a>
<a id='n24' href='#n24'>24</a>
<a id='n25' href='#n25'>25</a>
<a id='n26' href='#n26'>26</a>
<a id='n27' href='#n27'>27</a>
<a id='n28' href='#n28'>28</a>
<a id='n29' href='#n29'>29</a>
<a id='n30' href='#n30'>30</a>
<a id='n31' href='#n31'>31</a>
<a id='n32' href='#n32'>32</a>
<a id='n33' href='#n33'>33</a>
<a id='n34' href='#n34'>34</a>
<a id='n35' href='#n35'>35</a>
<a id='n36' href='#n36'>36</a>
<a id='n37' href='#n37'>37</a>
<a id='n38' href='#n38'>38</a>
<a id='n39' href='#n39'>39</a>
<a id='n40' href='#n40'>40</a>
<a id='n41' href='#n41'>41</a>
<a id='n42' href='#n42'>42</a>
<a id='n43' href='#n43'>43</a>
<a id='n44' href='#n44'>44</a>
<a id='n45' href='#n45'>45</a>
<a id='n46' href='#n46'>46</a>
<a id='n47' href='#n47'>47</a>
<a id='n48' href='#n48'>48</a>
<a id='n49' href='#n49'>49</a>
<a id='n50' href='#n50'>50</a>
<a id='n51' href='#n51'>51</a>
<a id='n52' href='#n52'>52</a>
<a id='n53' href='#n53'>53</a>
<a id='n54' href='#n54'>54</a>
<a id='n55' href='#n55'>55</a>
<a id='n56' href='#n56'>56</a>
<a id='n57' href='#n57'>57</a>
<a id='n58' href='#n58'>58</a>
<a id='n59' href='#n59'>59</a>
<a id='n60' href='#n60'>60</a>
<a id='n61' href='#n61'>61</a>
<a id='n62' href='#n62'>62</a>
<a id='n63' href='#n63'>63</a>
<a id='n64' href='#n64'>64</a>
</pre></td>
<td class='lines'><pre><code><span class="hl com">/* SPDX-License-Identifier: BSD-2-Clause */</span>
<span class="hl com">/*</span>
<span class="hl com"> * Copyright (C) 2019, Raspberry Pi (Trading) Limited</span>
<span class="hl com"> *</span>
<span class="hl com"> * histogram.cpp - histogram calculations</span>
<span class="hl com"> */</span>
<span class="hl ppc">#include <math.h></span>
<span class="hl ppc">#include <stdio.h></span>
<span class="hl ppc">#include</span> <span class="hl pps">"histogram.hpp"</span><span class="hl ppc"></span>
<span class="hl kwa">using namespace</span> RPiController<span class="hl opt">;</span>
<span class="hl kwb">uint64_t</span> <span class="hl kwc">Histogram</span><span class="hl opt">::</span><span class="hl kwd">CumulativeFreq</span><span class="hl opt">(</span><span class="hl kwb">double</span> bin<span class="hl opt">)</span> <span class="hl kwb">const</span>
<span class="hl opt">{</span>
<span class="hl kwa">if</span> <span class="hl opt">(</span>bin <span class="hl opt"><=</span> <span class="hl num">0</span><span class="hl opt">)</span>
<span class="hl kwa">return</span> <span class="hl num">0</span><span class="hl opt">;</span>
<span class="hl kwa">else if</span> <span class="hl opt">(</span>bin <span class="hl opt">>=</span> <span class="hl kwd">Bins</span><span class="hl opt">())</span>
<span class="hl kwa">return</span> <span class="hl kwd">Total</span><span class="hl opt">();</span>
<span class="hl kwb">int</span> b <span class="hl opt">= (</span><span class="hl kwb">int</span><span class="hl opt">)</span>bin<span class="hl opt">;</span>
<span class="hl kwa">return</span> cumulative_<span class="hl opt">[</span>b<span class="hl opt">] +</span>
<span class="hl opt">(</span>bin <span class="hl opt">-</span> b<span class="hl opt">) * (</span>cumulative_<span class="hl opt">[</span>b <span class="hl opt">+</span> <span class="hl num">1</span><span class="hl opt">] -</span> cumulative_<span class="hl opt">[</span>b<span class="hl opt">]);</span>
<span class="hl opt">}</span>
<span class="hl kwb">double</span> <span class="hl kwc">Histogram</span><span class="hl opt">::</span><span class="hl kwd">Quantile</span><span class="hl opt">(</span><span class="hl kwb">double</span> q<span class="hl opt">,</span> <span class="hl kwb">int</span> first<span class="hl opt">,</span> <span class="hl kwb">int</span> last<span class="hl opt">)</span> <span class="hl kwb">const</span>
<span class="hl opt">{</span>
<span class="hl kwa">if</span> <span class="hl opt">(</span>first <span class="hl opt">== -</span><span class="hl num">1</span><span class="hl opt">)</span>
first <span class="hl opt">=</span> <span class="hl num">0</span><span class="hl opt">;</span>
<span class="hl kwa">if</span> <span class="hl opt">(</span>last <span class="hl opt">== -</span><span class="hl num">1</span><span class="hl opt">)</span>
last <span class="hl opt">=</span> cumulative_<span class="hl opt">.</span><span class="hl kwd">size</span><span class="hl opt">() -</span> <span class="hl num">2</span><span class="hl opt">;</span>
<span class="hl kwa">assert</span><span class="hl opt">(</span>first <span class="hl opt"><=</span> last<span class="hl opt">);</span>
<span class="hl kwb">uint64_t</span> items <span class="hl opt">=</span> q <span class="hl opt">*</span> <span class="hl kwd">Total</span><span class="hl opt">();</span>
<span class="hl kwa">while</span> <span class="hl opt">(</span>first <span class="hl opt"><</span> last<span class="hl opt">)</span> <span class="hl slc">// binary search to find the right bin</span>
<span class="hl opt">{</span>
<span class="hl kwb">int</span> middle <span class="hl opt">= (</span>first <span class="hl opt">+</span> last<span class="hl opt">) /</span> <span class="hl num">2</span><span class="hl opt">;</span>
<span class="hl kwa">if</span> <span class="hl opt">(</span>cumulative_<span class="hl opt">[</span>middle <span class="hl opt">+</span> <span class="hl num">1</span><span class="hl opt">] ></span> items<span class="hl opt">)</span>
last <span class="hl opt">=</span> middle<span class="hl opt">;</span> <span class="hl slc">// between first and middle</span>
<span class="hl kwa">else</span>
first <span class="hl opt">=</span> middle <span class="hl opt">+</span> <span class="hl num">1</span><span class="hl opt">;</span> <span class="hl slc">// after middle</span>
<span class="hl opt">}</span>
<span class="hl kwa">assert</span><span class="hl opt">(</span>items <span class="hl opt">>=</span> cumulative_<span class="hl opt">[</span>first<span class="hl opt">] &&</span> items <span class="hl opt"><=</span> cumulative_<span class="hl opt">[</span>last <span class="hl opt">+</span> <span class="hl num">1</span><span class="hl opt">]);</span>
<span class="hl kwb">double</span> frac <span class="hl opt">=</span> cumulative_<span class="hl opt">[</span>first <span class="hl opt">+</span> <span class="hl num">1</span><span class="hl opt">] ==</span> cumulative_<span class="hl opt">[</span>first<span class="hl opt">] ?</span> <span class="hl num">0</span>
<span class="hl opt">: (</span><span class="hl kwb">double</span><span class="hl opt">)(</span>items <span class="hl opt">-</span> cumulative_<span class="hl opt">[</span>first<span class="hl opt">]) /</span>
<span class="hl opt">(</span>cumulative_<span class="hl opt">[</span>first <span class="hl opt">+</span> <span class="hl num">1</span><span class="hl opt">] -</span> cumulative_<span class="hl opt">[</span>first<span class="hl opt">]);</span>
<span class="hl kwa">return</span> first <span class="hl opt">+</span> frac<span class="hl opt">;</span>
<span class="hl opt">}</span>
<span class="hl kwb">double</span> <span class="hl kwc">Histogram</span><span class="hl opt">::</span><span class="hl kwd">InterQuantileMean</span><span class="hl opt">(</span><span class="hl kwb">double</span> q_lo<span class="hl opt">,</span> <span class="hl kwb">double</span> q_hi<span class="hl opt">)</span> <span class="hl kwb">const</span>
<span class="hl opt">{</span>
<span class="hl kwa">assert</span><span class="hl opt">(</span>q_hi <span class="hl opt">></span> q_lo<span class="hl opt">);</span>
<span class="hl kwb">double</span> p_lo <span class="hl opt">=</span> <span class="hl kwd">Quantile</span><span class="hl opt">(</span>q_lo<span class="hl opt">);</span>
<span class="hl kwb">double</span> p_hi <span class="hl opt">=</span> <span class="hl kwd">Quantile</span><span class="hl opt">(</span>q_hi<span class="hl opt">, (</span><span class="hl kwb">int</span><span class="hl opt">)</span>p_lo<span class="hl opt">);</span>
<span class="hl kwb">double</span> sum_bin_freq <span class="hl opt">=</span> <span class="hl num">0</span><span class="hl opt">,</span> cumul_freq <span class="hl opt">=</span> <span class="hl num">0</span><span class="hl opt">;</span>
<span class="hl kwa">for</span> <span class="hl opt">(</span><span class="hl kwb">double</span> p_next <span class="hl opt">=</span> <span class="hl kwd">floor</span><span class="hl opt">(</span>p_lo<span class="hl opt">) +</span> <span class="hl num">1.0</span><span class="hl opt">;</span> p_next <span class="hl opt"><=</span> <span class="hl kwd">ceil</span><span class="hl opt">(</span>p_hi<span class="hl opt">);</span>
p_lo <span class="hl opt">=</span> p_next<span class="hl opt">,</span> p_next <span class="hl opt">+=</span> <span class="hl num">1.0</span><span class="hl opt">) {</span>
<span class="hl kwb">int</span> bin <span class="hl opt">=</span> <span class="hl kwd">floor</span><span class="hl opt">(</span>p_lo<span class="hl opt">);</span>
<span class="hl kwb">double</span> freq <span class="hl opt">= (</span>cumulative_<span class="hl opt">[</span>bin <span class="hl opt">+</span> <span class="hl num">1</span><span class="hl opt">] -</span> cumulative_<span class="hl opt">[</span>bin<span class="hl opt">]) *</span>
<span class="hl opt">(</span><span class="hl kwc">std</span><span class="hl opt">::</span><span class="hl kwd">min</span><span class="hl opt">(</span>p_next<span class="hl opt">,</span> p_hi<span class="hl opt">) -</span> p_lo<span class="hl opt">);</span>
sum_bin_freq <span class="hl opt">+=</span> bin <span class="hl opt">*</span> freq<span class="hl opt">;</span>
cumul_freq <span class="hl opt">+=</span> freq<span class="hl opt">;</span>
<span class="hl opt">}</span>
<span class="hl slc">// add 0.5 to give an average for bin mid-points</span>
<span class="hl kwa">return</span> sum_bin_freq <span class="hl opt">/</span> cumul_freq <span class="hl opt">+</span> <span class="hl num">0.5</span><span class="hl opt">;</span>
<span class="hl opt">}</span>
</code></pre></td></tr></table>
</div> <!-- class=content -->
<div class='footer'>generated by <a href='https://git.zx2c4.com/cgit/about/'>cgit v1.2.1</a> (<a href='https://git-scm.com/'>git 2.18.0</a>) at 2025-03-11 02:23:43 +0000</div>
</div> <!-- id=cgit -->
</body>
</html>