536 lines
107 KiB
HTML
536 lines
107 KiB
HTML
<!-- HTML header for doxygen 1.8.18-->
|
|
<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "https://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd">
|
|
<html xmlns="http://www.w3.org/1999/xhtml">
|
|
<head>
|
|
<meta http-equiv="Content-Type" content="text/xhtml;charset=UTF-8"/>
|
|
<meta http-equiv="X-UA-Compatible" content="IE=9"/>
|
|
<meta name="generator" content="Doxygen 1.9.1"/>
|
|
<meta name="viewport" content="width=device-width, initial-scale=1"/>
|
|
<title>OR-Tools: lp_utils.h Source File</title>
|
|
<link href="tabs.css" rel="stylesheet" type="text/css"/>
|
|
<script type="text/javascript" src="jquery.js"></script>
|
|
<script type="text/javascript" src="dynsections.js"></script>
|
|
<link href="navtree.css" rel="stylesheet" type="text/css"/>
|
|
<script type="text/javascript" src="resize.js"></script>
|
|
<script type="text/javascript" src="navtreedata.js"></script>
|
|
<script type="text/javascript" src="navtree.js"></script>
|
|
<link href="search/search.css" rel="stylesheet" type="text/css"/>
|
|
<script type="text/javascript" src="search/searchdata.js"></script>
|
|
<script type="text/javascript" src="search/search.js"></script>
|
|
<link href="doxygen.css" rel="stylesheet" type="text/css" />
|
|
<link href="styleSheet.tmp.css" rel="stylesheet" type="text/css"/>
|
|
</head>
|
|
<body>
|
|
<div id="top"><!-- do not remove this div, it is closed by doxygen! -->
|
|
<div id="titlearea">
|
|
<table cellspacing="0" cellpadding="0">
|
|
<tbody>
|
|
<tr style="height: 56px;">
|
|
<td id="projectlogo"><img alt="Logo" src="orLogo.png"/></td>
|
|
<td id="projectalign" style="padding-left: 0.5em;">
|
|
<div id="projectname">OR-Tools
|
|
 <span id="projectnumber">8.2</span>
|
|
</div>
|
|
</td>
|
|
</tr>
|
|
</tbody>
|
|
</table>
|
|
</div>
|
|
<!-- end header part -->
|
|
<!-- Generated by Doxygen 1.9.1 -->
|
|
<script type="text/javascript">
|
|
/* @license magnet:?xt=urn:btih:cf05388f2679ee054f2beb29a391d25f4e673ac3&dn=gpl-2.0.txt GPL-v2 */
|
|
var searchBox = new SearchBox("searchBox", "search",false,'Search','.html');
|
|
/* @license-end */
|
|
</script>
|
|
<script type="text/javascript" src="menudata.js"></script>
|
|
<script type="text/javascript" src="menu.js"></script>
|
|
<script type="text/javascript">
|
|
/* @license magnet:?xt=urn:btih:cf05388f2679ee054f2beb29a391d25f4e673ac3&dn=gpl-2.0.txt GPL-v2 */
|
|
$(function() {
|
|
initMenu('',true,false,'search.php','Search');
|
|
$(document).ready(function() { init_search(); });
|
|
});
|
|
/* @license-end */</script>
|
|
<div id="main-nav"></div>
|
|
</div><!-- top -->
|
|
<div id="side-nav" class="ui-resizable side-nav-resizable">
|
|
<div id="nav-tree">
|
|
<div id="nav-tree-contents">
|
|
<div id="nav-sync" class="sync"></div>
|
|
</div>
|
|
</div>
|
|
<div id="splitbar" style="-moz-user-select:none;"
|
|
class="ui-resizable-handle">
|
|
</div>
|
|
</div>
|
|
<script type="text/javascript">
|
|
/* @license magnet:?xt=urn:btih:cf05388f2679ee054f2beb29a391d25f4e673ac3&dn=gpl-2.0.txt GPL-v2 */
|
|
$(document).ready(function(){initNavTree('lp__data_2lp__utils_8h_source.html',''); initResizable(); });
|
|
/* @license-end */
|
|
</script>
|
|
<div id="doc-content">
|
|
<!-- window showing the filter options -->
|
|
<div id="MSearchSelectWindow"
|
|
onmouseover="return searchBox.OnSearchSelectShow()"
|
|
onmouseout="return searchBox.OnSearchSelectHide()"
|
|
onkeydown="return searchBox.OnSearchSelectKey(event)">
|
|
</div>
|
|
|
|
<!-- iframe showing the search results (closed by default) -->
|
|
<div id="MSearchResultsWindow">
|
|
<iframe src="javascript:void(0)" frameborder="0"
|
|
name="MSearchResults" id="MSearchResults">
|
|
</iframe>
|
|
</div>
|
|
|
|
<div class="header">
|
|
<div class="headertitle">
|
|
<div class="title">lp_data/lp_utils.h</div> </div>
|
|
</div><!--header-->
|
|
<div class="contents">
|
|
<a href="lp__data_2lp__utils_8h.html">Go to the documentation of this file.</a><div class="fragment"><div class="line"><a name="l00001"></a><span class="lineno"> 1</span> <span class="comment">// Copyright 2010-2018 Google LLC</span></div>
|
|
<div class="line"><a name="l00002"></a><span class="lineno"> 2</span> <span class="comment">// Licensed under the Apache License, Version 2.0 (the "License");</span></div>
|
|
<div class="line"><a name="l00003"></a><span class="lineno"> 3</span> <span class="comment">// you may not use this file except in compliance with the License.</span></div>
|
|
<div class="line"><a name="l00004"></a><span class="lineno"> 4</span> <span class="comment">// You may obtain a copy of the License at</span></div>
|
|
<div class="line"><a name="l00005"></a><span class="lineno"> 5</span> <span class="comment">//</span></div>
|
|
<div class="line"><a name="l00006"></a><span class="lineno"> 6</span> <span class="comment">// http://www.apache.org/licenses/LICENSE-2.0</span></div>
|
|
<div class="line"><a name="l00007"></a><span class="lineno"> 7</span> <span class="comment">//</span></div>
|
|
<div class="line"><a name="l00008"></a><span class="lineno"> 8</span> <span class="comment">// Unless required by applicable law or agreed to in writing, software</span></div>
|
|
<div class="line"><a name="l00009"></a><span class="lineno"> 9</span> <span class="comment">// distributed under the License is distributed on an "AS IS" BASIS,</span></div>
|
|
<div class="line"><a name="l00010"></a><span class="lineno"> 10</span> <span class="comment">// WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.</span></div>
|
|
<div class="line"><a name="l00011"></a><span class="lineno"> 11</span> <span class="comment">// See the License for the specific language governing permissions and</span></div>
|
|
<div class="line"><a name="l00012"></a><span class="lineno"> 12</span> <span class="comment">// limitations under the License.</span></div>
|
|
<div class="line"><a name="l00013"></a><span class="lineno"> 13</span>  </div>
|
|
<div class="line"><a name="l00014"></a><span class="lineno"> 14</span> <span class="comment">// Basic utility functions on Fractional or row/column of Fractional.</span></div>
|
|
<div class="line"><a name="l00015"></a><span class="lineno"> 15</span>  </div>
|
|
<div class="line"><a name="l00016"></a><span class="lineno"> 16</span> <span class="preprocessor">#ifndef OR_TOOLS_LP_DATA_LP_UTILS_H_</span></div>
|
|
<div class="line"><a name="l00017"></a><span class="lineno"> 17</span> <span class="preprocessor">#define OR_TOOLS_LP_DATA_LP_UTILS_H_</span></div>
|
|
<div class="line"><a name="l00018"></a><span class="lineno"> 18</span>  </div>
|
|
<div class="line"><a name="l00019"></a><span class="lineno"> 19</span> <span class="preprocessor">#include "<a class="code" href="accurate__sum_8h.html">ortools/base/accurate_sum.h</a>"</span></div>
|
|
<div class="line"><a name="l00020"></a><span class="lineno"> 20</span> <span class="preprocessor">#include "<a class="code" href="lp__types_8h.html">ortools/lp_data/lp_types.h</a>"</span></div>
|
|
<div class="line"><a name="l00021"></a><span class="lineno"> 21</span> <span class="preprocessor">#include "<a class="code" href="scattered__vector_8h.html">ortools/lp_data/scattered_vector.h</a>"</span></div>
|
|
<div class="line"><a name="l00022"></a><span class="lineno"> 22</span> <span class="preprocessor">#include "<a class="code" href="sparse__column_8h.html">ortools/lp_data/sparse_column.h</a>"</span></div>
|
|
<div class="line"><a name="l00023"></a><span class="lineno"> 23</span>  </div>
|
|
<div class="line"><a name="l00024"></a><span class="lineno"> 24</span> <span class="keyword">namespace </span><a class="code" href="namespaceoperations__research.html">operations_research</a> {</div>
|
|
<div class="line"><a name="l00025"></a><span class="lineno"> 25</span> <span class="keyword">namespace </span>glop {</div>
|
|
<div class="line"><a name="l00026"></a><span class="lineno"> 26</span>  </div>
|
|
<div class="line"><a name="l00027"></a><span class="lineno"> 27</span> <span class="comment">// TODO(user): For some Fractional types, it may not gain much (or even nothing</span></div>
|
|
<div class="line"><a name="l00028"></a><span class="lineno"> 28</span> <span class="comment">// if we are in infinite precision) to use this sum. A solution is to templatize</span></div>
|
|
<div class="line"><a name="l00029"></a><span class="lineno"> 29</span> <span class="comment">// this class and specialize it to a normal sum for the Fractional type we want</span></div>
|
|
<div class="line"><a name="l00030"></a><span class="lineno"> 30</span> <span class="comment">// so in this case the PreciseXXX() functions below will become equivalent to</span></div>
|
|
<div class="line"><a name="l00031"></a><span class="lineno"> 31</span> <span class="comment">// their normal version.</span></div>
|
|
<div class="line"><a name="l00032"></a><span class="lineno"><a class="line" href="namespaceoperations__research_1_1glop.html#a394b55c9b64ac7ce7c9e8772ef9f7c73"> 32</a></span> <span class="keyword">typedef</span> <a class="code" href="classoperations__research_1_1_accurate_sum.html">AccurateSum<Fractional></a> <a class="code" href="namespaceoperations__research_1_1glop.html#a394b55c9b64ac7ce7c9e8772ef9f7c73">KahanSum</a>;</div>
|
|
<div class="line"><a name="l00033"></a><span class="lineno"> 33</span>  </div>
|
|
<div class="line"><a name="l00034"></a><span class="lineno"> 34</span> <span class="comment">// Returns the square of a Fractional.</span></div>
|
|
<div class="line"><a name="l00035"></a><span class="lineno"> 35</span> <span class="comment">// Useful to shorten the code when f is an expression or a long name.</span></div>
|
|
<div class="line"><a name="l00036"></a><span class="lineno"><a class="line" href="namespaceoperations__research_1_1glop.html#a1dcd08b0f6c19cd4a302bb5a3a6ea06e"> 36</a></span> <span class="keyword">inline</span> <a class="code" href="namespaceoperations__research_1_1glop.html#a733947145e3e1631165b618b05c9ccb7">Fractional</a> <a class="code" href="namespaceoperations__research_1_1glop.html#a1dcd08b0f6c19cd4a302bb5a3a6ea06e">Square</a>(<a class="code" href="namespaceoperations__research_1_1glop.html#a733947145e3e1631165b618b05c9ccb7">Fractional</a> f) { <span class="keywordflow">return</span> f * f; }</div>
|
|
<div class="line"><a name="l00037"></a><span class="lineno"> 37</span>  </div>
|
|
<div class="line"><a name="l00038"></a><span class="lineno"> 38</span> <span class="comment">// Returns distance from a given fractional number to the closest integer. It</span></div>
|
|
<div class="line"><a name="l00039"></a><span class="lineno"> 39</span> <span class="comment">// means that the result is always contained in range of [0.0, 0.5].</span></div>
|
|
<div class="line"><a name="l00040"></a><span class="lineno"><a class="line" href="namespaceoperations__research_1_1glop.html#ac51a45a6f706c77b25cea909b9d879db"> 40</a></span> <span class="keyword">static</span> <span class="keyword">inline</span> <a class="code" href="namespaceoperations__research_1_1glop.html#a733947145e3e1631165b618b05c9ccb7">Fractional</a> <a class="code" href="namespaceoperations__research_1_1glop.html#ac51a45a6f706c77b25cea909b9d879db">Fractionality</a>(<a class="code" href="namespaceoperations__research_1_1glop.html#a733947145e3e1631165b618b05c9ccb7">Fractional</a> f) {</div>
|
|
<div class="line"><a name="l00041"></a><span class="lineno"> 41</span>  <span class="keywordflow">return</span> std::abs(f - std::round(f));</div>
|
|
<div class="line"><a name="l00042"></a><span class="lineno"> 42</span> }</div>
|
|
<div class="line"><a name="l00043"></a><span class="lineno"> 43</span>  </div>
|
|
<div class="line"><a name="l00044"></a><span class="lineno"> 44</span> <span class="comment">// Returns the scalar product between u and v.</span></div>
|
|
<div class="line"><a name="l00045"></a><span class="lineno"> 45</span> <span class="comment">// The precise versions use KahanSum and are about two times slower.</span></div>
|
|
<div class="line"><a name="l00046"></a><span class="lineno"> 46</span> <span class="keyword">template</span> <<span class="keyword">class</span> DenseRowOrColumn1, <span class="keyword">class</span> DenseRowOrColumn2></div>
|
|
<div class="line"><a name="l00047"></a><span class="lineno"><a class="line" href="namespaceoperations__research_1_1glop.html#a3d9874910f9632900fe8677463b1c0be"> 47</a></span> <a class="code" href="namespaceoperations__research_1_1glop.html#a733947145e3e1631165b618b05c9ccb7">Fractional</a> <a class="code" href="namespaceoperations__research_1_1glop.html#a3d9874910f9632900fe8677463b1c0be">ScalarProduct</a>(<span class="keyword">const</span> DenseRowOrColumn1& u,</div>
|
|
<div class="line"><a name="l00048"></a><span class="lineno"> 48</span>  <span class="keyword">const</span> DenseRowOrColumn2& v) {</div>
|
|
<div class="line"><a name="l00049"></a><span class="lineno"> 49</span>  <a class="code" href="base_2logging_8h.html#ae89df3243bbb8341130c7b3f44145ea0">DCHECK_EQ</a>(u.size().value(), v.size().value());</div>
|
|
<div class="line"><a name="l00050"></a><span class="lineno"> 50</span>  <a class="code" href="namespaceoperations__research_1_1glop.html#a733947145e3e1631165b618b05c9ccb7">Fractional</a> sum(0.0);</div>
|
|
<div class="line"><a name="l00051"></a><span class="lineno"> 51</span>  <span class="keyword">typename</span> DenseRowOrColumn1::IndexType i(0);</div>
|
|
<div class="line"><a name="l00052"></a><span class="lineno"> 52</span>  <span class="keyword">typename</span> DenseRowOrColumn2::IndexType j(0);</div>
|
|
<div class="line"><a name="l00053"></a><span class="lineno"> 53</span>  <span class="keyword">const</span> <span class="keywordtype">size_t</span> num_blocks = u.size().value() / 4;</div>
|
|
<div class="line"><a name="l00054"></a><span class="lineno"> 54</span>  <span class="keywordflow">for</span> (<span class="keywordtype">size_t</span> block = 0; block < num_blocks; ++block) {</div>
|
|
<div class="line"><a name="l00055"></a><span class="lineno"> 55</span>  <span class="comment">// Computing the sum of 4 elements at once may allow the compiler to</span></div>
|
|
<div class="line"><a name="l00056"></a><span class="lineno"> 56</span>  <span class="comment">// generate more efficient code, e.g. using SIMD and checking the loop</span></div>
|
|
<div class="line"><a name="l00057"></a><span class="lineno"> 57</span>  <span class="comment">// condition much less frequently.</span></div>
|
|
<div class="line"><a name="l00058"></a><span class="lineno"> 58</span>  <span class="comment">//</span></div>
|
|
<div class="line"><a name="l00059"></a><span class="lineno"> 59</span>  <span class="comment">// This produces different results from the case where each multiplication</span></div>
|
|
<div class="line"><a name="l00060"></a><span class="lineno"> 60</span>  <span class="comment">// is added to sum separately. An extreme example of this can be derived</span></div>
|
|
<div class="line"><a name="l00061"></a><span class="lineno"> 61</span>  <span class="comment">// using the fact that 1e11 + 2e-6 == 1e11, but 1e11 + 8e-6 > 1e11.</span></div>
|
|
<div class="line"><a name="l00062"></a><span class="lineno"> 62</span>  <span class="comment">//</span></div>
|
|
<div class="line"><a name="l00063"></a><span class="lineno"> 63</span>  <span class="comment">// While the results are different, they aren't necessarily better or worse.</span></div>
|
|
<div class="line"><a name="l00064"></a><span class="lineno"> 64</span>  <span class="comment">// Typically, sum will be of larger magnitude than any individual</span></div>
|
|
<div class="line"><a name="l00065"></a><span class="lineno"> 65</span>  <span class="comment">// multiplication, so one might expect, in practice, this method to yield</span></div>
|
|
<div class="line"><a name="l00066"></a><span class="lineno"> 66</span>  <span class="comment">// more accurate results. However, if accuracy is vital, use the precise</span></div>
|
|
<div class="line"><a name="l00067"></a><span class="lineno"> 67</span>  <span class="comment">// version.</span></div>
|
|
<div class="line"><a name="l00068"></a><span class="lineno"> 68</span>  sum += (u[i++] * v[j++]) + (u[i++] * v[j++]) + (u[i++] * v[j++]) +</div>
|
|
<div class="line"><a name="l00069"></a><span class="lineno"> 69</span>  (u[i++] * v[j++]);</div>
|
|
<div class="line"><a name="l00070"></a><span class="lineno"> 70</span>  }</div>
|
|
<div class="line"><a name="l00071"></a><span class="lineno"> 71</span>  <span class="keywordflow">while</span> (i < u.size()) {</div>
|
|
<div class="line"><a name="l00072"></a><span class="lineno"> 72</span>  sum += u[i++] * v[j++];</div>
|
|
<div class="line"><a name="l00073"></a><span class="lineno"> 73</span>  }</div>
|
|
<div class="line"><a name="l00074"></a><span class="lineno"> 74</span>  <span class="keywordflow">return</span> sum;</div>
|
|
<div class="line"><a name="l00075"></a><span class="lineno"> 75</span> }</div>
|
|
<div class="line"><a name="l00076"></a><span class="lineno"> 76</span>  </div>
|
|
<div class="line"><a name="l00077"></a><span class="lineno"> 77</span> <span class="comment">// Note: This version is heavily used in the pricing.</span></div>
|
|
<div class="line"><a name="l00078"></a><span class="lineno"> 78</span> <span class="comment">// TODO(user): Optimize this more (SSE or unroll with two sums). Another</span></div>
|
|
<div class="line"><a name="l00079"></a><span class="lineno"> 79</span> <span class="comment">// option is to skip the u[col] that are 0.0 rather than fetching the coeff</span></div>
|
|
<div class="line"><a name="l00080"></a><span class="lineno"> 80</span> <span class="comment">// and doing a Fractional multiplication.</span></div>
|
|
<div class="line"><a name="l00081"></a><span class="lineno"> 81</span> <span class="keyword">template</span> <<span class="keyword">class</span> DenseRowOrColumn></div>
|
|
<div class="line"><a name="l00082"></a><span class="lineno"><a class="line" href="namespaceoperations__research_1_1glop.html#a46935a48ea28abfec1c2569fe4aea70f"> 82</a></span> <a class="code" href="namespaceoperations__research_1_1glop.html#a733947145e3e1631165b618b05c9ccb7">Fractional</a> <a class="code" href="namespaceoperations__research_1_1glop.html#a3d9874910f9632900fe8677463b1c0be">ScalarProduct</a>(<span class="keyword">const</span> DenseRowOrColumn& u, <span class="keyword">const</span> <a class="code" href="classoperations__research_1_1glop_1_1_sparse_column.html">SparseColumn</a>& v) {</div>
|
|
<div class="line"><a name="l00083"></a><span class="lineno"> 83</span>  <a class="code" href="namespaceoperations__research_1_1glop.html#a733947145e3e1631165b618b05c9ccb7">Fractional</a> sum(0.0);</div>
|
|
<div class="line"><a name="l00084"></a><span class="lineno"> 84</span>  <span class="keywordflow">for</span> (<span class="keyword">const</span> <a class="code" href="classoperations__research_1_1glop_1_1_sparse_vector.html#ab38326ea6cb6187267665dd8b2748f3d">SparseColumn::Entry</a> e : v) {</div>
|
|
<div class="line"><a name="l00085"></a><span class="lineno"> 85</span>  sum += u[<span class="keyword">typename</span> DenseRowOrColumn::IndexType(e.row().value())] *</div>
|
|
<div class="line"><a name="l00086"></a><span class="lineno"> 86</span>  e.coefficient();</div>
|
|
<div class="line"><a name="l00087"></a><span class="lineno"> 87</span>  }</div>
|
|
<div class="line"><a name="l00088"></a><span class="lineno"> 88</span>  <span class="keywordflow">return</span> sum;</div>
|
|
<div class="line"><a name="l00089"></a><span class="lineno"> 89</span> }</div>
|
|
<div class="line"><a name="l00090"></a><span class="lineno"> 90</span>  </div>
|
|
<div class="line"><a name="l00091"></a><span class="lineno"> 91</span> <span class="keyword">template</span> <<span class="keyword">class</span> DenseRowOrColumn, <span class="keyword">class</span> DenseRowOrColumn2></div>
|
|
<div class="line"><a name="l00092"></a><span class="lineno"><a class="line" href="namespaceoperations__research_1_1glop.html#a46fb729c0be27d1b97db15e0ce9c6067"> 92</a></span> <a class="code" href="namespaceoperations__research_1_1glop.html#a733947145e3e1631165b618b05c9ccb7">Fractional</a> <a class="code" href="namespaceoperations__research_1_1glop.html#a46fb729c0be27d1b97db15e0ce9c6067">PreciseScalarProduct</a>(<span class="keyword">const</span> DenseRowOrColumn& u,</div>
|
|
<div class="line"><a name="l00093"></a><span class="lineno"> 93</span>  <span class="keyword">const</span> DenseRowOrColumn2& v) {</div>
|
|
<div class="line"><a name="l00094"></a><span class="lineno"> 94</span>  <a class="code" href="base_2logging_8h.html#ae89df3243bbb8341130c7b3f44145ea0">DCHECK_EQ</a>(u.size().value(), v.size().value());</div>
|
|
<div class="line"><a name="l00095"></a><span class="lineno"> 95</span>  <a class="code" href="classoperations__research_1_1_accurate_sum.html">KahanSum</a> sum;</div>
|
|
<div class="line"><a name="l00096"></a><span class="lineno"> 96</span>  <span class="keywordflow">for</span> (<span class="keyword">typename</span> DenseRowOrColumn::IndexType i(0); i < u.size(); ++i) {</div>
|
|
<div class="line"><a name="l00097"></a><span class="lineno"> 97</span>  sum.<a class="code" href="classoperations__research_1_1_accurate_sum.html#acc8fad1f98b885fa7d0504f4dfc71f08">Add</a>(u[i] * v[<span class="keyword">typename</span> DenseRowOrColumn2::IndexType(i.value())]);</div>
|
|
<div class="line"><a name="l00098"></a><span class="lineno"> 98</span>  }</div>
|
|
<div class="line"><a name="l00099"></a><span class="lineno"> 99</span>  <span class="keywordflow">return</span> sum.<a class="code" href="classoperations__research_1_1_accurate_sum.html#a176a3e919acd979b67cea1ede094cdaa">Value</a>();</div>
|
|
<div class="line"><a name="l00100"></a><span class="lineno"> 100</span> }</div>
|
|
<div class="line"><a name="l00101"></a><span class="lineno"> 101</span>  </div>
|
|
<div class="line"><a name="l00102"></a><span class="lineno"> 102</span> <span class="keyword">template</span> <<span class="keyword">class</span> DenseRowOrColumn></div>
|
|
<div class="line"><a name="l00103"></a><span class="lineno"><a class="line" href="namespaceoperations__research_1_1glop.html#ab179616817239f2167055368df1e9f66"> 103</a></span> <a class="code" href="namespaceoperations__research_1_1glop.html#a733947145e3e1631165b618b05c9ccb7">Fractional</a> <a class="code" href="namespaceoperations__research_1_1glop.html#a46fb729c0be27d1b97db15e0ce9c6067">PreciseScalarProduct</a>(<span class="keyword">const</span> DenseRowOrColumn& u,</div>
|
|
<div class="line"><a name="l00104"></a><span class="lineno"> 104</span>  <span class="keyword">const</span> <a class="code" href="classoperations__research_1_1glop_1_1_sparse_column.html">SparseColumn</a>& v) {</div>
|
|
<div class="line"><a name="l00105"></a><span class="lineno"> 105</span>  <a class="code" href="classoperations__research_1_1_accurate_sum.html">KahanSum</a> sum;</div>
|
|
<div class="line"><a name="l00106"></a><span class="lineno"> 106</span>  <span class="keywordflow">for</span> (<span class="keyword">const</span> <a class="code" href="classoperations__research_1_1glop_1_1_sparse_vector.html#ab38326ea6cb6187267665dd8b2748f3d">SparseColumn::Entry</a> e : v) {</div>
|
|
<div class="line"><a name="l00107"></a><span class="lineno"> 107</span>  sum.<a class="code" href="classoperations__research_1_1_accurate_sum.html#acc8fad1f98b885fa7d0504f4dfc71f08">Add</a>(u[<span class="keyword">typename</span> DenseRowOrColumn::IndexType(e.row().value())] *</div>
|
|
<div class="line"><a name="l00108"></a><span class="lineno"> 108</span>  e.coefficient());</div>
|
|
<div class="line"><a name="l00109"></a><span class="lineno"> 109</span>  }</div>
|
|
<div class="line"><a name="l00110"></a><span class="lineno"> 110</span>  <span class="keywordflow">return</span> sum.<a class="code" href="classoperations__research_1_1_accurate_sum.html#a176a3e919acd979b67cea1ede094cdaa">Value</a>();</div>
|
|
<div class="line"><a name="l00111"></a><span class="lineno"> 111</span> }</div>
|
|
<div class="line"><a name="l00112"></a><span class="lineno"> 112</span>  </div>
|
|
<div class="line"><a name="l00113"></a><span class="lineno"> 113</span> <span class="keyword">template</span> <<span class="keyword">class</span> DenseRowOrColumn></div>
|
|
<div class="line"><a name="l00114"></a><span class="lineno"><a class="line" href="namespaceoperations__research_1_1glop.html#a434f75c61605b1ede60e834ee196660d"> 114</a></span> <a class="code" href="namespaceoperations__research_1_1glop.html#a733947145e3e1631165b618b05c9ccb7">Fractional</a> <a class="code" href="namespaceoperations__research_1_1glop.html#a46fb729c0be27d1b97db15e0ce9c6067">PreciseScalarProduct</a>(<span class="keyword">const</span> DenseRowOrColumn& u,</div>
|
|
<div class="line"><a name="l00115"></a><span class="lineno"> 115</span>  <span class="keyword">const</span> <a class="code" href="structoperations__research_1_1glop_1_1_scattered_column.html">ScatteredColumn</a>& v) {</div>
|
|
<div class="line"><a name="l00116"></a><span class="lineno"> 116</span>  <a class="code" href="base_2logging_8h.html#ae89df3243bbb8341130c7b3f44145ea0">DCHECK_EQ</a>(u.size().value(), v.<a class="code" href="structoperations__research_1_1glop_1_1_scattered_vector.html#ae0473e18a367af671dc3f08063c80da4">values</a>.<a class="code" href="classoperations__research_1_1glop_1_1_strict_i_t_i_vector.html#a967a5c081ad4195a30c78dc2c0bcabf5">size</a>().value());</div>
|
|
<div class="line"><a name="l00117"></a><span class="lineno"> 117</span>  <span class="keywordflow">if</span> (v.<a class="code" href="structoperations__research_1_1glop_1_1_scattered_vector.html#a73934bd40690bdc3cbe36d6fbb0ecad5">ShouldUseDenseIteration</a>()) {</div>
|
|
<div class="line"><a name="l00118"></a><span class="lineno"> 118</span>  <span class="keywordflow">return</span> <a class="code" href="namespaceoperations__research_1_1glop.html#a46fb729c0be27d1b97db15e0ce9c6067">PreciseScalarProduct</a>(u, v.<a class="code" href="structoperations__research_1_1glop_1_1_scattered_vector.html#ae0473e18a367af671dc3f08063c80da4">values</a>);</div>
|
|
<div class="line"><a name="l00119"></a><span class="lineno"> 119</span>  }</div>
|
|
<div class="line"><a name="l00120"></a><span class="lineno"> 120</span>  <a class="code" href="classoperations__research_1_1_accurate_sum.html">KahanSum</a> sum;</div>
|
|
<div class="line"><a name="l00121"></a><span class="lineno"> 121</span>  <span class="keywordflow">for</span> (<span class="keyword">const</span> <span class="keyword">auto</span> e : v) {</div>
|
|
<div class="line"><a name="l00122"></a><span class="lineno"> 122</span>  sum.<a class="code" href="classoperations__research_1_1_accurate_sum.html#acc8fad1f98b885fa7d0504f4dfc71f08">Add</a>(u[<span class="keyword">typename</span> DenseRowOrColumn::IndexType(e.row().value())] *</div>
|
|
<div class="line"><a name="l00123"></a><span class="lineno"> 123</span>  e.coefficient());</div>
|
|
<div class="line"><a name="l00124"></a><span class="lineno"> 124</span>  }</div>
|
|
<div class="line"><a name="l00125"></a><span class="lineno"> 125</span>  <span class="keywordflow">return</span> sum.<a class="code" href="classoperations__research_1_1_accurate_sum.html#a176a3e919acd979b67cea1ede094cdaa">Value</a>();</div>
|
|
<div class="line"><a name="l00126"></a><span class="lineno"> 126</span> }</div>
|
|
<div class="line"><a name="l00127"></a><span class="lineno"> 127</span>  </div>
|
|
<div class="line"><a name="l00128"></a><span class="lineno"> 128</span> <span class="comment">// Computes a scalar product for entries with index not greater than max_index.</span></div>
|
|
<div class="line"><a name="l00129"></a><span class="lineno"> 129</span> <span class="keyword">template</span> <<span class="keyword">class</span> DenseRowOrColumn></div>
|
|
<div class="line"><a name="l00130"></a><span class="lineno"><a class="line" href="namespaceoperations__research_1_1glop.html#ad2cb64ef0c88350c1ad4c731e0b7e183"> 130</a></span> <a class="code" href="namespaceoperations__research_1_1glop.html#a733947145e3e1631165b618b05c9ccb7">Fractional</a> <a class="code" href="namespaceoperations__research_1_1glop.html#ad2cb64ef0c88350c1ad4c731e0b7e183">PartialScalarProduct</a>(<span class="keyword">const</span> DenseRowOrColumn& u,</div>
|
|
<div class="line"><a name="l00131"></a><span class="lineno"> 131</span>  <span class="keyword">const</span> <a class="code" href="classoperations__research_1_1glop_1_1_sparse_column.html">SparseColumn</a>& v, <span class="keywordtype">int</span> max_index) {</div>
|
|
<div class="line"><a name="l00132"></a><span class="lineno"> 132</span>  <a class="code" href="namespaceoperations__research_1_1glop.html#a733947145e3e1631165b618b05c9ccb7">Fractional</a> sum(0.0);</div>
|
|
<div class="line"><a name="l00133"></a><span class="lineno"> 133</span>  <span class="keywordflow">for</span> (<span class="keyword">const</span> <a class="code" href="classoperations__research_1_1glop_1_1_sparse_vector.html#ab38326ea6cb6187267665dd8b2748f3d">SparseColumn::Entry</a> e : v) {</div>
|
|
<div class="line"><a name="l00134"></a><span class="lineno"> 134</span>  <span class="keywordflow">if</span> (e.row().value() >= max_index) {</div>
|
|
<div class="line"><a name="l00135"></a><span class="lineno"> 135</span>  <span class="keywordflow">return</span> sum;</div>
|
|
<div class="line"><a name="l00136"></a><span class="lineno"> 136</span>  }</div>
|
|
<div class="line"><a name="l00137"></a><span class="lineno"> 137</span>  sum += u[<span class="keyword">typename</span> DenseRowOrColumn::IndexType(e.row().value())] *</div>
|
|
<div class="line"><a name="l00138"></a><span class="lineno"> 138</span>  e.coefficient();</div>
|
|
<div class="line"><a name="l00139"></a><span class="lineno"> 139</span>  }</div>
|
|
<div class="line"><a name="l00140"></a><span class="lineno"> 140</span>  <span class="keywordflow">return</span> sum;</div>
|
|
<div class="line"><a name="l00141"></a><span class="lineno"> 141</span> }</div>
|
|
<div class="line"><a name="l00142"></a><span class="lineno"> 142</span>  </div>
|
|
<div class="line"><a name="l00143"></a><span class="lineno"> 143</span> <span class="comment">// Returns the norm^2 (sum of the square of the entries) of the given column.</span></div>
|
|
<div class="line"><a name="l00144"></a><span class="lineno"> 144</span> <span class="comment">// The precise version uses KahanSum and are about two times slower.</span></div>
|
|
<div class="line"><a name="l00145"></a><span class="lineno"> 145</span> <a class="code" href="namespaceoperations__research_1_1glop.html#a733947145e3e1631165b618b05c9ccb7">Fractional</a> <a class="code" href="namespaceoperations__research_1_1glop.html#a2d53948bf5e999d006e781105aa8bc77">SquaredNorm</a>(<span class="keyword">const</span> SparseColumn& v);</div>
|
|
<div class="line"><a name="l00146"></a><span class="lineno"> 146</span> <a class="code" href="namespaceoperations__research_1_1glop.html#a733947145e3e1631165b618b05c9ccb7">Fractional</a> <a class="code" href="namespaceoperations__research_1_1glop.html#a2d53948bf5e999d006e781105aa8bc77">SquaredNorm</a>(<span class="keyword">const</span> <a class="code" href="namespaceoperations__research_1_1glop.html#ae5fa9e57c2d31b294195ed8a9e17bfac">DenseColumn</a>& column);</div>
|
|
<div class="line"><a name="l00147"></a><span class="lineno"> 147</span> <a class="code" href="namespaceoperations__research_1_1glop.html#a733947145e3e1631165b618b05c9ccb7">Fractional</a> <a class="code" href="namespaceoperations__research_1_1glop.html#a2d53948bf5e999d006e781105aa8bc77">SquaredNorm</a>(<span class="keyword">const</span> ColumnView& v);</div>
|
|
<div class="line"><a name="l00148"></a><span class="lineno"> 148</span> <a class="code" href="namespaceoperations__research_1_1glop.html#a733947145e3e1631165b618b05c9ccb7">Fractional</a> <a class="code" href="namespaceoperations__research_1_1glop.html#a1faa927dd93b43b3dea3eb2a993e30a1">PreciseSquaredNorm</a>(<span class="keyword">const</span> SparseColumn& v);</div>
|
|
<div class="line"><a name="l00149"></a><span class="lineno"> 149</span> <a class="code" href="namespaceoperations__research_1_1glop.html#a733947145e3e1631165b618b05c9ccb7">Fractional</a> <a class="code" href="namespaceoperations__research_1_1glop.html#a1faa927dd93b43b3dea3eb2a993e30a1">PreciseSquaredNorm</a>(<span class="keyword">const</span> <a class="code" href="namespaceoperations__research_1_1glop.html#ae5fa9e57c2d31b294195ed8a9e17bfac">DenseColumn</a>& column);</div>
|
|
<div class="line"><a name="l00150"></a><span class="lineno"> 150</span> <a class="code" href="namespaceoperations__research_1_1glop.html#a733947145e3e1631165b618b05c9ccb7">Fractional</a> <a class="code" href="namespaceoperations__research_1_1glop.html#a1faa927dd93b43b3dea3eb2a993e30a1">PreciseSquaredNorm</a>(<span class="keyword">const</span> ScatteredColumn& v);</div>
|
|
<div class="line"><a name="l00151"></a><span class="lineno"> 151</span>  </div>
|
|
<div class="line"><a name="l00152"></a><span class="lineno"> 152</span> <span class="comment">// Returns the maximum of the |coefficients| of 'v'.</span></div>
|
|
<div class="line"><a name="l00153"></a><span class="lineno"> 153</span> <a class="code" href="namespaceoperations__research_1_1glop.html#a733947145e3e1631165b618b05c9ccb7">Fractional</a> <a class="code" href="namespaceoperations__research_1_1glop.html#a2ad7ea612bc859f5b637d5029fb875fb">InfinityNorm</a>(<span class="keyword">const</span> <a class="code" href="namespaceoperations__research_1_1glop.html#ae5fa9e57c2d31b294195ed8a9e17bfac">DenseColumn</a>& v);</div>
|
|
<div class="line"><a name="l00154"></a><span class="lineno"> 154</span> <a class="code" href="namespaceoperations__research_1_1glop.html#a733947145e3e1631165b618b05c9ccb7">Fractional</a> <a class="code" href="namespaceoperations__research_1_1glop.html#a2ad7ea612bc859f5b637d5029fb875fb">InfinityNorm</a>(<span class="keyword">const</span> SparseColumn& v);</div>
|
|
<div class="line"><a name="l00155"></a><span class="lineno"> 155</span> <a class="code" href="namespaceoperations__research_1_1glop.html#a733947145e3e1631165b618b05c9ccb7">Fractional</a> <a class="code" href="namespaceoperations__research_1_1glop.html#a2ad7ea612bc859f5b637d5029fb875fb">InfinityNorm</a>(<span class="keyword">const</span> ColumnView& v);</div>
|
|
<div class="line"><a name="l00156"></a><span class="lineno"> 156</span>  </div>
|
|
<div class="line"><a name="l00157"></a><span class="lineno"> 157</span> <span class="comment">// Returns the fraction of non-zero entries of the given row.</span></div>
|
|
<div class="line"><a name="l00158"></a><span class="lineno"> 158</span> <span class="comment">//</span></div>
|
|
<div class="line"><a name="l00159"></a><span class="lineno"> 159</span> <span class="comment">// TODO(user): Take a Scattered row/col instead. This is only used to report</span></div>
|
|
<div class="line"><a name="l00160"></a><span class="lineno"> 160</span> <span class="comment">// stats, but we should still have a sparse version to do it faster.</span></div>
|
|
<div class="line"><a name="l00161"></a><span class="lineno"> 161</span> <span class="keywordtype">double</span> <a class="code" href="namespaceoperations__research_1_1glop.html#a6b36bbbcd695191600dab9b88487ff18">Density</a>(<span class="keyword">const</span> <a class="code" href="namespaceoperations__research_1_1glop.html#a4dc8ae0f97c4110f3cc770715b2bbd12">DenseRow</a>& <a class="code" href="markowitz_8cc.html#aea35f36ba98d5bbd8d033382f50c9e52">row</a>);</div>
|
|
<div class="line"><a name="l00162"></a><span class="lineno"> 162</span>  </div>
|
|
<div class="line"><a name="l00163"></a><span class="lineno"> 163</span> <span class="comment">// Sets to 0.0 all entries of the given row whose fabs() is lower than the given</span></div>
|
|
<div class="line"><a name="l00164"></a><span class="lineno"> 164</span> <span class="comment">// threshold.</span></div>
|
|
<div class="line"><a name="l00165"></a><span class="lineno"> 165</span> <span class="keywordtype">void</span> <a class="code" href="namespaceoperations__research_1_1glop.html#a5e79e30b7239adc4fb2a27778335bca0">RemoveNearZeroEntries</a>(<a class="code" href="namespaceoperations__research_1_1glop.html#a733947145e3e1631165b618b05c9ccb7">Fractional</a> threshold, <a class="code" href="namespaceoperations__research_1_1glop.html#a4dc8ae0f97c4110f3cc770715b2bbd12">DenseRow</a>* <a class="code" href="markowitz_8cc.html#aea35f36ba98d5bbd8d033382f50c9e52">row</a>);</div>
|
|
<div class="line"><a name="l00166"></a><span class="lineno"> 166</span> <span class="keywordtype">void</span> <a class="code" href="namespaceoperations__research_1_1glop.html#a5e79e30b7239adc4fb2a27778335bca0">RemoveNearZeroEntries</a>(<a class="code" href="namespaceoperations__research_1_1glop.html#a733947145e3e1631165b618b05c9ccb7">Fractional</a> threshold, <a class="code" href="namespaceoperations__research_1_1glop.html#ae5fa9e57c2d31b294195ed8a9e17bfac">DenseColumn</a>* column);</div>
|
|
<div class="line"><a name="l00167"></a><span class="lineno"> 167</span>  </div>
|
|
<div class="line"><a name="l00168"></a><span class="lineno"> 168</span> <span class="comment">// Transposition functions implemented below with a cast so it should actually</span></div>
|
|
<div class="line"><a name="l00169"></a><span class="lineno"> 169</span> <span class="comment">// have no complexity cost.</span></div>
|
|
<div class="line"><a name="l00170"></a><span class="lineno"> 170</span> <span class="keyword">const</span> <a class="code" href="namespaceoperations__research_1_1glop.html#a4dc8ae0f97c4110f3cc770715b2bbd12">DenseRow</a>& <a class="code" href="namespaceoperations__research_1_1glop.html#aaa803ce9366dca251925e0bdde517430">Transpose</a>(<span class="keyword">const</span> <a class="code" href="namespaceoperations__research_1_1glop.html#ae5fa9e57c2d31b294195ed8a9e17bfac">DenseColumn</a>& <a class="code" href="markowitz_8cc.html#aa9d6c98fdf8d89b0e2321fda02adc82c">col</a>);</div>
|
|
<div class="line"><a name="l00171"></a><span class="lineno"> 171</span> <span class="keyword">const</span> <a class="code" href="namespaceoperations__research_1_1glop.html#ae5fa9e57c2d31b294195ed8a9e17bfac">DenseColumn</a>& <a class="code" href="namespaceoperations__research_1_1glop.html#aaa803ce9366dca251925e0bdde517430">Transpose</a>(<span class="keyword">const</span> <a class="code" href="namespaceoperations__research_1_1glop.html#a4dc8ae0f97c4110f3cc770715b2bbd12">DenseRow</a>& <a class="code" href="markowitz_8cc.html#aea35f36ba98d5bbd8d033382f50c9e52">row</a>);</div>
|
|
<div class="line"><a name="l00172"></a><span class="lineno"> 172</span>  </div>
|
|
<div class="line"><a name="l00173"></a><span class="lineno"> 173</span> <span class="comment">// Returns the maximum of the |coefficients| of the given column restricted</span></div>
|
|
<div class="line"><a name="l00174"></a><span class="lineno"> 174</span> <span class="comment">// to the rows_to_consider. Also returns the first RowIndex 'row' that attains</span></div>
|
|
<div class="line"><a name="l00175"></a><span class="lineno"> 175</span> <span class="comment">// this maximum. If the maximum is 0.0, then row_index is left untouched.</span></div>
|
|
<div class="line"><a name="l00176"></a><span class="lineno"> 176</span> <a class="code" href="namespaceoperations__research_1_1glop.html#a733947145e3e1631165b618b05c9ccb7">Fractional</a> <a class="code" href="namespaceoperations__research_1_1glop.html#ad8019bac1bde0ead6ff32980cd5bff52">RestrictedInfinityNorm</a>(<span class="keyword">const</span> ColumnView& column,</div>
|
|
<div class="line"><a name="l00177"></a><span class="lineno"> 177</span>  <span class="keyword">const</span> <a class="code" href="namespaceoperations__research_1_1glop.html#aec7177dd4fab33ada14a7689223fe89f">DenseBooleanColumn</a>& rows_to_consider,</div>
|
|
<div class="line"><a name="l00178"></a><span class="lineno"> 178</span>  RowIndex* row_index);</div>
|
|
<div class="line"><a name="l00179"></a><span class="lineno"> 179</span>  </div>
|
|
<div class="line"><a name="l00180"></a><span class="lineno"> 180</span> <span class="comment">// Sets to false the entry b[row] if column[row] is non null.</span></div>
|
|
<div class="line"><a name="l00181"></a><span class="lineno"> 181</span> <span class="comment">// Note that if 'b' was true only on the non-zero position of column, this can</span></div>
|
|
<div class="line"><a name="l00182"></a><span class="lineno"> 182</span> <span class="comment">// be used as a fast way to clear 'b'.</span></div>
|
|
<div class="line"><a name="l00183"></a><span class="lineno"> 183</span> <span class="keywordtype">void</span> <a class="code" href="namespaceoperations__research_1_1glop.html#a82b6c13b99cf0a0c7c92f4f1e44eda29">SetSupportToFalse</a>(<span class="keyword">const</span> ColumnView& column, <a class="code" href="namespaceoperations__research_1_1glop.html#aec7177dd4fab33ada14a7689223fe89f">DenseBooleanColumn</a>* <a class="code" href="constraint__solver_2table_8cc.html#a344010e26426d6a13411648d988bc9b6">b</a>);</div>
|
|
<div class="line"><a name="l00184"></a><span class="lineno"> 184</span>  </div>
|
|
<div class="line"><a name="l00185"></a><span class="lineno"> 185</span> <span class="comment">// Returns true iff for all 'row' we have '|column[row]| <= radius[row]'.</span></div>
|
|
<div class="line"><a name="l00186"></a><span class="lineno"> 186</span> <span class="keywordtype">bool</span> <a class="code" href="namespaceoperations__research_1_1glop.html#a9a9b90bb0105347953a41ae1d6f4dce3">IsDominated</a>(<span class="keyword">const</span> ColumnView& column, <span class="keyword">const</span> <a class="code" href="namespaceoperations__research_1_1glop.html#ae5fa9e57c2d31b294195ed8a9e17bfac">DenseColumn</a>& radius);</div>
|
|
<div class="line"><a name="l00187"></a><span class="lineno"> 187</span>  </div>
|
|
<div class="line"><a name="l00188"></a><span class="lineno"> 188</span> <span class="comment">// This cast based implementation should be safe, as long as DenseRow and</span></div>
|
|
<div class="line"><a name="l00189"></a><span class="lineno"> 189</span> <span class="comment">// DenseColumn are implemented by the same underlying type.</span></div>
|
|
<div class="line"><a name="l00190"></a><span class="lineno"> 190</span> <span class="comment">// We still do some DCHECK to be sure it works as expected in addition to the</span></div>
|
|
<div class="line"><a name="l00191"></a><span class="lineno"> 191</span> <span class="comment">// unit tests.</span></div>
|
|
<div class="line"><a name="l00192"></a><span class="lineno"><a class="line" href="namespaceoperations__research_1_1glop.html#aaa803ce9366dca251925e0bdde517430"> 192</a></span> <span class="keyword">inline</span> <span class="keyword">const</span> <a class="code" href="classoperations__research_1_1glop_1_1_strict_i_t_i_vector.html">DenseRow</a>& <a class="code" href="namespaceoperations__research_1_1glop.html#aaa803ce9366dca251925e0bdde517430">Transpose</a>(<span class="keyword">const</span> <a class="code" href="classoperations__research_1_1glop_1_1_strict_i_t_i_vector.html">DenseColumn</a>& <a class="code" href="markowitz_8cc.html#aa9d6c98fdf8d89b0e2321fda02adc82c">col</a>) {</div>
|
|
<div class="line"><a name="l00193"></a><span class="lineno"> 193</span>  <span class="keyword">const</span> <a class="code" href="classoperations__research_1_1glop_1_1_strict_i_t_i_vector.html">DenseRow</a>& <a class="code" href="markowitz_8cc.html#aea35f36ba98d5bbd8d033382f50c9e52">row</a> = <span class="keyword">reinterpret_cast<</span><span class="keyword">const </span><a class="code" href="classoperations__research_1_1glop_1_1_strict_i_t_i_vector.html">DenseRow</a>&<span class="keyword">></span>(<a class="code" href="markowitz_8cc.html#aa9d6c98fdf8d89b0e2321fda02adc82c">col</a>);</div>
|
|
<div class="line"><a name="l00194"></a><span class="lineno"> 194</span>  <a class="code" href="base_2logging_8h.html#ae89df3243bbb8341130c7b3f44145ea0">DCHECK_EQ</a>(<a class="code" href="markowitz_8cc.html#aa9d6c98fdf8d89b0e2321fda02adc82c">col</a>.size(), <a class="code" href="namespaceoperations__research_1_1glop.html#ab65a327cfc2a74c15fa26b91f19acc64">ColToRowIndex</a>(<a class="code" href="markowitz_8cc.html#aea35f36ba98d5bbd8d033382f50c9e52">row</a>.size()));</div>
|
|
<div class="line"><a name="l00195"></a><span class="lineno"> 195</span>  <a class="code" href="base_2logging_8h.html#ae17f8119c108cf3070bad3449c7e0006">DCHECK</a>(<a class="code" href="markowitz_8cc.html#aa9d6c98fdf8d89b0e2321fda02adc82c">col</a>.empty() || (&(<a class="code" href="markowitz_8cc.html#aa9d6c98fdf8d89b0e2321fda02adc82c">col</a>[RowIndex(0)]) == &(<a class="code" href="markowitz_8cc.html#aea35f36ba98d5bbd8d033382f50c9e52">row</a>[ColIndex(0)])));</div>
|
|
<div class="line"><a name="l00196"></a><span class="lineno"> 196</span>  <span class="keywordflow">return</span> <a class="code" href="markowitz_8cc.html#aea35f36ba98d5bbd8d033382f50c9e52">row</a>;</div>
|
|
<div class="line"><a name="l00197"></a><span class="lineno"> 197</span> }</div>
|
|
<div class="line"><a name="l00198"></a><span class="lineno"> 198</span>  </div>
|
|
<div class="line"><a name="l00199"></a><span class="lineno"> 199</span> <span class="comment">// Similar comment as the other Transpose() implementation above.</span></div>
|
|
<div class="line"><a name="l00200"></a><span class="lineno"><a class="line" href="namespaceoperations__research_1_1glop.html#a96eb7e615016e66686739537ebf5e1a4"> 200</a></span> <span class="keyword">inline</span> <span class="keyword">const</span> <a class="code" href="classoperations__research_1_1glop_1_1_strict_i_t_i_vector.html">DenseColumn</a>& <a class="code" href="namespaceoperations__research_1_1glop.html#aaa803ce9366dca251925e0bdde517430">Transpose</a>(<span class="keyword">const</span> <a class="code" href="classoperations__research_1_1glop_1_1_strict_i_t_i_vector.html">DenseRow</a>& <a class="code" href="markowitz_8cc.html#aea35f36ba98d5bbd8d033382f50c9e52">row</a>) {</div>
|
|
<div class="line"><a name="l00201"></a><span class="lineno"> 201</span>  <span class="keyword">const</span> <a class="code" href="classoperations__research_1_1glop_1_1_strict_i_t_i_vector.html">DenseColumn</a>& <a class="code" href="markowitz_8cc.html#aa9d6c98fdf8d89b0e2321fda02adc82c">col</a> = <span class="keyword">reinterpret_cast<</span><span class="keyword">const </span><a class="code" href="classoperations__research_1_1glop_1_1_strict_i_t_i_vector.html">DenseColumn</a>&<span class="keyword">></span>(<a class="code" href="markowitz_8cc.html#aea35f36ba98d5bbd8d033382f50c9e52">row</a>);</div>
|
|
<div class="line"><a name="l00202"></a><span class="lineno"> 202</span>  <a class="code" href="base_2logging_8h.html#ae89df3243bbb8341130c7b3f44145ea0">DCHECK_EQ</a>(<a class="code" href="markowitz_8cc.html#aa9d6c98fdf8d89b0e2321fda02adc82c">col</a>.size(), <a class="code" href="namespaceoperations__research_1_1glop.html#ab65a327cfc2a74c15fa26b91f19acc64">ColToRowIndex</a>(<a class="code" href="markowitz_8cc.html#aea35f36ba98d5bbd8d033382f50c9e52">row</a>.size()));</div>
|
|
<div class="line"><a name="l00203"></a><span class="lineno"> 203</span>  <a class="code" href="base_2logging_8h.html#ae17f8119c108cf3070bad3449c7e0006">DCHECK</a>(<a class="code" href="markowitz_8cc.html#aa9d6c98fdf8d89b0e2321fda02adc82c">col</a>.empty() || (&(<a class="code" href="markowitz_8cc.html#aa9d6c98fdf8d89b0e2321fda02adc82c">col</a>[RowIndex(0)]) == &(<a class="code" href="markowitz_8cc.html#aea35f36ba98d5bbd8d033382f50c9e52">row</a>[ColIndex(0)])));</div>
|
|
<div class="line"><a name="l00204"></a><span class="lineno"> 204</span>  <span class="keywordflow">return</span> <a class="code" href="markowitz_8cc.html#aa9d6c98fdf8d89b0e2321fda02adc82c">col</a>;</div>
|
|
<div class="line"><a name="l00205"></a><span class="lineno"> 205</span> }</div>
|
|
<div class="line"><a name="l00206"></a><span class="lineno"> 206</span>  </div>
|
|
<div class="line"><a name="l00207"></a><span class="lineno"> 207</span> <span class="comment">// Computes the positions of the non-zeros of a dense vector.</span></div>
|
|
<div class="line"><a name="l00208"></a><span class="lineno"> 208</span> <span class="keyword">template</span> <<span class="keyword">typename</span> IndexType></div>
|
|
<div class="line"><a name="l00209"></a><span class="lineno"><a class="line" href="namespaceoperations__research_1_1glop.html#a3e037ab543673629f84850a85c761132"> 209</a></span> <span class="keyword">inline</span> <span class="keywordtype">void</span> <a class="code" href="namespaceoperations__research_1_1glop.html#a3e037ab543673629f84850a85c761132">ComputeNonZeros</a>(<span class="keyword">const</span> <a class="code" href="classoperations__research_1_1glop_1_1_strict_i_t_i_vector.html">StrictITIVector<IndexType, Fractional></a>& <a class="code" href="parser_8yy_8cc.html#a5a634cf4429798b1c921a81de8250051">input</a>,</div>
|
|
<div class="line"><a name="l00210"></a><span class="lineno"> 210</span>  std::vector<IndexType>* non_zeros) {</div>
|
|
<div class="line"><a name="l00211"></a><span class="lineno"> 211</span>  non_zeros->clear();</div>
|
|
<div class="line"><a name="l00212"></a><span class="lineno"> 212</span>  <span class="keyword">const</span> IndexType end = <a class="code" href="parser_8yy_8cc.html#a5a634cf4429798b1c921a81de8250051">input</a>.size();</div>
|
|
<div class="line"><a name="l00213"></a><span class="lineno"> 213</span>  <span class="keywordflow">for</span> (IndexType <a class="code" href="pack_8cc.html#a750b5d744c39a06bfb13e6eb010e35d0">index</a>(0); <a class="code" href="pack_8cc.html#a750b5d744c39a06bfb13e6eb010e35d0">index</a> < end; ++<a class="code" href="pack_8cc.html#a750b5d744c39a06bfb13e6eb010e35d0">index</a>) {</div>
|
|
<div class="line"><a name="l00214"></a><span class="lineno"> 214</span>  <span class="keywordflow">if</span> (<a class="code" href="parser_8yy_8cc.html#a5a634cf4429798b1c921a81de8250051">input</a>[<a class="code" href="pack_8cc.html#a750b5d744c39a06bfb13e6eb010e35d0">index</a>] != 0.0) {</div>
|
|
<div class="line"><a name="l00215"></a><span class="lineno"> 215</span>  non_zeros->push_back(<a class="code" href="pack_8cc.html#a750b5d744c39a06bfb13e6eb010e35d0">index</a>);</div>
|
|
<div class="line"><a name="l00216"></a><span class="lineno"> 216</span>  }</div>
|
|
<div class="line"><a name="l00217"></a><span class="lineno"> 217</span>  }</div>
|
|
<div class="line"><a name="l00218"></a><span class="lineno"> 218</span> }</div>
|
|
<div class="line"><a name="l00219"></a><span class="lineno"> 219</span>  </div>
|
|
<div class="line"><a name="l00220"></a><span class="lineno"> 220</span> <span class="comment">// Returns true if the given Fractional container is all zeros.</span></div>
|
|
<div class="line"><a name="l00221"></a><span class="lineno"> 221</span> <span class="keyword">template</span> <<span class="keyword">typename</span> Container></div>
|
|
<div class="line"><a name="l00222"></a><span class="lineno"><a class="line" href="namespaceoperations__research_1_1glop.html#a386d851d7c8b409205dd82d87e44f2fe"> 222</a></span> <span class="keyword">inline</span> <span class="keywordtype">bool</span> <a class="code" href="namespaceoperations__research_1_1glop.html#a386d851d7c8b409205dd82d87e44f2fe">IsAllZero</a>(<span class="keyword">const</span> Container& <a class="code" href="parser_8yy_8cc.html#a5a634cf4429798b1c921a81de8250051">input</a>) {</div>
|
|
<div class="line"><a name="l00223"></a><span class="lineno"> 223</span>  <span class="keywordflow">for</span> (<a class="code" href="namespaceoperations__research_1_1glop.html#a733947145e3e1631165b618b05c9ccb7">Fractional</a> <a class="code" href="demon__profiler_8cc.html#a21edc7ca4cc5802c8779d68556bc09cf">value</a> : <a class="code" href="parser_8yy_8cc.html#a5a634cf4429798b1c921a81de8250051">input</a>) {</div>
|
|
<div class="line"><a name="l00224"></a><span class="lineno"> 224</span>  <span class="keywordflow">if</span> (<a class="code" href="demon__profiler_8cc.html#a21edc7ca4cc5802c8779d68556bc09cf">value</a> != 0.0) <span class="keywordflow">return</span> <span class="keyword">false</span>;</div>
|
|
<div class="line"><a name="l00225"></a><span class="lineno"> 225</span>  }</div>
|
|
<div class="line"><a name="l00226"></a><span class="lineno"> 226</span>  <span class="keywordflow">return</span> <span class="keyword">true</span>;</div>
|
|
<div class="line"><a name="l00227"></a><span class="lineno"> 227</span> }</div>
|
|
<div class="line"><a name="l00228"></a><span class="lineno"> 228</span>  </div>
|
|
<div class="line"><a name="l00229"></a><span class="lineno"> 229</span> <span class="comment">// Returns true if the given vector of bool is all false.</span></div>
|
|
<div class="line"><a name="l00230"></a><span class="lineno"> 230</span> <span class="keyword">template</span> <<span class="keyword">typename</span> BoolVector></div>
|
|
<div class="line"><a name="l00231"></a><span class="lineno"><a class="line" href="namespaceoperations__research_1_1glop.html#a66f88d7a4bcc601c81b7d694bcfae840"> 231</a></span> <span class="keywordtype">bool</span> <a class="code" href="namespaceoperations__research_1_1glop.html#a66f88d7a4bcc601c81b7d694bcfae840">IsAllFalse</a>(<span class="keyword">const</span> BoolVector& v) {</div>
|
|
<div class="line"><a name="l00232"></a><span class="lineno"> 232</span>  <span class="keywordflow">return</span> std::all_of(v.begin(), v.end(), [](<span class="keywordtype">bool</span> <a class="code" href="demon__profiler_8cc.html#a21edc7ca4cc5802c8779d68556bc09cf">value</a>) { return !value; });</div>
|
|
<div class="line"><a name="l00233"></a><span class="lineno"> 233</span> }</div>
|
|
<div class="line"><a name="l00234"></a><span class="lineno"> 234</span>  </div>
|
|
<div class="line"><a name="l00235"></a><span class="lineno"> 235</span> <span class="comment">// Permutes the given dense vector. It uses for this an all zero scratchpad.</span></div>
|
|
<div class="line"><a name="l00236"></a><span class="lineno"> 236</span> <span class="keyword">template</span> <<span class="keyword">typename</span> IndexType, <span class="keyword">typename</span> PermutationIndexType></div>
|
|
<div class="line"><a name="l00237"></a><span class="lineno"><a class="line" href="namespaceoperations__research_1_1glop.html#a08d7a83791c6677d1008336cacf3d591"> 237</a></span> <span class="keyword">inline</span> <span class="keywordtype">void</span> <a class="code" href="namespaceoperations__research_1_1glop.html#a08d7a83791c6677d1008336cacf3d591">PermuteWithScratchpad</a>(</div>
|
|
<div class="line"><a name="l00238"></a><span class="lineno"> 238</span>  <span class="keyword">const</span> <a class="code" href="classoperations__research_1_1glop_1_1_permutation.html">Permutation<PermutationIndexType></a>& permutation,</div>
|
|
<div class="line"><a name="l00239"></a><span class="lineno"> 239</span>  <a class="code" href="classoperations__research_1_1glop_1_1_strict_i_t_i_vector.html">StrictITIVector<IndexType, Fractional></a>* zero_scratchpad,</div>
|
|
<div class="line"><a name="l00240"></a><span class="lineno"> 240</span>  <a class="code" href="classoperations__research_1_1glop_1_1_strict_i_t_i_vector.html">StrictITIVector<IndexType, Fractional></a>* input_output) {</div>
|
|
<div class="line"><a name="l00241"></a><span class="lineno"> 241</span>  <a class="code" href="base_2logging_8h.html#ae17f8119c108cf3070bad3449c7e0006">DCHECK</a>(<a class="code" href="namespaceoperations__research_1_1glop.html#a386d851d7c8b409205dd82d87e44f2fe">IsAllZero</a>(*zero_scratchpad));</div>
|
|
<div class="line"><a name="l00242"></a><span class="lineno"> 242</span>  <span class="keyword">const</span> IndexType size = input_output-><a class="code" href="classoperations__research_1_1glop_1_1_strict_i_t_i_vector.html#a967a5c081ad4195a30c78dc2c0bcabf5">size</a>();</div>
|
|
<div class="line"><a name="l00243"></a><span class="lineno"> 243</span>  zero_scratchpad-><a class="code" href="classabsl_1_1_strong_vector.html#aa5d16d85614c5d518ae10f882e6fb981">swap</a>(*input_output);</div>
|
|
<div class="line"><a name="l00244"></a><span class="lineno"> 244</span>  input_output-><a class="code" href="classoperations__research_1_1glop_1_1_strict_i_t_i_vector.html#a64b6b04f3a519d2c61d49daaa88bf06e">resize</a>(size, 0.0);</div>
|
|
<div class="line"><a name="l00245"></a><span class="lineno"> 245</span>  <span class="keywordflow">for</span> (IndexType <a class="code" href="pack_8cc.html#a750b5d744c39a06bfb13e6eb010e35d0">index</a>(0); <a class="code" href="pack_8cc.html#a750b5d744c39a06bfb13e6eb010e35d0">index</a> < size; ++<a class="code" href="pack_8cc.html#a750b5d744c39a06bfb13e6eb010e35d0">index</a>) {</div>
|
|
<div class="line"><a name="l00246"></a><span class="lineno"> 246</span>  <span class="keyword">const</span> <a class="code" href="namespaceoperations__research_1_1glop.html#a733947145e3e1631165b618b05c9ccb7">Fractional</a> <a class="code" href="demon__profiler_8cc.html#a21edc7ca4cc5802c8779d68556bc09cf">value</a> = (*zero_scratchpad)[<a class="code" href="pack_8cc.html#a750b5d744c39a06bfb13e6eb010e35d0">index</a>];</div>
|
|
<div class="line"><a name="l00247"></a><span class="lineno"> 247</span>  <span class="keywordflow">if</span> (<a class="code" href="demon__profiler_8cc.html#a21edc7ca4cc5802c8779d68556bc09cf">value</a> != 0.0) {</div>
|
|
<div class="line"><a name="l00248"></a><span class="lineno"> 248</span>  <span class="keyword">const</span> IndexType permuted_index(</div>
|
|
<div class="line"><a name="l00249"></a><span class="lineno"> 249</span>  permutation[PermutationIndexType(<a class="code" href="pack_8cc.html#a750b5d744c39a06bfb13e6eb010e35d0">index</a>.value())].value());</div>
|
|
<div class="line"><a name="l00250"></a><span class="lineno"> 250</span>  (*input_output)[permuted_index] = <a class="code" href="demon__profiler_8cc.html#a21edc7ca4cc5802c8779d68556bc09cf">value</a>;</div>
|
|
<div class="line"><a name="l00251"></a><span class="lineno"> 251</span>  }</div>
|
|
<div class="line"><a name="l00252"></a><span class="lineno"> 252</span>  }</div>
|
|
<div class="line"><a name="l00253"></a><span class="lineno"> 253</span>  zero_scratchpad-><a class="code" href="classoperations__research_1_1glop_1_1_strict_i_t_i_vector.html#af8d7048738ceb4c753b040e6d29db79c">assign</a>(size, 0.0);</div>
|
|
<div class="line"><a name="l00254"></a><span class="lineno"> 254</span> }</div>
|
|
<div class="line"><a name="l00255"></a><span class="lineno"> 255</span>  </div>
|
|
<div class="line"><a name="l00256"></a><span class="lineno"> 256</span> <span class="comment">// Same as PermuteAndComputeNonZeros() except that we assume that the given</span></div>
|
|
<div class="line"><a name="l00257"></a><span class="lineno"> 257</span> <span class="comment">// non-zeros are the initial non-zeros positions of output.</span></div>
|
|
<div class="line"><a name="l00258"></a><span class="lineno"> 258</span> <span class="keyword">template</span> <<span class="keyword">typename</span> IndexType></div>
|
|
<div class="line"><a name="l00259"></a><span class="lineno"><a class="line" href="namespaceoperations__research_1_1glop.html#a6a2019fc6c15a0413896d3f35057a070"> 259</a></span> <span class="keyword">inline</span> <span class="keywordtype">void</span> <a class="code" href="namespaceoperations__research_1_1glop.html#a6a2019fc6c15a0413896d3f35057a070">PermuteWithKnownNonZeros</a>(</div>
|
|
<div class="line"><a name="l00260"></a><span class="lineno"> 260</span>  <span class="keyword">const</span> <a class="code" href="classoperations__research_1_1glop_1_1_permutation.html">Permutation<IndexType></a>& permutation,</div>
|
|
<div class="line"><a name="l00261"></a><span class="lineno"> 261</span>  <a class="code" href="classoperations__research_1_1glop_1_1_strict_i_t_i_vector.html">StrictITIVector<IndexType, Fractional></a>* zero_scratchpad,</div>
|
|
<div class="line"><a name="l00262"></a><span class="lineno"> 262</span>  <a class="code" href="classoperations__research_1_1glop_1_1_strict_i_t_i_vector.html">StrictITIVector<IndexType, Fractional></a>* output,</div>
|
|
<div class="line"><a name="l00263"></a><span class="lineno"> 263</span>  std::vector<IndexType>* non_zeros) {</div>
|
|
<div class="line"><a name="l00264"></a><span class="lineno"> 264</span>  <a class="code" href="base_2logging_8h.html#ae17f8119c108cf3070bad3449c7e0006">DCHECK</a>(<a class="code" href="namespaceoperations__research_1_1glop.html#a386d851d7c8b409205dd82d87e44f2fe">IsAllZero</a>(*zero_scratchpad));</div>
|
|
<div class="line"><a name="l00265"></a><span class="lineno"> 265</span>  zero_scratchpad-><a class="code" href="classabsl_1_1_strong_vector.html#aa5d16d85614c5d518ae10f882e6fb981">swap</a>(*output);</div>
|
|
<div class="line"><a name="l00266"></a><span class="lineno"> 266</span>  output-><a class="code" href="classoperations__research_1_1glop_1_1_strict_i_t_i_vector.html#a64b6b04f3a519d2c61d49daaa88bf06e">resize</a>(zero_scratchpad-><a class="code" href="classoperations__research_1_1glop_1_1_strict_i_t_i_vector.html#a967a5c081ad4195a30c78dc2c0bcabf5">size</a>(), 0.0);</div>
|
|
<div class="line"><a name="l00267"></a><span class="lineno"> 267</span>  <span class="keywordflow">for</span> (IndexType& index_ref : *non_zeros) {</div>
|
|
<div class="line"><a name="l00268"></a><span class="lineno"> 268</span>  <span class="keyword">const</span> <a class="code" href="namespaceoperations__research_1_1glop.html#a733947145e3e1631165b618b05c9ccb7">Fractional</a> <a class="code" href="demon__profiler_8cc.html#a21edc7ca4cc5802c8779d68556bc09cf">value</a> = (*zero_scratchpad)[index_ref];</div>
|
|
<div class="line"><a name="l00269"></a><span class="lineno"> 269</span>  (*zero_scratchpad)[index_ref] = 0.0;</div>
|
|
<div class="line"><a name="l00270"></a><span class="lineno"> 270</span>  <span class="keyword">const</span> IndexType permuted_index(permutation[index_ref]);</div>
|
|
<div class="line"><a name="l00271"></a><span class="lineno"> 271</span>  (*output)[permuted_index] = <a class="code" href="demon__profiler_8cc.html#a21edc7ca4cc5802c8779d68556bc09cf">value</a>;</div>
|
|
<div class="line"><a name="l00272"></a><span class="lineno"> 272</span>  index_ref = permuted_index;</div>
|
|
<div class="line"><a name="l00273"></a><span class="lineno"> 273</span>  }</div>
|
|
<div class="line"><a name="l00274"></a><span class="lineno"> 274</span> }</div>
|
|
<div class="line"><a name="l00275"></a><span class="lineno"> 275</span>  </div>
|
|
<div class="line"><a name="l00276"></a><span class="lineno"> 276</span> <span class="comment">// Sets a dense vector for which the non zeros are known to be non_zeros.</span></div>
|
|
<div class="line"><a name="l00277"></a><span class="lineno"> 277</span> <span class="keyword">template</span> <<span class="keyword">typename</span> IndexType, <span class="keyword">typename</span> ScatteredRowOrCol></div>
|
|
<div class="line"><a name="l00278"></a><span class="lineno"><a class="line" href="namespaceoperations__research_1_1glop.html#aa6c552b94fa80def1d4d1ea64697afb1"> 278</a></span> <span class="keyword">inline</span> <span class="keywordtype">void</span> <a class="code" href="namespaceoperations__research_1_1glop.html#aa6c552b94fa80def1d4d1ea64697afb1">ClearAndResizeVectorWithNonZeros</a>(IndexType size,</div>
|
|
<div class="line"><a name="l00279"></a><span class="lineno"> 279</span>  ScatteredRowOrCol* v) {</div>
|
|
<div class="line"><a name="l00280"></a><span class="lineno"> 280</span>  <span class="comment">// Only use the sparse version if there is less than 5% non-zeros positions</span></div>
|
|
<div class="line"><a name="l00281"></a><span class="lineno"> 281</span>  <span class="comment">// compared to the wanted size. Note that in most cases the vector will</span></div>
|
|
<div class="line"><a name="l00282"></a><span class="lineno"> 282</span>  <span class="comment">// already be of the correct size.</span></div>
|
|
<div class="line"><a name="l00283"></a><span class="lineno"> 283</span>  <span class="keyword">const</span> <span class="keywordtype">double</span> kSparseThreshold = 0.05;</div>
|
|
<div class="line"><a name="l00284"></a><span class="lineno"> 284</span>  <span class="keywordflow">if</span> (!v->non_zeros.empty() &&</div>
|
|
<div class="line"><a name="l00285"></a><span class="lineno"> 285</span>  v->non_zeros.size() < kSparseThreshold * size.value()) {</div>
|
|
<div class="line"><a name="l00286"></a><span class="lineno"> 286</span>  <span class="keywordflow">for</span> (<span class="keyword">const</span> IndexType <a class="code" href="pack_8cc.html#a750b5d744c39a06bfb13e6eb010e35d0">index</a> : v->non_zeros) {</div>
|
|
<div class="line"><a name="l00287"></a><span class="lineno"> 287</span>  <a class="code" href="base_2logging_8h.html#ab62f5ed8f2d48e29802be0cbbcd1359a">DCHECK_LT</a>(<a class="code" href="pack_8cc.html#a750b5d744c39a06bfb13e6eb010e35d0">index</a>, v->values.size());</div>
|
|
<div class="line"><a name="l00288"></a><span class="lineno"> 288</span>  (*v)[<a class="code" href="pack_8cc.html#a750b5d744c39a06bfb13e6eb010e35d0">index</a>] = 0.0;</div>
|
|
<div class="line"><a name="l00289"></a><span class="lineno"> 289</span>  }</div>
|
|
<div class="line"><a name="l00290"></a><span class="lineno"> 290</span>  v->values.resize(size, 0.0);</div>
|
|
<div class="line"><a name="l00291"></a><span class="lineno"> 291</span>  <a class="code" href="base_2logging_8h.html#ae17f8119c108cf3070bad3449c7e0006">DCHECK</a>(<a class="code" href="namespaceoperations__research_1_1glop.html#a386d851d7c8b409205dd82d87e44f2fe">IsAllZero</a>(v->values));</div>
|
|
<div class="line"><a name="l00292"></a><span class="lineno"> 292</span>  } <span class="keywordflow">else</span> {</div>
|
|
<div class="line"><a name="l00293"></a><span class="lineno"> 293</span>  v->values.AssignToZero(size);</div>
|
|
<div class="line"><a name="l00294"></a><span class="lineno"> 294</span>  }</div>
|
|
<div class="line"><a name="l00295"></a><span class="lineno"> 295</span>  v->non_zeros.clear();</div>
|
|
<div class="line"><a name="l00296"></a><span class="lineno"> 296</span> }</div>
|
|
<div class="line"><a name="l00297"></a><span class="lineno"> 297</span>  </div>
|
|
<div class="line"><a name="l00298"></a><span class="lineno"> 298</span> <span class="comment">// Changes the sign of all the entries in the given vector.</span></div>
|
|
<div class="line"><a name="l00299"></a><span class="lineno"> 299</span> <span class="keyword">template</span> <<span class="keyword">typename</span> IndexType></div>
|
|
<div class="line"><a name="l00300"></a><span class="lineno"><a class="line" href="namespaceoperations__research_1_1glop.html#ab7ff1f9383648fd132667b1e601229c9"> 300</a></span> <span class="keyword">inline</span> <span class="keywordtype">void</span> <a class="code" href="namespaceoperations__research_1_1glop.html#ab7ff1f9383648fd132667b1e601229c9">ChangeSign</a>(<a class="code" href="classoperations__research_1_1glop_1_1_strict_i_t_i_vector.html">StrictITIVector<IndexType, Fractional></a>* data) {</div>
|
|
<div class="line"><a name="l00301"></a><span class="lineno"> 301</span>  <span class="keyword">const</span> IndexType end = data-><a class="code" href="classoperations__research_1_1glop_1_1_strict_i_t_i_vector.html#a967a5c081ad4195a30c78dc2c0bcabf5">size</a>();</div>
|
|
<div class="line"><a name="l00302"></a><span class="lineno"> 302</span>  <span class="keywordflow">for</span> (IndexType i(0); i < end; ++i) {</div>
|
|
<div class="line"><a name="l00303"></a><span class="lineno"> 303</span>  (*data)[i] = -(*data)[i];</div>
|
|
<div class="line"><a name="l00304"></a><span class="lineno"> 304</span>  }</div>
|
|
<div class="line"><a name="l00305"></a><span class="lineno"> 305</span> }</div>
|
|
<div class="line"><a name="l00306"></a><span class="lineno"> 306</span>  </div>
|
|
<div class="line"><a name="l00307"></a><span class="lineno"> 307</span> <span class="comment">// Given N Fractional elements, this class maintains their sum and can</span></div>
|
|
<div class="line"><a name="l00308"></a><span class="lineno"> 308</span> <span class="comment">// provide, for each element X, the sum of all elements except X.</span></div>
|
|
<div class="line"><a name="l00309"></a><span class="lineno"> 309</span> <span class="comment">// The subtelty is that it works well with infinities: for example, if there is</span></div>
|
|
<div class="line"><a name="l00310"></a><span class="lineno"> 310</span> <span class="comment">// exactly one infinite element X, then SumWithout(X) will be finite.</span></div>
|
|
<div class="line"><a name="l00311"></a><span class="lineno"> 311</span> <span class="comment">//</span></div>
|
|
<div class="line"><a name="l00312"></a><span class="lineno"> 312</span> <span class="comment">// Two flavors of this class are provided: SumWithPositiveInfiniteAndOneMissing</span></div>
|
|
<div class="line"><a name="l00313"></a><span class="lineno"> 313</span> <span class="comment">// supports calling Add() with normal numbers and positive infinities (and will</span></div>
|
|
<div class="line"><a name="l00314"></a><span class="lineno"> 314</span> <span class="comment">// DCHECK() that), and SumWithNegativeInfiniteAndOneMissing does the same with</span></div>
|
|
<div class="line"><a name="l00315"></a><span class="lineno"> 315</span> <span class="comment">// negative infinities.</span></div>
|
|
<div class="line"><a name="l00316"></a><span class="lineno"> 316</span> <span class="comment">//</span></div>
|
|
<div class="line"><a name="l00317"></a><span class="lineno"> 317</span> <span class="comment">// The numerical accuracy suffers however. If X is 1e100 and SumWithout(X)</span></div>
|
|
<div class="line"><a name="l00318"></a><span class="lineno"> 318</span> <span class="comment">// should be 1e-100, then the value actually returned by SumWithout(X) is likely</span></div>
|
|
<div class="line"><a name="l00319"></a><span class="lineno"> 319</span> <span class="comment">// to be wrong (by up to std::numeric_limits<Fractional>::epsilon() ^ 2).</span></div>
|
|
<div class="line"><a name="l00320"></a><span class="lineno"> 320</span> <span class="keyword">template</span> <<span class="keywordtype">bool</span> supported_infinity_is_positive></div>
|
|
<div class="line"><a name="l00321"></a><span class="lineno"><a class="line" href="classoperations__research_1_1glop_1_1_sum_with_one_missing.html"> 321</a></span> <span class="keyword">class </span><a class="code" href="classoperations__research_1_1glop_1_1_sum_with_one_missing.html">SumWithOneMissing</a> {</div>
|
|
<div class="line"><a name="l00322"></a><span class="lineno"> 322</span>  <span class="keyword">public</span>:</div>
|
|
<div class="line"><a name="l00323"></a><span class="lineno"><a class="line" href="classoperations__research_1_1glop_1_1_sum_with_one_missing.html#ac49fbdf48898a72ccf672da218def05b"> 323</a></span>  <a class="code" href="classoperations__research_1_1glop_1_1_sum_with_one_missing.html#ac49fbdf48898a72ccf672da218def05b">SumWithOneMissing</a>() : num_infinities_(0), sum_() {}</div>
|
|
<div class="line"><a name="l00324"></a><span class="lineno"> 324</span>  </div>
|
|
<div class="line"><a name="l00325"></a><span class="lineno"><a class="line" href="classoperations__research_1_1glop_1_1_sum_with_one_missing.html#a8e16e2adc11d40ca12b2735ae21623b8"> 325</a></span>  <span class="keywordtype">void</span> <a class="code" href="classoperations__research_1_1glop_1_1_sum_with_one_missing.html#a8e16e2adc11d40ca12b2735ae21623b8">Add</a>(<a class="code" href="namespaceoperations__research_1_1glop.html#a733947145e3e1631165b618b05c9ccb7">Fractional</a> x) {</div>
|
|
<div class="line"><a name="l00326"></a><span class="lineno"> 326</span>  <span class="keywordflow">if</span> (num_infinities_ > 1) <span class="keywordflow">return</span>;</div>
|
|
<div class="line"><a name="l00327"></a><span class="lineno"> 327</span>  <span class="keywordflow">if</span> (<a class="code" href="namespaceoperations__research_1_1glop.html#a95879916d90daeba91c40399ae5ddcc6">IsFinite</a>(x)) {</div>
|
|
<div class="line"><a name="l00328"></a><span class="lineno"> 328</span>  sum_.<a class="code" href="classoperations__research_1_1_accurate_sum.html#acc8fad1f98b885fa7d0504f4dfc71f08">Add</a>(x);</div>
|
|
<div class="line"><a name="l00329"></a><span class="lineno"> 329</span>  <span class="keywordflow">return</span>;</div>
|
|
<div class="line"><a name="l00330"></a><span class="lineno"> 330</span>  }</div>
|
|
<div class="line"><a name="l00331"></a><span class="lineno"> 331</span>  <a class="code" href="base_2logging_8h.html#ae89df3243bbb8341130c7b3f44145ea0">DCHECK_EQ</a>(Infinity(), x);</div>
|
|
<div class="line"><a name="l00332"></a><span class="lineno"> 332</span>  ++num_infinities_;</div>
|
|
<div class="line"><a name="l00333"></a><span class="lineno"> 333</span>  }</div>
|
|
<div class="line"><a name="l00334"></a><span class="lineno"> 334</span>  </div>
|
|
<div class="line"><a name="l00335"></a><span class="lineno"><a class="line" href="classoperations__research_1_1glop_1_1_sum_with_one_missing.html#af7c8b5d68bfc286dbb52a547917de1af"> 335</a></span>  <a class="code" href="namespaceoperations__research_1_1glop.html#a733947145e3e1631165b618b05c9ccb7">Fractional</a> <a class="code" href="classoperations__research_1_1glop_1_1_sum_with_one_missing.html#af7c8b5d68bfc286dbb52a547917de1af">Sum</a>()<span class="keyword"> const </span>{</div>
|
|
<div class="line"><a name="l00336"></a><span class="lineno"> 336</span>  <span class="keywordflow">if</span> (num_infinities_ > 0) <span class="keywordflow">return</span> Infinity();</div>
|
|
<div class="line"><a name="l00337"></a><span class="lineno"> 337</span>  <span class="keywordflow">return</span> sum_.<a class="code" href="classoperations__research_1_1_accurate_sum.html#a176a3e919acd979b67cea1ede094cdaa">Value</a>();</div>
|
|
<div class="line"><a name="l00338"></a><span class="lineno"> 338</span>  }</div>
|
|
<div class="line"><a name="l00339"></a><span class="lineno"> 339</span>  </div>
|
|
<div class="line"><a name="l00340"></a><span class="lineno"><a class="line" href="classoperations__research_1_1glop_1_1_sum_with_one_missing.html#af6e4d9c2405447b44b4cf701e43e2200"> 340</a></span>  <a class="code" href="namespaceoperations__research_1_1glop.html#a733947145e3e1631165b618b05c9ccb7">Fractional</a> <a class="code" href="classoperations__research_1_1glop_1_1_sum_with_one_missing.html#af6e4d9c2405447b44b4cf701e43e2200">SumWithout</a>(<a class="code" href="namespaceoperations__research_1_1glop.html#a733947145e3e1631165b618b05c9ccb7">Fractional</a> x)<span class="keyword"> const </span>{</div>
|
|
<div class="line"><a name="l00341"></a><span class="lineno"> 341</span>  <span class="keywordflow">if</span> (<a class="code" href="namespaceoperations__research_1_1glop.html#a95879916d90daeba91c40399ae5ddcc6">IsFinite</a>(x)) {</div>
|
|
<div class="line"><a name="l00342"></a><span class="lineno"> 342</span>  <span class="keywordflow">if</span> (num_infinities_ > 0) <span class="keywordflow">return</span> Infinity();</div>
|
|
<div class="line"><a name="l00343"></a><span class="lineno"> 343</span>  <span class="keywordflow">return</span> sum_.<a class="code" href="classoperations__research_1_1_accurate_sum.html#a176a3e919acd979b67cea1ede094cdaa">Value</a>() - x;</div>
|
|
<div class="line"><a name="l00344"></a><span class="lineno"> 344</span>  }</div>
|
|
<div class="line"><a name="l00345"></a><span class="lineno"> 345</span>  <a class="code" href="base_2logging_8h.html#ae89df3243bbb8341130c7b3f44145ea0">DCHECK_EQ</a>(Infinity(), x);</div>
|
|
<div class="line"><a name="l00346"></a><span class="lineno"> 346</span>  <span class="keywordflow">if</span> (num_infinities_ > 1) <span class="keywordflow">return</span> Infinity();</div>
|
|
<div class="line"><a name="l00347"></a><span class="lineno"> 347</span>  <span class="keywordflow">return</span> sum_.<a class="code" href="classoperations__research_1_1_accurate_sum.html#a176a3e919acd979b67cea1ede094cdaa">Value</a>();</div>
|
|
<div class="line"><a name="l00348"></a><span class="lineno"> 348</span>  }</div>
|
|
<div class="line"><a name="l00349"></a><span class="lineno"> 349</span>  </div>
|
|
<div class="line"><a name="l00350"></a><span class="lineno"> 350</span>  <span class="keyword">private</span>:</div>
|
|
<div class="line"><a name="l00351"></a><span class="lineno"> 351</span>  <a class="code" href="namespaceoperations__research_1_1glop.html#a733947145e3e1631165b618b05c9ccb7">Fractional</a> Infinity()<span class="keyword"> const </span>{</div>
|
|
<div class="line"><a name="l00352"></a><span class="lineno"> 352</span>  <span class="keywordflow">return</span> supported_infinity_is_positive ? <a class="code" href="namespaceoperations__research_1_1glop.html#af9a790b7e8c5b0c6d55b336177378e78">kInfinity</a> : -<a class="code" href="namespaceoperations__research_1_1glop.html#af9a790b7e8c5b0c6d55b336177378e78">kInfinity</a>;</div>
|
|
<div class="line"><a name="l00353"></a><span class="lineno"> 353</span>  }</div>
|
|
<div class="line"><a name="l00354"></a><span class="lineno"> 354</span>  </div>
|
|
<div class="line"><a name="l00355"></a><span class="lineno"> 355</span>  <span class="comment">// Count how many times Add() was called with an infinite value. The count is</span></div>
|
|
<div class="line"><a name="l00356"></a><span class="lineno"> 356</span>  <span class="comment">// stopped at 2 to be a bit faster.</span></div>
|
|
<div class="line"><a name="l00357"></a><span class="lineno"> 357</span>  <span class="keywordtype">int</span> num_infinities_;</div>
|
|
<div class="line"><a name="l00358"></a><span class="lineno"> 358</span>  <a class="code" href="namespaceoperations__research_1_1glop.html#a394b55c9b64ac7ce7c9e8772ef9f7c73">KahanSum</a> sum_; <span class="comment">// stripped of all the infinite values.</span></div>
|
|
<div class="line"><a name="l00359"></a><span class="lineno"> 359</span> };</div>
|
|
<div class="line"><a name="l00360"></a><span class="lineno"><a class="line" href="namespaceoperations__research_1_1glop.html#aeb4b2cc773e71eeb6b9d2b6f4c05a858"> 360</a></span> <span class="keyword">typedef</span> <a class="code" href="classoperations__research_1_1glop_1_1_sum_with_one_missing.html">SumWithOneMissing<true></a> <a class="code" href="namespaceoperations__research_1_1glop.html#aeb4b2cc773e71eeb6b9d2b6f4c05a858">SumWithPositiveInfiniteAndOneMissing</a>;</div>
|
|
<div class="line"><a name="l00361"></a><span class="lineno"><a class="line" href="namespaceoperations__research_1_1glop.html#a64c3eaa146467633bb8fdd8fbc0f9482"> 361</a></span> <span class="keyword">typedef</span> <a class="code" href="classoperations__research_1_1glop_1_1_sum_with_one_missing.html">SumWithOneMissing<false></a> <a class="code" href="namespaceoperations__research_1_1glop.html#a64c3eaa146467633bb8fdd8fbc0f9482">SumWithNegativeInfiniteAndOneMissing</a>;</div>
|
|
<div class="line"><a name="l00362"></a><span class="lineno"> 362</span>  </div>
|
|
<div class="line"><a name="l00363"></a><span class="lineno"> 363</span> } <span class="comment">// namespace glop</span></div>
|
|
<div class="line"><a name="l00364"></a><span class="lineno"> 364</span> } <span class="comment">// namespace operations_research</span></div>
|
|
<div class="line"><a name="l00365"></a><span class="lineno"> 365</span>  </div>
|
|
<div class="line"><a name="l00366"></a><span class="lineno"> 366</span> <span class="preprocessor">#endif </span><span class="comment">// OR_TOOLS_LP_DATA_LP_UTILS_H_</span></div>
|
|
<div class="ttc" id="aaccurate__sum_8h_html"><div class="ttname"><a href="accurate__sum_8h.html">accurate_sum.h</a></div></div>
|
|
<div class="ttc" id="abase_2logging_8h_html_ab62f5ed8f2d48e29802be0cbbcd1359a"><div class="ttname"><a href="base_2logging_8h.html#ab62f5ed8f2d48e29802be0cbbcd1359a">DCHECK_LT</a></div><div class="ttdeci">#define DCHECK_LT(val1, val2)</div><div class="ttdef"><b>Definition:</b> <a href="base_2logging_8h_source.html#l00888">base/logging.h:888</a></div></div>
|
|
<div class="ttc" id="abase_2logging_8h_html_ae17f8119c108cf3070bad3449c7e0006"><div class="ttname"><a href="base_2logging_8h.html#ae17f8119c108cf3070bad3449c7e0006">DCHECK</a></div><div class="ttdeci">#define DCHECK(condition)</div><div class="ttdef"><b>Definition:</b> <a href="base_2logging_8h_source.html#l00884">base/logging.h:884</a></div></div>
|
|
<div class="ttc" id="abase_2logging_8h_html_ae89df3243bbb8341130c7b3f44145ea0"><div class="ttname"><a href="base_2logging_8h.html#ae89df3243bbb8341130c7b3f44145ea0">DCHECK_EQ</a></div><div class="ttdeci">#define DCHECK_EQ(val1, val2)</div><div class="ttdef"><b>Definition:</b> <a href="base_2logging_8h_source.html#l00885">base/logging.h:885</a></div></div>
|
|
<div class="ttc" id="aclassabsl_1_1_strong_vector_html_aa5d16d85614c5d518ae10f882e6fb981"><div class="ttname"><a href="classabsl_1_1_strong_vector.html#aa5d16d85614c5d518ae10f882e6fb981">absl::StrongVector::swap</a></div><div class="ttdeci">void swap(StrongVector &x)</div><div class="ttdef"><b>Definition:</b> <a href="strong__vector_8h_source.html#l00169">strong_vector.h:169</a></div></div>
|
|
<div class="ttc" id="aclassoperations__research_1_1_accurate_sum_html"><div class="ttname"><a href="classoperations__research_1_1_accurate_sum.html">operations_research::AccurateSum< Fractional ></a></div></div>
|
|
<div class="ttc" id="aclassoperations__research_1_1_accurate_sum_html_a176a3e919acd979b67cea1ede094cdaa"><div class="ttname"><a href="classoperations__research_1_1_accurate_sum.html#a176a3e919acd979b67cea1ede094cdaa">operations_research::AccurateSum::Value</a></div><div class="ttdeci">FpNumber Value() const</div><div class="ttdef"><b>Definition:</b> <a href="accurate__sum_8h_source.html#l00037">accurate_sum.h:37</a></div></div>
|
|
<div class="ttc" id="aclassoperations__research_1_1_accurate_sum_html_acc8fad1f98b885fa7d0504f4dfc71f08"><div class="ttname"><a href="classoperations__research_1_1_accurate_sum.html#acc8fad1f98b885fa7d0504f4dfc71f08">operations_research::AccurateSum::Add</a></div><div class="ttdeci">void Add(const FpNumber &value)</div><div class="ttdef"><b>Definition:</b> <a href="accurate__sum_8h_source.html#l00029">accurate_sum.h:29</a></div></div>
|
|
<div class="ttc" id="aclassoperations__research_1_1glop_1_1_permutation_html"><div class="ttname"><a href="classoperations__research_1_1glop_1_1_permutation.html">operations_research::glop::Permutation</a></div><div class="ttdef"><b>Definition:</b> <a href="lp__data_2permutation_8h_source.html#l00044">lp_data/permutation.h:44</a></div></div>
|
|
<div class="ttc" id="aclassoperations__research_1_1glop_1_1_sparse_column_html"><div class="ttname"><a href="classoperations__research_1_1glop_1_1_sparse_column.html">operations_research::glop::SparseColumn</a></div><div class="ttdef"><b>Definition:</b> <a href="sparse__column_8h_source.html#l00044">sparse_column.h:44</a></div></div>
|
|
<div class="ttc" id="aclassoperations__research_1_1glop_1_1_sparse_vector_html_ab38326ea6cb6187267665dd8b2748f3d"><div class="ttname"><a href="classoperations__research_1_1glop_1_1_sparse_vector.html#ab38326ea6cb6187267665dd8b2748f3d">operations_research::glop::SparseVector< RowIndex, SparseColumnIterator >::Entry</a></div><div class="ttdeci">typename Iterator::Entry Entry</div><div class="ttdef"><b>Definition:</b> <a href="sparse__vector_8h_source.html#l00091">sparse_vector.h:91</a></div></div>
|
|
<div class="ttc" id="aclassoperations__research_1_1glop_1_1_strict_i_t_i_vector_html"><div class="ttname"><a href="classoperations__research_1_1glop_1_1_strict_i_t_i_vector.html">operations_research::glop::StrictITIVector< ColIndex, Fractional ></a></div></div>
|
|
<div class="ttc" id="aclassoperations__research_1_1glop_1_1_strict_i_t_i_vector_html_a64b6b04f3a519d2c61d49daaa88bf06e"><div class="ttname"><a href="classoperations__research_1_1glop_1_1_strict_i_t_i_vector.html#a64b6b04f3a519d2c61d49daaa88bf06e">operations_research::glop::StrictITIVector::resize</a></div><div class="ttdeci">void resize(IntType size)</div><div class="ttdef"><b>Definition:</b> <a href="lp__types_8h_source.html#l00269">lp_types.h:269</a></div></div>
|
|
<div class="ttc" id="aclassoperations__research_1_1glop_1_1_strict_i_t_i_vector_html_a967a5c081ad4195a30c78dc2c0bcabf5"><div class="ttname"><a href="classoperations__research_1_1glop_1_1_strict_i_t_i_vector.html#a967a5c081ad4195a30c78dc2c0bcabf5">operations_research::glop::StrictITIVector::size</a></div><div class="ttdeci">IntType size() const</div><div class="ttdef"><b>Definition:</b> <a href="lp__types_8h_source.html#l00276">lp_types.h:276</a></div></div>
|
|
<div class="ttc" id="aclassoperations__research_1_1glop_1_1_strict_i_t_i_vector_html_af8d7048738ceb4c753b040e6d29db79c"><div class="ttname"><a href="classoperations__research_1_1glop_1_1_strict_i_t_i_vector.html#af8d7048738ceb4c753b040e6d29db79c">operations_research::glop::StrictITIVector::assign</a></div><div class="ttdeci">void assign(IntType size, const T &v)</div><div class="ttdef"><b>Definition:</b> <a href="lp__types_8h_source.html#l00274">lp_types.h:274</a></div></div>
|
|
<div class="ttc" id="aclassoperations__research_1_1glop_1_1_sum_with_one_missing_html"><div class="ttname"><a href="classoperations__research_1_1glop_1_1_sum_with_one_missing.html">operations_research::glop::SumWithOneMissing</a></div><div class="ttdef"><b>Definition:</b> <a href="lp__data_2lp__utils_8h_source.html#l00321">lp_data/lp_utils.h:321</a></div></div>
|
|
<div class="ttc" id="aclassoperations__research_1_1glop_1_1_sum_with_one_missing_html_a8e16e2adc11d40ca12b2735ae21623b8"><div class="ttname"><a href="classoperations__research_1_1glop_1_1_sum_with_one_missing.html#a8e16e2adc11d40ca12b2735ae21623b8">operations_research::glop::SumWithOneMissing::Add</a></div><div class="ttdeci">void Add(Fractional x)</div><div class="ttdef"><b>Definition:</b> <a href="lp__data_2lp__utils_8h_source.html#l00325">lp_data/lp_utils.h:325</a></div></div>
|
|
<div class="ttc" id="aclassoperations__research_1_1glop_1_1_sum_with_one_missing_html_ac49fbdf48898a72ccf672da218def05b"><div class="ttname"><a href="classoperations__research_1_1glop_1_1_sum_with_one_missing.html#ac49fbdf48898a72ccf672da218def05b">operations_research::glop::SumWithOneMissing::SumWithOneMissing</a></div><div class="ttdeci">SumWithOneMissing()</div><div class="ttdef"><b>Definition:</b> <a href="lp__data_2lp__utils_8h_source.html#l00323">lp_data/lp_utils.h:323</a></div></div>
|
|
<div class="ttc" id="aclassoperations__research_1_1glop_1_1_sum_with_one_missing_html_af6e4d9c2405447b44b4cf701e43e2200"><div class="ttname"><a href="classoperations__research_1_1glop_1_1_sum_with_one_missing.html#af6e4d9c2405447b44b4cf701e43e2200">operations_research::glop::SumWithOneMissing::SumWithout</a></div><div class="ttdeci">Fractional SumWithout(Fractional x) const</div><div class="ttdef"><b>Definition:</b> <a href="lp__data_2lp__utils_8h_source.html#l00340">lp_data/lp_utils.h:340</a></div></div>
|
|
<div class="ttc" id="aclassoperations__research_1_1glop_1_1_sum_with_one_missing_html_af7c8b5d68bfc286dbb52a547917de1af"><div class="ttname"><a href="classoperations__research_1_1glop_1_1_sum_with_one_missing.html#af7c8b5d68bfc286dbb52a547917de1af">operations_research::glop::SumWithOneMissing::Sum</a></div><div class="ttdeci">Fractional Sum() const</div><div class="ttdef"><b>Definition:</b> <a href="lp__data_2lp__utils_8h_source.html#l00335">lp_data/lp_utils.h:335</a></div></div>
|
|
<div class="ttc" id="aconstraint__solver_2table_8cc_html_a344010e26426d6a13411648d988bc9b6"><div class="ttname"><a href="constraint__solver_2table_8cc.html#a344010e26426d6a13411648d988bc9b6">b</a></div><div class="ttdeci">int64 b</div><div class="ttdef"><b>Definition:</b> <a href="constraint__solver_2table_8cc_source.html#l00043">constraint_solver/table.cc:43</a></div></div>
|
|
<div class="ttc" id="ademon__profiler_8cc_html_a21edc7ca4cc5802c8779d68556bc09cf"><div class="ttname"><a href="demon__profiler_8cc.html#a21edc7ca4cc5802c8779d68556bc09cf">value</a></div><div class="ttdeci">int64 value</div><div class="ttdef"><b>Definition:</b> <a href="demon__profiler_8cc_source.html#l00043">demon_profiler.cc:43</a></div></div>
|
|
<div class="ttc" id="alp__types_8h_html"><div class="ttname"><a href="lp__types_8h.html">lp_types.h</a></div></div>
|
|
<div class="ttc" id="amarkowitz_8cc_html_aa9d6c98fdf8d89b0e2321fda02adc82c"><div class="ttname"><a href="markowitz_8cc.html#aa9d6c98fdf8d89b0e2321fda02adc82c">col</a></div><div class="ttdeci">ColIndex col</div><div class="ttdef"><b>Definition:</b> <a href="markowitz_8cc_source.html#l00176">markowitz.cc:176</a></div></div>
|
|
<div class="ttc" id="amarkowitz_8cc_html_aea35f36ba98d5bbd8d033382f50c9e52"><div class="ttname"><a href="markowitz_8cc.html#aea35f36ba98d5bbd8d033382f50c9e52">row</a></div><div class="ttdeci">RowIndex row</div><div class="ttdef"><b>Definition:</b> <a href="markowitz_8cc_source.html#l00175">markowitz.cc:175</a></div></div>
|
|
<div class="ttc" id="anamespaceoperations__research_1_1glop_html_a08d7a83791c6677d1008336cacf3d591"><div class="ttname"><a href="namespaceoperations__research_1_1glop.html#a08d7a83791c6677d1008336cacf3d591">operations_research::glop::PermuteWithScratchpad</a></div><div class="ttdeci">void PermuteWithScratchpad(const Permutation< PermutationIndexType > &permutation, StrictITIVector< IndexType, Fractional > *zero_scratchpad, StrictITIVector< IndexType, Fractional > *input_output)</div><div class="ttdef"><b>Definition:</b> <a href="lp__data_2lp__utils_8h_source.html#l00237">lp_data/lp_utils.h:237</a></div></div>
|
|
<div class="ttc" id="anamespaceoperations__research_1_1glop_html_a1dcd08b0f6c19cd4a302bb5a3a6ea06e"><div class="ttname"><a href="namespaceoperations__research_1_1glop.html#a1dcd08b0f6c19cd4a302bb5a3a6ea06e">operations_research::glop::Square</a></div><div class="ttdeci">Fractional Square(Fractional f)</div><div class="ttdef"><b>Definition:</b> <a href="lp__data_2lp__utils_8h_source.html#l00036">lp_data/lp_utils.h:36</a></div></div>
|
|
<div class="ttc" id="anamespaceoperations__research_1_1glop_html_a1faa927dd93b43b3dea3eb2a993e30a1"><div class="ttname"><a href="namespaceoperations__research_1_1glop.html#a1faa927dd93b43b3dea3eb2a993e30a1">operations_research::glop::PreciseSquaredNorm</a></div><div class="ttdeci">Fractional PreciseSquaredNorm(const SparseColumn &v)</div><div class="ttdef"><b>Definition:</b> <a href="lp__data_2lp__utils_8cc_source.html#l00038">lp_data/lp_utils.cc:38</a></div></div>
|
|
<div class="ttc" id="anamespaceoperations__research_1_1glop_html_a2ad7ea612bc859f5b637d5029fb875fb"><div class="ttname"><a href="namespaceoperations__research_1_1glop.html#a2ad7ea612bc859f5b637d5029fb875fb">operations_research::glop::InfinityNorm</a></div><div class="ttdeci">Fractional InfinityNorm(const DenseColumn &v)</div><div class="ttdef"><b>Definition:</b> <a href="lp__data_2lp__utils_8cc_source.html#l00081">lp_data/lp_utils.cc:81</a></div></div>
|
|
<div class="ttc" id="anamespaceoperations__research_1_1glop_html_a2d53948bf5e999d006e781105aa8bc77"><div class="ttname"><a href="namespaceoperations__research_1_1glop.html#a2d53948bf5e999d006e781105aa8bc77">operations_research::glop::SquaredNorm</a></div><div class="ttdeci">Fractional SquaredNorm(const SparseColumn &v)</div><div class="ttdef"><b>Definition:</b> <a href="lp__data_2lp__utils_8cc_source.html#l00030">lp_data/lp_utils.cc:30</a></div></div>
|
|
<div class="ttc" id="anamespaceoperations__research_1_1glop_html_a386d851d7c8b409205dd82d87e44f2fe"><div class="ttname"><a href="namespaceoperations__research_1_1glop.html#a386d851d7c8b409205dd82d87e44f2fe">operations_research::glop::IsAllZero</a></div><div class="ttdeci">bool IsAllZero(const Container &input)</div><div class="ttdef"><b>Definition:</b> <a href="lp__data_2lp__utils_8h_source.html#l00222">lp_data/lp_utils.h:222</a></div></div>
|
|
<div class="ttc" id="anamespaceoperations__research_1_1glop_html_a394b55c9b64ac7ce7c9e8772ef9f7c73"><div class="ttname"><a href="namespaceoperations__research_1_1glop.html#a394b55c9b64ac7ce7c9e8772ef9f7c73">operations_research::glop::KahanSum</a></div><div class="ttdeci">AccurateSum< Fractional > KahanSum</div><div class="ttdef"><b>Definition:</b> <a href="lp__data_2lp__utils_8h_source.html#l00032">lp_data/lp_utils.h:32</a></div></div>
|
|
<div class="ttc" id="anamespaceoperations__research_1_1glop_html_a3d9874910f9632900fe8677463b1c0be"><div class="ttname"><a href="namespaceoperations__research_1_1glop.html#a3d9874910f9632900fe8677463b1c0be">operations_research::glop::ScalarProduct</a></div><div class="ttdeci">Fractional ScalarProduct(const DenseRowOrColumn1 &u, const DenseRowOrColumn2 &v)</div><div class="ttdef"><b>Definition:</b> <a href="lp__data_2lp__utils_8h_source.html#l00047">lp_data/lp_utils.h:47</a></div></div>
|
|
<div class="ttc" id="anamespaceoperations__research_1_1glop_html_a3e037ab543673629f84850a85c761132"><div class="ttname"><a href="namespaceoperations__research_1_1glop.html#a3e037ab543673629f84850a85c761132">operations_research::glop::ComputeNonZeros</a></div><div class="ttdeci">void ComputeNonZeros(const StrictITIVector< IndexType, Fractional > &input, std::vector< IndexType > *non_zeros)</div><div class="ttdef"><b>Definition:</b> <a href="lp__data_2lp__utils_8h_source.html#l00209">lp_data/lp_utils.h:209</a></div></div>
|
|
<div class="ttc" id="anamespaceoperations__research_1_1glop_html_a46fb729c0be27d1b97db15e0ce9c6067"><div class="ttname"><a href="namespaceoperations__research_1_1glop.html#a46fb729c0be27d1b97db15e0ce9c6067">operations_research::glop::PreciseScalarProduct</a></div><div class="ttdeci">Fractional PreciseScalarProduct(const DenseRowOrColumn &u, const DenseRowOrColumn2 &v)</div><div class="ttdef"><b>Definition:</b> <a href="lp__data_2lp__utils_8h_source.html#l00092">lp_data/lp_utils.h:92</a></div></div>
|
|
<div class="ttc" id="anamespaceoperations__research_1_1glop_html_a4dc8ae0f97c4110f3cc770715b2bbd12"><div class="ttname"><a href="namespaceoperations__research_1_1glop.html#a4dc8ae0f97c4110f3cc770715b2bbd12">operations_research::glop::DenseRow</a></div><div class="ttdeci">StrictITIVector< ColIndex, Fractional > DenseRow</div><div class="ttdef"><b>Definition:</b> <a href="lp__types_8h_source.html#l00299">lp_types.h:299</a></div></div>
|
|
<div class="ttc" id="anamespaceoperations__research_1_1glop_html_a5e79e30b7239adc4fb2a27778335bca0"><div class="ttname"><a href="namespaceoperations__research_1_1glop.html#a5e79e30b7239adc4fb2a27778335bca0">operations_research::glop::RemoveNearZeroEntries</a></div><div class="ttdeci">void RemoveNearZeroEntries(Fractional threshold, DenseRow *row)</div><div class="ttdef"><b>Definition:</b> <a href="lp__data_2lp__utils_8cc_source.html#l00115">lp_data/lp_utils.cc:115</a></div></div>
|
|
<div class="ttc" id="anamespaceoperations__research_1_1glop_html_a64c3eaa146467633bb8fdd8fbc0f9482"><div class="ttname"><a href="namespaceoperations__research_1_1glop.html#a64c3eaa146467633bb8fdd8fbc0f9482">operations_research::glop::SumWithNegativeInfiniteAndOneMissing</a></div><div class="ttdeci">SumWithOneMissing< false > SumWithNegativeInfiniteAndOneMissing</div><div class="ttdef"><b>Definition:</b> <a href="lp__data_2lp__utils_8h_source.html#l00361">lp_data/lp_utils.h:361</a></div></div>
|
|
<div class="ttc" id="anamespaceoperations__research_1_1glop_html_a66f88d7a4bcc601c81b7d694bcfae840"><div class="ttname"><a href="namespaceoperations__research_1_1glop.html#a66f88d7a4bcc601c81b7d694bcfae840">operations_research::glop::IsAllFalse</a></div><div class="ttdeci">bool IsAllFalse(const BoolVector &v)</div><div class="ttdef"><b>Definition:</b> <a href="lp__data_2lp__utils_8h_source.html#l00231">lp_data/lp_utils.h:231</a></div></div>
|
|
<div class="ttc" id="anamespaceoperations__research_1_1glop_html_a6a2019fc6c15a0413896d3f35057a070"><div class="ttname"><a href="namespaceoperations__research_1_1glop.html#a6a2019fc6c15a0413896d3f35057a070">operations_research::glop::PermuteWithKnownNonZeros</a></div><div class="ttdeci">void PermuteWithKnownNonZeros(const Permutation< IndexType > &permutation, StrictITIVector< IndexType, Fractional > *zero_scratchpad, StrictITIVector< IndexType, Fractional > *output, std::vector< IndexType > *non_zeros)</div><div class="ttdef"><b>Definition:</b> <a href="lp__data_2lp__utils_8h_source.html#l00259">lp_data/lp_utils.h:259</a></div></div>
|
|
<div class="ttc" id="anamespaceoperations__research_1_1glop_html_a6b36bbbcd695191600dab9b88487ff18"><div class="ttname"><a href="namespaceoperations__research_1_1glop.html#a6b36bbbcd695191600dab9b88487ff18">operations_research::glop::Density</a></div><div class="ttdeci">double Density(const DenseRow &row)</div><div class="ttdef"><b>Definition:</b> <a href="lp__data_2lp__utils_8cc_source.html#l00106">lp_data/lp_utils.cc:106</a></div></div>
|
|
<div class="ttc" id="anamespaceoperations__research_1_1glop_html_a733947145e3e1631165b618b05c9ccb7"><div class="ttname"><a href="namespaceoperations__research_1_1glop.html#a733947145e3e1631165b618b05c9ccb7">operations_research::glop::Fractional</a></div><div class="ttdeci">double Fractional</div><div class="ttdef"><b>Definition:</b> <a href="lp__types_8h_source.html#l00077">lp_types.h:77</a></div></div>
|
|
<div class="ttc" id="anamespaceoperations__research_1_1glop_html_a82b6c13b99cf0a0c7c92f4f1e44eda29"><div class="ttname"><a href="namespaceoperations__research_1_1glop.html#a82b6c13b99cf0a0c7c92f4f1e44eda29">operations_research::glop::SetSupportToFalse</a></div><div class="ttdeci">void SetSupportToFalse(const ColumnView &column, DenseBooleanColumn *b)</div><div class="ttdef"><b>Definition:</b> <a href="lp__data_2lp__utils_8cc_source.html#l00146">lp_data/lp_utils.cc:146</a></div></div>
|
|
<div class="ttc" id="anamespaceoperations__research_1_1glop_html_a95879916d90daeba91c40399ae5ddcc6"><div class="ttname"><a href="namespaceoperations__research_1_1glop.html#a95879916d90daeba91c40399ae5ddcc6">operations_research::glop::IsFinite</a></div><div class="ttdeci">bool IsFinite(Fractional value)</div><div class="ttdef"><b>Definition:</b> <a href="lp__types_8h_source.html#l00090">lp_types.h:90</a></div></div>
|
|
<div class="ttc" id="anamespaceoperations__research_1_1glop_html_a9a9b90bb0105347953a41ae1d6f4dce3"><div class="ttname"><a href="namespaceoperations__research_1_1glop.html#a9a9b90bb0105347953a41ae1d6f4dce3">operations_research::glop::IsDominated</a></div><div class="ttdeci">bool IsDominated(const ColumnView &column, const DenseColumn &radius)</div><div class="ttdef"><b>Definition:</b> <a href="lp__data_2lp__utils_8cc_source.html#l00154">lp_data/lp_utils.cc:154</a></div></div>
|
|
<div class="ttc" id="anamespaceoperations__research_1_1glop_html_aa6c552b94fa80def1d4d1ea64697afb1"><div class="ttname"><a href="namespaceoperations__research_1_1glop.html#aa6c552b94fa80def1d4d1ea64697afb1">operations_research::glop::ClearAndResizeVectorWithNonZeros</a></div><div class="ttdeci">void ClearAndResizeVectorWithNonZeros(IndexType size, ScatteredRowOrCol *v)</div><div class="ttdef"><b>Definition:</b> <a href="lp__data_2lp__utils_8h_source.html#l00278">lp_data/lp_utils.h:278</a></div></div>
|
|
<div class="ttc" id="anamespaceoperations__research_1_1glop_html_aaa803ce9366dca251925e0bdde517430"><div class="ttname"><a href="namespaceoperations__research_1_1glop.html#aaa803ce9366dca251925e0bdde517430">operations_research::glop::Transpose</a></div><div class="ttdeci">const DenseRow & Transpose(const DenseColumn &col)</div><div class="ttdef"><b>Definition:</b> <a href="lp__data_2lp__utils_8h_source.html#l00192">lp_data/lp_utils.h:192</a></div></div>
|
|
<div class="ttc" id="anamespaceoperations__research_1_1glop_html_ab65a327cfc2a74c15fa26b91f19acc64"><div class="ttname"><a href="namespaceoperations__research_1_1glop.html#ab65a327cfc2a74c15fa26b91f19acc64">operations_research::glop::ColToRowIndex</a></div><div class="ttdeci">RowIndex ColToRowIndex(ColIndex col)</div><div class="ttdef"><b>Definition:</b> <a href="lp__types_8h_source.html#l00051">lp_types.h:51</a></div></div>
|
|
<div class="ttc" id="anamespaceoperations__research_1_1glop_html_ab7ff1f9383648fd132667b1e601229c9"><div class="ttname"><a href="namespaceoperations__research_1_1glop.html#ab7ff1f9383648fd132667b1e601229c9">operations_research::glop::ChangeSign</a></div><div class="ttdeci">void ChangeSign(StrictITIVector< IndexType, Fractional > *data)</div><div class="ttdef"><b>Definition:</b> <a href="lp__data_2lp__utils_8h_source.html#l00300">lp_data/lp_utils.h:300</a></div></div>
|
|
<div class="ttc" id="anamespaceoperations__research_1_1glop_html_ac51a45a6f706c77b25cea909b9d879db"><div class="ttname"><a href="namespaceoperations__research_1_1glop.html#ac51a45a6f706c77b25cea909b9d879db">operations_research::glop::Fractionality</a></div><div class="ttdeci">static Fractional Fractionality(Fractional f)</div><div class="ttdef"><b>Definition:</b> <a href="lp__data_2lp__utils_8h_source.html#l00040">lp_data/lp_utils.h:40</a></div></div>
|
|
<div class="ttc" id="anamespaceoperations__research_1_1glop_html_ad2cb64ef0c88350c1ad4c731e0b7e183"><div class="ttname"><a href="namespaceoperations__research_1_1glop.html#ad2cb64ef0c88350c1ad4c731e0b7e183">operations_research::glop::PartialScalarProduct</a></div><div class="ttdeci">Fractional PartialScalarProduct(const DenseRowOrColumn &u, const SparseColumn &v, int max_index)</div><div class="ttdef"><b>Definition:</b> <a href="lp__data_2lp__utils_8h_source.html#l00130">lp_data/lp_utils.h:130</a></div></div>
|
|
<div class="ttc" id="anamespaceoperations__research_1_1glop_html_ad8019bac1bde0ead6ff32980cd5bff52"><div class="ttname"><a href="namespaceoperations__research_1_1glop.html#ad8019bac1bde0ead6ff32980cd5bff52">operations_research::glop::RestrictedInfinityNorm</a></div><div class="ttdeci">Fractional RestrictedInfinityNorm(const ColumnView &column, const DenseBooleanColumn &rows_to_consider, RowIndex *row_index)</div><div class="ttdef"><b>Definition:</b> <a href="lp__data_2lp__utils_8cc_source.html#l00133">lp_data/lp_utils.cc:133</a></div></div>
|
|
<div class="ttc" id="anamespaceoperations__research_1_1glop_html_ae5fa9e57c2d31b294195ed8a9e17bfac"><div class="ttname"><a href="namespaceoperations__research_1_1glop.html#ae5fa9e57c2d31b294195ed8a9e17bfac">operations_research::glop::DenseColumn</a></div><div class="ttdeci">StrictITIVector< RowIndex, Fractional > DenseColumn</div><div class="ttdef"><b>Definition:</b> <a href="lp__types_8h_source.html#l00328">lp_types.h:328</a></div></div>
|
|
<div class="ttc" id="anamespaceoperations__research_1_1glop_html_aeb4b2cc773e71eeb6b9d2b6f4c05a858"><div class="ttname"><a href="namespaceoperations__research_1_1glop.html#aeb4b2cc773e71eeb6b9d2b6f4c05a858">operations_research::glop::SumWithPositiveInfiniteAndOneMissing</a></div><div class="ttdeci">SumWithOneMissing< true > SumWithPositiveInfiniteAndOneMissing</div><div class="ttdef"><b>Definition:</b> <a href="lp__data_2lp__utils_8h_source.html#l00360">lp_data/lp_utils.h:360</a></div></div>
|
|
<div class="ttc" id="anamespaceoperations__research_1_1glop_html_aec7177dd4fab33ada14a7689223fe89f"><div class="ttname"><a href="namespaceoperations__research_1_1glop.html#aec7177dd4fab33ada14a7689223fe89f">operations_research::glop::DenseBooleanColumn</a></div><div class="ttdeci">StrictITIVector< RowIndex, bool > DenseBooleanColumn</div><div class="ttdef"><b>Definition:</b> <a href="lp__types_8h_source.html#l00331">lp_types.h:331</a></div></div>
|
|
<div class="ttc" id="anamespaceoperations__research_1_1glop_html_af9a790b7e8c5b0c6d55b336177378e78"><div class="ttname"><a href="namespaceoperations__research_1_1glop.html#af9a790b7e8c5b0c6d55b336177378e78">operations_research::glop::kInfinity</a></div><div class="ttdeci">const double kInfinity</div><div class="ttdef"><b>Definition:</b> <a href="lp__types_8h_source.html#l00083">lp_types.h:83</a></div></div>
|
|
<div class="ttc" id="anamespaceoperations__research_html"><div class="ttname"><a href="namespaceoperations__research.html">operations_research</a></div><div class="ttdoc">The vehicle routing library lets one model and solve generic vehicle routing problems ranging from th...</div><div class="ttdef"><b>Definition:</b> <a href="dense__doubly__linked__list_8h_source.html#l00021">dense_doubly_linked_list.h:21</a></div></div>
|
|
<div class="ttc" id="apack_8cc_html_a750b5d744c39a06bfb13e6eb010e35d0"><div class="ttname"><a href="pack_8cc.html#a750b5d744c39a06bfb13e6eb010e35d0">index</a></div><div class="ttdeci">int index</div><div class="ttdef"><b>Definition:</b> <a href="pack_8cc_source.html#l00508">pack.cc:508</a></div></div>
|
|
<div class="ttc" id="aparser_8yy_8cc_html_a5a634cf4429798b1c921a81de8250051"><div class="ttname"><a href="parser_8yy_8cc.html#a5a634cf4429798b1c921a81de8250051">input</a></div><div class="ttdeci">static int input(yyscan_t yyscanner)</div></div>
|
|
<div class="ttc" id="ascattered__vector_8h_html"><div class="ttname"><a href="scattered__vector_8h.html">scattered_vector.h</a></div></div>
|
|
<div class="ttc" id="asparse__column_8h_html"><div class="ttname"><a href="sparse__column_8h.html">sparse_column.h</a></div></div>
|
|
<div class="ttc" id="astructoperations__research_1_1glop_1_1_scattered_column_html"><div class="ttname"><a href="structoperations__research_1_1glop_1_1_scattered_column.html">operations_research::glop::ScatteredColumn</a></div><div class="ttdef"><b>Definition:</b> <a href="scattered__vector_8h_source.html#l00190">scattered_vector.h:191</a></div></div>
|
|
<div class="ttc" id="astructoperations__research_1_1glop_1_1_scattered_vector_html_a73934bd40690bdc3cbe36d6fbb0ecad5"><div class="ttname"><a href="structoperations__research_1_1glop_1_1_scattered_vector.html#a73934bd40690bdc3cbe36d6fbb0ecad5">operations_research::glop::ScatteredVector::ShouldUseDenseIteration</a></div><div class="ttdeci">bool ShouldUseDenseIteration(double ratio_for_using_dense_representation) const</div><div class="ttdef"><b>Definition:</b> <a href="scattered__vector_8h_source.html#l00119">scattered_vector.h:119</a></div></div>
|
|
<div class="ttc" id="astructoperations__research_1_1glop_1_1_scattered_vector_html_ae0473e18a367af671dc3f08063c80da4"><div class="ttname"><a href="structoperations__research_1_1glop_1_1_scattered_vector.html#ae0473e18a367af671dc3f08063c80da4">operations_research::glop::ScatteredVector::values</a></div><div class="ttdeci">StrictITIVector< Index, Fractional > values</div><div class="ttdef"><b>Definition:</b> <a href="scattered__vector_8h_source.html#l00057">scattered_vector.h:57</a></div></div>
|
|
</div><!-- fragment --></div><!-- contents -->
|
|
</div><!-- doc-content -->
|
|
<!-- HTML footer for doxygen 1.8.18-->
|
|
<!-- start footer part -->
|
|
<div id="nav-path" class="navpath"><!-- id is needed for treeview function! -->
|
|
<ul>
|
|
<li class="navelem"><a class="el" href="dir_a7cc1eeded8f693d0da6c729bc88c45a.html">ortools</a></li><li class="navelem"><a class="el" href="dir_903a7ff076ee7f01be86992fe717fc82.html">lp_data</a></li><li class="navelem"><a class="el" href="lp__data_2lp__utils_8h.html">lp_utils.h</a></li>
|
|
<li class="footer">Generated by
|
|
<a href="http://www.doxygen.org/index.html">
|
|
<img class="footer" src="doxygen.png" alt="doxygen"/></a> 1.9.1 </li>
|
|
</ul>
|
|
</div>
|
|
</body>
|
|
</html>
|