<!-- HTML header for doxygen 1.8.8-->
<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd">
<html xmlns="http://www.w3.org/1999/xhtml">
<head>
<meta http-equiv="X-UA-Compatible" content="IE=edge">
<!-- For Mobile Devices -->
<meta name="viewport" content="width=device-width, initial-scale=1">
<meta http-equiv="Content-Type" content="text/xhtml;charset=UTF-8"/>
<meta name="generator" content="Doxygen 1.8.17"/>
<script type="text/javascript" src="https://code.jquery.com/jquery-2.1.1.min.js"></script>
<title>QuEST: QuEST_common.c Source File</title>
<!--<link href="tabs.css" rel="stylesheet" type="text/css"/>-->
<script type="text/javascript" src="dynsections.js"></script>
<link href="doxygen.css" rel="stylesheet" type="text/css" />
<link href="style.css" rel="stylesheet" type="text/css"/>
<link href='https://fonts.googleapis.com/css?family=Roboto+Slab' rel='stylesheet' type='text/css'>
<link rel="stylesheet" href="https://maxcdn.bootstrapcdn.com/bootstrap/3.3.6/css/bootstrap.min.css">
<script src="https://maxcdn.bootstrapcdn.com/bootstrap/3.3.6/js/bootstrap.min.js"></script>
<link href="jquery.smartmenus.bootstrap.css" rel="stylesheet">
<script type="text/javascript" src="jquery.smartmenus.js"></script>
<!-- SmartMenus jQuery Bootstrap Addon -->
<script type="text/javascript" src="jquery.smartmenus.bootstrap.js"></script>
<!-- SmartMenus jQuery plugin -->
</head>
<body>
<nav class="navbar navbar-default" role="navigation">
<div class="container">
<div class="navbar-header">
<a class="navbar-brand"><img alt="Logo" src="logo.png"/></a>
<a class="navbar-brand"><b>QuEST</b> v3.4.0<br>The Quantum Exact Simulation Toolkit</a>
</div>
</div>
</nav>
<div id="top"><!-- do not remove this div, it is closed by doxygen! -->
<div class="content" id="content">
<div class="container">
<div class="row">
<div class="col-sm-12 panel " style="padding-bottom: 15px;">
<div style="margin-bottom: 15px;">
<!-- end header part -->
<!-- Generated by Doxygen 1.8.17 -->
<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('',false,false,'search.php','Search');
});
/* @license-end */</script>
<div id="main-nav"></div>
<div id="nav-path" class="navpath">
<ul>
<li class="navelem"><a class="el" href="dir_e145b6391efe0b6e575fd5bb5f76dbf8.html">QuEST</a></li><li class="navelem"><a class="el" href="dir_aaf4e0c86438000dc030ca371e863336.html">QuEST</a></li><li class="navelem"><a class="el" href="dir_288db1ec25fc0ec29a6dc915511dd336.html">src</a></li> </ul>
</div>
</div><!-- top -->
<div class="header">
<div class="headertitle">
<div class="title">QuEST_common.c</div> </div>
</div><!--header-->
<div class="contents">
<a href="QuEST__common_8c.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">// Distributed under MIT licence. See https://github.com/QuEST-Kit/QuEST/blob/master/LICENCE.txt for details</span></div>
<div class="line"><a name="l00002"></a><span class="lineno"> 2</span>  </div>
<div class="line"><a name="l00018"></a><span class="lineno"> 18</span> <span class="preprocessor"># include "<a class="code" href="QuEST_8h.html">QuEST.h</a>"</span></div>
<div class="line"><a name="l00019"></a><span class="lineno"> 19</span> <span class="preprocessor"># include "<a class="code" href="QuEST__internal_8h.html">QuEST_internal.h</a>"</span></div>
<div class="line"><a name="l00020"></a><span class="lineno"> 20</span> <span class="preprocessor"># include "<a class="code" href="QuEST__precision_8h.html">QuEST_precision.h</a>"</span></div>
<div class="line"><a name="l00021"></a><span class="lineno"> 21</span> <span class="preprocessor"># include "<a class="code" href="QuEST__validation_8h.html">QuEST_validation.h</a>"</span></div>
<div class="line"><a name="l00022"></a><span class="lineno"> 22</span> <span class="preprocessor"># include "<a class="code" href="QuEST__qasm_8h.html">QuEST_qasm.h</a>"</span></div>
<div class="line"><a name="l00023"></a><span class="lineno"> 23</span> <span class="preprocessor"># include "<a class="code" href="mt19937ar_8h.html">mt19937ar.h</a>"</span></div>
<div class="line"><a name="l00024"></a><span class="lineno"> 24</span>  </div>
<div class="line"><a name="l00025"></a><span class="lineno"> 25</span> <span class="preprocessor">#if defined(_WIN32) && ! defined(__MINGW32__)</span></div>
<div class="line"><a name="l00026"></a><span class="lineno"> 26</span> <span class="preprocessor"> #include <Windows.h></span></div>
<div class="line"><a name="l00027"></a><span class="lineno"> 27</span> <span class="preprocessor"> #include <io.h></span></div>
<div class="line"><a name="l00028"></a><span class="lineno"> 28</span> <span class="preprocessor"> #include <process.h></span></div>
<div class="line"><a name="l00029"></a><span class="lineno"> 29</span> <span class="preprocessor">#else</span></div>
<div class="line"><a name="l00030"></a><span class="lineno"> 30</span> <span class="preprocessor"> #include <unistd.h></span></div>
<div class="line"><a name="l00031"></a><span class="lineno"> 31</span> <span class="preprocessor"> #include <sys/time.h></span></div>
<div class="line"><a name="l00032"></a><span class="lineno"> 32</span> <span class="preprocessor">#endif</span></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="preprocessor"># include <sys/types.h></span> </div>
<div class="line"><a name="l00035"></a><span class="lineno"> 35</span> <span class="preprocessor"># include <stdio.h></span></div>
<div class="line"><a name="l00036"></a><span class="lineno"> 36</span> <span class="preprocessor"># include <stdlib.h></span></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>  </div>
<div class="line"><a name="l00039"></a><span class="lineno"> 39</span> <span class="comment">// expose PI on GPU build</span></div>
<div class="line"><a name="l00040"></a><span class="lineno"> 40</span> <span class="preprocessor">#ifndef M_PI </span></div>
<div class="line"><a name="l00041"></a><span class="lineno"><a class="line" href="QuEST__common_8c.html#ae71449b1cc6e6250b91f539153a7a0d3"> 41</a></span> <span class="preprocessor">#define M_PI 3.1415926535897932384626433832795028841971</span></div>
<div class="line"><a name="l00042"></a><span class="lineno"> 42</span> <span class="preprocessor">#endif</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>  </div>
<div class="line"><a name="l00045"></a><span class="lineno"> 45</span> <span class="preprocessor">#ifdef __cplusplus</span></div>
<div class="line"><a name="l00046"></a><span class="lineno"> 46</span> <span class="keyword">extern</span> <span class="stringliteral">"C"</span> {</div>
<div class="line"><a name="l00047"></a><span class="lineno"> 47</span> <span class="preprocessor">#endif</span></div>
<div class="line"><a name="l00048"></a><span class="lineno"> 48</span>  </div>
<div class="line"><a name="l00049"></a><span class="lineno"> 49</span> <span class="comment">/* builds a bit-string where 1 indicates a qubit is present in this list */</span></div>
<div class="line"><a name="l00050"></a><span class="lineno"><a class="line" href="QuEST__internal_8h.html#a0182eedcff44469771c971d8e2d73dc5"> 50</a></span> <span class="keywordtype">long</span> <span class="keywordtype">long</span> <span class="keywordtype">int</span> <a class="code" href="QuEST__common_8c.html#acc45bb8300922ea5933e294921e25ea2">getQubitBitMask</a>(<span class="keywordtype">int</span>* qubits, <span class="keywordtype">int</span> numQubits) {</div>
<div class="line"><a name="l00051"></a><span class="lineno"> 51</span>  </div>
<div class="line"><a name="l00052"></a><span class="lineno"> 52</span>  <span class="keywordtype">long</span> <span class="keywordtype">long</span> <span class="keywordtype">int</span> mask=0; </div>
<div class="line"><a name="l00053"></a><span class="lineno"> 53</span>  <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i=0; i<numQubits; i++)</div>
<div class="line"><a name="l00054"></a><span class="lineno"> 54</span>  mask = mask | (1LL << qubits[i]);</div>
<div class="line"><a name="l00055"></a><span class="lineno"> 55</span>  </div>
<div class="line"><a name="l00056"></a><span class="lineno"> 56</span>  <span class="keywordflow">return</span> mask;</div>
<div class="line"><a name="l00057"></a><span class="lineno"> 57</span> }</div>
<div class="line"><a name="l00058"></a><span class="lineno"> 58</span>  </div>
<div class="line"><a name="l00059"></a><span class="lineno"> 59</span> <span class="comment">/* builds a bit-string where 1 indicates control qubits conditioned on 0 ('flipped') */</span></div>
<div class="line"><a name="l00060"></a><span class="lineno"><a class="line" href="QuEST__internal_8h.html#a6fa6f1bc617e9fef3cd87ecc3677a02b"> 60</a></span> <span class="keywordtype">long</span> <span class="keywordtype">long</span> <span class="keywordtype">int</span> <a class="code" href="QuEST__common_8c.html#a6fa6f1bc617e9fef3cd87ecc3677a02b">getControlFlipMask</a>(<span class="keywordtype">int</span>* controlQubits, <span class="keywordtype">int</span>* controlState, <span class="keywordtype">int</span> numControlQubits) {</div>
<div class="line"><a name="l00061"></a><span class="lineno"> 61</span>  </div>
<div class="line"><a name="l00062"></a><span class="lineno"> 62</span>  <span class="keywordtype">long</span> <span class="keywordtype">long</span> <span class="keywordtype">int</span> mask=0;</div>
<div class="line"><a name="l00063"></a><span class="lineno"> 63</span>  <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i=0; i<numControlQubits; i++)</div>
<div class="line"><a name="l00064"></a><span class="lineno"> 64</span>  <span class="keywordflow">if</span> (controlState[i] == 0)</div>
<div class="line"><a name="l00065"></a><span class="lineno"> 65</span>  mask = mask | (1LL << controlQubits[i]);</div>
<div class="line"><a name="l00066"></a><span class="lineno"> 66</span>  </div>
<div class="line"><a name="l00067"></a><span class="lineno"> 67</span>  <span class="keywordflow">return</span> mask;</div>
<div class="line"><a name="l00068"></a><span class="lineno"> 68</span> }</div>
<div class="line"><a name="l00069"></a><span class="lineno"> 69</span>  </div>
<div class="line"><a name="l00070"></a><span class="lineno"><a class="line" href="QuEST__internal_8h.html#ac8ea84dd85eca525daeeb30b3c96b673"> 70</a></span> <span class="keywordtype">void</span> <a class="code" href="QuEST__common_8c.html#ac8ea84dd85eca525daeeb30b3c96b673">ensureIndsIncrease</a>(<span class="keywordtype">int</span>* ind1, <span class="keywordtype">int</span>* ind2) {</div>
<div class="line"><a name="l00071"></a><span class="lineno"> 71</span>  </div>
<div class="line"><a name="l00072"></a><span class="lineno"> 72</span>  <span class="keywordflow">if</span> (*ind1 > *ind2) {</div>
<div class="line"><a name="l00073"></a><span class="lineno"> 73</span>  <span class="keywordtype">int</span> copy = *ind1;</div>
<div class="line"><a name="l00074"></a><span class="lineno"> 74</span>  *ind1 = *ind2;</div>
<div class="line"><a name="l00075"></a><span class="lineno"> 75</span>  *ind2 = copy;</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> }</div>
<div class="line"><a name="l00078"></a><span class="lineno"> 78</span>  </div>
<div class="line"><a name="l00079"></a><span class="lineno"><a class="line" href="QuEST__internal_8h.html#a3cf54148e2162b04aeb634badb90826c"> 79</a></span> <a class="code" href="group__type.html#ga7740e349b4f8bae6451547680f0ce2d6">qreal</a> <a class="code" href="QuEST__common_8c.html#a3cf54148e2162b04aeb634badb90826c">getVectorMagnitude</a>(<a class="code" href="structVector.html">Vector</a> vec) {</div>
<div class="line"><a name="l00080"></a><span class="lineno"> 80</span>  </div>
<div class="line"><a name="l00081"></a><span class="lineno"> 81</span>  <span class="keywordflow">return</span> sqrt(vec.<a class="code" href="structVector.html#a0d0459a5b731ac5fb3dff2c6134e6e88">x</a>*vec.<a class="code" href="structVector.html#a0d0459a5b731ac5fb3dff2c6134e6e88">x</a> + vec.<a class="code" href="structVector.html#abc3cfc6c8e297b1153a373d84b981e6b">y</a>*vec.<a class="code" href="structVector.html#abc3cfc6c8e297b1153a373d84b981e6b">y</a> + vec.<a class="code" href="structVector.html#ab33976d8fda2e44b770f71269d698a25">z</a>*vec.<a class="code" href="structVector.html#ab33976d8fda2e44b770f71269d698a25">z</a>);</div>
<div class="line"><a name="l00082"></a><span class="lineno"> 82</span> }</div>
<div class="line"><a name="l00083"></a><span class="lineno"> 83</span>  </div>
<div class="line"><a name="l00084"></a><span class="lineno"><a class="line" href="QuEST__common_8c.html#a041bb75200ee316de7b90667a4fe04f3"> 84</a></span> <a class="code" href="structVector.html">Vector</a> <a class="code" href="QuEST__common_8c.html#a041bb75200ee316de7b90667a4fe04f3">getUnitVector</a>(<a class="code" href="structVector.html">Vector</a> vec) {</div>
<div class="line"><a name="l00085"></a><span class="lineno"> 85</span>  </div>
<div class="line"><a name="l00086"></a><span class="lineno"> 86</span>  <a class="code" href="group__type.html#ga7740e349b4f8bae6451547680f0ce2d6">qreal</a> mag = <a class="code" href="QuEST__common_8c.html#a3cf54148e2162b04aeb634badb90826c">getVectorMagnitude</a>(vec);</div>
<div class="line"><a name="l00087"></a><span class="lineno"> 87</span>  <a class="code" href="structVector.html">Vector</a> unitVec = (<a class="code" href="structVector.html">Vector</a>) {.<a class="code" href="structVector.html#a0d0459a5b731ac5fb3dff2c6134e6e88">x</a>=vec.<a class="code" href="structVector.html#a0d0459a5b731ac5fb3dff2c6134e6e88">x</a>/mag, .y=vec.<a class="code" href="structVector.html#abc3cfc6c8e297b1153a373d84b981e6b">y</a>/mag, .z=vec.<a class="code" href="structVector.html#ab33976d8fda2e44b770f71269d698a25">z</a>/mag};</div>
<div class="line"><a name="l00088"></a><span class="lineno"> 88</span>  <span class="keywordflow">return</span> unitVec;</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"><a class="line" href="QuEST__internal_8h.html#af92fa737f4491324408b3d15543dced4"> 91</a></span> <a class="code" href="structComplex.html">Complex</a> <a class="code" href="QuEST__common_8c.html#af92fa737f4491324408b3d15543dced4">getConjugateScalar</a>(<a class="code" href="structComplex.html">Complex</a> scalar) {</div>
<div class="line"><a name="l00092"></a><span class="lineno"> 92</span>  </div>
<div class="line"><a name="l00093"></a><span class="lineno"> 93</span>  <a class="code" href="structComplex.html">Complex</a> conjScalar;</div>
<div class="line"><a name="l00094"></a><span class="lineno"> 94</span>  conjScalar.<a class="code" href="structComplex.html#ab5b2e2eca02131fc74b289a83636cbe3">real</a> = scalar.<a class="code" href="structComplex.html#ab5b2e2eca02131fc74b289a83636cbe3">real</a>;</div>
<div class="line"><a name="l00095"></a><span class="lineno"> 95</span>  conjScalar.<a class="code" href="structComplex.html#a84f5439aad0ef495efdd3b4c1c02d27e">imag</a> = - scalar.<a class="code" href="structComplex.html#a84f5439aad0ef495efdd3b4c1c02d27e">imag</a>;</div>
<div class="line"><a name="l00096"></a><span class="lineno"> 96</span>  <span class="keywordflow">return</span> conjScalar;</div>
<div class="line"><a name="l00097"></a><span class="lineno"> 97</span> }</div>
<div class="line"><a name="l00098"></a><span class="lineno"> 98</span>  </div>
<div class="line"><a name="l00099"></a><span class="lineno"><a class="line" href="QuEST__common_8c.html#ae2b222fda52c10fb420cb3ac78ca31e4"> 99</a></span> <span class="preprocessor">#define macro_setConjugateMatrix(dest, src, dim) \</span></div>
<div class="line"><a name="l00100"></a><span class="lineno"> 100</span> <span class="preprocessor"> for (int i=0; i<dim; i++) \</span></div>
<div class="line"><a name="l00101"></a><span class="lineno"> 101</span> <span class="preprocessor"> for (int j=0; j<dim; j++) { \</span></div>
<div class="line"><a name="l00102"></a><span class="lineno"> 102</span> <span class="preprocessor"> dest.real[i][j] = src.real[i][j]; \</span></div>
<div class="line"><a name="l00103"></a><span class="lineno"> 103</span> <span class="preprocessor"> dest.imag[i][j] = - src.imag[i][j]; </span><span class="comment">/* negative for conjugate */</span><span class="preprocessor"> \</span></div>
<div class="line"><a name="l00104"></a><span class="lineno"> 104</span> <span class="preprocessor"> } </span></div>
<div class="line"><a name="l00105"></a><span class="lineno"><a class="line" href="QuEST__internal_8h.html#ad34d98eeaa6a3cb45e92c696cfd2d1c8"> 105</a></span> <a class="code" href="structComplexMatrix2.html">ComplexMatrix2</a> <a class="code" href="QuEST__common_8c.html#ad34d98eeaa6a3cb45e92c696cfd2d1c8">getConjugateMatrix2</a>(<a class="code" href="structComplexMatrix2.html">ComplexMatrix2</a> src) {</div>
<div class="line"><a name="l00106"></a><span class="lineno"> 106</span>  <a class="code" href="structComplexMatrix2.html">ComplexMatrix2</a> conj;</div>
<div class="line"><a name="l00107"></a><span class="lineno"> 107</span>  <a class="code" href="QuEST__common_8c.html#ae2b222fda52c10fb420cb3ac78ca31e4">macro_setConjugateMatrix</a>(conj, src, 2);</div>
<div class="line"><a name="l00108"></a><span class="lineno"> 108</span>  <span class="keywordflow">return</span> conj;</div>
<div class="line"><a name="l00109"></a><span class="lineno"> 109</span> }</div>
<div class="line"><a name="l00110"></a><span class="lineno"><a class="line" href="QuEST__internal_8h.html#a20d45e85dd49dd3bf6045990b590b24b"> 110</a></span> <a class="code" href="structComplexMatrix4.html">ComplexMatrix4</a> <a class="code" href="QuEST__common_8c.html#a20d45e85dd49dd3bf6045990b590b24b">getConjugateMatrix4</a>(<a class="code" href="structComplexMatrix4.html">ComplexMatrix4</a> src) {</div>
<div class="line"><a name="l00111"></a><span class="lineno"> 111</span>  <a class="code" href="structComplexMatrix4.html">ComplexMatrix4</a> conj;</div>
<div class="line"><a name="l00112"></a><span class="lineno"> 112</span>  <a class="code" href="QuEST__common_8c.html#ae2b222fda52c10fb420cb3ac78ca31e4">macro_setConjugateMatrix</a>(conj, src, 4);</div>
<div class="line"><a name="l00113"></a><span class="lineno"> 113</span>  <span class="keywordflow">return</span> conj;</div>
<div class="line"><a name="l00114"></a><span class="lineno"> 114</span> }</div>
<div class="line"><a name="l00115"></a><span class="lineno"><a class="line" href="QuEST__internal_8h.html#a9fd917537204e074119d4f283dea78e3"> 115</a></span> <span class="keywordtype">void</span> <a class="code" href="QuEST__common_8c.html#a9fd917537204e074119d4f283dea78e3">setConjugateMatrixN</a>(<a class="code" href="structComplexMatrixN.html">ComplexMatrixN</a> m) {</div>
<div class="line"><a name="l00116"></a><span class="lineno"> 116</span>  <span class="keywordtype">int</span> len = 1 << m.<a class="code" href="structComplexMatrixN.html#a606f7afbc3f65a596bdca21cb43be6a8">numQubits</a>;</div>
<div class="line"><a name="l00117"></a><span class="lineno"> 117</span>  <a class="code" href="QuEST__common_8c.html#ae2b222fda52c10fb420cb3ac78ca31e4">macro_setConjugateMatrix</a>(m, m, len);</div>
<div class="line"><a name="l00118"></a><span class="lineno"> 118</span> }</div>
<div class="line"><a name="l00119"></a><span class="lineno"> 119</span>  </div>
<div class="line"><a name="l00120"></a><span class="lineno"><a class="line" href="QuEST__internal_8h.html#a71b3ae4402c1c5039d3086ca0e47a522"> 120</a></span> <span class="keywordtype">void</span> <a class="code" href="QuEST__common_8c.html#a71b3ae4402c1c5039d3086ca0e47a522">getComplexPairFromRotation</a>(<a class="code" href="group__type.html#ga7740e349b4f8bae6451547680f0ce2d6">qreal</a> angle, <a class="code" href="structVector.html">Vector</a> axis, <a class="code" href="structComplex.html">Complex</a>* alpha, <a class="code" href="structComplex.html">Complex</a>* beta) {</div>
<div class="line"><a name="l00121"></a><span class="lineno"> 121</span>  </div>
<div class="line"><a name="l00122"></a><span class="lineno"> 122</span>  <a class="code" href="structVector.html">Vector</a> unitAxis = <a class="code" href="QuEST__common_8c.html#a041bb75200ee316de7b90667a4fe04f3">getUnitVector</a>(axis);</div>
<div class="line"><a name="l00123"></a><span class="lineno"> 123</span>  alpha-><a class="code" href="structComplex.html#ab5b2e2eca02131fc74b289a83636cbe3">real</a> = cos(angle/2.0);</div>
<div class="line"><a name="l00124"></a><span class="lineno"> 124</span>  alpha-><a class="code" href="structComplex.html#a84f5439aad0ef495efdd3b4c1c02d27e">imag</a> = - sin(angle/2.0)*unitAxis.<a class="code" href="structVector.html#ab33976d8fda2e44b770f71269d698a25">z</a>; </div>
<div class="line"><a name="l00125"></a><span class="lineno"> 125</span>  beta-><a class="code" href="structComplex.html#ab5b2e2eca02131fc74b289a83636cbe3">real</a> = sin(angle/2.0)*unitAxis.<a class="code" href="structVector.html#abc3cfc6c8e297b1153a373d84b981e6b">y</a>;</div>
<div class="line"><a name="l00126"></a><span class="lineno"> 126</span>  beta-><a class="code" href="structComplex.html#a84f5439aad0ef495efdd3b4c1c02d27e">imag</a> = - sin(angle/2.0)*unitAxis.<a class="code" href="structVector.html#a0d0459a5b731ac5fb3dff2c6134e6e88">x</a>;</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>  </div>
<div class="line"><a name="l00130"></a><span class="lineno"><a class="line" href="QuEST__internal_8h.html#a404d177b894ce008289e44a599ad77d6"> 130</a></span> <span class="keywordtype">void</span> <a class="code" href="QuEST__common_8c.html#a404d177b894ce008289e44a599ad77d6">getZYZRotAnglesFromComplexPair</a>(<a class="code" href="structComplex.html">Complex</a> alpha, <a class="code" href="structComplex.html">Complex</a> beta, <a class="code" href="group__type.html#ga7740e349b4f8bae6451547680f0ce2d6">qreal</a>* rz2, <a class="code" href="group__type.html#ga7740e349b4f8bae6451547680f0ce2d6">qreal</a>* ry, <a class="code" href="group__type.html#ga7740e349b4f8bae6451547680f0ce2d6">qreal</a>* rz1) {</div>
<div class="line"><a name="l00131"></a><span class="lineno"> 131</span>  </div>
<div class="line"><a name="l00132"></a><span class="lineno"> 132</span>  <a class="code" href="group__type.html#ga7740e349b4f8bae6451547680f0ce2d6">qreal</a> alphaMag = sqrt(alpha.<a class="code" href="structComplex.html#ab5b2e2eca02131fc74b289a83636cbe3">real</a>*alpha.<a class="code" href="structComplex.html#ab5b2e2eca02131fc74b289a83636cbe3">real</a> + alpha.<a class="code" href="structComplex.html#a84f5439aad0ef495efdd3b4c1c02d27e">imag</a>*alpha.<a class="code" href="structComplex.html#a84f5439aad0ef495efdd3b4c1c02d27e">imag</a>);</div>
<div class="line"><a name="l00133"></a><span class="lineno"> 133</span>  *ry = 2.0 * acos(alphaMag);</div>
<div class="line"><a name="l00134"></a><span class="lineno"> 134</span>  </div>
<div class="line"><a name="l00135"></a><span class="lineno"> 135</span>  <a class="code" href="group__type.html#ga7740e349b4f8bae6451547680f0ce2d6">qreal</a> alphaPhase = atan2(alpha.<a class="code" href="structComplex.html#a84f5439aad0ef495efdd3b4c1c02d27e">imag</a>, alpha.<a class="code" href="structComplex.html#ab5b2e2eca02131fc74b289a83636cbe3">real</a>);</div>
<div class="line"><a name="l00136"></a><span class="lineno"> 136</span>  <a class="code" href="group__type.html#ga7740e349b4f8bae6451547680f0ce2d6">qreal</a> betaPhase = atan2(beta.<a class="code" href="structComplex.html#a84f5439aad0ef495efdd3b4c1c02d27e">imag</a>, beta.<a class="code" href="structComplex.html#ab5b2e2eca02131fc74b289a83636cbe3">real</a>);</div>
<div class="line"><a name="l00137"></a><span class="lineno"> 137</span>  *rz2 = - alphaPhase + betaPhase;</div>
<div class="line"><a name="l00138"></a><span class="lineno"> 138</span>  *rz1 = - alphaPhase - betaPhase;</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>  </div>
<div class="line"><a name="l00142"></a><span class="lineno"><a class="line" href="QuEST__internal_8h.html#a73300e3cdca56a7615980274905ba097"> 142</a></span> <span class="keywordtype">void</span> <a class="code" href="QuEST__common_8c.html#a73300e3cdca56a7615980274905ba097">getComplexPairAndPhaseFromUnitary</a>(<a class="code" href="structComplexMatrix2.html">ComplexMatrix2</a> u, <a class="code" href="structComplex.html">Complex</a>* alpha, <a class="code" href="structComplex.html">Complex</a>* beta, <a class="code" href="group__type.html#ga7740e349b4f8bae6451547680f0ce2d6">qreal</a>* globalPhase) {</div>
<div class="line"><a name="l00143"></a><span class="lineno"> 143</span>  </div>
<div class="line"><a name="l00144"></a><span class="lineno"> 144</span>  <a class="code" href="group__type.html#ga7740e349b4f8bae6451547680f0ce2d6">qreal</a> r0c0Phase = atan2(u.<a class="code" href="structComplexMatrix2.html#a3a2b60a9e00ccef054a6de3b9307c1be">imag</a>[0][0], u.<a class="code" href="structComplexMatrix2.html#a3ea98d545b67a190ee6abe936c72076f">real</a>[0][0]);</div>
<div class="line"><a name="l00145"></a><span class="lineno"> 145</span>  <a class="code" href="group__type.html#ga7740e349b4f8bae6451547680f0ce2d6">qreal</a> r1c1Phase = atan2(u.<a class="code" href="structComplexMatrix2.html#a3a2b60a9e00ccef054a6de3b9307c1be">imag</a>[1][1], u.<a class="code" href="structComplexMatrix2.html#a3ea98d545b67a190ee6abe936c72076f">real</a>[1][1]);</div>
<div class="line"><a name="l00146"></a><span class="lineno"> 146</span>  *globalPhase = (r0c0Phase + r1c1Phase)/2.0;</div>
<div class="line"><a name="l00147"></a><span class="lineno"> 147</span>  </div>
<div class="line"><a name="l00148"></a><span class="lineno"> 148</span>  <a class="code" href="group__type.html#ga7740e349b4f8bae6451547680f0ce2d6">qreal</a> cosPhase = cos(*globalPhase);</div>
<div class="line"><a name="l00149"></a><span class="lineno"> 149</span>  <a class="code" href="group__type.html#ga7740e349b4f8bae6451547680f0ce2d6">qreal</a> sinPhase = sin(*globalPhase);</div>
<div class="line"><a name="l00150"></a><span class="lineno"> 150</span>  alpha-><a class="code" href="structComplex.html#ab5b2e2eca02131fc74b289a83636cbe3">real</a> = u.<a class="code" href="structComplexMatrix2.html#a3ea98d545b67a190ee6abe936c72076f">real</a>[0][0]*cosPhase + u.<a class="code" href="structComplexMatrix2.html#a3a2b60a9e00ccef054a6de3b9307c1be">imag</a>[0][0]*sinPhase;</div>
<div class="line"><a name="l00151"></a><span class="lineno"> 151</span>  alpha-><a class="code" href="structComplex.html#a84f5439aad0ef495efdd3b4c1c02d27e">imag</a> = u.<a class="code" href="structComplexMatrix2.html#a3a2b60a9e00ccef054a6de3b9307c1be">imag</a>[0][0]*cosPhase - u.<a class="code" href="structComplexMatrix2.html#a3ea98d545b67a190ee6abe936c72076f">real</a>[0][0]*sinPhase;</div>
<div class="line"><a name="l00152"></a><span class="lineno"> 152</span>  beta-><a class="code" href="structComplex.html#ab5b2e2eca02131fc74b289a83636cbe3">real</a> = u.<a class="code" href="structComplexMatrix2.html#a3ea98d545b67a190ee6abe936c72076f">real</a>[1][0]*cosPhase + u.<a class="code" href="structComplexMatrix2.html#a3a2b60a9e00ccef054a6de3b9307c1be">imag</a>[1][0]*sinPhase;</div>
<div class="line"><a name="l00153"></a><span class="lineno"> 153</span>  beta-><a class="code" href="structComplex.html#a84f5439aad0ef495efdd3b4c1c02d27e">imag</a> = u.<a class="code" href="structComplexMatrix2.html#a3a2b60a9e00ccef054a6de3b9307c1be">imag</a>[1][0]*cosPhase - u.<a class="code" href="structComplexMatrix2.html#a3ea98d545b67a190ee6abe936c72076f">real</a>[1][0]*sinPhase;</div>
<div class="line"><a name="l00154"></a><span class="lineno"> 154</span> }</div>
<div class="line"><a name="l00155"></a><span class="lineno"> 155</span>  </div>
<div class="line"><a name="l00156"></a><span class="lineno"><a class="line" href="QuEST__internal_8h.html#acf1c651a4359a370a2d9b3e1e1ebd430"> 156</a></span> <span class="keywordtype">void</span> <a class="code" href="QuEST__common_8c.html#acf1c651a4359a370a2d9b3e1e1ebd430">shiftIndices</a>(<span class="keywordtype">int</span>* indices, <span class="keywordtype">int</span> numIndices, <span class="keywordtype">int</span> shift) {</div>
<div class="line"><a name="l00157"></a><span class="lineno"> 157</span>  <span class="keywordflow">for</span> (<span class="keywordtype">int</span> j=0; j < numIndices; j++)</div>
<div class="line"><a name="l00158"></a><span class="lineno"> 158</span>  indices[j] += shift;</div>
<div class="line"><a name="l00159"></a><span class="lineno"> 159</span> }</div>
<div class="line"><a name="l00160"></a><span class="lineno"> 160</span>  </div>
<div class="line"><a name="l00161"></a><span class="lineno"><a class="line" href="QuEST__internal_8h.html#aecc3c36b275c53b321de1611c2f06f78"> 161</a></span> <span class="keywordtype">void</span> <a class="code" href="QuEST__common_8c.html#aecc3c36b275c53b321de1611c2f06f78">shiftSubregIndices</a>(<span class="keywordtype">int</span>* allInds, <span class="keywordtype">int</span>* numIndsPerReg, <span class="keywordtype">int</span> numRegs, <span class="keywordtype">int</span> shift) {</div>
<div class="line"><a name="l00162"></a><span class="lineno"> 162</span>  <span class="keywordtype">int</span> i=0;</div>
<div class="line"><a name="l00163"></a><span class="lineno"> 163</span>  <span class="keywordflow">for</span> (<span class="keywordtype">int</span> r=0; r<numRegs; r++)</div>
<div class="line"><a name="l00164"></a><span class="lineno"> 164</span>  <span class="keywordflow">for</span> (<span class="keywordtype">int</span> j=0; j<numIndsPerReg[r]; j++)</div>
<div class="line"><a name="l00165"></a><span class="lineno"> 165</span>  allInds[i++] += shift;</div>
<div class="line"><a name="l00166"></a><span class="lineno"> 166</span> }</div>
<div class="line"><a name="l00167"></a><span class="lineno"> 167</span>  </div>
<div class="line"><a name="l00168"></a><span class="lineno"><a class="line" href="QuEST__common_8c.html#a9d19fc690afb425d0df5feb09be567a3"> 168</a></span> <span class="keywordtype">int</span> <a class="code" href="QuEST__common_8c.html#a9d19fc690afb425d0df5feb09be567a3">generateMeasurementOutcome</a>(<a class="code" href="group__type.html#ga7740e349b4f8bae6451547680f0ce2d6">qreal</a> zeroProb, <a class="code" href="group__type.html#ga7740e349b4f8bae6451547680f0ce2d6">qreal</a> *outcomeProb) {</div>
<div class="line"><a name="l00169"></a><span class="lineno"> 169</span>  </div>
<div class="line"><a name="l00170"></a><span class="lineno"> 170</span>  <span class="comment">// randomly choose outcome</span></div>
<div class="line"><a name="l00171"></a><span class="lineno"> 171</span>  <span class="keywordtype">int</span> outcome;</div>
<div class="line"><a name="l00172"></a><span class="lineno"> 172</span>  <span class="keywordflow">if</span> (zeroProb < REAL_EPS) </div>
<div class="line"><a name="l00173"></a><span class="lineno"> 173</span>  outcome = 1;</div>
<div class="line"><a name="l00174"></a><span class="lineno"> 174</span>  <span class="keywordflow">else</span> <span class="keywordflow">if</span> (1-zeroProb < REAL_EPS) </div>
<div class="line"><a name="l00175"></a><span class="lineno"> 175</span>  outcome = 0;</div>
<div class="line"><a name="l00176"></a><span class="lineno"> 176</span>  <span class="keywordflow">else</span></div>
<div class="line"><a name="l00177"></a><span class="lineno"> 177</span>  outcome = (<a class="code" href="mt19937ar_8c.html#ac94ab75771800274ed1a2bedeca86f04">genrand_real1</a>() > zeroProb);</div>
<div class="line"><a name="l00178"></a><span class="lineno"> 178</span>  </div>
<div class="line"><a name="l00179"></a><span class="lineno"> 179</span>  <span class="comment">// set probability of outcome</span></div>
<div class="line"><a name="l00180"></a><span class="lineno"> 180</span>  *outcomeProb = (outcome==0)? zeroProb : 1-zeroProb;</div>
<div class="line"><a name="l00181"></a><span class="lineno"> 181</span>  </div>
<div class="line"><a name="l00182"></a><span class="lineno"> 182</span>  <span class="keywordflow">return</span> outcome;</div>
<div class="line"><a name="l00183"></a><span class="lineno"> 183</span> }</div>
<div class="line"><a name="l00184"></a><span class="lineno"> 184</span>  </div>
<div class="line"><a name="l00185"></a><span class="lineno"><a class="line" href="QuEST__internal_8h.html#ab76254cfde16f0808476649507a1a2fc"> 185</a></span> <span class="keywordtype">unsigned</span> <span class="keywordtype">long</span> <span class="keywordtype">int</span> <a class="code" href="QuEST__common_8c.html#ab76254cfde16f0808476649507a1a2fc">hashString</a>(<span class="keywordtype">char</span> *str){</div>
<div class="line"><a name="l00186"></a><span class="lineno"> 186</span>  <span class="keywordtype">unsigned</span> <span class="keywordtype">long</span> <span class="keywordtype">int</span> hash = 5381;</div>
<div class="line"><a name="l00187"></a><span class="lineno"> 187</span>  <span class="keywordtype">int</span> c;</div>
<div class="line"><a name="l00188"></a><span class="lineno"> 188</span>  </div>
<div class="line"><a name="l00189"></a><span class="lineno"> 189</span>  <span class="keywordflow">while</span> ((c = *str++))</div>
<div class="line"><a name="l00190"></a><span class="lineno"> 190</span>  hash = ((hash << 5) + hash) + c; <span class="comment">/* hash * 33 + c */</span></div>
<div class="line"><a name="l00191"></a><span class="lineno"> 191</span>  </div>
<div class="line"><a name="l00192"></a><span class="lineno"> 192</span>  <span class="keywordflow">return</span> hash; </div>
<div class="line"><a name="l00193"></a><span class="lineno"> 193</span> }</div>
<div class="line"><a name="l00194"></a><span class="lineno"> 194</span>  </div>
<div class="line"><a name="l00195"></a><span class="lineno"><a class="line" href="QuEST__internal_8h.html#a2b323b65cd6512720e6aa27c182fcea0"> 195</a></span> <span class="keywordtype">void</span> <a class="code" href="QuEST__common_8c.html#a2b323b65cd6512720e6aa27c182fcea0">getQuESTDefaultSeedKey</a>(<span class="keywordtype">unsigned</span> <span class="keywordtype">long</span> <span class="keywordtype">int</span> *key){</div>
<div class="line"><a name="l00196"></a><span class="lineno"> 196</span>  <span class="comment">// init MT random number generator with two keys -- time and pid</span></div>
<div class="line"><a name="l00197"></a><span class="lineno"> 197</span>  <span class="comment">// for the MPI version, it is ok that all procs will get the same seed as random numbers will only be </span></div>
<div class="line"><a name="l00198"></a><span class="lineno"> 198</span>  <span class="comment">// used by the master process</span></div>
<div class="line"><a name="l00199"></a><span class="lineno"> 199</span> <span class="preprocessor">#if defined(_WIN32) && ! defined(__MINGW32__)</span></div>
<div class="line"><a name="l00200"></a><span class="lineno"> 200</span>  </div>
<div class="line"><a name="l00201"></a><span class="lineno"> 201</span>  <span class="keywordtype">unsigned</span> <span class="keywordtype">long</span> <span class="keywordtype">int</span> pid = (<span class="keywordtype">unsigned</span> <span class="keywordtype">long</span> int) _getpid();</div>
<div class="line"><a name="l00202"></a><span class="lineno"> 202</span>  <span class="keywordtype">unsigned</span> <span class="keywordtype">long</span> <span class="keywordtype">int</span> msecs = (<span class="keywordtype">unsigned</span> <span class="keywordtype">long</span> int) GetTickCount64();</div>
<div class="line"><a name="l00203"></a><span class="lineno"> 203</span>  </div>
<div class="line"><a name="l00204"></a><span class="lineno"> 204</span>  key[0] = msecs; key[1] = pid;</div>
<div class="line"><a name="l00205"></a><span class="lineno"> 205</span> <span class="preprocessor">#else</span></div>
<div class="line"><a name="l00206"></a><span class="lineno"> 206</span>  <span class="keyword">struct </span>timeval tv;</div>
<div class="line"><a name="l00207"></a><span class="lineno"> 207</span>  gettimeofday(&tv, NULL);</div>
<div class="line"><a name="l00208"></a><span class="lineno"> 208</span>  </div>
<div class="line"><a name="l00209"></a><span class="lineno"> 209</span>  <span class="keywordtype">double</span> time_in_mill =</div>
<div class="line"><a name="l00210"></a><span class="lineno"> 210</span>  (tv.tv_sec) * 1000 + (tv.tv_usec) / 1000 ; <span class="comment">// convert tv_sec & tv_usec to millisecond</span></div>
<div class="line"><a name="l00211"></a><span class="lineno"> 211</span>  </div>
<div class="line"><a name="l00212"></a><span class="lineno"> 212</span>  <span class="keywordtype">unsigned</span> <span class="keywordtype">long</span> <span class="keywordtype">int</span> pid = getpid();</div>
<div class="line"><a name="l00213"></a><span class="lineno"> 213</span>  <span class="keywordtype">unsigned</span> <span class="keywordtype">long</span> <span class="keywordtype">int</span> msecs = (<span class="keywordtype">unsigned</span> <span class="keywordtype">long</span> int) time_in_mill;</div>
<div class="line"><a name="l00214"></a><span class="lineno"> 214</span>  </div>
<div class="line"><a name="l00215"></a><span class="lineno"> 215</span>  key[0] = msecs; key[1] = pid;</div>
<div class="line"><a name="l00216"></a><span class="lineno"> 216</span> <span class="preprocessor">#endif </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"><a class="line" href="group__debug.html#gacd5f0ee488157837a8483c1515817f98"> 219</a></span> <span class="keywordtype">void</span> <a class="code" href="group__debug.html#gacd5f0ee488157837a8483c1515817f98">reportState</a>(<a class="code" href="structQureg.html">Qureg</a> qureg){</div>
<div class="line"><a name="l00220"></a><span class="lineno"> 220</span>  FILE *state;</div>
<div class="line"><a name="l00221"></a><span class="lineno"> 221</span>  <span class="keywordtype">char</span> filename[100];</div>
<div class="line"><a name="l00222"></a><span class="lineno"> 222</span>  <span class="keywordtype">long</span> <span class="keywordtype">long</span> <span class="keywordtype">int</span> index;</div>
<div class="line"><a name="l00223"></a><span class="lineno"> 223</span>  sprintf(filename, <span class="stringliteral">"state_rank_%d.csv"</span>, qureg.<a class="code" href="structQureg.html#ac2929e681d3d95591c18cf168dbbe4f0">chunkId</a>);</div>
<div class="line"><a name="l00224"></a><span class="lineno"> 224</span>  state = fopen(filename, <span class="stringliteral">"w"</span>);</div>
<div class="line"><a name="l00225"></a><span class="lineno"> 225</span>  <span class="keywordflow">if</span> (qureg.<a class="code" href="structQureg.html#ac2929e681d3d95591c18cf168dbbe4f0">chunkId</a>==0) fprintf(state, <span class="stringliteral">"real, imag\n"</span>);</div>
<div class="line"><a name="l00226"></a><span class="lineno"> 226</span>  </div>
<div class="line"><a name="l00227"></a><span class="lineno"> 227</span>  <span class="keywordflow">for</span>(index=0; index<qureg.<a class="code" href="structQureg.html#ab0ea0358482b62f43fdd781469607d97">numAmpsPerChunk</a>; index++){</div>
<div class="line"><a name="l00228"></a><span class="lineno"> 228</span> <span class="preprocessor"> # if QuEST_PREC==1 || QuEST_PREC==2</span></div>
<div class="line"><a name="l00229"></a><span class="lineno"> 229</span>  fprintf(state, <span class="stringliteral">"%.12f, %.12f\n"</span>, qureg.<a class="code" href="structQureg.html#a441e4cacef6bd17adb9813c7442d42fe">stateVec</a>.real[index], qureg.<a class="code" href="structQureg.html#a441e4cacef6bd17adb9813c7442d42fe">stateVec</a>.imag[index]);</div>
<div class="line"><a name="l00230"></a><span class="lineno"> 230</span> <span class="preprocessor"> # elif QuEST_PREC == 4</span></div>
<div class="line"><a name="l00231"></a><span class="lineno"> 231</span>  fprintf(state, <span class="stringliteral">"%.12Lf, %.12Lf\n"</span>, qureg.<a class="code" href="structQureg.html#a441e4cacef6bd17adb9813c7442d42fe">stateVec</a>.real[index], qureg.<a class="code" href="structQureg.html#a441e4cacef6bd17adb9813c7442d42fe">stateVec</a>.imag[index]);</div>
<div class="line"><a name="l00232"></a><span class="lineno"> 232</span> <span class="preprocessor"> #endif</span></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>  fclose(state);</div>
<div class="line"><a name="l00235"></a><span class="lineno"> 235</span> }</div>
<div class="line"><a name="l00236"></a><span class="lineno"> 236</span>  </div>
<div class="line"><a name="l00237"></a><span class="lineno"><a class="line" href="group__debug.html#ga97d96af7c7ea7b31e32cbe3b25377e09"> 237</a></span> <span class="keywordtype">void</span> <a class="code" href="group__debug.html#ga97d96af7c7ea7b31e32cbe3b25377e09">reportQuregParams</a>(<a class="code" href="structQureg.html">Qureg</a> qureg){</div>
<div class="line"><a name="l00238"></a><span class="lineno"> 238</span>  <span class="keywordtype">long</span> <span class="keywordtype">long</span> <span class="keywordtype">int</span> numAmps = 1LL << qureg.<a class="code" href="structQureg.html#a0a86e8a50ea8c998b4a0a6640c5a1218">numQubitsInStateVec</a>;</div>
<div class="line"><a name="l00239"></a><span class="lineno"> 239</span>  <span class="keywordtype">long</span> <span class="keywordtype">long</span> <span class="keywordtype">int</span> numAmpsPerRank = numAmps/qureg.<a class="code" href="structQureg.html#aa3deb7163aab90a61e755a1f200f5413">numChunks</a>;</div>
<div class="line"><a name="l00240"></a><span class="lineno"> 240</span>  <span class="keywordflow">if</span> (qureg.<a class="code" href="structQureg.html#ac2929e681d3d95591c18cf168dbbe4f0">chunkId</a>==0){</div>
<div class="line"><a name="l00241"></a><span class="lineno"> 241</span>  printf(<span class="stringliteral">"QUBITS:\n"</span>);</div>
<div class="line"><a name="l00242"></a><span class="lineno"> 242</span>  printf(<span class="stringliteral">"Number of qubits is %d.\n"</span>, qureg.<a class="code" href="structQureg.html#a0a86e8a50ea8c998b4a0a6640c5a1218">numQubitsInStateVec</a>);</div>
<div class="line"><a name="l00243"></a><span class="lineno"> 243</span>  printf(<span class="stringliteral">"Number of amps is %lld.\n"</span>, numAmps);</div>
<div class="line"><a name="l00244"></a><span class="lineno"> 244</span>  printf(<span class="stringliteral">"Number of amps per rank is %lld.\n"</span>, numAmpsPerRank);</div>
<div class="line"><a name="l00245"></a><span class="lineno"> 245</span>  }</div>
<div class="line"><a name="l00246"></a><span class="lineno"> 246</span> }</div>
<div class="line"><a name="l00247"></a><span class="lineno"> 247</span>  </div>
<div class="line"><a name="l00248"></a><span class="lineno"><a class="line" href="QuEST__internal_8h.html#a3afe58e49dce5536bc4c1b7caddf44a6"> 248</a></span> <a class="code" href="group__type.html#ga7740e349b4f8bae6451547680f0ce2d6">qreal</a> <a class="code" href="QuEST__common_8c.html#a3afe58e49dce5536bc4c1b7caddf44a6">statevec_getProbAmp</a>(<a class="code" href="structQureg.html">Qureg</a> qureg, <span class="keywordtype">long</span> <span class="keywordtype">long</span> <span class="keywordtype">int</span> index){</div>
<div class="line"><a name="l00249"></a><span class="lineno"> 249</span>  <a class="code" href="group__type.html#ga7740e349b4f8bae6451547680f0ce2d6">qreal</a> real = <a class="code" href="QuEST__internal_8h.html#abc9a9ef4344c7faaaf28ac25c76649b9">statevec_getRealAmp</a>(qureg, index);</div>
<div class="line"><a name="l00250"></a><span class="lineno"> 250</span>  <a class="code" href="group__type.html#ga7740e349b4f8bae6451547680f0ce2d6">qreal</a> imag = <a class="code" href="QuEST__internal_8h.html#abd509244d57657e148e4084c5ab5d28f">statevec_getImagAmp</a>(qureg, index);</div>
<div class="line"><a name="l00251"></a><span class="lineno"> 251</span>  <span class="keywordflow">return</span> real*real + imag*imag;</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>  </div>
<div class="line"><a name="l00254"></a><span class="lineno"><a class="line" href="QuEST__internal_8h.html#a5ef139196852f92fe8271c82964eadff"> 254</a></span> <span class="keywordtype">void</span> <a class="code" href="QuEST__common_8c.html#a5ef139196852f92fe8271c82964eadff">statevec_phaseShift</a>(<a class="code" href="structQureg.html">Qureg</a> qureg, <span class="keywordtype">int</span> targetQubit, <a class="code" href="group__type.html#ga7740e349b4f8bae6451547680f0ce2d6">qreal</a> angle) {</div>
<div class="line"><a name="l00255"></a><span class="lineno"> 255</span>  <a class="code" href="structComplex.html">Complex</a> term; </div>
<div class="line"><a name="l00256"></a><span class="lineno"> 256</span>  term.<a class="code" href="structComplex.html#ab5b2e2eca02131fc74b289a83636cbe3">real</a> = cos(angle); </div>
<div class="line"><a name="l00257"></a><span class="lineno"> 257</span>  term.<a class="code" href="structComplex.html#a84f5439aad0ef495efdd3b4c1c02d27e">imag</a> = sin(angle);</div>
<div class="line"><a name="l00258"></a><span class="lineno"> 258</span>  <a class="code" href="QuEST__internal_8h.html#a535cc3c9962d511078691e4d64b4e272">statevec_phaseShiftByTerm</a>(qureg, targetQubit, term);</div>
<div class="line"><a name="l00259"></a><span class="lineno"> 259</span> }</div>
<div class="line"><a name="l00260"></a><span class="lineno"> 260</span>  </div>
<div class="line"><a name="l00261"></a><span class="lineno"><a class="line" href="QuEST__internal_8h.html#ab14df201e5c05527dc77cfeb6f583042"> 261</a></span> <span class="keywordtype">void</span> <a class="code" href="QuEST__common_8c.html#ab14df201e5c05527dc77cfeb6f583042">statevec_pauliZ</a>(<a class="code" href="structQureg.html">Qureg</a> qureg, <span class="keywordtype">int</span> targetQubit) {</div>
<div class="line"><a name="l00262"></a><span class="lineno"> 262</span>  <a class="code" href="structComplex.html">Complex</a> term; </div>
<div class="line"><a name="l00263"></a><span class="lineno"> 263</span>  term.<a class="code" href="structComplex.html#ab5b2e2eca02131fc74b289a83636cbe3">real</a> = -1;</div>
<div class="line"><a name="l00264"></a><span class="lineno"> 264</span>  term.<a class="code" href="structComplex.html#a84f5439aad0ef495efdd3b4c1c02d27e">imag</a> = 0;</div>
<div class="line"><a name="l00265"></a><span class="lineno"> 265</span>  <a class="code" href="QuEST__internal_8h.html#a535cc3c9962d511078691e4d64b4e272">statevec_phaseShiftByTerm</a>(qureg, targetQubit, term);</div>
<div class="line"><a name="l00266"></a><span class="lineno"> 266</span> }</div>
<div class="line"><a name="l00267"></a><span class="lineno"> 267</span>  </div>
<div class="line"><a name="l00268"></a><span class="lineno"><a class="line" href="QuEST__internal_8h.html#a559d175d31cec552fc09aec3e8a78d11"> 268</a></span> <span class="keywordtype">void</span> <a class="code" href="QuEST__common_8c.html#a559d175d31cec552fc09aec3e8a78d11">statevec_sGate</a>(<a class="code" href="structQureg.html">Qureg</a> qureg, <span class="keywordtype">int</span> targetQubit) {</div>
<div class="line"><a name="l00269"></a><span class="lineno"> 269</span>  <a class="code" href="structComplex.html">Complex</a> term; </div>
<div class="line"><a name="l00270"></a><span class="lineno"> 270</span>  term.<a class="code" href="structComplex.html#ab5b2e2eca02131fc74b289a83636cbe3">real</a> = 0;</div>
<div class="line"><a name="l00271"></a><span class="lineno"> 271</span>  term.<a class="code" href="structComplex.html#a84f5439aad0ef495efdd3b4c1c02d27e">imag</a> = 1;</div>
<div class="line"><a name="l00272"></a><span class="lineno"> 272</span>  <a class="code" href="QuEST__internal_8h.html#a535cc3c9962d511078691e4d64b4e272">statevec_phaseShiftByTerm</a>(qureg, targetQubit, term);</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"><a class="line" href="QuEST__internal_8h.html#a2d96b4601cc3a9e1221e11a366fe1699"> 275</a></span> <span class="keywordtype">void</span> <a class="code" href="QuEST__common_8c.html#a2d96b4601cc3a9e1221e11a366fe1699">statevec_tGate</a>(<a class="code" href="structQureg.html">Qureg</a> qureg, <span class="keywordtype">int</span> targetQubit) {</div>
<div class="line"><a name="l00276"></a><span class="lineno"> 276</span>  <a class="code" href="structComplex.html">Complex</a> term; </div>
<div class="line"><a name="l00277"></a><span class="lineno"> 277</span>  term.<a class="code" href="structComplex.html#ab5b2e2eca02131fc74b289a83636cbe3">real</a> = 1/sqrt(2);</div>
<div class="line"><a name="l00278"></a><span class="lineno"> 278</span>  term.<a class="code" href="structComplex.html#a84f5439aad0ef495efdd3b4c1c02d27e">imag</a> = 1/sqrt(2);</div>
<div class="line"><a name="l00279"></a><span class="lineno"> 279</span>  <a class="code" href="QuEST__internal_8h.html#a535cc3c9962d511078691e4d64b4e272">statevec_phaseShiftByTerm</a>(qureg, targetQubit, term);</div>
<div class="line"><a name="l00280"></a><span class="lineno"> 280</span> }</div>
<div class="line"><a name="l00281"></a><span class="lineno"> 281</span>  </div>
<div class="line"><a name="l00282"></a><span class="lineno"><a class="line" href="QuEST__internal_8h.html#a142b424d903dfaca7a7ac166076f543a"> 282</a></span> <span class="keywordtype">void</span> <a class="code" href="QuEST__common_8c.html#a142b424d903dfaca7a7ac166076f543a">statevec_sGateConj</a>(<a class="code" href="structQureg.html">Qureg</a> qureg, <span class="keywordtype">int</span> targetQubit) {</div>
<div class="line"><a name="l00283"></a><span class="lineno"> 283</span>  <a class="code" href="structComplex.html">Complex</a> term; </div>
<div class="line"><a name="l00284"></a><span class="lineno"> 284</span>  term.<a class="code" href="structComplex.html#ab5b2e2eca02131fc74b289a83636cbe3">real</a> = 0;</div>
<div class="line"><a name="l00285"></a><span class="lineno"> 285</span>  term.<a class="code" href="structComplex.html#a84f5439aad0ef495efdd3b4c1c02d27e">imag</a> = -1;</div>
<div class="line"><a name="l00286"></a><span class="lineno"> 286</span>  <a class="code" href="QuEST__internal_8h.html#a535cc3c9962d511078691e4d64b4e272">statevec_phaseShiftByTerm</a>(qureg, targetQubit, term);</div>
<div class="line"><a name="l00287"></a><span class="lineno"> 287</span> } </div>
<div class="line"><a name="l00288"></a><span class="lineno"> 288</span>  </div>
<div class="line"><a name="l00289"></a><span class="lineno"><a class="line" href="QuEST__internal_8h.html#a086fffa88b1aeec1405aa2dc41eab8a3"> 289</a></span> <span class="keywordtype">void</span> <a class="code" href="QuEST__common_8c.html#a086fffa88b1aeec1405aa2dc41eab8a3">statevec_tGateConj</a>(<a class="code" href="structQureg.html">Qureg</a> qureg, <span class="keywordtype">int</span> targetQubit) {</div>
<div class="line"><a name="l00290"></a><span class="lineno"> 290</span>  <a class="code" href="structComplex.html">Complex</a> term; </div>
<div class="line"><a name="l00291"></a><span class="lineno"> 291</span>  term.<a class="code" href="structComplex.html#ab5b2e2eca02131fc74b289a83636cbe3">real</a> = 1/sqrt(2);</div>
<div class="line"><a name="l00292"></a><span class="lineno"> 292</span>  term.<a class="code" href="structComplex.html#a84f5439aad0ef495efdd3b4c1c02d27e">imag</a> = -1/sqrt(2);</div>
<div class="line"><a name="l00293"></a><span class="lineno"> 293</span>  <a class="code" href="QuEST__internal_8h.html#a535cc3c9962d511078691e4d64b4e272">statevec_phaseShiftByTerm</a>(qureg, targetQubit, term);</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>  </div>
<div class="line"><a name="l00296"></a><span class="lineno"><a class="line" href="QuEST__internal_8h.html#ab8687e6ea5cf35503c223051837f3b94"> 296</a></span> <span class="keywordtype">void</span> <a class="code" href="QuEST__common_8c.html#ab8687e6ea5cf35503c223051837f3b94">statevec_rotateX</a>(<a class="code" href="structQureg.html">Qureg</a> qureg, <span class="keywordtype">int</span> rotQubit, <a class="code" href="group__type.html#ga7740e349b4f8bae6451547680f0ce2d6">qreal</a> angle){</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>  <a class="code" href="structVector.html">Vector</a> unitAxis = {1, 0, 0};</div>
<div class="line"><a name="l00299"></a><span class="lineno"> 299</span>  <a class="code" href="QuEST__common_8c.html#adf5fbea75826cde7d86fd7c31bcf3dec">statevec_rotateAroundAxis</a>(qureg, rotQubit, angle, unitAxis);</div>
<div class="line"><a name="l00300"></a><span class="lineno"> 300</span> }</div>
<div class="line"><a name="l00301"></a><span class="lineno"> 301</span>  </div>
<div class="line"><a name="l00302"></a><span class="lineno"><a class="line" href="QuEST__internal_8h.html#a1f0db2800d13b6c2b120e05129cd13cb"> 302</a></span> <span class="keywordtype">void</span> <a class="code" href="QuEST__common_8c.html#a1f0db2800d13b6c2b120e05129cd13cb">statevec_rotateY</a>(<a class="code" href="structQureg.html">Qureg</a> qureg, <span class="keywordtype">int</span> rotQubit, <a class="code" href="group__type.html#ga7740e349b4f8bae6451547680f0ce2d6">qreal</a> angle){</div>
<div class="line"><a name="l00303"></a><span class="lineno"> 303</span>  </div>
<div class="line"><a name="l00304"></a><span class="lineno"> 304</span>  <a class="code" href="structVector.html">Vector</a> unitAxis = {0, 1, 0};</div>
<div class="line"><a name="l00305"></a><span class="lineno"> 305</span>  <a class="code" href="QuEST__common_8c.html#adf5fbea75826cde7d86fd7c31bcf3dec">statevec_rotateAroundAxis</a>(qureg, rotQubit, angle, unitAxis);</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>  </div>
<div class="line"><a name="l00308"></a><span class="lineno"><a class="line" href="QuEST__internal_8h.html#a2b2b8c69566a49cdf4a3fa2cb4f5e03c"> 308</a></span> <span class="keywordtype">void</span> <a class="code" href="QuEST__common_8c.html#a2b2b8c69566a49cdf4a3fa2cb4f5e03c">statevec_rotateZ</a>(<a class="code" href="structQureg.html">Qureg</a> qureg, <span class="keywordtype">int</span> rotQubit, <a class="code" href="group__type.html#ga7740e349b4f8bae6451547680f0ce2d6">qreal</a> angle){</div>
<div class="line"><a name="l00309"></a><span class="lineno"> 309</span>  </div>
<div class="line"><a name="l00310"></a><span class="lineno"> 310</span>  <a class="code" href="structVector.html">Vector</a> unitAxis = {0, 0, 1};</div>
<div class="line"><a name="l00311"></a><span class="lineno"> 311</span>  <a class="code" href="QuEST__common_8c.html#adf5fbea75826cde7d86fd7c31bcf3dec">statevec_rotateAroundAxis</a>(qureg, rotQubit, angle, unitAxis);</div>
<div class="line"><a name="l00312"></a><span class="lineno"> 312</span> }</div>
<div class="line"><a name="l00313"></a><span class="lineno"> 313</span>  </div>
<div class="line"><a name="l00314"></a><span class="lineno"><a class="line" href="QuEST__internal_8h.html#adf5fbea75826cde7d86fd7c31bcf3dec"> 314</a></span> <span class="keywordtype">void</span> <a class="code" href="QuEST__common_8c.html#adf5fbea75826cde7d86fd7c31bcf3dec">statevec_rotateAroundAxis</a>(<a class="code" href="structQureg.html">Qureg</a> qureg, <span class="keywordtype">int</span> rotQubit, <a class="code" href="group__type.html#ga7740e349b4f8bae6451547680f0ce2d6">qreal</a> angle, <a class="code" href="structVector.html">Vector</a> axis){</div>
<div class="line"><a name="l00315"></a><span class="lineno"> 315</span>  </div>
<div class="line"><a name="l00316"></a><span class="lineno"> 316</span>  <a class="code" href="structComplex.html">Complex</a> alpha, beta;</div>
<div class="line"><a name="l00317"></a><span class="lineno"> 317</span>  <a class="code" href="QuEST__common_8c.html#a71b3ae4402c1c5039d3086ca0e47a522">getComplexPairFromRotation</a>(angle, axis, &alpha, &beta);</div>
<div class="line"><a name="l00318"></a><span class="lineno"> 318</span>  <a class="code" href="QuEST__internal_8h.html#aec63559e8b9ab0e17efe18d64d2ceca6">statevec_compactUnitary</a>(qureg, rotQubit, alpha, beta);</div>
<div class="line"><a name="l00319"></a><span class="lineno"> 319</span> }</div>
<div class="line"><a name="l00320"></a><span class="lineno"> 320</span>  </div>
<div class="line"><a name="l00321"></a><span class="lineno"><a class="line" href="QuEST__internal_8h.html#a2cbfe9f416d4d8acb56213841f62cb53"> 321</a></span> <span class="keywordtype">void</span> <a class="code" href="QuEST__common_8c.html#a2cbfe9f416d4d8acb56213841f62cb53">statevec_rotateAroundAxisConj</a>(<a class="code" href="structQureg.html">Qureg</a> qureg, <span class="keywordtype">int</span> rotQubit, <a class="code" href="group__type.html#ga7740e349b4f8bae6451547680f0ce2d6">qreal</a> angle, <a class="code" href="structVector.html">Vector</a> axis){</div>
<div class="line"><a name="l00322"></a><span class="lineno"> 322</span>  </div>
<div class="line"><a name="l00323"></a><span class="lineno"> 323</span>  <a class="code" href="structComplex.html">Complex</a> alpha, beta;</div>
<div class="line"><a name="l00324"></a><span class="lineno"> 324</span>  <a class="code" href="QuEST__common_8c.html#a71b3ae4402c1c5039d3086ca0e47a522">getComplexPairFromRotation</a>(angle, axis, &alpha, &beta);</div>
<div class="line"><a name="l00325"></a><span class="lineno"> 325</span>  alpha.<a class="code" href="structComplex.html#a84f5439aad0ef495efdd3b4c1c02d27e">imag</a> *= -1; </div>
<div class="line"><a name="l00326"></a><span class="lineno"> 326</span>  beta.<a class="code" href="structComplex.html#a84f5439aad0ef495efdd3b4c1c02d27e">imag</a> *= -1;</div>
<div class="line"><a name="l00327"></a><span class="lineno"> 327</span>  <a class="code" href="QuEST__internal_8h.html#aec63559e8b9ab0e17efe18d64d2ceca6">statevec_compactUnitary</a>(qureg, rotQubit, alpha, beta);</div>
<div class="line"><a name="l00328"></a><span class="lineno"> 328</span> }</div>
<div class="line"><a name="l00329"></a><span class="lineno"> 329</span>  </div>
<div class="line"><a name="l00330"></a><span class="lineno"><a class="line" href="QuEST__internal_8h.html#a336c3286d2434fc892ebe20e9fa22a24"> 330</a></span> <span class="keywordtype">void</span> <a class="code" href="QuEST__common_8c.html#a336c3286d2434fc892ebe20e9fa22a24">statevec_controlledRotateAroundAxis</a>(<a class="code" href="structQureg.html">Qureg</a> qureg, <span class="keywordtype">int</span> controlQubit, <span class="keywordtype">int</span> targetQubit, <a class="code" href="group__type.html#ga7740e349b4f8bae6451547680f0ce2d6">qreal</a> angle, <a class="code" href="structVector.html">Vector</a> axis){</div>
<div class="line"><a name="l00331"></a><span class="lineno"> 331</span>  </div>
<div class="line"><a name="l00332"></a><span class="lineno"> 332</span>  <a class="code" href="structComplex.html">Complex</a> alpha, beta;</div>
<div class="line"><a name="l00333"></a><span class="lineno"> 333</span>  <a class="code" href="QuEST__common_8c.html#a71b3ae4402c1c5039d3086ca0e47a522">getComplexPairFromRotation</a>(angle, axis, &alpha, &beta);</div>
<div class="line"><a name="l00334"></a><span class="lineno"> 334</span>  <a class="code" href="QuEST__internal_8h.html#a1275d604674224d87a173fb5bac78835">statevec_controlledCompactUnitary</a>(qureg, controlQubit, targetQubit, alpha, beta);</div>
<div class="line"><a name="l00335"></a><span class="lineno"> 335</span> }</div>
<div class="line"><a name="l00336"></a><span class="lineno"> 336</span>  </div>
<div class="line"><a name="l00337"></a><span class="lineno"><a class="line" href="QuEST__internal_8h.html#a5ea2384a9b33c94eac891aacbfa5da6a"> 337</a></span> <span class="keywordtype">void</span> <a class="code" href="QuEST__common_8c.html#a5ea2384a9b33c94eac891aacbfa5da6a">statevec_controlledRotateAroundAxisConj</a>(<a class="code" href="structQureg.html">Qureg</a> qureg, <span class="keywordtype">int</span> controlQubit, <span class="keywordtype">int</span> targetQubit, <a class="code" href="group__type.html#ga7740e349b4f8bae6451547680f0ce2d6">qreal</a> angle, <a class="code" href="structVector.html">Vector</a> axis){</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>  <a class="code" href="structComplex.html">Complex</a> alpha, beta;</div>
<div class="line"><a name="l00340"></a><span class="lineno"> 340</span>  <a class="code" href="QuEST__common_8c.html#a71b3ae4402c1c5039d3086ca0e47a522">getComplexPairFromRotation</a>(angle, axis, &alpha, &beta);</div>
<div class="line"><a name="l00341"></a><span class="lineno"> 341</span>  alpha.<a class="code" href="structComplex.html#a84f5439aad0ef495efdd3b4c1c02d27e">imag</a> *= -1; </div>
<div class="line"><a name="l00342"></a><span class="lineno"> 342</span>  beta.<a class="code" href="structComplex.html#a84f5439aad0ef495efdd3b4c1c02d27e">imag</a> *= -1;</div>
<div class="line"><a name="l00343"></a><span class="lineno"> 343</span>  <a class="code" href="QuEST__internal_8h.html#a1275d604674224d87a173fb5bac78835">statevec_controlledCompactUnitary</a>(qureg, controlQubit, targetQubit, alpha, beta);</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>  </div>
<div class="line"><a name="l00346"></a><span class="lineno"><a class="line" href="QuEST__internal_8h.html#aa32eed5bc38e85193e1aa23b74e8541b"> 346</a></span> <span class="keywordtype">void</span> <a class="code" href="QuEST__common_8c.html#aa32eed5bc38e85193e1aa23b74e8541b">statevec_controlledRotateX</a>(<a class="code" href="structQureg.html">Qureg</a> qureg, <span class="keywordtype">int</span> controlQubit, <span class="keywordtype">int</span> targetQubit, <a class="code" href="group__type.html#ga7740e349b4f8bae6451547680f0ce2d6">qreal</a> angle){</div>
<div class="line"><a name="l00347"></a><span class="lineno"> 347</span>  </div>
<div class="line"><a name="l00348"></a><span class="lineno"> 348</span>  <a class="code" href="structVector.html">Vector</a> unitAxis = {1, 0, 0};</div>
<div class="line"><a name="l00349"></a><span class="lineno"> 349</span>  <a class="code" href="QuEST__common_8c.html#a336c3286d2434fc892ebe20e9fa22a24">statevec_controlledRotateAroundAxis</a>(qureg, controlQubit, targetQubit, angle, unitAxis);</div>
<div class="line"><a name="l00350"></a><span class="lineno"> 350</span> }</div>
<div class="line"><a name="l00351"></a><span class="lineno"> 351</span>  </div>
<div class="line"><a name="l00352"></a><span class="lineno"><a class="line" href="QuEST__internal_8h.html#ade67874c89179541f616b3685595d3ae"> 352</a></span> <span class="keywordtype">void</span> <a class="code" href="QuEST__common_8c.html#ade67874c89179541f616b3685595d3ae">statevec_controlledRotateY</a>(<a class="code" href="structQureg.html">Qureg</a> qureg, <span class="keywordtype">int</span> controlQubit, <span class="keywordtype">int</span> targetQubit, <a class="code" href="group__type.html#ga7740e349b4f8bae6451547680f0ce2d6">qreal</a> angle){</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>  <a class="code" href="structVector.html">Vector</a> unitAxis = {0, 1, 0};</div>
<div class="line"><a name="l00355"></a><span class="lineno"> 355</span>  <a class="code" href="QuEST__common_8c.html#a336c3286d2434fc892ebe20e9fa22a24">statevec_controlledRotateAroundAxis</a>(qureg, controlQubit, targetQubit, angle, unitAxis);</div>
<div class="line"><a name="l00356"></a><span class="lineno"> 356</span> }</div>
<div class="line"><a name="l00357"></a><span class="lineno"> 357</span>  </div>
<div class="line"><a name="l00358"></a><span class="lineno"><a class="line" href="QuEST__internal_8h.html#a3e11a9b9fa5ee16d4738a97f4590ee61"> 358</a></span> <span class="keywordtype">void</span> <a class="code" href="QuEST__common_8c.html#a3e11a9b9fa5ee16d4738a97f4590ee61">statevec_controlledRotateZ</a>(<a class="code" href="structQureg.html">Qureg</a> qureg, <span class="keywordtype">int</span> controlQubit, <span class="keywordtype">int</span> targetQubit, <a class="code" href="group__type.html#ga7740e349b4f8bae6451547680f0ce2d6">qreal</a> angle){</div>
<div class="line"><a name="l00359"></a><span class="lineno"> 359</span>  </div>
<div class="line"><a name="l00360"></a><span class="lineno"> 360</span>  <a class="code" href="structVector.html">Vector</a> unitAxis = {0, 0, 1};</div>
<div class="line"><a name="l00361"></a><span class="lineno"> 361</span>  <a class="code" href="QuEST__common_8c.html#a336c3286d2434fc892ebe20e9fa22a24">statevec_controlledRotateAroundAxis</a>(qureg, controlQubit, targetQubit, angle, unitAxis);</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>  </div>
<div class="line"><a name="l00364"></a><span class="lineno"><a class="line" href="QuEST__internal_8h.html#ac2eb81c71d0e0fbec7be5b45b1db1448"> 364</a></span> <span class="keywordtype">int</span> <a class="code" href="QuEST__common_8c.html#ac2eb81c71d0e0fbec7be5b45b1db1448">statevec_measureWithStats</a>(<a class="code" href="structQureg.html">Qureg</a> qureg, <span class="keywordtype">int</span> measureQubit, <a class="code" href="group__type.html#ga7740e349b4f8bae6451547680f0ce2d6">qreal</a> *outcomeProb) {</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>  <a class="code" href="group__type.html#ga7740e349b4f8bae6451547680f0ce2d6">qreal</a> zeroProb = <a class="code" href="QuEST__internal_8h.html#a63cb88a5f4eb12fef087435b8db3f8eb">statevec_calcProbOfOutcome</a>(qureg, measureQubit, 0);</div>
<div class="line"><a name="l00367"></a><span class="lineno"> 367</span>  <span class="keywordtype">int</span> outcome = <a class="code" href="QuEST__common_8c.html#a9d19fc690afb425d0df5feb09be567a3">generateMeasurementOutcome</a>(zeroProb, outcomeProb);</div>
<div class="line"><a name="l00368"></a><span class="lineno"> 368</span>  <a class="code" href="QuEST__internal_8h.html#a008becc4a18f868340836e0ab9fd6df6">statevec_collapseToKnownProbOutcome</a>(qureg, measureQubit, outcome, *outcomeProb);</div>
<div class="line"><a name="l00369"></a><span class="lineno"> 369</span>  <span class="keywordflow">return</span> outcome;</div>
<div class="line"><a name="l00370"></a><span class="lineno"> 370</span> }</div>
<div class="line"><a name="l00371"></a><span class="lineno"> 371</span>  </div>
<div class="line"><a name="l00372"></a><span class="lineno"><a class="line" href="QuEST__internal_8h.html#ae4ebc1f4dd48922bfb30ad1692310269"> 372</a></span> <span class="keywordtype">int</span> <a class="code" href="QuEST__common_8c.html#ae4ebc1f4dd48922bfb30ad1692310269">densmatr_measureWithStats</a>(<a class="code" href="structQureg.html">Qureg</a> qureg, <span class="keywordtype">int</span> measureQubit, <a class="code" href="group__type.html#ga7740e349b4f8bae6451547680f0ce2d6">qreal</a> *outcomeProb) {</div>
<div class="line"><a name="l00373"></a><span class="lineno"> 373</span>  </div>
<div class="line"><a name="l00374"></a><span class="lineno"> 374</span>  <a class="code" href="group__type.html#ga7740e349b4f8bae6451547680f0ce2d6">qreal</a> zeroProb = <a class="code" href="QuEST__internal_8h.html#ae37b3eff5d3affbc1edfe216cabc1c53">densmatr_calcProbOfOutcome</a>(qureg, measureQubit, 0);</div>
<div class="line"><a name="l00375"></a><span class="lineno"> 375</span>  <span class="keywordtype">int</span> outcome = <a class="code" href="QuEST__common_8c.html#a9d19fc690afb425d0df5feb09be567a3">generateMeasurementOutcome</a>(zeroProb, outcomeProb);</div>
<div class="line"><a name="l00376"></a><span class="lineno"> 376</span>  <a class="code" href="QuEST__internal_8h.html#a5658f2ecbbd4ea425db5f7dc7fc4ba92">densmatr_collapseToKnownProbOutcome</a>(qureg, measureQubit, outcome, *outcomeProb);</div>
<div class="line"><a name="l00377"></a><span class="lineno"> 377</span>  <span class="keywordflow">return</span> outcome;</div>
<div class="line"><a name="l00378"></a><span class="lineno"> 378</span> }</div>
<div class="line"><a name="l00379"></a><span class="lineno"> 379</span>  </div>
<div class="line"><a name="l00380"></a><span class="lineno"><a class="line" href="QuEST__internal_8h.html#aeeec9fe1501e25930adf6c38763a158e"> 380</a></span> <a class="code" href="group__type.html#ga7740e349b4f8bae6451547680f0ce2d6">qreal</a> <a class="code" href="QuEST__common_8c.html#aeeec9fe1501e25930adf6c38763a158e">statevec_calcFidelity</a>(<a class="code" href="structQureg.html">Qureg</a> qureg, <a class="code" href="structQureg.html">Qureg</a> pureState) {</div>
<div class="line"><a name="l00381"></a><span class="lineno"> 381</span>  </div>
<div class="line"><a name="l00382"></a><span class="lineno"> 382</span>  <a class="code" href="structComplex.html">Complex</a> innerProd = <a class="code" href="QuEST__internal_8h.html#a7ebd3198a198f4cd20840f64fd8b84d0">statevec_calcInnerProduct</a>(qureg, pureState);</div>
<div class="line"><a name="l00383"></a><span class="lineno"> 383</span>  <a class="code" href="group__type.html#ga7740e349b4f8bae6451547680f0ce2d6">qreal</a> innerProdMag = innerProd.<a class="code" href="structComplex.html#ab5b2e2eca02131fc74b289a83636cbe3">real</a>*innerProd.<a class="code" href="structComplex.html#ab5b2e2eca02131fc74b289a83636cbe3">real</a> + innerProd.<a class="code" href="structComplex.html#a84f5439aad0ef495efdd3b4c1c02d27e">imag</a>*innerProd.<a class="code" href="structComplex.html#a84f5439aad0ef495efdd3b4c1c02d27e">imag</a>;</div>
<div class="line"><a name="l00384"></a><span class="lineno"> 384</span>  <span class="keywordflow">return</span> innerProdMag;</div>
<div class="line"><a name="l00385"></a><span class="lineno"> 385</span> }</div>
<div class="line"><a name="l00386"></a><span class="lineno"> 386</span>  </div>
<div class="line"><a name="l00387"></a><span class="lineno"><a class="line" href="QuEST__internal_8h.html#a6550462400942be61b171c8840b36aa1"> 387</a></span> <span class="keywordtype">void</span> <a class="code" href="QuEST__common_8c.html#a6550462400942be61b171c8840b36aa1">statevec_sqrtSwapGate</a>(<a class="code" href="structQureg.html">Qureg</a> qureg, <span class="keywordtype">int</span> qb1, <span class="keywordtype">int</span> qb2) {</div>
<div class="line"><a name="l00388"></a><span class="lineno"> 388</span>  </div>
<div class="line"><a name="l00389"></a><span class="lineno"> 389</span>  <a class="code" href="structComplexMatrix4.html">ComplexMatrix4</a> u = (<a class="code" href="structComplexMatrix4.html">ComplexMatrix4</a>) {.<a class="code" href="structComplexMatrix4.html#aeadc9893e1b16cdd585dac523e73926c">real</a>={{0}}, .imag={{0}}};</div>
<div class="line"><a name="l00390"></a><span class="lineno"> 390</span>  u.<a class="code" href="structComplexMatrix4.html#aeadc9893e1b16cdd585dac523e73926c">real</a>[0][0]=1;</div>
<div class="line"><a name="l00391"></a><span class="lineno"> 391</span>  u.<a class="code" href="structComplexMatrix4.html#aeadc9893e1b16cdd585dac523e73926c">real</a>[3][3]=1;</div>
<div class="line"><a name="l00392"></a><span class="lineno"> 392</span>  u.<a class="code" href="structComplexMatrix4.html#aeadc9893e1b16cdd585dac523e73926c">real</a>[1][1] = .5; u.<a class="code" href="structComplexMatrix4.html#a1f6bb840b5a9e6252980cab7a14f1a35">imag</a>[1][1] = .5;</div>
<div class="line"><a name="l00393"></a><span class="lineno"> 393</span>  u.<a class="code" href="structComplexMatrix4.html#aeadc9893e1b16cdd585dac523e73926c">real</a>[1][2] = .5; u.<a class="code" href="structComplexMatrix4.html#a1f6bb840b5a9e6252980cab7a14f1a35">imag</a>[1][2] =-.5;</div>
<div class="line"><a name="l00394"></a><span class="lineno"> 394</span>  u.<a class="code" href="structComplexMatrix4.html#aeadc9893e1b16cdd585dac523e73926c">real</a>[2][1] = .5; u.<a class="code" href="structComplexMatrix4.html#a1f6bb840b5a9e6252980cab7a14f1a35">imag</a>[2][1] =-.5;</div>
<div class="line"><a name="l00395"></a><span class="lineno"> 395</span>  u.<a class="code" href="structComplexMatrix4.html#aeadc9893e1b16cdd585dac523e73926c">real</a>[2][2] = .5; u.<a class="code" href="structComplexMatrix4.html#a1f6bb840b5a9e6252980cab7a14f1a35">imag</a>[2][2] = .5;</div>
<div class="line"><a name="l00396"></a><span class="lineno"> 396</span>  </div>
<div class="line"><a name="l00397"></a><span class="lineno"> 397</span>  <a class="code" href="QuEST__common_8c.html#a95abbe0cedf45b1ff2f070dba219fd7f">statevec_twoQubitUnitary</a>(qureg, qb1, qb2, u);</div>
<div class="line"><a name="l00398"></a><span class="lineno"> 398</span> }</div>
<div class="line"><a name="l00399"></a><span class="lineno"> 399</span>  </div>
<div class="line"><a name="l00400"></a><span class="lineno"><a class="line" href="QuEST__internal_8h.html#a33cde9b6518d55a544eff20301a5d844"> 400</a></span> <span class="keywordtype">void</span> <a class="code" href="QuEST__common_8c.html#a33cde9b6518d55a544eff20301a5d844">statevec_sqrtSwapGateConj</a>(<a class="code" href="structQureg.html">Qureg</a> qureg, <span class="keywordtype">int</span> qb1, <span class="keywordtype">int</span> qb2) {</div>
<div class="line"><a name="l00401"></a><span class="lineno"> 401</span>  </div>
<div class="line"><a name="l00402"></a><span class="lineno"> 402</span>  <a class="code" href="structComplexMatrix4.html">ComplexMatrix4</a> u = (<a class="code" href="structComplexMatrix4.html">ComplexMatrix4</a>) {.<a class="code" href="structComplexMatrix4.html#aeadc9893e1b16cdd585dac523e73926c">real</a>={{0}}, .imag={{0}}};</div>
<div class="line"><a name="l00403"></a><span class="lineno"> 403</span>  u.<a class="code" href="structComplexMatrix4.html#aeadc9893e1b16cdd585dac523e73926c">real</a>[0][0]=1;</div>
<div class="line"><a name="l00404"></a><span class="lineno"> 404</span>  u.<a class="code" href="structComplexMatrix4.html#aeadc9893e1b16cdd585dac523e73926c">real</a>[3][3]=1;</div>
<div class="line"><a name="l00405"></a><span class="lineno"> 405</span>  u.<a class="code" href="structComplexMatrix4.html#aeadc9893e1b16cdd585dac523e73926c">real</a>[1][1] = .5; u.<a class="code" href="structComplexMatrix4.html#a1f6bb840b5a9e6252980cab7a14f1a35">imag</a>[1][1] =-.5;</div>
<div class="line"><a name="l00406"></a><span class="lineno"> 406</span>  u.<a class="code" href="structComplexMatrix4.html#aeadc9893e1b16cdd585dac523e73926c">real</a>[1][2] = .5; u.<a class="code" href="structComplexMatrix4.html#a1f6bb840b5a9e6252980cab7a14f1a35">imag</a>[1][2] = .5;</div>
<div class="line"><a name="l00407"></a><span class="lineno"> 407</span>  u.<a class="code" href="structComplexMatrix4.html#aeadc9893e1b16cdd585dac523e73926c">real</a>[2][1] = .5; u.<a class="code" href="structComplexMatrix4.html#a1f6bb840b5a9e6252980cab7a14f1a35">imag</a>[2][1] = .5;</div>
<div class="line"><a name="l00408"></a><span class="lineno"> 408</span>  u.<a class="code" href="structComplexMatrix4.html#aeadc9893e1b16cdd585dac523e73926c">real</a>[2][2] = .5; u.<a class="code" href="structComplexMatrix4.html#a1f6bb840b5a9e6252980cab7a14f1a35">imag</a>[2][2] =-.5;</div>
<div class="line"><a name="l00409"></a><span class="lineno"> 409</span>  </div>
<div class="line"><a name="l00410"></a><span class="lineno"> 410</span>  <a class="code" href="QuEST__common_8c.html#a95abbe0cedf45b1ff2f070dba219fd7f">statevec_twoQubitUnitary</a>(qureg, qb1, qb2, u);</div>
<div class="line"><a name="l00411"></a><span class="lineno"> 411</span> }</div>
<div class="line"><a name="l00412"></a><span class="lineno"> 412</span>  </div>
<div class="line"><a name="l00414"></a><span class="lineno"><a class="line" href="QuEST__internal_8h.html#a2153d8611247dd3a2ef483355265cc99"> 414</a></span> <span class="keywordtype">void</span> <a class="code" href="QuEST__common_8c.html#a2153d8611247dd3a2ef483355265cc99">statevec_multiRotatePauli</a>(</div>
<div class="line"><a name="l00415"></a><span class="lineno"> 415</span>  <a class="code" href="structQureg.html">Qureg</a> qureg, <span class="keywordtype">int</span>* targetQubits, <span class="keyword">enum</span> <a class="code" href="group__type.html#ga1c703cf89629e4e9c7023cd402d67028">pauliOpType</a>* targetPaulis, <span class="keywordtype">int</span> numTargets, <a class="code" href="group__type.html#ga7740e349b4f8bae6451547680f0ce2d6">qreal</a> angle,</div>
<div class="line"><a name="l00416"></a><span class="lineno"> 416</span>  <span class="keywordtype">int</span> applyConj</div>
<div class="line"><a name="l00417"></a><span class="lineno"> 417</span> ) {</div>
<div class="line"><a name="l00418"></a><span class="lineno"> 418</span>  <a class="code" href="group__type.html#ga7740e349b4f8bae6451547680f0ce2d6">qreal</a> fac = 1/sqrt(2);</div>
<div class="line"><a name="l00419"></a><span class="lineno"> 419</span>  <a class="code" href="structComplex.html">Complex</a> uRxAlpha = {.<a class="code" href="structComplex.html#ab5b2e2eca02131fc74b289a83636cbe3">real</a> = fac, .imag = 0}; <span class="comment">// Rx(pi/2)* rotates Z -> Y</span></div>
<div class="line"><a name="l00420"></a><span class="lineno"> 420</span>  <a class="code" href="structComplex.html">Complex</a> uRxBeta = {.<a class="code" href="structComplex.html#ab5b2e2eca02131fc74b289a83636cbe3">real</a> = 0, .imag = (applyConj)? fac : -fac};</div>
<div class="line"><a name="l00421"></a><span class="lineno"> 421</span>  <a class="code" href="structComplex.html">Complex</a> uRyAlpha = {.<a class="code" href="structComplex.html#ab5b2e2eca02131fc74b289a83636cbe3">real</a> = fac, .imag = 0}; <span class="comment">// Ry(-pi/2) rotates Z -> X</span></div>
<div class="line"><a name="l00422"></a><span class="lineno"> 422</span>  <a class="code" href="structComplex.html">Complex</a> uRyBeta = {.<a class="code" href="structComplex.html#ab5b2e2eca02131fc74b289a83636cbe3">real</a> = -fac, .imag = 0};</div>
<div class="line"><a name="l00423"></a><span class="lineno"> 423</span>  </div>
<div class="line"><a name="l00424"></a><span class="lineno"> 424</span>  <span class="comment">// mask may be modified to remove superfluous Identity ops</span></div>
<div class="line"><a name="l00425"></a><span class="lineno"> 425</span>  <span class="keywordtype">long</span> <span class="keywordtype">long</span> <span class="keywordtype">int</span> mask = <a class="code" href="QuEST__common_8c.html#acc45bb8300922ea5933e294921e25ea2">getQubitBitMask</a>(targetQubits, numTargets);</div>
<div class="line"><a name="l00426"></a><span class="lineno"> 426</span>  </div>
<div class="line"><a name="l00427"></a><span class="lineno"> 427</span>  <span class="comment">// rotate basis so that exp(Z) will effect exp(Y) and exp(X)</span></div>
<div class="line"><a name="l00428"></a><span class="lineno"> 428</span>  <span class="keywordflow">for</span> (<span class="keywordtype">int</span> t=0; t < numTargets; t++) {</div>
<div class="line"><a name="l00429"></a><span class="lineno"> 429</span>  <span class="keywordflow">if</span> (targetPaulis[t] == <a class="code" href="group__type.html#gga1c703cf89629e4e9c7023cd402d67028a63800f45b01ccc2608f754aee850bf4a">PAULI_I</a>)</div>
<div class="line"><a name="l00430"></a><span class="lineno"> 430</span>  mask -= 1LL << targetQubits[t]; <span class="comment">// remove target from mask</span></div>
<div class="line"><a name="l00431"></a><span class="lineno"> 431</span>  <span class="keywordflow">if</span> (targetPaulis[t] == <a class="code" href="group__type.html#gga1c703cf89629e4e9c7023cd402d67028a7abac7bb5b71e17382014c443244ad5c">PAULI_X</a>)</div>
<div class="line"><a name="l00432"></a><span class="lineno"> 432</span>  <a class="code" href="QuEST__internal_8h.html#aec63559e8b9ab0e17efe18d64d2ceca6">statevec_compactUnitary</a>(qureg, targetQubits[t], uRyAlpha, uRyBeta);</div>
<div class="line"><a name="l00433"></a><span class="lineno"> 433</span>  <span class="keywordflow">if</span> (targetPaulis[t] == <a class="code" href="group__type.html#gga1c703cf89629e4e9c7023cd402d67028a2c64a9aaee27f6bcfef894d159e759bd">PAULI_Y</a>)</div>
<div class="line"><a name="l00434"></a><span class="lineno"> 434</span>  <a class="code" href="QuEST__internal_8h.html#aec63559e8b9ab0e17efe18d64d2ceca6">statevec_compactUnitary</a>(qureg, targetQubits[t], uRxAlpha, uRxBeta);</div>
<div class="line"><a name="l00435"></a><span class="lineno"> 435</span>  <span class="comment">// (targetPaulis[t] == 3) is Z basis</span></div>
<div class="line"><a name="l00436"></a><span class="lineno"> 436</span>  }</div>
<div class="line"><a name="l00437"></a><span class="lineno"> 437</span>  </div>
<div class="line"><a name="l00438"></a><span class="lineno"> 438</span>  <span class="comment">// does nothing if there are no qubits to 'rotate'</span></div>
<div class="line"><a name="l00439"></a><span class="lineno"> 439</span>  <span class="keywordflow">if</span> (mask != 0)</div>
<div class="line"><a name="l00440"></a><span class="lineno"> 440</span>  <a class="code" href="QuEST__internal_8h.html#a8d88b4ed94ab024e5e7b2f1f4f505b44">statevec_multiRotateZ</a>(qureg, mask, (applyConj)? -angle : angle);</div>
<div class="line"><a name="l00441"></a><span class="lineno"> 441</span>  </div>
<div class="line"><a name="l00442"></a><span class="lineno"> 442</span>  <span class="comment">// undo X and Y basis rotations</span></div>
<div class="line"><a name="l00443"></a><span class="lineno"> 443</span>  uRxBeta.<a class="code" href="structComplex.html#a84f5439aad0ef495efdd3b4c1c02d27e">imag</a> *= -1;</div>
<div class="line"><a name="l00444"></a><span class="lineno"> 444</span>  uRyBeta.<a class="code" href="structComplex.html#ab5b2e2eca02131fc74b289a83636cbe3">real</a> *= -1;</div>
<div class="line"><a name="l00445"></a><span class="lineno"> 445</span>  <span class="keywordflow">for</span> (<span class="keywordtype">int</span> t=0; t < numTargets; t++) {</div>
<div class="line"><a name="l00446"></a><span class="lineno"> 446</span>  <span class="keywordflow">if</span> (targetPaulis[t] == <a class="code" href="group__type.html#gga1c703cf89629e4e9c7023cd402d67028a7abac7bb5b71e17382014c443244ad5c">PAULI_X</a>)</div>
<div class="line"><a name="l00447"></a><span class="lineno"> 447</span>  <a class="code" href="QuEST__internal_8h.html#aec63559e8b9ab0e17efe18d64d2ceca6">statevec_compactUnitary</a>(qureg, targetQubits[t], uRyAlpha, uRyBeta);</div>
<div class="line"><a name="l00448"></a><span class="lineno"> 448</span>  <span class="keywordflow">if</span> (targetPaulis[t] == <a class="code" href="group__type.html#gga1c703cf89629e4e9c7023cd402d67028a2c64a9aaee27f6bcfef894d159e759bd">PAULI_Y</a>)</div>
<div class="line"><a name="l00449"></a><span class="lineno"> 449</span>  <a class="code" href="QuEST__internal_8h.html#aec63559e8b9ab0e17efe18d64d2ceca6">statevec_compactUnitary</a>(qureg, targetQubits[t], uRxAlpha, uRxBeta);</div>
<div class="line"><a name="l00450"></a><span class="lineno"> 450</span>  }</div>
<div class="line"><a name="l00451"></a><span class="lineno"> 451</span> }</div>
<div class="line"><a name="l00452"></a><span class="lineno"> 452</span>  </div>
<div class="line"><a name="l00453"></a><span class="lineno"><a class="line" href="QuEST__internal_8h.html#a08d2e22870f3ef0289e355d634ea17a2"> 453</a></span> <span class="keywordtype">void</span> <a class="code" href="QuEST__common_8c.html#a08d2e22870f3ef0289e355d634ea17a2">statevec_multiControlledMultiRotatePauli</a>(</div>
<div class="line"><a name="l00454"></a><span class="lineno"> 454</span>  <a class="code" href="structQureg.html">Qureg</a> qureg, <span class="keywordtype">long</span> <span class="keywordtype">long</span> <span class="keywordtype">int</span> ctrlMask, <span class="keywordtype">int</span>* targetQubits, <span class="keyword">enum</span> <a class="code" href="group__type.html#ga1c703cf89629e4e9c7023cd402d67028">pauliOpType</a>* targetPaulis, <span class="keywordtype">int</span> numTargets, <a class="code" href="group__type.html#ga7740e349b4f8bae6451547680f0ce2d6">qreal</a> angle,</div>
<div class="line"><a name="l00455"></a><span class="lineno"> 455</span>  <span class="keywordtype">int</span> applyConj</div>
<div class="line"><a name="l00456"></a><span class="lineno"> 456</span> ) {</div>
<div class="line"><a name="l00457"></a><span class="lineno"> 457</span>  <a class="code" href="group__type.html#ga7740e349b4f8bae6451547680f0ce2d6">qreal</a> fac = 1/sqrt(2);</div>
<div class="line"><a name="l00458"></a><span class="lineno"> 458</span>  <a class="code" href="group__type.html#ga7740e349b4f8bae6451547680f0ce2d6">qreal</a> sgn = (applyConj)? 1 : -1;</div>
<div class="line"><a name="l00459"></a><span class="lineno"> 459</span>  <a class="code" href="structComplexMatrix2.html">ComplexMatrix2</a> uRx = {.<a class="code" href="structComplexMatrix2.html#a3ea98d545b67a190ee6abe936c72076f">real</a>={{fac,0},{0,fac}}, .imag={{0,sgn*fac},{sgn*fac,0}}}; <span class="comment">// Rx(pi/2)* rotates Z -> Y</span></div>
<div class="line"><a name="l00460"></a><span class="lineno"> 460</span>  <a class="code" href="structComplexMatrix2.html">ComplexMatrix2</a> uRy = {.<a class="code" href="structComplexMatrix2.html#a3ea98d545b67a190ee6abe936c72076f">real</a>={{fac,fac},{-fac,fac}}, .imag={{0,0},{0,0}}}; <span class="comment">// Ry(-pi/2) rotates Z -> X</span></div>
<div class="line"><a name="l00461"></a><span class="lineno"> 461</span>  </div>
<div class="line"><a name="l00462"></a><span class="lineno"> 462</span>  <span class="comment">// this function is controlled on the all-one state, so no ctrl flips</span></div>
<div class="line"><a name="l00463"></a><span class="lineno"> 463</span>  <span class="keywordtype">long</span> <span class="keywordtype">long</span> <span class="keywordtype">int</span> ctrlFlipMask = 0;</div>
<div class="line"><a name="l00464"></a><span class="lineno"> 464</span>  </div>
<div class="line"><a name="l00465"></a><span class="lineno"> 465</span>  <span class="comment">// mask may be modified to remove superfluous Identity ops</span></div>
<div class="line"><a name="l00466"></a><span class="lineno"> 466</span>  <span class="keywordtype">long</span> <span class="keywordtype">long</span> <span class="keywordtype">int</span> targMask = <a class="code" href="QuEST__common_8c.html#acc45bb8300922ea5933e294921e25ea2">getQubitBitMask</a>(targetQubits, numTargets);</div>
<div class="line"><a name="l00467"></a><span class="lineno"> 467</span>  </div>
<div class="line"><a name="l00468"></a><span class="lineno"> 468</span>  <span class="comment">// rotate basis so that exp(Z) will effect exp(Y) and exp(X)</span></div>
<div class="line"><a name="l00469"></a><span class="lineno"> 469</span>  <span class="keywordflow">for</span> (<span class="keywordtype">int</span> t=0; t < numTargets; t++) {</div>
<div class="line"><a name="l00470"></a><span class="lineno"> 470</span>  <span class="keywordflow">if</span> (targetPaulis[t] == <a class="code" href="group__type.html#gga1c703cf89629e4e9c7023cd402d67028a63800f45b01ccc2608f754aee850bf4a">PAULI_I</a>)</div>
<div class="line"><a name="l00471"></a><span class="lineno"> 471</span>  targMask -= 1LL << targetQubits[t]; <span class="comment">// remove target from mask</span></div>
<div class="line"><a name="l00472"></a><span class="lineno"> 472</span>  <span class="keywordflow">if</span> (targetPaulis[t] == <a class="code" href="group__type.html#gga1c703cf89629e4e9c7023cd402d67028a7abac7bb5b71e17382014c443244ad5c">PAULI_X</a>)</div>
<div class="line"><a name="l00473"></a><span class="lineno"> 473</span>  <a class="code" href="QuEST__internal_8h.html#a600bda10463e2078ce746ec28f2d39c9">statevec_multiControlledUnitary</a>(qureg, ctrlMask, ctrlFlipMask, targetQubits[t], uRy);</div>
<div class="line"><a name="l00474"></a><span class="lineno"> 474</span>  <span class="keywordflow">if</span> (targetPaulis[t] == <a class="code" href="group__type.html#gga1c703cf89629e4e9c7023cd402d67028a2c64a9aaee27f6bcfef894d159e759bd">PAULI_Y</a>)</div>
<div class="line"><a name="l00475"></a><span class="lineno"> 475</span>  <a class="code" href="QuEST__internal_8h.html#a600bda10463e2078ce746ec28f2d39c9">statevec_multiControlledUnitary</a>(qureg, ctrlMask, ctrlFlipMask, targetQubits[t], uRx);</div>
<div class="line"><a name="l00476"></a><span class="lineno"> 476</span>  <span class="comment">// (targetPaulis[t] == 3) is Z basis</span></div>
<div class="line"><a name="l00477"></a><span class="lineno"> 477</span>  }</div>
<div class="line"><a name="l00478"></a><span class="lineno"> 478</span>  </div>
<div class="line"><a name="l00479"></a><span class="lineno"> 479</span>  <span class="comment">// does nothing if there are no qubits to 'rotate'</span></div>
<div class="line"><a name="l00480"></a><span class="lineno"> 480</span>  <span class="keywordflow">if</span> (targMask != 0)</div>
<div class="line"><a name="l00481"></a><span class="lineno"> 481</span>  <a class="code" href="QuEST__internal_8h.html#a3ed30f823a59890be17441d97ce90c5c">statevec_multiControlledMultiRotateZ</a>(qureg, ctrlMask, targMask, (applyConj)? -angle : angle);</div>
<div class="line"><a name="l00482"></a><span class="lineno"> 482</span>  </div>
<div class="line"><a name="l00483"></a><span class="lineno"> 483</span>  <span class="comment">// undo X and Y basis rotations</span></div>
<div class="line"><a name="l00484"></a><span class="lineno"> 484</span>  uRx.<a class="code" href="structComplexMatrix2.html#a3a2b60a9e00ccef054a6de3b9307c1be">imag</a>[0][1] *= -1; uRx.<a class="code" href="structComplexMatrix2.html#a3a2b60a9e00ccef054a6de3b9307c1be">imag</a>[1][0] *= -1;</div>
<div class="line"><a name="l00485"></a><span class="lineno"> 485</span>  uRy.<a class="code" href="structComplexMatrix2.html#a3ea98d545b67a190ee6abe936c72076f">real</a>[0][1] *= -1; uRy.<a class="code" href="structComplexMatrix2.html#a3ea98d545b67a190ee6abe936c72076f">real</a>[1][0] *= -1;</div>
<div class="line"><a name="l00486"></a><span class="lineno"> 486</span>  <span class="keywordflow">for</span> (<span class="keywordtype">int</span> t=0; t < numTargets; t++) {</div>
<div class="line"><a name="l00487"></a><span class="lineno"> 487</span>  <span class="keywordflow">if</span> (targetPaulis[t] == <a class="code" href="group__type.html#gga1c703cf89629e4e9c7023cd402d67028a7abac7bb5b71e17382014c443244ad5c">PAULI_X</a>)</div>
<div class="line"><a name="l00488"></a><span class="lineno"> 488</span>  <a class="code" href="QuEST__internal_8h.html#a600bda10463e2078ce746ec28f2d39c9">statevec_multiControlledUnitary</a>(qureg, ctrlMask, ctrlFlipMask, targetQubits[t], uRy);</div>
<div class="line"><a name="l00489"></a><span class="lineno"> 489</span>  <span class="keywordflow">if</span> (targetPaulis[t] == <a class="code" href="group__type.html#gga1c703cf89629e4e9c7023cd402d67028a2c64a9aaee27f6bcfef894d159e759bd">PAULI_Y</a>)</div>
<div class="line"><a name="l00490"></a><span class="lineno"> 490</span>  <a class="code" href="QuEST__internal_8h.html#a600bda10463e2078ce746ec28f2d39c9">statevec_multiControlledUnitary</a>(qureg, ctrlMask, ctrlFlipMask, targetQubits[t], uRx);</div>
<div class="line"><a name="l00491"></a><span class="lineno"> 491</span>  }</div>
<div class="line"><a name="l00492"></a><span class="lineno"> 492</span> }</div>
<div class="line"><a name="l00493"></a><span class="lineno"> 493</span>  </div>
<div class="line"><a name="l00494"></a><span class="lineno"> 494</span> <span class="comment">/* produces both pauli|qureg> or pauli * qureg (as a density matrix) */</span></div>
<div class="line"><a name="l00495"></a><span class="lineno"><a class="line" href="QuEST__common_8c.html#accae504dfcd6846159baa16a1e8263a0"> 495</a></span> <span class="keywordtype">void</span> <a class="code" href="QuEST__common_8c.html#accae504dfcd6846159baa16a1e8263a0">statevec_applyPauliProd</a>(<a class="code" href="structQureg.html">Qureg</a> workspace, <span class="keywordtype">int</span>* targetQubits, <span class="keyword">enum</span> <a class="code" href="group__type.html#ga1c703cf89629e4e9c7023cd402d67028">pauliOpType</a>* pauliCodes, <span class="keywordtype">int</span> numTargets) {</div>
<div class="line"><a name="l00496"></a><span class="lineno"> 496</span>  </div>
<div class="line"><a name="l00497"></a><span class="lineno"> 497</span>  <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i=0; i < numTargets; i++) {</div>
<div class="line"><a name="l00498"></a><span class="lineno"> 498</span>  <span class="comment">// (pauliCodes[i] == PAULI_I) applies no operation</span></div>
<div class="line"><a name="l00499"></a><span class="lineno"> 499</span>  <span class="keywordflow">if</span> (pauliCodes[i] == <a class="code" href="group__type.html#gga1c703cf89629e4e9c7023cd402d67028a7abac7bb5b71e17382014c443244ad5c">PAULI_X</a>)</div>
<div class="line"><a name="l00500"></a><span class="lineno"> 500</span>  <a class="code" href="QuEST__internal_8h.html#a1c82b5e914a8a39c5c7b782dbc16fa76">statevec_pauliX</a>(workspace, targetQubits[i]);</div>
<div class="line"><a name="l00501"></a><span class="lineno"> 501</span>  <span class="keywordflow">if</span> (pauliCodes[i] == <a class="code" href="group__type.html#gga1c703cf89629e4e9c7023cd402d67028a2c64a9aaee27f6bcfef894d159e759bd">PAULI_Y</a>)</div>
<div class="line"><a name="l00502"></a><span class="lineno"> 502</span>  <a class="code" href="QuEST__internal_8h.html#af1cca787f3b0e0c91ffb2495f9ad8d5a">statevec_pauliY</a>(workspace, targetQubits[i]);</div>
<div class="line"><a name="l00503"></a><span class="lineno"> 503</span>  <span class="keywordflow">if</span> (pauliCodes[i] == <a class="code" href="group__type.html#gga1c703cf89629e4e9c7023cd402d67028a085c85cff6ba79c52d8b1e61c42ddddf">PAULI_Z</a>)</div>
<div class="line"><a name="l00504"></a><span class="lineno"> 504</span>  <a class="code" href="QuEST__common_8c.html#ab14df201e5c05527dc77cfeb6f583042">statevec_pauliZ</a>(workspace, targetQubits[i]);</div>
<div class="line"><a name="l00505"></a><span class="lineno"> 505</span>  }</div>
<div class="line"><a name="l00506"></a><span class="lineno"> 506</span> }</div>
<div class="line"><a name="l00507"></a><span class="lineno"> 507</span>  </div>
<div class="line"><a name="l00508"></a><span class="lineno"> 508</span> <span class="comment">// <pauli> = <qureg|pauli|qureg> = qureg . pauli(qureg)</span></div>
<div class="line"><a name="l00509"></a><span class="lineno"><a class="line" href="QuEST__internal_8h.html#aef380928d1b0003411d6159ca6594cd5"> 509</a></span> <a class="code" href="group__type.html#ga7740e349b4f8bae6451547680f0ce2d6">qreal</a> <a class="code" href="QuEST__common_8c.html#aef380928d1b0003411d6159ca6594cd5">statevec_calcExpecPauliProd</a>(<a class="code" href="structQureg.html">Qureg</a> qureg, <span class="keywordtype">int</span>* targetQubits, <span class="keyword">enum</span> <a class="code" href="group__type.html#ga1c703cf89629e4e9c7023cd402d67028">pauliOpType</a>* pauliCodes, <span class="keywordtype">int</span> numTargets, <a class="code" href="structQureg.html">Qureg</a> workspace) {</div>
<div class="line"><a name="l00510"></a><span class="lineno"> 510</span>  </div>
<div class="line"><a name="l00511"></a><span class="lineno"> 511</span>  <a class="code" href="QuEST__internal_8h.html#aa4b77089145184d990848cd66993aa8b">statevec_cloneQureg</a>(workspace, qureg);</div>
<div class="line"><a name="l00512"></a><span class="lineno"> 512</span>  <a class="code" href="QuEST__common_8c.html#accae504dfcd6846159baa16a1e8263a0">statevec_applyPauliProd</a>(workspace, targetQubits, pauliCodes, numTargets);</div>
<div class="line"><a name="l00513"></a><span class="lineno"> 513</span>  </div>
<div class="line"><a name="l00514"></a><span class="lineno"> 514</span>  <span class="comment">// compute the expected value</span></div>
<div class="line"><a name="l00515"></a><span class="lineno"> 515</span>  <a class="code" href="group__type.html#ga7740e349b4f8bae6451547680f0ce2d6">qreal</a> value;</div>
<div class="line"><a name="l00516"></a><span class="lineno"> 516</span>  <span class="keywordflow">if</span> (qureg.<a class="code" href="structQureg.html#acf78445e9435d09f44f0cc832c6aee79">isDensityMatrix</a>)</div>
<div class="line"><a name="l00517"></a><span class="lineno"> 517</span>  value = <a class="code" href="QuEST__internal_8h.html#a6153547f245c05874161a105e9a2f02c">densmatr_calcTotalProb</a>(workspace); <span class="comment">// Trace(ops qureg)</span></div>
<div class="line"><a name="l00518"></a><span class="lineno"> 518</span>  <span class="keywordflow">else</span></div>
<div class="line"><a name="l00519"></a><span class="lineno"> 519</span>  value = <a class="code" href="QuEST__internal_8h.html#a7ebd3198a198f4cd20840f64fd8b84d0">statevec_calcInnerProduct</a>(workspace, qureg).<a class="code" href="structComplex.html#ab5b2e2eca02131fc74b289a83636cbe3">real</a>; <span class="comment">// <qureg|ops|qureg></span></div>
<div class="line"><a name="l00520"></a><span class="lineno"> 520</span>  </div>
<div class="line"><a name="l00521"></a><span class="lineno"> 521</span>  <span class="keywordflow">return</span> value;</div>
<div class="line"><a name="l00522"></a><span class="lineno"> 522</span> }</div>
<div class="line"><a name="l00523"></a><span class="lineno"> 523</span>  </div>
<div class="line"><a name="l00524"></a><span class="lineno"><a class="line" href="QuEST__internal_8h.html#ae0cfea24091c8336146aa84c2253992e"> 524</a></span> <a class="code" href="group__type.html#ga7740e349b4f8bae6451547680f0ce2d6">qreal</a> <a class="code" href="QuEST__common_8c.html#ae0cfea24091c8336146aa84c2253992e">statevec_calcExpecPauliSum</a>(<a class="code" href="structQureg.html">Qureg</a> qureg, <span class="keyword">enum</span> <a class="code" href="group__type.html#ga1c703cf89629e4e9c7023cd402d67028">pauliOpType</a>* allCodes, <a class="code" href="group__type.html#ga7740e349b4f8bae6451547680f0ce2d6">qreal</a>* termCoeffs, <span class="keywordtype">int</span> numSumTerms, <a class="code" href="structQureg.html">Qureg</a> workspace) {</div>
<div class="line"><a name="l00525"></a><span class="lineno"> 525</span>  </div>
<div class="line"><a name="l00526"></a><span class="lineno"> 526</span>  <span class="keywordtype">int</span> numQb = qureg.<a class="code" href="structQureg.html#ad08dff5316b8937f4b2a1417591543dc">numQubitsRepresented</a>;</div>
<div class="line"><a name="l00527"></a><span class="lineno"> 527</span>  <span class="keywordtype">int</span> targs[100]; <span class="comment">// [numQb];</span></div>
<div class="line"><a name="l00528"></a><span class="lineno"> 528</span>  <span class="keywordflow">for</span> (<span class="keywordtype">int</span> q=0; q < numQb; q++)</div>
<div class="line"><a name="l00529"></a><span class="lineno"> 529</span>  targs[q] = q;</div>
<div class="line"><a name="l00530"></a><span class="lineno"> 530</span>  </div>
<div class="line"><a name="l00531"></a><span class="lineno"> 531</span>  <a class="code" href="group__type.html#ga7740e349b4f8bae6451547680f0ce2d6">qreal</a> value = 0;</div>
<div class="line"><a name="l00532"></a><span class="lineno"> 532</span>  <span class="keywordflow">for</span> (<span class="keywordtype">int</span> t=0; t < numSumTerms; t++)</div>
<div class="line"><a name="l00533"></a><span class="lineno"> 533</span>  value += termCoeffs[t] * <a class="code" href="QuEST__common_8c.html#aef380928d1b0003411d6159ca6594cd5">statevec_calcExpecPauliProd</a>(qureg, targs, &allCodes[t*numQb], numQb, workspace);</div>
<div class="line"><a name="l00534"></a><span class="lineno"> 534</span>  </div>
<div class="line"><a name="l00535"></a><span class="lineno"> 535</span>  <span class="keywordflow">return</span> value;</div>
<div class="line"><a name="l00536"></a><span class="lineno"> 536</span> }</div>
<div class="line"><a name="l00537"></a><span class="lineno"> 537</span>  </div>
<div class="line"><a name="l00538"></a><span class="lineno"><a class="line" href="QuEST__internal_8h.html#af71c2470694795a8d639be8ff9c34ed5"> 538</a></span> <span class="keywordtype">void</span> <a class="code" href="QuEST__common_8c.html#af71c2470694795a8d639be8ff9c34ed5">statevec_applyPauliSum</a>(<a class="code" href="structQureg.html">Qureg</a> inQureg, <span class="keyword">enum</span> <a class="code" href="group__type.html#ga1c703cf89629e4e9c7023cd402d67028">pauliOpType</a>* allCodes, <a class="code" href="group__type.html#ga7740e349b4f8bae6451547680f0ce2d6">qreal</a>* termCoeffs, <span class="keywordtype">int</span> numSumTerms, <a class="code" href="structQureg.html">Qureg</a> outQureg) {</div>
<div class="line"><a name="l00539"></a><span class="lineno"> 539</span>  </div>
<div class="line"><a name="l00540"></a><span class="lineno"> 540</span>  <span class="keywordtype">int</span> numQb = inQureg.<a class="code" href="structQureg.html#ad08dff5316b8937f4b2a1417591543dc">numQubitsRepresented</a>;</div>
<div class="line"><a name="l00541"></a><span class="lineno"> 541</span>  <span class="keywordtype">int</span> targs[100]; <span class="comment">// [numQb];</span></div>
<div class="line"><a name="l00542"></a><span class="lineno"> 542</span>  <span class="keywordflow">for</span> (<span class="keywordtype">int</span> q=0; q < numQb; q++)</div>
<div class="line"><a name="l00543"></a><span class="lineno"> 543</span>  targs[q] = q;</div>
<div class="line"><a name="l00544"></a><span class="lineno"> 544</span>  </div>
<div class="line"><a name="l00545"></a><span class="lineno"> 545</span>  <a class="code" href="QuEST__internal_8h.html#a758bad4237ff0bf3b4ff5be626a982ae">statevec_initBlankState</a>(outQureg);</div>
<div class="line"><a name="l00546"></a><span class="lineno"> 546</span>  </div>
<div class="line"><a name="l00547"></a><span class="lineno"> 547</span>  <span class="keywordflow">for</span> (<span class="keywordtype">int</span> t=0; t < numSumTerms; t++) {</div>
<div class="line"><a name="l00548"></a><span class="lineno"> 548</span>  <a class="code" href="structComplex.html">Complex</a> coef = (<a class="code" href="structComplex.html">Complex</a>) {.<a class="code" href="structComplex.html#ab5b2e2eca02131fc74b289a83636cbe3">real</a>=termCoeffs[t], .imag=0};</div>
<div class="line"><a name="l00549"></a><span class="lineno"> 549</span>  <a class="code" href="structComplex.html">Complex</a> iden = (<a class="code" href="structComplex.html">Complex</a>) {.<a class="code" href="structComplex.html#ab5b2e2eca02131fc74b289a83636cbe3">real</a>=1, .imag=0};</div>
<div class="line"><a name="l00550"></a><span class="lineno"> 550</span>  <a class="code" href="structComplex.html">Complex</a> zero = (<a class="code" href="structComplex.html">Complex</a>) {.<a class="code" href="structComplex.html#ab5b2e2eca02131fc74b289a83636cbe3">real</a>=0, .imag=0};</div>
<div class="line"><a name="l00551"></a><span class="lineno"> 551</span>  </div>
<div class="line"><a name="l00552"></a><span class="lineno"> 552</span>  <span class="comment">// outQureg += coef paulis(inQureg)</span></div>
<div class="line"><a name="l00553"></a><span class="lineno"> 553</span>  <a class="code" href="QuEST__common_8c.html#accae504dfcd6846159baa16a1e8263a0">statevec_applyPauliProd</a>(inQureg, targs, &allCodes[t*numQb], numQb);</div>
<div class="line"><a name="l00554"></a><span class="lineno"> 554</span>  <a class="code" href="QuEST__internal_8h.html#ad64a23933361fafc1b6b1787759c2663">statevec_setWeightedQureg</a>(coef, inQureg, iden, outQureg, zero, outQureg); </div>
<div class="line"><a name="l00555"></a><span class="lineno"> 555</span>  </div>
<div class="line"><a name="l00556"></a><span class="lineno"> 556</span>  <span class="comment">// undero paulis(inQureg), exploiting XX=YY=ZZ=I</span></div>
<div class="line"><a name="l00557"></a><span class="lineno"> 557</span>  <a class="code" href="QuEST__common_8c.html#accae504dfcd6846159baa16a1e8263a0">statevec_applyPauliProd</a>(inQureg, targs, &allCodes[t*numQb], numQb);</div>
<div class="line"><a name="l00558"></a><span class="lineno"> 558</span>  }</div>
<div class="line"><a name="l00559"></a><span class="lineno"> 559</span> }</div>
<div class="line"><a name="l00560"></a><span class="lineno"> 560</span>  </div>
<div class="line"><a name="l00561"></a><span class="lineno"><a class="line" href="QuEST__internal_8h.html#a95abbe0cedf45b1ff2f070dba219fd7f"> 561</a></span> <span class="keywordtype">void</span> <a class="code" href="QuEST__common_8c.html#a95abbe0cedf45b1ff2f070dba219fd7f">statevec_twoQubitUnitary</a>(<a class="code" href="structQureg.html">Qureg</a> qureg, <span class="keywordtype">int</span> targetQubit1, <span class="keywordtype">int</span> targetQubit2, <a class="code" href="structComplexMatrix4.html">ComplexMatrix4</a> u) {</div>
<div class="line"><a name="l00562"></a><span class="lineno"> 562</span>  </div>
<div class="line"><a name="l00563"></a><span class="lineno"> 563</span>  <span class="keywordtype">long</span> <span class="keywordtype">long</span> <span class="keywordtype">int</span> ctrlMask = 0;</div>
<div class="line"><a name="l00564"></a><span class="lineno"> 564</span>  <a class="code" href="QuEST__internal_8h.html#a3eadcb581f2623e21e4967f791d0439c">statevec_multiControlledTwoQubitUnitary</a>(qureg, ctrlMask, targetQubit1, targetQubit2, u);</div>
<div class="line"><a name="l00565"></a><span class="lineno"> 565</span> }</div>
<div class="line"><a name="l00566"></a><span class="lineno"> 566</span>  </div>
<div class="line"><a name="l00567"></a><span class="lineno"><a class="line" href="QuEST__internal_8h.html#ad7ac4934371d72e24c6d1cea3d3da391"> 567</a></span> <span class="keywordtype">void</span> <a class="code" href="QuEST__common_8c.html#ad7ac4934371d72e24c6d1cea3d3da391">statevec_controlledTwoQubitUnitary</a>(<a class="code" href="structQureg.html">Qureg</a> qureg, <span class="keywordtype">int</span> controlQubit, <span class="keywordtype">int</span> targetQubit1, <span class="keywordtype">int</span> targetQubit2, <a class="code" href="structComplexMatrix4.html">ComplexMatrix4</a> u) {</div>
<div class="line"><a name="l00568"></a><span class="lineno"> 568</span>  </div>
<div class="line"><a name="l00569"></a><span class="lineno"> 569</span>  <span class="keywordtype">long</span> <span class="keywordtype">long</span> <span class="keywordtype">int</span> ctrlMask = 1LL << controlQubit;</div>
<div class="line"><a name="l00570"></a><span class="lineno"> 570</span>  <a class="code" href="QuEST__internal_8h.html#a3eadcb581f2623e21e4967f791d0439c">statevec_multiControlledTwoQubitUnitary</a>(qureg, ctrlMask, targetQubit1, targetQubit2, u);</div>
<div class="line"><a name="l00571"></a><span class="lineno"> 571</span> }</div>
<div class="line"><a name="l00572"></a><span class="lineno"> 572</span>  </div>
<div class="line"><a name="l00573"></a><span class="lineno"><a class="line" href="QuEST__internal_8h.html#a7a76950824f17e821415070a4bfdad5b"> 573</a></span> <span class="keywordtype">void</span> <a class="code" href="QuEST__common_8c.html#a7a76950824f17e821415070a4bfdad5b">statevec_multiQubitUnitary</a>(<a class="code" href="structQureg.html">Qureg</a> qureg, <span class="keywordtype">int</span>* targets, <span class="keywordtype">int</span> numTargets, <a class="code" href="structComplexMatrixN.html">ComplexMatrixN</a> u) {</div>
<div class="line"><a name="l00574"></a><span class="lineno"> 574</span>  </div>
<div class="line"><a name="l00575"></a><span class="lineno"> 575</span>  <span class="keywordtype">long</span> <span class="keywordtype">long</span> <span class="keywordtype">int</span> ctrlMask = 0;</div>
<div class="line"><a name="l00576"></a><span class="lineno"> 576</span>  <a class="code" href="QuEST__internal_8h.html#a0acfe30083d0bc2a84da196378a36122">statevec_multiControlledMultiQubitUnitary</a>(qureg, ctrlMask, targets, numTargets, u);</div>
<div class="line"><a name="l00577"></a><span class="lineno"> 577</span> }</div>
<div class="line"><a name="l00578"></a><span class="lineno"> 578</span>  </div>
<div class="line"><a name="l00579"></a><span class="lineno"><a class="line" href="QuEST__internal_8h.html#a06b8474aec51579ff680fb66439e7fc6"> 579</a></span> <span class="keywordtype">void</span> <a class="code" href="QuEST__common_8c.html#a06b8474aec51579ff680fb66439e7fc6">statevec_controlledMultiQubitUnitary</a>(<a class="code" href="structQureg.html">Qureg</a> qureg, <span class="keywordtype">int</span> ctrl, <span class="keywordtype">int</span>* targets, <span class="keywordtype">int</span> numTargets, <a class="code" href="structComplexMatrixN.html">ComplexMatrixN</a> u) {</div>
<div class="line"><a name="l00580"></a><span class="lineno"> 580</span>  </div>
<div class="line"><a name="l00581"></a><span class="lineno"> 581</span>  <span class="keywordtype">long</span> <span class="keywordtype">long</span> <span class="keywordtype">int</span> ctrlMask = 1LL << ctrl;</div>
<div class="line"><a name="l00582"></a><span class="lineno"> 582</span>  <a class="code" href="QuEST__internal_8h.html#a0acfe30083d0bc2a84da196378a36122">statevec_multiControlledMultiQubitUnitary</a>(qureg, ctrlMask, targets, numTargets, u);</div>
<div class="line"><a name="l00583"></a><span class="lineno"> 583</span> }</div>
<div class="line"><a name="l00584"></a><span class="lineno"> 584</span>  </div>
<div class="line"><a name="l00585"></a><span class="lineno"><a class="line" href="QuEST__common_8c.html#ab9babafb4a8cc2af97e3d815eadec9db"> 585</a></span> <span class="preprocessor">#define macro_populateKrausOperator(superOp, ops, numOps, opDim) \</span></div>
<div class="line"><a name="l00586"></a><span class="lineno"> 586</span> <span class="preprocessor"> </span><span class="comment">/* clear the superop */</span><span class="preprocessor"> \</span></div>
<div class="line"><a name="l00587"></a><span class="lineno"> 587</span> <span class="preprocessor"> for (int r=0; r < (opDim)*(opDim); r++) \</span></div>
<div class="line"><a name="l00588"></a><span class="lineno"> 588</span> <span class="preprocessor"> for (int c=0; c < (opDim)*(opDim); c++) { \</span></div>
<div class="line"><a name="l00589"></a><span class="lineno"> 589</span> <span class="preprocessor"> superOp->real[r][c] = 0; \</span></div>
<div class="line"><a name="l00590"></a><span class="lineno"> 590</span> <span class="preprocessor"> superOp->imag[r][c] = 0; \</span></div>
<div class="line"><a name="l00591"></a><span class="lineno"> 591</span> <span class="preprocessor"> } \</span></div>
<div class="line"><a name="l00592"></a><span class="lineno"> 592</span> <span class="preprocessor"> </span><span class="comment">/* add each op's contribution to the superop */</span><span class="preprocessor"> \</span></div>
<div class="line"><a name="l00593"></a><span class="lineno"> 593</span> <span class="preprocessor"> for (int n=0; n<(numOps); n++) \</span></div>
<div class="line"><a name="l00594"></a><span class="lineno"> 594</span> <span class="preprocessor"> </span><span class="comment">/* superop += conjugate(op) (x) op, where (x) is a tensor product */</span><span class="preprocessor"> \</span></div>
<div class="line"><a name="l00595"></a><span class="lineno"> 595</span> <span class="preprocessor"> for (int i = 0; i < (opDim); i++) \</span></div>
<div class="line"><a name="l00596"></a><span class="lineno"> 596</span> <span class="preprocessor"> for (int j = 0; j < (opDim); j++) \</span></div>
<div class="line"><a name="l00597"></a><span class="lineno"> 597</span> <span class="preprocessor"> for (int k = 0; k < (opDim); k++) \</span></div>
<div class="line"><a name="l00598"></a><span class="lineno"> 598</span> <span class="preprocessor"> for (int l = 0; l < (opDim); l++) { \</span></div>
<div class="line"><a name="l00599"></a><span class="lineno"> 599</span> <span class="preprocessor"> superOp->real[i*(opDim) + k][j*(opDim) + l] += \</span></div>
<div class="line"><a name="l00600"></a><span class="lineno"> 600</span> <span class="preprocessor"> ops[n].real[i][j]*ops[n].real[k][l] + \</span></div>
<div class="line"><a name="l00601"></a><span class="lineno"> 601</span> <span class="preprocessor"> ops[n].imag[i][j]*ops[n].imag[k][l]; \</span></div>
<div class="line"><a name="l00602"></a><span class="lineno"> 602</span> <span class="preprocessor"> superOp->imag[i*(opDim) + k][j*(opDim) + l] += \</span></div>
<div class="line"><a name="l00603"></a><span class="lineno"> 603</span> <span class="preprocessor"> ops[n].real[i][j]*ops[n].imag[k][l] - \</span></div>
<div class="line"><a name="l00604"></a><span class="lineno"> 604</span> <span class="preprocessor"> ops[n].imag[i][j]*ops[n].real[k][l]; \</span></div>
<div class="line"><a name="l00605"></a><span class="lineno"> 605</span> <span class="preprocessor"> } </span></div>
<div class="line"><a name="l00606"></a><span class="lineno"> 606</span>  </div>
<div class="line"><a name="l00607"></a><span class="lineno"><a class="line" href="QuEST__common_8c.html#abd00ba1321c9bd43c3a9119f5b3ed1be"> 607</a></span> <span class="keywordtype">void</span> <a class="code" href="QuEST__common_8c.html#abd00ba1321c9bd43c3a9119f5b3ed1be">populateKrausSuperOperator2</a>(<a class="code" href="structComplexMatrix4.html">ComplexMatrix4</a>* superOp, <a class="code" href="structComplexMatrix2.html">ComplexMatrix2</a>* ops, <span class="keywordtype">int</span> numOps) {</div>
<div class="line"><a name="l00608"></a><span class="lineno"> 608</span>  <span class="keywordtype">int</span> opDim = 2;</div>
<div class="line"><a name="l00609"></a><span class="lineno"> 609</span>  <a class="code" href="QuEST__common_8c.html#ab9babafb4a8cc2af97e3d815eadec9db">macro_populateKrausOperator</a>(superOp, ops, numOps, opDim);</div>
<div class="line"><a name="l00610"></a><span class="lineno"> 610</span> }</div>
<div class="line"><a name="l00611"></a><span class="lineno"><a class="line" href="QuEST__common_8c.html#a2a6999e47bed978eb66f589c8abc7766"> 611</a></span> <span class="keywordtype">void</span> <a class="code" href="QuEST__common_8c.html#a2a6999e47bed978eb66f589c8abc7766">populateKrausSuperOperator4</a>(<a class="code" href="structComplexMatrixN.html">ComplexMatrixN</a>* superOp, <a class="code" href="structComplexMatrix4.html">ComplexMatrix4</a>* ops, <span class="keywordtype">int</span> numOps) {</div>
<div class="line"><a name="l00612"></a><span class="lineno"> 612</span>  <span class="keywordtype">int</span> opDim = 4;</div>
<div class="line"><a name="l00613"></a><span class="lineno"> 613</span>  <a class="code" href="QuEST__common_8c.html#ab9babafb4a8cc2af97e3d815eadec9db">macro_populateKrausOperator</a>(superOp, ops, numOps, opDim);</div>
<div class="line"><a name="l00614"></a><span class="lineno"> 614</span> }</div>
<div class="line"><a name="l00615"></a><span class="lineno"><a class="line" href="QuEST__common_8c.html#a6dcb85ef8e1257813af0dc306082caf8"> 615</a></span> <span class="keywordtype">void</span> <a class="code" href="QuEST__common_8c.html#a6dcb85ef8e1257813af0dc306082caf8">populateKrausSuperOperatorN</a>(<a class="code" href="structComplexMatrixN.html">ComplexMatrixN</a>* superOp, <a class="code" href="structComplexMatrixN.html">ComplexMatrixN</a>* ops, <span class="keywordtype">int</span> numOps) {</div>
<div class="line"><a name="l00616"></a><span class="lineno"> 616</span>  <span class="keywordtype">int</span> opDim = 1 << ops[0].<a class="code" href="structComplexMatrixN.html#a606f7afbc3f65a596bdca21cb43be6a8">numQubits</a>;</div>
<div class="line"><a name="l00617"></a><span class="lineno"> 617</span>  <a class="code" href="QuEST__common_8c.html#ab9babafb4a8cc2af97e3d815eadec9db">macro_populateKrausOperator</a>(superOp, ops, numOps, opDim);</div>
<div class="line"><a name="l00618"></a><span class="lineno"> 618</span> }</div>
<div class="line"><a name="l00619"></a><span class="lineno"> 619</span>  </div>
<div class="line"><a name="l00620"></a><span class="lineno"><a class="line" href="QuEST__common_8c.html#af3db21114efc431ae7a7ae06fc2d7b5b"> 620</a></span> <span class="keywordtype">void</span> <a class="code" href="QuEST__common_8c.html#af3db21114efc431ae7a7ae06fc2d7b5b">densmatr_applyKrausSuperoperator</a>(<a class="code" href="structQureg.html">Qureg</a> qureg, <span class="keywordtype">int</span> target, <a class="code" href="structComplexMatrix4.html">ComplexMatrix4</a> superOp) {</div>
<div class="line"><a name="l00621"></a><span class="lineno"> 621</span>  </div>
<div class="line"><a name="l00622"></a><span class="lineno"> 622</span>  <span class="keywordtype">long</span> <span class="keywordtype">long</span> <span class="keywordtype">int</span> ctrlMask = 0;</div>
<div class="line"><a name="l00623"></a><span class="lineno"> 623</span>  <a class="code" href="QuEST__internal_8h.html#a3eadcb581f2623e21e4967f791d0439c">statevec_multiControlledTwoQubitUnitary</a>(qureg, ctrlMask, target, target + qureg.<a class="code" href="structQureg.html#ad08dff5316b8937f4b2a1417591543dc">numQubitsRepresented</a>, superOp);</div>
<div class="line"><a name="l00624"></a><span class="lineno"> 624</span> }</div>
<div class="line"><a name="l00625"></a><span class="lineno"> 625</span>  </div>
<div class="line"><a name="l00626"></a><span class="lineno"><a class="line" href="QuEST__common_8c.html#a3d39314af5d88822d4aefb5e9cc4e2f6"> 626</a></span> <span class="keywordtype">void</span> <a class="code" href="QuEST__common_8c.html#a3d39314af5d88822d4aefb5e9cc4e2f6">densmatr_applyTwoQubitKrausSuperoperator</a>(<a class="code" href="structQureg.html">Qureg</a> qureg, <span class="keywordtype">int</span> target1, <span class="keywordtype">int</span> target2, <a class="code" href="structComplexMatrixN.html">ComplexMatrixN</a> superOp) {</div>
<div class="line"><a name="l00627"></a><span class="lineno"> 627</span>  </div>
<div class="line"><a name="l00628"></a><span class="lineno"> 628</span>  <span class="keywordtype">long</span> <span class="keywordtype">long</span> <span class="keywordtype">int</span> ctrlMask = 0;</div>
<div class="line"><a name="l00629"></a><span class="lineno"> 629</span>  <span class="keywordtype">int</span> numQb = qureg.<a class="code" href="structQureg.html#ad08dff5316b8937f4b2a1417591543dc">numQubitsRepresented</a>;</div>
<div class="line"><a name="l00630"></a><span class="lineno"> 630</span>  <span class="keywordtype">int</span> allTargets[4] = {target1, target2, target1+numQb, target2+numQb};</div>
<div class="line"><a name="l00631"></a><span class="lineno"> 631</span>  <a class="code" href="QuEST__internal_8h.html#a0acfe30083d0bc2a84da196378a36122">statevec_multiControlledMultiQubitUnitary</a>(qureg, ctrlMask, allTargets, 4, superOp);</div>
<div class="line"><a name="l00632"></a><span class="lineno"> 632</span> }</div>
<div class="line"><a name="l00633"></a><span class="lineno"> 633</span>  </div>
<div class="line"><a name="l00634"></a><span class="lineno"><a class="line" href="QuEST__common_8c.html#a514ca44595676abf4bc2e10fc48ef970"> 634</a></span> <span class="keywordtype">void</span> <a class="code" href="QuEST__common_8c.html#a514ca44595676abf4bc2e10fc48ef970">densmatr_applyMultiQubitKrausSuperoperator</a>(<a class="code" href="structQureg.html">Qureg</a> qureg, <span class="keywordtype">int</span> *targets, <span class="keywordtype">int</span> numTargets, <a class="code" href="structComplexMatrixN.html">ComplexMatrixN</a> superOp) {</div>
<div class="line"><a name="l00635"></a><span class="lineno"> 635</span>  <span class="keywordtype">long</span> <span class="keywordtype">long</span> <span class="keywordtype">int</span> ctrlMask = 0;</div>
<div class="line"><a name="l00636"></a><span class="lineno"> 636</span>  <span class="keywordtype">int</span> allTargets[200]; <span class="comment">// [2*numTargets];</span></div>
<div class="line"><a name="l00637"></a><span class="lineno"> 637</span>  <span class="keywordflow">for</span> (<span class="keywordtype">int</span> t=0; t < numTargets; t++) {</div>
<div class="line"><a name="l00638"></a><span class="lineno"> 638</span>  allTargets[t] = targets[t];</div>
<div class="line"><a name="l00639"></a><span class="lineno"> 639</span>  allTargets[t+numTargets] = targets[t] + qureg.<a class="code" href="structQureg.html#ad08dff5316b8937f4b2a1417591543dc">numQubitsRepresented</a>;</div>
<div class="line"><a name="l00640"></a><span class="lineno"> 640</span>  }</div>
<div class="line"><a name="l00641"></a><span class="lineno"> 641</span>  <a class="code" href="QuEST__internal_8h.html#a0acfe30083d0bc2a84da196378a36122">statevec_multiControlledMultiQubitUnitary</a>(qureg, ctrlMask, allTargets, 2*numTargets, superOp);</div>
<div class="line"><a name="l00642"></a><span class="lineno"> 642</span> }</div>
<div class="line"><a name="l00643"></a><span class="lineno"> 643</span>  </div>
<div class="line"><a name="l00644"></a><span class="lineno"><a class="line" href="QuEST__internal_8h.html#ab987904f721746d1176aba7ac6b566fa"> 644</a></span> <span class="keywordtype">void</span> <a class="code" href="QuEST__common_8c.html#ab987904f721746d1176aba7ac6b566fa">densmatr_mixKrausMap</a>(<a class="code" href="structQureg.html">Qureg</a> qureg, <span class="keywordtype">int</span> target, <a class="code" href="structComplexMatrix2.html">ComplexMatrix2</a> *ops, <span class="keywordtype">int</span> numOps) {</div>
<div class="line"><a name="l00645"></a><span class="lineno"> 645</span>  </div>
<div class="line"><a name="l00646"></a><span class="lineno"> 646</span>  <a class="code" href="structComplexMatrix4.html">ComplexMatrix4</a> superOp; </div>
<div class="line"><a name="l00647"></a><span class="lineno"> 647</span>  <a class="code" href="QuEST__common_8c.html#abd00ba1321c9bd43c3a9119f5b3ed1be">populateKrausSuperOperator2</a>(&superOp, ops, numOps);</div>
<div class="line"><a name="l00648"></a><span class="lineno"> 648</span>  <a class="code" href="QuEST__common_8c.html#af3db21114efc431ae7a7ae06fc2d7b5b">densmatr_applyKrausSuperoperator</a>(qureg, target, superOp);</div>
<div class="line"><a name="l00649"></a><span class="lineno"> 649</span> }</div>
<div class="line"><a name="l00650"></a><span class="lineno"> 650</span>  </div>
<div class="line"><a name="l00651"></a><span class="lineno"> 651</span> <span class="preprocessor">#ifndef _WIN32</span></div>
<div class="line"><a name="l00652"></a><span class="lineno"><a class="line" href="QuEST__common_8c.html#ad67cb67877c941189b0de2c2e67e2026"> 652</a></span> <a class="code" href="structComplexMatrixN.html">ComplexMatrixN</a> <a class="code" href="QuEST__common_8c.html#ad67cb67877c941189b0de2c2e67e2026">bindArraysToStackComplexMatrixN</a>(</div>
<div class="line"><a name="l00653"></a><span class="lineno"> 653</span>  <span class="keywordtype">int</span> numQubits, <a class="code" href="group__type.html#ga7740e349b4f8bae6451547680f0ce2d6">qreal</a> re[][1<<numQubits], <a class="code" href="group__type.html#ga7740e349b4f8bae6451547680f0ce2d6">qreal</a> im[][1<<numQubits], </div>
<div class="line"><a name="l00654"></a><span class="lineno"> 654</span>  <a class="code" href="group__type.html#ga7740e349b4f8bae6451547680f0ce2d6">qreal</a>** reStorage, <a class="code" href="group__type.html#ga7740e349b4f8bae6451547680f0ce2d6">qreal</a>** imStorage</div>
<div class="line"><a name="l00655"></a><span class="lineno"> 655</span> ) {</div>
<div class="line"><a name="l00656"></a><span class="lineno"> 656</span>  <a class="code" href="structComplexMatrixN.html">ComplexMatrixN</a> m;</div>
<div class="line"><a name="l00657"></a><span class="lineno"> 657</span>  m.<a class="code" href="structComplexMatrixN.html#a606f7afbc3f65a596bdca21cb43be6a8">numQubits</a> = numQubits;</div>
<div class="line"><a name="l00658"></a><span class="lineno"> 658</span>  m.<a class="code" href="structComplexMatrixN.html#a446e8a077597e4f76049d547df0c9738">real</a> = reStorage;</div>
<div class="line"><a name="l00659"></a><span class="lineno"> 659</span>  m.<a class="code" href="structComplexMatrixN.html#afde13e469a3b4f179f45e53d1ec0a7cf">imag</a> = imStorage;</div>
<div class="line"><a name="l00660"></a><span class="lineno"> 660</span>  </div>
<div class="line"><a name="l00661"></a><span class="lineno"> 661</span>  <span class="keywordtype">int</span> len = 1<<numQubits;</div>
<div class="line"><a name="l00662"></a><span class="lineno"> 662</span>  <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i=0; i<len; i++) {</div>
<div class="line"><a name="l00663"></a><span class="lineno"> 663</span>  m.<a class="code" href="structComplexMatrixN.html#a446e8a077597e4f76049d547df0c9738">real</a>[i] = re[i];</div>
<div class="line"><a name="l00664"></a><span class="lineno"> 664</span>  m.<a class="code" href="structComplexMatrixN.html#afde13e469a3b4f179f45e53d1ec0a7cf">imag</a>[i] = im[i];</div>
<div class="line"><a name="l00665"></a><span class="lineno"> 665</span>  }</div>
<div class="line"><a name="l00666"></a><span class="lineno"> 666</span>  <span class="keywordflow">return</span> m;</div>
<div class="line"><a name="l00667"></a><span class="lineno"> 667</span> }</div>
<div class="line"><a name="l00668"></a><span class="lineno"> 668</span>  </div>
<div class="line"><a name="l00669"></a><span class="lineno"><a class="line" href="QuEST__common_8c.html#a65b0614efc2e0237d56f86ccec182294"> 669</a></span> <span class="preprocessor">#define macro_initialiseStackComplexMatrixN(matrix, numQubits, real, imag) \</span></div>
<div class="line"><a name="l00670"></a><span class="lineno"> 670</span> <span class="preprocessor"> </span><span class="comment">/* reStorage_ and imStorage_ must not exist in calling scope */</span><span class="preprocessor"> \</span></div>
<div class="line"><a name="l00671"></a><span class="lineno"> 671</span> <span class="preprocessor"> qreal* reStorage_[1<<(numQubits)]; \</span></div>
<div class="line"><a name="l00672"></a><span class="lineno"> 672</span> <span class="preprocessor"> qreal* imStorage_[1<<(numQubits)]; \</span></div>
<div class="line"><a name="l00673"></a><span class="lineno"> 673</span> <span class="preprocessor"> matrix = bindArraysToStackComplexMatrixN((numQubits), real, imag, reStorage_, imStorage_);</span></div>
<div class="line"><a name="l00674"></a><span class="lineno"> 674</span>  </div>
<div class="line"><a name="l00675"></a><span class="lineno"><a class="line" href="QuEST__common_8c.html#a4da103b16292ffd99197d05c5fa83712"> 675</a></span> <span class="preprocessor">#define macro_allocStackComplexMatrixN(matrix, numQubits) \</span></div>
<div class="line"><a name="l00676"></a><span class="lineno"> 676</span> <span class="preprocessor"> </span><span class="comment">/* reArr_, imArr_, reStorage_, and imStorage_ must not exist in calling scope */</span><span class="preprocessor"> \</span></div>
<div class="line"><a name="l00677"></a><span class="lineno"> 677</span> <span class="preprocessor"> qreal reArr_[1<<(numQubits)][1<<(numQubits)]; \</span></div>
<div class="line"><a name="l00678"></a><span class="lineno"> 678</span> <span class="preprocessor"> qreal imArr_[1<<(numQubits)][1<<(numQubits)]; \</span></div>
<div class="line"><a name="l00679"></a><span class="lineno"> 679</span> <span class="preprocessor"> macro_initialiseStackComplexMatrixN(matrix, (numQubits), reArr_, imArr_);</span></div>
<div class="line"><a name="l00680"></a><span class="lineno"> 680</span> <span class="preprocessor">#endif</span></div>
<div class="line"><a name="l00681"></a><span class="lineno"> 681</span>  </div>
<div class="line"><a name="l00682"></a><span class="lineno"><a class="line" href="QuEST__internal_8h.html#ab004038b5f906ea897621c88bf213129"> 682</a></span> <span class="keywordtype">void</span> <a class="code" href="QuEST__common_8c.html#ab004038b5f906ea897621c88bf213129">densmatr_mixTwoQubitKrausMap</a>(<a class="code" href="structQureg.html">Qureg</a> qureg, <span class="keywordtype">int</span> target1, <span class="keywordtype">int</span> target2, <a class="code" href="structComplexMatrix4.html">ComplexMatrix4</a> *ops, <span class="keywordtype">int</span> numOps) {</div>
<div class="line"><a name="l00683"></a><span class="lineno"> 683</span>  </div>
<div class="line"><a name="l00684"></a><span class="lineno"> 684</span>  <span class="comment">// if NOT on Windows, allocate ComplexN on stack</span></div>
<div class="line"><a name="l00685"></a><span class="lineno"> 685</span> <span class="preprocessor"> #ifndef _WIN32</span></div>
<div class="line"><a name="l00686"></a><span class="lineno"> 686</span>  <a class="code" href="structComplexMatrixN.html">ComplexMatrixN</a> superOp;</div>
<div class="line"><a name="l00687"></a><span class="lineno"> 687</span>  <a class="code" href="QuEST__common_8c.html#a4da103b16292ffd99197d05c5fa83712">macro_allocStackComplexMatrixN</a>(superOp, 4);</div>
<div class="line"><a name="l00688"></a><span class="lineno"> 688</span>  <a class="code" href="QuEST__common_8c.html#a2a6999e47bed978eb66f589c8abc7766">populateKrausSuperOperator4</a>(&superOp, ops, numOps);</div>
<div class="line"><a name="l00689"></a><span class="lineno"> 689</span>  <a class="code" href="QuEST__common_8c.html#a3d39314af5d88822d4aefb5e9cc4e2f6">densmatr_applyTwoQubitKrausSuperoperator</a>(qureg, target1, target2, superOp);</div>
<div class="line"><a name="l00690"></a><span class="lineno"> 690</span>  </div>
<div class="line"><a name="l00691"></a><span class="lineno"> 691</span>  <span class="comment">// but on Windows, we MUST allocated dynamically</span></div>
<div class="line"><a name="l00692"></a><span class="lineno"> 692</span> <span class="preprocessor"> #else</span></div>
<div class="line"><a name="l00693"></a><span class="lineno"> 693</span>  <a class="code" href="structComplexMatrixN.html">ComplexMatrixN</a> superOp = <a class="code" href="group__type.html#ga815103261fb22ea9690e1427571df00e">createComplexMatrixN</a>(4);</div>
<div class="line"><a name="l00694"></a><span class="lineno"> 694</span>  <a class="code" href="QuEST__common_8c.html#a2a6999e47bed978eb66f589c8abc7766">populateKrausSuperOperator4</a>(&superOp, ops, numOps);</div>
<div class="line"><a name="l00695"></a><span class="lineno"> 695</span>  <a class="code" href="QuEST__common_8c.html#a3d39314af5d88822d4aefb5e9cc4e2f6">densmatr_applyTwoQubitKrausSuperoperator</a>(qureg, target1, target2, superOp);</div>
<div class="line"><a name="l00696"></a><span class="lineno"> 696</span>  <a class="code" href="group__type.html#ga9df2f3d86be4a6e9aad481665e5e6753">destroyComplexMatrixN</a>(superOp);</div>
<div class="line"><a name="l00697"></a><span class="lineno"> 697</span>  </div>
<div class="line"><a name="l00698"></a><span class="lineno"> 698</span> <span class="preprocessor"> #endif</span></div>
<div class="line"><a name="l00699"></a><span class="lineno"> 699</span> }</div>
<div class="line"><a name="l00700"></a><span class="lineno"> 700</span>  </div>
<div class="line"><a name="l00701"></a><span class="lineno"><a class="line" href="QuEST__internal_8h.html#a99702f4667c429969d7e7b4ecd45c2f3"> 701</a></span> <span class="keywordtype">void</span> <a class="code" href="QuEST__common_8c.html#a99702f4667c429969d7e7b4ecd45c2f3">densmatr_mixMultiQubitKrausMap</a>(<a class="code" href="structQureg.html">Qureg</a> qureg, <span class="keywordtype">int</span>* targets, <span class="keywordtype">int</span> numTargets, <a class="code" href="structComplexMatrixN.html">ComplexMatrixN</a>* ops, <span class="keywordtype">int</span> numOps) {</div>
<div class="line"><a name="l00702"></a><span class="lineno"> 702</span>  </div>
<div class="line"><a name="l00703"></a><span class="lineno"> 703</span>  <a class="code" href="structComplexMatrixN.html">ComplexMatrixN</a> superOp;</div>
<div class="line"><a name="l00704"></a><span class="lineno"> 704</span>  </div>
<div class="line"><a name="l00705"></a><span class="lineno"> 705</span>  <span class="comment">/* superOp will contain 2^(4 numTargets) complex numbers.</span></div>
<div class="line"><a name="l00706"></a><span class="lineno"> 706</span> <span class="comment"> * At double precision, superOp will cost additional memory:</span></div>
<div class="line"><a name="l00707"></a><span class="lineno"> 707</span> <span class="comment"> * numTargs=1 -> 0.25 KiB</span></div>
<div class="line"><a name="l00708"></a><span class="lineno"> 708</span> <span class="comment"> * numTargs=2 -> 4 KiB </span></div>
<div class="line"><a name="l00709"></a><span class="lineno"> 709</span> <span class="comment"> * numTargs=3 -> 64 KiB</span></div>
<div class="line"><a name="l00710"></a><span class="lineno"> 710</span> <span class="comment"> * numTargs=4 -> 1 MiB</span></div>
<div class="line"><a name="l00711"></a><span class="lineno"> 711</span> <span class="comment"> * numTargs=5 -> 16 MiB.</span></div>
<div class="line"><a name="l00712"></a><span class="lineno"> 712</span> <span class="comment"> * At quad precision (usually 10 B per number, but possibly 16 B due to alignment),</span></div>
<div class="line"><a name="l00713"></a><span class="lineno"> 713</span> <span class="comment"> * this costs at most double.</span></div>
<div class="line"><a name="l00714"></a><span class="lineno"> 714</span> <span class="comment"> *</span></div>
<div class="line"><a name="l00715"></a><span class="lineno"> 715</span> <span class="comment"> * Hence, if superOp is kept in the stack, numTargs >= 4 would exceed Windows' 1 MB </span></div>
<div class="line"><a name="l00716"></a><span class="lineno"> 716</span> <span class="comment"> * maximum stack-space allocation (numTargs >= 5 exceeding Linux' 8 MB). Therefore, </span></div>
<div class="line"><a name="l00717"></a><span class="lineno"> 717</span> <span class="comment"> * for numTargets < 4, superOp will be kept in the stack, else in the heap</span></div>
<div class="line"><a name="l00718"></a><span class="lineno"> 718</span> <span class="comment"> */</span></div>
<div class="line"><a name="l00719"></a><span class="lineno"> 719</span>  </div>
<div class="line"><a name="l00720"></a><span class="lineno"> 720</span>  <span class="comment">// if NOT on Windows, allocate ComplexN on stack depending on size</span></div>
<div class="line"><a name="l00721"></a><span class="lineno"> 721</span> <span class="preprocessor"> #ifndef _WIN32</span></div>
<div class="line"><a name="l00722"></a><span class="lineno"> 722</span>  <span class="keywordflow">if</span> (numTargets < 4) {</div>
<div class="line"><a name="l00723"></a><span class="lineno"> 723</span>  <span class="comment">// everything must live in 'if' since this macro declares local vars</span></div>
<div class="line"><a name="l00724"></a><span class="lineno"> 724</span>  <a class="code" href="QuEST__common_8c.html#a4da103b16292ffd99197d05c5fa83712">macro_allocStackComplexMatrixN</a>(superOp, 2*numTargets);</div>
<div class="line"><a name="l00725"></a><span class="lineno"> 725</span>  <a class="code" href="QuEST__common_8c.html#a6dcb85ef8e1257813af0dc306082caf8">populateKrausSuperOperatorN</a>(&superOp, ops, numOps);</div>
<div class="line"><a name="l00726"></a><span class="lineno"> 726</span>  <a class="code" href="QuEST__common_8c.html#a514ca44595676abf4bc2e10fc48ef970">densmatr_applyMultiQubitKrausSuperoperator</a>(qureg, targets, numTargets, superOp);</div>
<div class="line"><a name="l00727"></a><span class="lineno"> 727</span>  }</div>
<div class="line"><a name="l00728"></a><span class="lineno"> 728</span>  <span class="keywordflow">else</span> {</div>
<div class="line"><a name="l00729"></a><span class="lineno"> 729</span>  superOp = <a class="code" href="group__type.html#ga815103261fb22ea9690e1427571df00e">createComplexMatrixN</a>(2*numTargets);</div>
<div class="line"><a name="l00730"></a><span class="lineno"> 730</span>  <a class="code" href="QuEST__common_8c.html#a6dcb85ef8e1257813af0dc306082caf8">populateKrausSuperOperatorN</a>(&superOp, ops, numOps);</div>
<div class="line"><a name="l00731"></a><span class="lineno"> 731</span>  <a class="code" href="QuEST__common_8c.html#a514ca44595676abf4bc2e10fc48ef970">densmatr_applyMultiQubitKrausSuperoperator</a>(qureg, targets, numTargets, superOp);</div>
<div class="line"><a name="l00732"></a><span class="lineno"> 732</span>  <a class="code" href="group__type.html#ga9df2f3d86be4a6e9aad481665e5e6753">destroyComplexMatrixN</a>(superOp);</div>
<div class="line"><a name="l00733"></a><span class="lineno"> 733</span>  }</div>
<div class="line"><a name="l00734"></a><span class="lineno"> 734</span>  <span class="comment">// on Windows, we must always create in heap</span></div>
<div class="line"><a name="l00735"></a><span class="lineno"> 735</span> <span class="preprocessor"> #else</span></div>
<div class="line"><a name="l00736"></a><span class="lineno"> 736</span>  superOp = <a class="code" href="group__type.html#ga815103261fb22ea9690e1427571df00e">createComplexMatrixN</a>(2*numTargets);</div>
<div class="line"><a name="l00737"></a><span class="lineno"> 737</span>  <a class="code" href="QuEST__common_8c.html#a6dcb85ef8e1257813af0dc306082caf8">populateKrausSuperOperatorN</a>(&superOp, ops, numOps);</div>
<div class="line"><a name="l00738"></a><span class="lineno"> 738</span>  <a class="code" href="QuEST__common_8c.html#a514ca44595676abf4bc2e10fc48ef970">densmatr_applyMultiQubitKrausSuperoperator</a>(qureg, targets, numTargets, superOp);</div>
<div class="line"><a name="l00739"></a><span class="lineno"> 739</span>  <a class="code" href="group__type.html#ga9df2f3d86be4a6e9aad481665e5e6753">destroyComplexMatrixN</a>(superOp);</div>
<div class="line"><a name="l00740"></a><span class="lineno"> 740</span> <span class="preprocessor"> #endif</span></div>
<div class="line"><a name="l00741"></a><span class="lineno"> 741</span> }</div>
<div class="line"><a name="l00742"></a><span class="lineno"> 742</span>  </div>
<div class="line"><a name="l00743"></a><span class="lineno"><a class="line" href="QuEST__internal_8h.html#a74c14b1c46d838f5bf2bed1d5cf63271"> 743</a></span> <span class="keywordtype">void</span> <a class="code" href="QuEST__common_8c.html#a63bf6575f12d7113fb854a2b6cc1a756">densmatr_mixPauli</a>(<a class="code" href="structQureg.html">Qureg</a> qureg, <span class="keywordtype">int</span> qubit, <a class="code" href="group__type.html#ga7740e349b4f8bae6451547680f0ce2d6">qreal</a> probX, <a class="code" href="group__type.html#ga7740e349b4f8bae6451547680f0ce2d6">qreal</a> probY, <a class="code" href="group__type.html#ga7740e349b4f8bae6451547680f0ce2d6">qreal</a> probZ) {</div>
<div class="line"><a name="l00744"></a><span class="lineno"> 744</span>  </div>
<div class="line"><a name="l00745"></a><span class="lineno"> 745</span>  <span class="comment">// convert pauli probabilities into Kraus map</span></div>
<div class="line"><a name="l00746"></a><span class="lineno"> 746</span>  <span class="keyword">const</span> <span class="keywordtype">int</span> numOps = 4;</div>
<div class="line"><a name="l00747"></a><span class="lineno"> 747</span>  <a class="code" href="structComplexMatrix2.html">ComplexMatrix2</a> ops[4]; <span class="comment">// [numOps];</span></div>
<div class="line"><a name="l00748"></a><span class="lineno"> 748</span>  <span class="keywordflow">for</span> (<span class="keywordtype">int</span> n=0; n < numOps; n++)</div>
<div class="line"><a name="l00749"></a><span class="lineno"> 749</span>  ops[n] = (<a class="code" href="structComplexMatrix2.html">ComplexMatrix2</a>) {.<a class="code" href="structComplexMatrix2.html#a3ea98d545b67a190ee6abe936c72076f">real</a>={{0}}, .imag={{0}}};</div>
<div class="line"><a name="l00750"></a><span class="lineno"> 750</span>  </div>
<div class="line"><a name="l00751"></a><span class="lineno"> 751</span>  <a class="code" href="group__type.html#ga7740e349b4f8bae6451547680f0ce2d6">qreal</a> facs[4] = { <span class="comment">// literal numOps=4 for valid initialisation</span></div>
<div class="line"><a name="l00752"></a><span class="lineno"> 752</span>  sqrt(1-(probX + probY + probZ)),</div>
<div class="line"><a name="l00753"></a><span class="lineno"> 753</span>  sqrt(probX),</div>
<div class="line"><a name="l00754"></a><span class="lineno"> 754</span>  sqrt(probY),</div>
<div class="line"><a name="l00755"></a><span class="lineno"> 755</span>  sqrt(probZ)</div>
<div class="line"><a name="l00756"></a><span class="lineno"> 756</span>  };</div>
<div class="line"><a name="l00757"></a><span class="lineno"> 757</span>  ops[0].<a class="code" href="structComplexMatrix2.html#a3ea98d545b67a190ee6abe936c72076f">real</a>[0][0] = facs[0]; ops[0].<a class="code" href="structComplexMatrix2.html#a3ea98d545b67a190ee6abe936c72076f">real</a>[1][1] = facs[0];</div>
<div class="line"><a name="l00758"></a><span class="lineno"> 758</span>  ops[1].<a class="code" href="structComplexMatrix2.html#a3ea98d545b67a190ee6abe936c72076f">real</a>[0][1] = facs[1]; ops[1].<a class="code" href="structComplexMatrix2.html#a3ea98d545b67a190ee6abe936c72076f">real</a>[1][0] = facs[1];</div>
<div class="line"><a name="l00759"></a><span class="lineno"> 759</span>  ops[2].<a class="code" href="structComplexMatrix2.html#a3a2b60a9e00ccef054a6de3b9307c1be">imag</a>[0][1] = -facs[2]; ops[2].<a class="code" href="structComplexMatrix2.html#a3a2b60a9e00ccef054a6de3b9307c1be">imag</a>[1][0] = facs[2];</div>
<div class="line"><a name="l00760"></a><span class="lineno"> 760</span>  ops[3].<a class="code" href="structComplexMatrix2.html#a3ea98d545b67a190ee6abe936c72076f">real</a>[0][0] = facs[3]; ops[3].<a class="code" href="structComplexMatrix2.html#a3ea98d545b67a190ee6abe936c72076f">real</a>[1][1] = -facs[3];</div>
<div class="line"><a name="l00761"></a><span class="lineno"> 761</span>  </div>
<div class="line"><a name="l00762"></a><span class="lineno"> 762</span>  <a class="code" href="QuEST__common_8c.html#ab987904f721746d1176aba7ac6b566fa">densmatr_mixKrausMap</a>(qureg, qubit, ops, numOps);</div>
<div class="line"><a name="l00763"></a><span class="lineno"> 763</span> }</div>
<div class="line"><a name="l00764"></a><span class="lineno"> 764</span>  </div>
<div class="line"><a name="l00765"></a><span class="lineno"><a class="line" href="QuEST__common_8c.html#a9a28932c97151b4499450326f501164b"> 765</a></span> <span class="keywordtype">void</span> <a class="code" href="QuEST__common_8c.html#a9a28932c97151b4499450326f501164b">applyExponentiatedPauliHamil</a>(<a class="code" href="structQureg.html">Qureg</a> qureg, <a class="code" href="structPauliHamil.html">PauliHamil</a> hamil, <a class="code" href="group__type.html#ga7740e349b4f8bae6451547680f0ce2d6">qreal</a> fac, <span class="keywordtype">int</span> reverse) {</div>
<div class="line"><a name="l00766"></a><span class="lineno"> 766</span>  </div>
<div class="line"><a name="l00767"></a><span class="lineno"> 767</span>  <span class="comment">/* applies a first-order one-repetition approximation of exp(-i fac H)</span></div>
<div class="line"><a name="l00768"></a><span class="lineno"> 768</span> <span class="comment"> * to qureg. Letting H = sum_j c_j h_j, it does this via </span></div>
<div class="line"><a name="l00769"></a><span class="lineno"> 769</span> <span class="comment"> * exp(-i fac H) ~ prod_j exp(-i fac c_j h_j), where each inner exp </span></div>
<div class="line"><a name="l00770"></a><span class="lineno"> 770</span> <span class="comment"> * is performed with multiRotatePauli (with pre-factor 2).</span></div>
<div class="line"><a name="l00771"></a><span class="lineno"> 771</span> <span class="comment"> */</span></div>
<div class="line"><a name="l00772"></a><span class="lineno"> 772</span>  </div>
<div class="line"><a name="l00773"></a><span class="lineno"> 773</span>  <span class="comment">// prepare targets for multiRotatePauli </span></div>
<div class="line"><a name="l00774"></a><span class="lineno"> 774</span>  <span class="comment">// (all qubits; actual targets are determined by Pauli codes)</span></div>
<div class="line"><a name="l00775"></a><span class="lineno"> 775</span>  <span class="keywordtype">int</span> vecTargs[100]; <span class="comment">// [hamil.numQubits];</span></div>
<div class="line"><a name="l00776"></a><span class="lineno"> 776</span>  <span class="keywordtype">int</span> densTargs[100]; <span class="comment">// [hamil.numQubits];</span></div>
<div class="line"><a name="l00777"></a><span class="lineno"> 777</span>  <span class="keywordflow">for</span> (<span class="keywordtype">int</span> q=0; q<hamil.<a class="code" href="structPauliHamil.html#acb49fff0e3160ea1904c9ca34408b4e1">numQubits</a>; q++) {</div>
<div class="line"><a name="l00778"></a><span class="lineno"> 778</span>  vecTargs[q] = q;</div>
<div class="line"><a name="l00779"></a><span class="lineno"> 779</span>  densTargs[q] = q + hamil.<a class="code" href="structPauliHamil.html#acb49fff0e3160ea1904c9ca34408b4e1">numQubits</a>;</div>
<div class="line"><a name="l00780"></a><span class="lineno"> 780</span>  }</div>
<div class="line"><a name="l00781"></a><span class="lineno"> 781</span>  </div>
<div class="line"><a name="l00782"></a><span class="lineno"> 782</span>  <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i=0; i<hamil.<a class="code" href="structPauliHamil.html#a737e46b40ef8990cf004d22eb9dbe2e4">numSumTerms</a>; i++) {</div>
<div class="line"><a name="l00783"></a><span class="lineno"> 783</span>  </div>
<div class="line"><a name="l00784"></a><span class="lineno"> 784</span>  <span class="keywordtype">int</span> t=i;</div>
<div class="line"><a name="l00785"></a><span class="lineno"> 785</span>  <span class="keywordflow">if</span> (reverse)</div>
<div class="line"><a name="l00786"></a><span class="lineno"> 786</span>  t=hamil.<a class="code" href="structPauliHamil.html#a737e46b40ef8990cf004d22eb9dbe2e4">numSumTerms</a>-1-i;</div>
<div class="line"><a name="l00787"></a><span class="lineno"> 787</span>  </div>
<div class="line"><a name="l00788"></a><span class="lineno"> 788</span>  <a class="code" href="group__type.html#ga7740e349b4f8bae6451547680f0ce2d6">qreal</a> angle = 2*fac*hamil.<a class="code" href="structPauliHamil.html#ac1c6c9f6299cf83bbba88123216c3cac">termCoeffs</a>[t];</div>
<div class="line"><a name="l00789"></a><span class="lineno"> 789</span>  </div>
<div class="line"><a name="l00790"></a><span class="lineno"> 790</span>  <a class="code" href="QuEST__common_8c.html#a2153d8611247dd3a2ef483355265cc99">statevec_multiRotatePauli</a>(</div>
<div class="line"><a name="l00791"></a><span class="lineno"> 791</span>  qureg, vecTargs, &(hamil.<a class="code" href="structPauliHamil.html#a4b8d38c403553c07ceddcd40c46ce6c8">pauliCodes</a>[t*hamil.<a class="code" href="structPauliHamil.html#acb49fff0e3160ea1904c9ca34408b4e1">numQubits</a>]), </div>
<div class="line"><a name="l00792"></a><span class="lineno"> 792</span>  hamil.<a class="code" href="structPauliHamil.html#acb49fff0e3160ea1904c9ca34408b4e1">numQubits</a>, angle, 0);</div>
<div class="line"><a name="l00793"></a><span class="lineno"> 793</span>  </div>
<div class="line"><a name="l00794"></a><span class="lineno"> 794</span>  <span class="keywordflow">if</span> (qureg.<a class="code" href="structQureg.html#acf78445e9435d09f44f0cc832c6aee79">isDensityMatrix</a>)</div>
<div class="line"><a name="l00795"></a><span class="lineno"> 795</span>  <a class="code" href="QuEST__common_8c.html#a2153d8611247dd3a2ef483355265cc99">statevec_multiRotatePauli</a>(</div>
<div class="line"><a name="l00796"></a><span class="lineno"> 796</span>  qureg, densTargs, &(hamil.<a class="code" href="structPauliHamil.html#a4b8d38c403553c07ceddcd40c46ce6c8">pauliCodes</a>[t*hamil.<a class="code" href="structPauliHamil.html#acb49fff0e3160ea1904c9ca34408b4e1">numQubits</a>]), </div>
<div class="line"><a name="l00797"></a><span class="lineno"> 797</span>  hamil.<a class="code" href="structPauliHamil.html#acb49fff0e3160ea1904c9ca34408b4e1">numQubits</a>, angle, 1);</div>
<div class="line"><a name="l00798"></a><span class="lineno"> 798</span>  </div>
<div class="line"><a name="l00799"></a><span class="lineno"> 799</span>  <span class="comment">// record qasm </span></div>
<div class="line"><a name="l00800"></a><span class="lineno"> 800</span>  <span class="keywordtype">char</span> buff[1024];</div>
<div class="line"><a name="l00801"></a><span class="lineno"> 801</span>  <span class="keywordtype">int</span> b=0;</div>
<div class="line"><a name="l00802"></a><span class="lineno"> 802</span>  <span class="keywordflow">for</span> (<span class="keywordtype">int</span> q=0; q<hamil.<a class="code" href="structPauliHamil.html#acb49fff0e3160ea1904c9ca34408b4e1">numQubits</a>; q++) {</div>
<div class="line"><a name="l00803"></a><span class="lineno"> 803</span>  <span class="keyword">enum</span> <a class="code" href="group__type.html#ga1c703cf89629e4e9c7023cd402d67028">pauliOpType</a> op = hamil.<a class="code" href="structPauliHamil.html#a4b8d38c403553c07ceddcd40c46ce6c8">pauliCodes</a>[q + t*hamil.<a class="code" href="structPauliHamil.html#acb49fff0e3160ea1904c9ca34408b4e1">numQubits</a>];</div>
<div class="line"><a name="l00804"></a><span class="lineno"> 804</span>  </div>
<div class="line"><a name="l00805"></a><span class="lineno"> 805</span>  <span class="keywordtype">char</span> p = <span class="charliteral">'I'</span>;</div>
<div class="line"><a name="l00806"></a><span class="lineno"> 806</span>  <span class="keywordflow">if</span> (op == <a class="code" href="group__type.html#gga1c703cf89629e4e9c7023cd402d67028a7abac7bb5b71e17382014c443244ad5c">PAULI_X</a>) p = <span class="charliteral">'X'</span>;</div>
<div class="line"><a name="l00807"></a><span class="lineno"> 807</span>  <span class="keywordflow">if</span> (op == <a class="code" href="group__type.html#gga1c703cf89629e4e9c7023cd402d67028a2c64a9aaee27f6bcfef894d159e759bd">PAULI_Y</a>) p = <span class="charliteral">'Y'</span>;</div>
<div class="line"><a name="l00808"></a><span class="lineno"> 808</span>  <span class="keywordflow">if</span> (op == <a class="code" href="group__type.html#gga1c703cf89629e4e9c7023cd402d67028a085c85cff6ba79c52d8b1e61c42ddddf">PAULI_Z</a>) p = <span class="charliteral">'Z'</span>;</div>
<div class="line"><a name="l00809"></a><span class="lineno"> 809</span>  buff[b++] = p;</div>
<div class="line"><a name="l00810"></a><span class="lineno"> 810</span>  buff[b++] = <span class="charliteral">' '</span>;</div>
<div class="line"><a name="l00811"></a><span class="lineno"> 811</span>  }</div>
<div class="line"><a name="l00812"></a><span class="lineno"> 812</span>  buff[b] = <span class="charliteral">'\0'</span>;</div>
<div class="line"><a name="l00813"></a><span class="lineno"> 813</span>  </div>
<div class="line"><a name="l00814"></a><span class="lineno"> 814</span>  <a class="code" href="QuEST__qasm_8c.html#a9a385fde708958617bbeb599064b8dcb">qasm_recordComment</a>(qureg, </div>
<div class="line"><a name="l00815"></a><span class="lineno"> 815</span>  <span class="stringliteral">"Here, a multiRotatePauli with angle %g and paulis %s was applied."</span>,</div>
<div class="line"><a name="l00816"></a><span class="lineno"> 816</span>  angle, buff);</div>
<div class="line"><a name="l00817"></a><span class="lineno"> 817</span>  }</div>
<div class="line"><a name="l00818"></a><span class="lineno"> 818</span> }</div>
<div class="line"><a name="l00819"></a><span class="lineno"> 819</span>  </div>
<div class="line"><a name="l00820"></a><span class="lineno"><a class="line" href="QuEST__common_8c.html#a4f419c557e4505da481879407a920541"> 820</a></span> <span class="keywordtype">void</span> <a class="code" href="QuEST__common_8c.html#a4f419c557e4505da481879407a920541">applySymmetrizedTrotterCircuit</a>(<a class="code" href="structQureg.html">Qureg</a> qureg, <a class="code" href="structPauliHamil.html">PauliHamil</a> hamil, <a class="code" href="group__type.html#ga7740e349b4f8bae6451547680f0ce2d6">qreal</a> time, <span class="keywordtype">int</span> order) {</div>
<div class="line"><a name="l00821"></a><span class="lineno"> 821</span>  </div>
<div class="line"><a name="l00822"></a><span class="lineno"> 822</span>  <span class="keywordflow">if</span> (order == 1) {</div>
<div class="line"><a name="l00823"></a><span class="lineno"> 823</span>  <a class="code" href="QuEST__common_8c.html#a9a28932c97151b4499450326f501164b">applyExponentiatedPauliHamil</a>(qureg, hamil, time, 0);</div>
<div class="line"><a name="l00824"></a><span class="lineno"> 824</span>  }</div>
<div class="line"><a name="l00825"></a><span class="lineno"> 825</span>  <span class="keywordflow">else</span> <span class="keywordflow">if</span> (order == 2) {</div>
<div class="line"><a name="l00826"></a><span class="lineno"> 826</span>  <a class="code" href="QuEST__common_8c.html#a9a28932c97151b4499450326f501164b">applyExponentiatedPauliHamil</a>(qureg, hamil, time/2., 0);</div>
<div class="line"><a name="l00827"></a><span class="lineno"> 827</span>  <a class="code" href="QuEST__common_8c.html#a9a28932c97151b4499450326f501164b">applyExponentiatedPauliHamil</a>(qureg, hamil, time/2., 1);</div>
<div class="line"><a name="l00828"></a><span class="lineno"> 828</span>  }</div>
<div class="line"><a name="l00829"></a><span class="lineno"> 829</span>  <span class="keywordflow">else</span> {</div>
<div class="line"><a name="l00830"></a><span class="lineno"> 830</span>  <a class="code" href="group__type.html#ga7740e349b4f8bae6451547680f0ce2d6">qreal</a> p = 1. / (4 - pow(4, 1./(order-1)));</div>
<div class="line"><a name="l00831"></a><span class="lineno"> 831</span>  <span class="keywordtype">int</span> lower = order-2;</div>
<div class="line"><a name="l00832"></a><span class="lineno"> 832</span>  <a class="code" href="QuEST__common_8c.html#a4f419c557e4505da481879407a920541">applySymmetrizedTrotterCircuit</a>(qureg, hamil, p*time, lower);</div>
<div class="line"><a name="l00833"></a><span class="lineno"> 833</span>  <a class="code" href="QuEST__common_8c.html#a4f419c557e4505da481879407a920541">applySymmetrizedTrotterCircuit</a>(qureg, hamil, p*time, lower);</div>
<div class="line"><a name="l00834"></a><span class="lineno"> 834</span>  <a class="code" href="QuEST__common_8c.html#a4f419c557e4505da481879407a920541">applySymmetrizedTrotterCircuit</a>(qureg, hamil, (1-4*p)*time, lower);</div>
<div class="line"><a name="l00835"></a><span class="lineno"> 835</span>  <a class="code" href="QuEST__common_8c.html#a4f419c557e4505da481879407a920541">applySymmetrizedTrotterCircuit</a>(qureg, hamil, p*time, lower);</div>
<div class="line"><a name="l00836"></a><span class="lineno"> 836</span>  <a class="code" href="QuEST__common_8c.html#a4f419c557e4505da481879407a920541">applySymmetrizedTrotterCircuit</a>(qureg, hamil, p*time, lower);</div>
<div class="line"><a name="l00837"></a><span class="lineno"> 837</span>  }</div>
<div class="line"><a name="l00838"></a><span class="lineno"> 838</span> }</div>
<div class="line"><a name="l00839"></a><span class="lineno"> 839</span>  </div>
<div class="line"><a name="l00840"></a><span class="lineno"><a class="line" href="QuEST__internal_8h.html#a8a8e789585651f90fcc587042d4cdfb8"> 840</a></span> <span class="keywordtype">void</span> <a class="code" href="QuEST__common_8c.html#a8a8e789585651f90fcc587042d4cdfb8">agnostic_applyTrotterCircuit</a>(<a class="code" href="structQureg.html">Qureg</a> qureg, <a class="code" href="structPauliHamil.html">PauliHamil</a> hamil, <a class="code" href="group__type.html#ga7740e349b4f8bae6451547680f0ce2d6">qreal</a> time, <span class="keywordtype">int</span> order, <span class="keywordtype">int</span> reps) {</div>
<div class="line"><a name="l00841"></a><span class="lineno"> 841</span>  </div>
<div class="line"><a name="l00842"></a><span class="lineno"> 842</span>  <span class="keywordflow">if</span> (time == 0)</div>
<div class="line"><a name="l00843"></a><span class="lineno"> 843</span>  <span class="keywordflow">return</span>;</div>
<div class="line"><a name="l00844"></a><span class="lineno"> 844</span>  </div>
<div class="line"><a name="l00845"></a><span class="lineno"> 845</span>  <span class="keywordflow">for</span> (<span class="keywordtype">int</span> r=0; r<reps; r++)</div>
<div class="line"><a name="l00846"></a><span class="lineno"> 846</span>  <a class="code" href="QuEST__common_8c.html#a4f419c557e4505da481879407a920541">applySymmetrizedTrotterCircuit</a>(qureg, hamil, time/reps, order);</div>
<div class="line"><a name="l00847"></a><span class="lineno"> 847</span> }</div>
<div class="line"><a name="l00848"></a><span class="lineno"> 848</span>  </div>
<div class="line"><a name="l00849"></a><span class="lineno"><a class="line" href="QuEST__internal_8h.html#a2b78ab7dd04eab615a0d7c796e73ae8c"> 849</a></span> <span class="keywordtype">void</span> <a class="code" href="QuEST__common_8c.html#a2b78ab7dd04eab615a0d7c796e73ae8c">agnostic_applyQFT</a>(<a class="code" href="structQureg.html">Qureg</a> qureg, <span class="keywordtype">int</span>* qubits, <span class="keywordtype">int</span> numQubits) {</div>
<div class="line"><a name="l00850"></a><span class="lineno"> 850</span>  </div>
<div class="line"><a name="l00851"></a><span class="lineno"> 851</span>  <span class="keywordtype">int</span> densShift = qureg.<a class="code" href="structQureg.html#ad08dff5316b8937f4b2a1417591543dc">numQubitsRepresented</a>;</div>
<div class="line"><a name="l00852"></a><span class="lineno"> 852</span>  </div>
<div class="line"><a name="l00853"></a><span class="lineno"> 853</span>  <span class="comment">// start with top/left-most qubit, work down</span></div>
<div class="line"><a name="l00854"></a><span class="lineno"> 854</span>  <span class="keywordflow">for</span> (<span class="keywordtype">int</span> q=numQubits-1; q >= 0; q--) {</div>
<div class="line"><a name="l00855"></a><span class="lineno"> 855</span>  </div>
<div class="line"><a name="l00856"></a><span class="lineno"> 856</span>  <span class="comment">// H</span></div>
<div class="line"><a name="l00857"></a><span class="lineno"> 857</span>  <a class="code" href="QuEST__internal_8h.html#af3fc7bd149c70f083b149050242cd55e">statevec_hadamard</a>(qureg, qubits[q]);</div>
<div class="line"><a name="l00858"></a><span class="lineno"> 858</span>  <span class="keywordflow">if</span> (qureg.<a class="code" href="structQureg.html#acf78445e9435d09f44f0cc832c6aee79">isDensityMatrix</a>)</div>
<div class="line"><a name="l00859"></a><span class="lineno"> 859</span>  <a class="code" href="QuEST__internal_8h.html#af3fc7bd149c70f083b149050242cd55e">statevec_hadamard</a>(qureg, qubits[q] + densShift);</div>
<div class="line"><a name="l00860"></a><span class="lineno"> 860</span>  <a class="code" href="QuEST__qasm_8c.html#a81749d23b8c3c1acf87c9db675872c24">qasm_recordGate</a>(qureg, <a class="code" href="QuEST__qasm_8h.html#a32bbd0276b80668812ac248cd81c725fa8417e56d80add49216437e0d451e8dd9">GATE_HADAMARD</a>, qubits[q]);</div>
<div class="line"><a name="l00861"></a><span class="lineno"> 861</span>  </div>
<div class="line"><a name="l00862"></a><span class="lineno"> 862</span>  <span class="keywordflow">if</span> (q == 0)</div>
<div class="line"><a name="l00863"></a><span class="lineno"> 863</span>  <span class="keywordflow">break</span>;</div>
<div class="line"><a name="l00864"></a><span class="lineno"> 864</span>  </div>
<div class="line"><a name="l00865"></a><span class="lineno"> 865</span>  <span class="comment">// succession of C-phases, control on qubits[q], targeting each of </span></div>
<div class="line"><a name="l00866"></a><span class="lineno"> 866</span>  <span class="comment">// qubits[q-1], qubits[q-1], ... qubits[0]. This can be expressed by </span></div>
<div class="line"><a name="l00867"></a><span class="lineno"> 867</span>  <span class="comment">// a single invocation of applyNamedPhaseFunc product</span></div>
<div class="line"><a name="l00868"></a><span class="lineno"> 868</span>  </div>
<div class="line"><a name="l00869"></a><span class="lineno"> 869</span>  <span class="keywordtype">int</span> numRegs = 2;</div>
<div class="line"><a name="l00870"></a><span class="lineno"> 870</span>  <span class="keywordtype">int</span> numQubitsPerReg[2] = {q, 1};</div>
<div class="line"><a name="l00871"></a><span class="lineno"> 871</span>  <span class="keywordtype">int</span> regs[100]; <span class="comment">// [q+1];</span></div>
<div class="line"><a name="l00872"></a><span class="lineno"> 872</span>  <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i=0; i<q+1; i++)</div>
<div class="line"><a name="l00873"></a><span class="lineno"> 873</span>  regs[i] = qubits[i]; <span class="comment">// qubits[q] is in own register</span></div>
<div class="line"><a name="l00874"></a><span class="lineno"> 874</span>  </div>
<div class="line"><a name="l00875"></a><span class="lineno"> 875</span>  <span class="keywordtype">int</span> numParams = 1;</div>
<div class="line"><a name="l00876"></a><span class="lineno"> 876</span>  <a class="code" href="group__type.html#ga7740e349b4f8bae6451547680f0ce2d6">qreal</a> params[1] = { <a class="code" href="QuEST__common_8c.html#ae71449b1cc6e6250b91f539153a7a0d3">M_PI</a> / (1 << q) };</div>
<div class="line"><a name="l00877"></a><span class="lineno"> 877</span>  </div>
<div class="line"><a name="l00878"></a><span class="lineno"> 878</span>  <span class="keywordtype">int</span> conj = 0;</div>
<div class="line"><a name="l00879"></a><span class="lineno"> 879</span>  <a class="code" href="QuEST__internal_8h.html#a7d8a817207847d2646a52044a63fd469">statevec_applyParamNamedPhaseFuncOverrides</a>(</div>
<div class="line"><a name="l00880"></a><span class="lineno"> 880</span>  qureg, regs, numQubitsPerReg, numRegs, </div>
<div class="line"><a name="l00881"></a><span class="lineno"> 881</span>  <a class="code" href="group__type.html#gga9c91c41aa42ae537995e89ffc616afe1a7165f9a47792f47c718ca128556fb3ae">UNSIGNED</a>, <a class="code" href="group__type.html#ggaa7d869b117ba5024d6b84938e8cdfc65a31a29a30f92fe34c35098f941adbbc93">SCALED_PRODUCT</a>, params, numParams, </div>
<div class="line"><a name="l00882"></a><span class="lineno"> 882</span>  NULL, NULL, 0, </div>
<div class="line"><a name="l00883"></a><span class="lineno"> 883</span>  conj);</div>
<div class="line"><a name="l00884"></a><span class="lineno"> 884</span>  <span class="keywordflow">if</span> (qureg.<a class="code" href="structQureg.html#acf78445e9435d09f44f0cc832c6aee79">isDensityMatrix</a>) {</div>
<div class="line"><a name="l00885"></a><span class="lineno"> 885</span>  conj = 1;</div>
<div class="line"><a name="l00886"></a><span class="lineno"> 886</span>  <a class="code" href="QuEST__common_8c.html#aecc3c36b275c53b321de1611c2f06f78">shiftSubregIndices</a>(regs, numQubitsPerReg, numRegs, densShift);</div>
<div class="line"><a name="l00887"></a><span class="lineno"> 887</span>  <a class="code" href="QuEST__internal_8h.html#a7d8a817207847d2646a52044a63fd469">statevec_applyParamNamedPhaseFuncOverrides</a>(</div>
<div class="line"><a name="l00888"></a><span class="lineno"> 888</span>  qureg, regs, numQubitsPerReg, numRegs, </div>
<div class="line"><a name="l00889"></a><span class="lineno"> 889</span>  <a class="code" href="group__type.html#gga9c91c41aa42ae537995e89ffc616afe1a7165f9a47792f47c718ca128556fb3ae">UNSIGNED</a>, <a class="code" href="group__type.html#ggaa7d869b117ba5024d6b84938e8cdfc65a31a29a30f92fe34c35098f941adbbc93">SCALED_PRODUCT</a>, params, numParams, </div>
<div class="line"><a name="l00890"></a><span class="lineno"> 890</span>  NULL, NULL, 0, </div>
<div class="line"><a name="l00891"></a><span class="lineno"> 891</span>  conj);</div>
<div class="line"><a name="l00892"></a><span class="lineno"> 892</span>  <a class="code" href="QuEST__common_8c.html#aecc3c36b275c53b321de1611c2f06f78">shiftSubregIndices</a>(regs, numQubitsPerReg, numRegs, - densShift);</div>
<div class="line"><a name="l00893"></a><span class="lineno"> 893</span>  }</div>
<div class="line"><a name="l00894"></a><span class="lineno"> 894</span>  <a class="code" href="QuEST__qasm_8c.html#ab60a4a8af0530cbfa8cba1b02692d448">qasm_recordNamedPhaseFunc</a>(</div>
<div class="line"><a name="l00895"></a><span class="lineno"> 895</span>  qureg, regs, numQubitsPerReg, numRegs, </div>
<div class="line"><a name="l00896"></a><span class="lineno"> 896</span>  <a class="code" href="group__type.html#gga9c91c41aa42ae537995e89ffc616afe1a7165f9a47792f47c718ca128556fb3ae">UNSIGNED</a>, <a class="code" href="group__type.html#ggaa7d869b117ba5024d6b84938e8cdfc65a31a29a30f92fe34c35098f941adbbc93">SCALED_PRODUCT</a>, params, numParams, </div>
<div class="line"><a name="l00897"></a><span class="lineno"> 897</span>  NULL, NULL, 0);</div>
<div class="line"><a name="l00898"></a><span class="lineno"> 898</span>  }</div>
<div class="line"><a name="l00899"></a><span class="lineno"> 899</span>  </div>
<div class="line"><a name="l00900"></a><span class="lineno"> 900</span>  <span class="comment">// final swaps</span></div>
<div class="line"><a name="l00901"></a><span class="lineno"> 901</span>  <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i=0; i<(numQubits/2); i++) {</div>
<div class="line"><a name="l00902"></a><span class="lineno"> 902</span>  </div>
<div class="line"><a name="l00903"></a><span class="lineno"> 903</span>  <span class="keywordtype">int</span> qb1 = qubits[i];</div>
<div class="line"><a name="l00904"></a><span class="lineno"> 904</span>  <span class="keywordtype">int</span> qb2 = qubits[numQubits-i-1];</div>
<div class="line"><a name="l00905"></a><span class="lineno"> 905</span>  </div>
<div class="line"><a name="l00906"></a><span class="lineno"> 906</span>  <a class="code" href="QuEST__internal_8h.html#ab7478fcaf6357ab202d72853b672e375">statevec_swapQubitAmps</a>(qureg, qb1, qb2);</div>
<div class="line"><a name="l00907"></a><span class="lineno"> 907</span>  <span class="keywordflow">if</span> (qureg.<a class="code" href="structQureg.html#acf78445e9435d09f44f0cc832c6aee79">isDensityMatrix</a>)</div>
<div class="line"><a name="l00908"></a><span class="lineno"> 908</span>  <a class="code" href="QuEST__internal_8h.html#ab7478fcaf6357ab202d72853b672e375">statevec_swapQubitAmps</a>(qureg, qb1 + densShift, qb2 + densShift);</div>
<div class="line"><a name="l00909"></a><span class="lineno"> 909</span>  <a class="code" href="QuEST__qasm_8c.html#aa4c380db9eda18892d5a86825fa4b165">qasm_recordControlledGate</a>(qureg, <a class="code" href="QuEST__qasm_8h.html#a32bbd0276b80668812ac248cd81c725fab12d052de2b7c59c6ab87710fd7d3c78">GATE_SWAP</a>, qb1, qb2);</div>
<div class="line"><a name="l00910"></a><span class="lineno"> 910</span>  }</div>
<div class="line"><a name="l00911"></a><span class="lineno"> 911</span> }</div>
<div class="line"><a name="l00912"></a><span class="lineno"> 912</span>  </div>
<div class="line"><a name="l00913"></a><span class="lineno"> 913</span> <span class="preprocessor">#ifdef __cplusplus</span></div>
<div class="line"><a name="l00914"></a><span class="lineno"> 914</span> }</div>
<div class="line"><a name="l00915"></a><span class="lineno"> 915</span> <span class="preprocessor">#endif</span></div>
</div><!-- fragment --></div><!-- contents -->
<div class="ttc" id="aQuEST__common_8c_html_a5ef139196852f92fe8271c82964eadff"><div class="ttname"><a href="QuEST__common_8c.html#a5ef139196852f92fe8271c82964eadff">statevec_phaseShift</a></div><div class="ttdeci">void statevec_phaseShift(Qureg qureg, int targetQubit, qreal angle)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00254">QuEST_common.c:254</a></div></div>
<div class="ttc" id="astructVector_html"><div class="ttname"><a href="structVector.html">Vector</a></div><div class="ttdoc">Represents a 3-vector of real numbers.</div><div class="ttdef"><b>Definition:</b> <a href="QuEST_8h_source.html#l00198">QuEST.h:198</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_a559d175d31cec552fc09aec3e8a78d11"><div class="ttname"><a href="QuEST__common_8c.html#a559d175d31cec552fc09aec3e8a78d11">statevec_sGate</a></div><div class="ttdeci">void statevec_sGate(Qureg qureg, int targetQubit)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00268">QuEST_common.c:268</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_ab14df201e5c05527dc77cfeb6f583042"><div class="ttname"><a href="QuEST__common_8c.html#ab14df201e5c05527dc77cfeb6f583042">statevec_pauliZ</a></div><div class="ttdeci">void statevec_pauliZ(Qureg qureg, int targetQubit)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00261">QuEST_common.c:261</a></div></div>
<div class="ttc" id="agroup__type_html_ga1c703cf89629e4e9c7023cd402d67028"><div class="ttname"><a href="group__type.html#ga1c703cf89629e4e9c7023cd402d67028">pauliOpType</a></div><div class="ttdeci">pauliOpType</div><div class="ttdoc">Codes for specifying Pauli operators.</div><div class="ttdef"><b>Definition:</b> <a href="QuEST_8h_source.html#l00096">QuEST.h:96</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_ab8687e6ea5cf35503c223051837f3b94"><div class="ttname"><a href="QuEST__common_8c.html#ab8687e6ea5cf35503c223051837f3b94">statevec_rotateX</a></div><div class="ttdeci">void statevec_rotateX(Qureg qureg, int rotQubit, qreal angle)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00296">QuEST_common.c:296</a></div></div>
<div class="ttc" id="agroup__debug_html_ga97d96af7c7ea7b31e32cbe3b25377e09"><div class="ttname"><a href="group__debug.html#ga97d96af7c7ea7b31e32cbe3b25377e09">reportQuregParams</a></div><div class="ttdeci">void reportQuregParams(Qureg qureg)</div><div class="ttdoc">Report metainformation about a set of qubits: number of qubits, number of probability amplitudes.</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00237">QuEST_common.c:237</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_a2b78ab7dd04eab615a0d7c796e73ae8c"><div class="ttname"><a href="QuEST__common_8c.html#a2b78ab7dd04eab615a0d7c796e73ae8c">agnostic_applyQFT</a></div><div class="ttdeci">void agnostic_applyQFT(Qureg qureg, int *qubits, int numQubits)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00849">QuEST_common.c:849</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_ae2b222fda52c10fb420cb3ac78ca31e4"><div class="ttname"><a href="QuEST__common_8c.html#ae2b222fda52c10fb420cb3ac78ca31e4">macro_setConjugateMatrix</a></div><div class="ttdeci">#define macro_setConjugateMatrix(dest, src, dim)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00099">QuEST_common.c:99</a></div></div>
<div class="ttc" id="astructComplexMatrix4_html_aeadc9893e1b16cdd585dac523e73926c"><div class="ttname"><a href="structComplexMatrix4.html#aeadc9893e1b16cdd585dac523e73926c">ComplexMatrix4::real</a></div><div class="ttdeci">qreal real[4][4]</div><div class="ttdef"><b>Definition:</b> <a href="QuEST_8h_source.html#l00177">QuEST.h:177</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_ab987904f721746d1176aba7ac6b566fa"><div class="ttname"><a href="QuEST__common_8c.html#ab987904f721746d1176aba7ac6b566fa">densmatr_mixKrausMap</a></div><div class="ttdeci">void densmatr_mixKrausMap(Qureg qureg, int target, ComplexMatrix2 *ops, int numOps)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00644">QuEST_common.c:644</a></div></div>
<div class="ttc" id="agroup__type_html_gga1c703cf89629e4e9c7023cd402d67028a085c85cff6ba79c52d8b1e61c42ddddf"><div class="ttname"><a href="group__type.html#gga1c703cf89629e4e9c7023cd402d67028a085c85cff6ba79c52d8b1e61c42ddddf">PAULI_Z</a></div><div class="ttdeci">@ PAULI_Z</div><div class="ttdef"><b>Definition:</b> <a href="QuEST_8h_source.html#l00096">QuEST.h:96</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_a9a28932c97151b4499450326f501164b"><div class="ttname"><a href="QuEST__common_8c.html#a9a28932c97151b4499450326f501164b">applyExponentiatedPauliHamil</a></div><div class="ttdeci">void applyExponentiatedPauliHamil(Qureg qureg, PauliHamil hamil, qreal fac, int reverse)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00765">QuEST_common.c:765</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_a2a6999e47bed978eb66f589c8abc7766"><div class="ttname"><a href="QuEST__common_8c.html#a2a6999e47bed978eb66f589c8abc7766">populateKrausSuperOperator4</a></div><div class="ttdeci">void populateKrausSuperOperator4(ComplexMatrixN *superOp, ComplexMatrix4 *ops, int numOps)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00611">QuEST_common.c:611</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_a6dcb85ef8e1257813af0dc306082caf8"><div class="ttname"><a href="QuEST__common_8c.html#a6dcb85ef8e1257813af0dc306082caf8">populateKrausSuperOperatorN</a></div><div class="ttdeci">void populateKrausSuperOperatorN(ComplexMatrixN *superOp, ComplexMatrixN *ops, int numOps)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00615">QuEST_common.c:615</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_a08d2e22870f3ef0289e355d634ea17a2"><div class="ttname"><a href="QuEST__common_8c.html#a08d2e22870f3ef0289e355d634ea17a2">statevec_multiControlledMultiRotatePauli</a></div><div class="ttdeci">void statevec_multiControlledMultiRotatePauli(Qureg qureg, long long int ctrlMask, int *targetQubits, enum pauliOpType *targetPaulis, int numTargets, qreal angle, int applyConj)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00453">QuEST_common.c:453</a></div></div>
<div class="ttc" id="agroup__type_html_ga9df2f3d86be4a6e9aad481665e5e6753"><div class="ttname"><a href="group__type.html#ga9df2f3d86be4a6e9aad481665e5e6753">destroyComplexMatrixN</a></div><div class="ttdeci">void destroyComplexMatrixN(ComplexMatrixN m)</div><div class="ttdoc">Destroy a ComplexMatrixN instance created with createComplexMatrixN()</div><div class="ttdef"><b>Definition:</b> <a href="QuEST_8c_source.html#l01369">QuEST.c:1369</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_acf1c651a4359a370a2d9b3e1e1ebd430"><div class="ttname"><a href="QuEST__common_8c.html#acf1c651a4359a370a2d9b3e1e1ebd430">shiftIndices</a></div><div class="ttdeci">void shiftIndices(int *indices, int numIndices, int shift)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00156">QuEST_common.c:156</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_aef380928d1b0003411d6159ca6594cd5"><div class="ttname"><a href="QuEST__common_8c.html#aef380928d1b0003411d6159ca6594cd5">statevec_calcExpecPauliProd</a></div><div class="ttdeci">qreal statevec_calcExpecPauliProd(Qureg qureg, int *targetQubits, enum pauliOpType *pauliCodes, int numTargets, Qureg workspace)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00509">QuEST_common.c:509</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_ad7ac4934371d72e24c6d1cea3d3da391"><div class="ttname"><a href="QuEST__common_8c.html#ad7ac4934371d72e24c6d1cea3d3da391">statevec_controlledTwoQubitUnitary</a></div><div class="ttdeci">void statevec_controlledTwoQubitUnitary(Qureg qureg, int controlQubit, int targetQubit1, int targetQubit2, ComplexMatrix4 u)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00567">QuEST_common.c:567</a></div></div>
<div class="ttc" id="agroup__type_html_gga1c703cf89629e4e9c7023cd402d67028a63800f45b01ccc2608f754aee850bf4a"><div class="ttname"><a href="group__type.html#gga1c703cf89629e4e9c7023cd402d67028a63800f45b01ccc2608f754aee850bf4a">PAULI_I</a></div><div class="ttdeci">@ PAULI_I</div><div class="ttdef"><b>Definition:</b> <a href="QuEST_8h_source.html#l00096">QuEST.h:96</a></div></div>
<div class="ttc" id="aQuEST__internal_8h_html"><div class="ttname"><a href="QuEST__internal_8h.html">QuEST_internal.h</a></div></div>
<div class="ttc" id="agroup__type_html_ga815103261fb22ea9690e1427571df00e"><div class="ttname"><a href="group__type.html#ga815103261fb22ea9690e1427571df00e">createComplexMatrixN</a></div><div class="ttdeci">ComplexMatrixN createComplexMatrixN(int numQubits)</div><div class="ttdoc">Allocate dynamic memory for a square complex matrix of any size, which can be passed to functions lik...</div><div class="ttdef"><b>Definition:</b> <a href="QuEST_8c_source.html#l01348">QuEST.c:1348</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_ae0cfea24091c8336146aa84c2253992e"><div class="ttname"><a href="QuEST__common_8c.html#ae0cfea24091c8336146aa84c2253992e">statevec_calcExpecPauliSum</a></div><div class="ttdeci">qreal statevec_calcExpecPauliSum(Qureg qureg, enum pauliOpType *allCodes, qreal *termCoeffs, int numSumTerms, Qureg workspace)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00524">QuEST_common.c:524</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_a95abbe0cedf45b1ff2f070dba219fd7f"><div class="ttname"><a href="QuEST__common_8c.html#a95abbe0cedf45b1ff2f070dba219fd7f">statevec_twoQubitUnitary</a></div><div class="ttdeci">void statevec_twoQubitUnitary(Qureg qureg, int targetQubit1, int targetQubit2, ComplexMatrix4 u)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00561">QuEST_common.c:561</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_a086fffa88b1aeec1405aa2dc41eab8a3"><div class="ttname"><a href="QuEST__common_8c.html#a086fffa88b1aeec1405aa2dc41eab8a3">statevec_tGateConj</a></div><div class="ttdeci">void statevec_tGateConj(Qureg qureg, int targetQubit)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00289">QuEST_common.c:289</a></div></div>
<div class="ttc" id="astructVector_html_ab33976d8fda2e44b770f71269d698a25"><div class="ttname"><a href="structVector.html#ab33976d8fda2e44b770f71269d698a25">Vector::z</a></div><div class="ttdeci">qreal z</div><div class="ttdef"><b>Definition:</b> <a href="QuEST_8h_source.html#l00200">QuEST.h:200</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_a20d45e85dd49dd3bf6045990b590b24b"><div class="ttname"><a href="QuEST__common_8c.html#a20d45e85dd49dd3bf6045990b590b24b">getConjugateMatrix4</a></div><div class="ttdeci">ComplexMatrix4 getConjugateMatrix4(ComplexMatrix4 src)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00110">QuEST_common.c:110</a></div></div>
<div class="ttc" id="aQuEST_8h_html"><div class="ttname"><a href="QuEST_8h.html">QuEST.h</a></div></div>
<div class="ttc" id="astructQureg_html_aa3deb7163aab90a61e755a1f200f5413"><div class="ttname"><a href="structQureg.html#aa3deb7163aab90a61e755a1f200f5413">Qureg::numChunks</a></div><div class="ttdeci">int numChunks</div><div class="ttdoc">Number of chunks the state vector is broken up into – the number of MPI processes used.</div><div class="ttdef"><b>Definition:</b> <a href="QuEST_8h_source.html#l00338">QuEST.h:338</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_accae504dfcd6846159baa16a1e8263a0"><div class="ttname"><a href="QuEST__common_8c.html#accae504dfcd6846159baa16a1e8263a0">statevec_applyPauliProd</a></div><div class="ttdeci">void statevec_applyPauliProd(Qureg workspace, int *targetQubits, enum pauliOpType *pauliCodes, int numTargets)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00495">QuEST_common.c:495</a></div></div>
<div class="ttc" id="aQuEST__internal_8h_html_a0acfe30083d0bc2a84da196378a36122"><div class="ttname"><a href="QuEST__internal_8h.html#a0acfe30083d0bc2a84da196378a36122">statevec_multiControlledMultiQubitUnitary</a></div><div class="ttdeci">void statevec_multiControlledMultiQubitUnitary(Qureg qureg, long long int ctrlMask, int *targs, int numTargs, ComplexMatrixN u)</div><div class="ttdoc">This calls swapQubitAmps only when it would involve a distributed communication; if the qubit chunks ...</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__cpu__distributed_8c_source.html#l01514">QuEST_cpu_distributed.c:1514</a></div></div>
<div class="ttc" id="aQuEST__internal_8h_html_a8d88b4ed94ab024e5e7b2f1f4f505b44"><div class="ttname"><a href="QuEST__internal_8h.html#a8d88b4ed94ab024e5e7b2f1f4f505b44">statevec_multiRotateZ</a></div><div class="ttdeci">void statevec_multiRotateZ(Qureg qureg, long long int mask, qreal angle)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__cpu_8c_source.html#l03316">QuEST_cpu.c:3316</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_a2b323b65cd6512720e6aa27c182fcea0"><div class="ttname"><a href="QuEST__common_8c.html#a2b323b65cd6512720e6aa27c182fcea0">getQuESTDefaultSeedKey</a></div><div class="ttdeci">void getQuESTDefaultSeedKey(unsigned long int *key)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00195">QuEST_common.c:195</a></div></div>
<div class="ttc" id="aQuEST__internal_8h_html_a63cb88a5f4eb12fef087435b8db3f8eb"><div class="ttname"><a href="QuEST__internal_8h.html#a63cb88a5f4eb12fef087435b8db3f8eb">statevec_calcProbOfOutcome</a></div><div class="ttdeci">qreal statevec_calcProbOfOutcome(Qureg qureg, int measureQubit, int outcome)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__cpu__distributed_8c_source.html#l01312">QuEST_cpu_distributed.c:1312</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_a99702f4667c429969d7e7b4ecd45c2f3"><div class="ttname"><a href="QuEST__common_8c.html#a99702f4667c429969d7e7b4ecd45c2f3">densmatr_mixMultiQubitKrausMap</a></div><div class="ttdeci">void densmatr_mixMultiQubitKrausMap(Qureg qureg, int *targets, int numTargets, ComplexMatrixN *ops, int numOps)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00701">QuEST_common.c:701</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_af92fa737f4491324408b3d15543dced4"><div class="ttname"><a href="QuEST__common_8c.html#af92fa737f4491324408b3d15543dced4">getConjugateScalar</a></div><div class="ttdeci">Complex getConjugateScalar(Complex scalar)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00091">QuEST_common.c:91</a></div></div>
<div class="ttc" id="aQuEST__qasm_8h_html_a32bbd0276b80668812ac248cd81c725fa8417e56d80add49216437e0d451e8dd9"><div class="ttname"><a href="QuEST__qasm_8h.html#a32bbd0276b80668812ac248cd81c725fa8417e56d80add49216437e0d451e8dd9">GATE_HADAMARD</a></div><div class="ttdeci">@ GATE_HADAMARD</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__qasm_8h_source.html#l00026">QuEST_qasm.h:26</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_ad34d98eeaa6a3cb45e92c696cfd2d1c8"><div class="ttname"><a href="QuEST__common_8c.html#ad34d98eeaa6a3cb45e92c696cfd2d1c8">getConjugateMatrix2</a></div><div class="ttdeci">ComplexMatrix2 getConjugateMatrix2(ComplexMatrix2 src)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00105">QuEST_common.c:105</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_a041bb75200ee316de7b90667a4fe04f3"><div class="ttname"><a href="QuEST__common_8c.html#a041bb75200ee316de7b90667a4fe04f3">getUnitVector</a></div><div class="ttdeci">Vector getUnitVector(Vector vec)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00084">QuEST_common.c:84</a></div></div>
<div class="ttc" id="astructComplexMatrix4_html"><div class="ttname"><a href="structComplexMatrix4.html">ComplexMatrix4</a></div><div class="ttdoc">Represents a 4x4 matrix of complex numbers.</div><div class="ttdef"><b>Definition:</b> <a href="QuEST_8h_source.html#l00175">QuEST.h:175</a></div></div>
<div class="ttc" id="agroup__type_html_gga9c91c41aa42ae537995e89ffc616afe1a7165f9a47792f47c718ca128556fb3ae"><div class="ttname"><a href="group__type.html#gga9c91c41aa42ae537995e89ffc616afe1a7165f9a47792f47c718ca128556fb3ae">UNSIGNED</a></div><div class="ttdeci">@ UNSIGNED</div><div class="ttdef"><b>Definition:</b> <a href="QuEST_8h_source.html#l00269">QuEST.h:269</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_a71b3ae4402c1c5039d3086ca0e47a522"><div class="ttname"><a href="QuEST__common_8c.html#a71b3ae4402c1c5039d3086ca0e47a522">getComplexPairFromRotation</a></div><div class="ttdeci">void getComplexPairFromRotation(qreal angle, Vector axis, Complex *alpha, Complex *beta)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00120">QuEST_common.c:120</a></div></div>
<div class="ttc" id="aQuEST__internal_8h_html_af1cca787f3b0e0c91ffb2495f9ad8d5a"><div class="ttname"><a href="QuEST__internal_8h.html#af1cca787f3b0e0c91ffb2495f9ad8d5a">statevec_pauliY</a></div><div class="ttdeci">void statevec_pauliY(Qureg qureg, int targetQubit)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__cpu__distributed_8c_source.html#l01142">QuEST_cpu_distributed.c:1142</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_ad67cb67877c941189b0de2c2e67e2026"><div class="ttname"><a href="QuEST__common_8c.html#ad67cb67877c941189b0de2c2e67e2026">bindArraysToStackComplexMatrixN</a></div><div class="ttdeci">ComplexMatrixN bindArraysToStackComplexMatrixN(int numQubits, qreal re[][1<< numQubits], qreal im[][1<< numQubits], qreal **reStorage, qreal **imStorage)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00652">QuEST_common.c:652</a></div></div>
<div class="ttc" id="astructComplexMatrixN_html"><div class="ttname"><a href="structComplexMatrixN.html">ComplexMatrixN</a></div><div class="ttdoc">Represents a general 2^N by 2^N matrix of complex numbers.</div><div class="ttdef"><b>Definition:</b> <a href="QuEST_8h_source.html#l00186">QuEST.h:186</a></div></div>
<div class="ttc" id="amt19937ar_8h_html"><div class="ttname"><a href="mt19937ar_8h.html">mt19937ar.h</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_a06b8474aec51579ff680fb66439e7fc6"><div class="ttname"><a href="QuEST__common_8c.html#a06b8474aec51579ff680fb66439e7fc6">statevec_controlledMultiQubitUnitary</a></div><div class="ttdeci">void statevec_controlledMultiQubitUnitary(Qureg qureg, int ctrl, int *targets, int numTargets, ComplexMatrixN u)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00579">QuEST_common.c:579</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_aeeec9fe1501e25930adf6c38763a158e"><div class="ttname"><a href="QuEST__common_8c.html#aeeec9fe1501e25930adf6c38763a158e">statevec_calcFidelity</a></div><div class="ttdeci">qreal statevec_calcFidelity(Qureg qureg, Qureg pureState)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00380">QuEST_common.c:380</a></div></div>
<div class="ttc" id="agroup__type_html_ga7740e349b4f8bae6451547680f0ce2d6"><div class="ttname"><a href="group__type.html#ga7740e349b4f8bae6451547680f0ce2d6">qreal</a></div><div class="ttdeci">#define qreal</div></div>
<div class="ttc" id="aQuEST__common_8c_html_a4da103b16292ffd99197d05c5fa83712"><div class="ttname"><a href="QuEST__common_8c.html#a4da103b16292ffd99197d05c5fa83712">macro_allocStackComplexMatrixN</a></div><div class="ttdeci">#define macro_allocStackComplexMatrixN(matrix, numQubits)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00675">QuEST_common.c:675</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_a6550462400942be61b171c8840b36aa1"><div class="ttname"><a href="QuEST__common_8c.html#a6550462400942be61b171c8840b36aa1">statevec_sqrtSwapGate</a></div><div class="ttdeci">void statevec_sqrtSwapGate(Qureg qureg, int qb1, int qb2)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00387">QuEST_common.c:387</a></div></div>
<div class="ttc" id="agroup__type_html_gga1c703cf89629e4e9c7023cd402d67028a7abac7bb5b71e17382014c443244ad5c"><div class="ttname"><a href="group__type.html#gga1c703cf89629e4e9c7023cd402d67028a7abac7bb5b71e17382014c443244ad5c">PAULI_X</a></div><div class="ttdeci">@ PAULI_X</div><div class="ttdef"><b>Definition:</b> <a href="QuEST_8h_source.html#l00096">QuEST.h:96</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_a3afe58e49dce5536bc4c1b7caddf44a6"><div class="ttname"><a href="QuEST__common_8c.html#a3afe58e49dce5536bc4c1b7caddf44a6">statevec_getProbAmp</a></div><div class="ttdeci">qreal statevec_getProbAmp(Qureg qureg, long long int index)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00248">QuEST_common.c:248</a></div></div>
<div class="ttc" id="agroup__debug_html_gacd5f0ee488157837a8483c1515817f98"><div class="ttname"><a href="group__debug.html#gacd5f0ee488157837a8483c1515817f98">reportState</a></div><div class="ttdeci">void reportState(Qureg qureg)</div><div class="ttdoc">Print the current state vector of probability amplitudes for a set of qubits to file.</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00219">QuEST_common.c:219</a></div></div>
<div class="ttc" id="aQuEST__internal_8h_html_ab7478fcaf6357ab202d72853b672e375"><div class="ttname"><a href="QuEST__internal_8h.html#ab7478fcaf6357ab202d72853b672e375">statevec_swapQubitAmps</a></div><div class="ttdeci">void statevec_swapQubitAmps(Qureg qureg, int qb1, int qb2)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__cpu__distributed_8c_source.html#l01431">QuEST_cpu_distributed.c:1431</a></div></div>
<div class="ttc" id="astructQureg_html_a0a86e8a50ea8c998b4a0a6640c5a1218"><div class="ttname"><a href="structQureg.html#a0a86e8a50ea8c998b4a0a6640c5a1218">Qureg::numQubitsInStateVec</a></div><div class="ttdeci">int numQubitsInStateVec</div><div class="ttdoc">Number of qubits in the state-vector - this is double the number represented for mixed states.</div><div class="ttdef"><b>Definition:</b> <a href="QuEST_8h_source.html#l00329">QuEST.h:329</a></div></div>
<div class="ttc" id="aQuEST__internal_8h_html_a6153547f245c05874161a105e9a2f02c"><div class="ttname"><a href="QuEST__internal_8h.html#a6153547f245c05874161a105e9a2f02c">densmatr_calcTotalProb</a></div><div class="ttdeci">qreal densmatr_calcTotalProb(Qureg qureg)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__cpu__distributed_8c_source.html#l00053">QuEST_cpu_distributed.c:53</a></div></div>
<div class="ttc" id="aQuEST__internal_8h_html_a5658f2ecbbd4ea425db5f7dc7fc4ba92"><div class="ttname"><a href="QuEST__internal_8h.html#a5658f2ecbbd4ea425db5f7dc7fc4ba92">densmatr_collapseToKnownProbOutcome</a></div><div class="ttdeci">void densmatr_collapseToKnownProbOutcome(Qureg qureg, int measureQubit, int outcome, qreal outcomeProb)</div><div class="ttdoc">Renorms (/prob) every | * outcome * >< * outcome * | state, setting all others to zero.</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__cpu_8c_source.html#l00791">QuEST_cpu.c:791</a></div></div>
<div class="ttc" id="astructQureg_html_ac2929e681d3d95591c18cf168dbbe4f0"><div class="ttname"><a href="structQureg.html#ac2929e681d3d95591c18cf168dbbe4f0">Qureg::chunkId</a></div><div class="ttdeci">int chunkId</div><div class="ttdoc">The position of the chunk of the state vector held by this process in the full state vector.</div><div class="ttdef"><b>Definition:</b> <a href="QuEST_8h_source.html#l00336">QuEST.h:336</a></div></div>
<div class="ttc" id="aQuEST__qasm_8c_html_aa4c380db9eda18892d5a86825fa4b165"><div class="ttname"><a href="QuEST__qasm_8c.html#aa4c380db9eda18892d5a86825fa4b165">qasm_recordControlledGate</a></div><div class="ttdeci">void qasm_recordControlledGate(Qureg qureg, TargetGate gate, int controlQubit, int targetQubit)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__qasm_8c_source.html#l00239">QuEST_qasm.c:239</a></div></div>
<div class="ttc" id="astructVector_html_abc3cfc6c8e297b1153a373d84b981e6b"><div class="ttname"><a href="structVector.html#abc3cfc6c8e297b1153a373d84b981e6b">Vector::y</a></div><div class="ttdeci">qreal y</div><div class="ttdef"><b>Definition:</b> <a href="QuEST_8h_source.html#l00200">QuEST.h:200</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_ab76254cfde16f0808476649507a1a2fc"><div class="ttname"><a href="QuEST__common_8c.html#ab76254cfde16f0808476649507a1a2fc">hashString</a></div><div class="ttdeci">unsigned long int hashString(char *str)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00185">QuEST_common.c:185</a></div></div>
<div class="ttc" id="astructComplexMatrix2_html_a3a2b60a9e00ccef054a6de3b9307c1be"><div class="ttname"><a href="structComplexMatrix2.html#a3a2b60a9e00ccef054a6de3b9307c1be">ComplexMatrix2::imag</a></div><div class="ttdeci">qreal imag[2][2]</div><div class="ttdef"><b>Definition:</b> <a href="QuEST_8h_source.html#l00140">QuEST.h:140</a></div></div>
<div class="ttc" id="aQuEST__internal_8h_html_a7d8a817207847d2646a52044a63fd469"><div class="ttname"><a href="QuEST__internal_8h.html#a7d8a817207847d2646a52044a63fd469">statevec_applyParamNamedPhaseFuncOverrides</a></div><div class="ttdeci">void statevec_applyParamNamedPhaseFuncOverrides(Qureg qureg, int *qubits, int *numQubitsPerReg, int numRegs, enum bitEncoding encoding, enum phaseFunc functionNameCode, qreal *params, int numParams, long long int *overrideInds, qreal *overridePhases, int numOverrides, int conj)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__cpu_8c_source.html#l04446">QuEST_cpu.c:4446</a></div></div>
<div class="ttc" id="aQuEST__internal_8h_html_a3ed30f823a59890be17441d97ce90c5c"><div class="ttname"><a href="QuEST__internal_8h.html#a3ed30f823a59890be17441d97ce90c5c">statevec_multiControlledMultiRotateZ</a></div><div class="ttdeci">void statevec_multiControlledMultiRotateZ(Qureg qureg, long long int ctrlMask, long long int targMask, qreal angle)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__cpu_8c_source.html#l03358">QuEST_cpu.c:3358</a></div></div>
<div class="ttc" id="astructVector_html_a0d0459a5b731ac5fb3dff2c6134e6e88"><div class="ttname"><a href="structVector.html#a0d0459a5b731ac5fb3dff2c6134e6e88">Vector::x</a></div><div class="ttdeci">qreal x</div><div class="ttdef"><b>Definition:</b> <a href="QuEST_8h_source.html#l00200">QuEST.h:200</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_a9d19fc690afb425d0df5feb09be567a3"><div class="ttname"><a href="QuEST__common_8c.html#a9d19fc690afb425d0df5feb09be567a3">generateMeasurementOutcome</a></div><div class="ttdeci">int generateMeasurementOutcome(qreal zeroProb, qreal *outcomeProb)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00168">QuEST_common.c:168</a></div></div>
<div class="ttc" id="astructQureg_html_ab0ea0358482b62f43fdd781469607d97"><div class="ttname"><a href="structQureg.html#ab0ea0358482b62f43fdd781469607d97">Qureg::numAmpsPerChunk</a></div><div class="ttdeci">long long int numAmpsPerChunk</div><div class="ttdoc">Number of probability amplitudes held in stateVec by this process In the non-MPI version,...</div><div class="ttdef"><b>Definition:</b> <a href="QuEST_8h_source.html#l00332">QuEST.h:332</a></div></div>
<div class="ttc" id="astructPauliHamil_html_ac1c6c9f6299cf83bbba88123216c3cac"><div class="ttname"><a href="structPauliHamil.html#ac1c6c9f6299cf83bbba88123216c3cac">PauliHamil::termCoeffs</a></div><div class="ttdeci">qreal * termCoeffs</div><div class="ttdoc">The real coefficient of each Pauli product. This is an array of length PauliHamil....</div><div class="ttdef"><b>Definition:</b> <a href="QuEST_8h_source.html#l00283">QuEST.h:283</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_a8a8e789585651f90fcc587042d4cdfb8"><div class="ttname"><a href="QuEST__common_8c.html#a8a8e789585651f90fcc587042d4cdfb8">agnostic_applyTrotterCircuit</a></div><div class="ttdeci">void agnostic_applyTrotterCircuit(Qureg qureg, PauliHamil hamil, qreal time, int order, int reps)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00840">QuEST_common.c:840</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_af3db21114efc431ae7a7ae06fc2d7b5b"><div class="ttname"><a href="QuEST__common_8c.html#af3db21114efc431ae7a7ae06fc2d7b5b">densmatr_applyKrausSuperoperator</a></div><div class="ttdeci">void densmatr_applyKrausSuperoperator(Qureg qureg, int target, ComplexMatrix4 superOp)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00620">QuEST_common.c:620</a></div></div>
<div class="ttc" id="astructPauliHamil_html_a4b8d38c403553c07ceddcd40c46ce6c8"><div class="ttname"><a href="structPauliHamil.html#a4b8d38c403553c07ceddcd40c46ce6c8">PauliHamil::pauliCodes</a></div><div class="ttdeci">enum pauliOpType * pauliCodes</div><div class="ttdoc">The Pauli operators acting on each qubit, flattened over every operator.</div><div class="ttdef"><b>Definition:</b> <a href="QuEST_8h_source.html#l00281">QuEST.h:281</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_a2d96b4601cc3a9e1221e11a366fe1699"><div class="ttname"><a href="QuEST__common_8c.html#a2d96b4601cc3a9e1221e11a366fe1699">statevec_tGate</a></div><div class="ttdeci">void statevec_tGate(Qureg qureg, int targetQubit)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00275">QuEST_common.c:275</a></div></div>
<div class="ttc" id="aQuEST__internal_8h_html_a758bad4237ff0bf3b4ff5be626a982ae"><div class="ttname"><a href="QuEST__internal_8h.html#a758bad4237ff0bf3b4ff5be626a982ae">statevec_initBlankState</a></div><div class="ttdeci">void statevec_initBlankState(Qureg qureg)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__cpu_8c_source.html#l01464">QuEST_cpu.c:1464</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_abd00ba1321c9bd43c3a9119f5b3ed1be"><div class="ttname"><a href="QuEST__common_8c.html#abd00ba1321c9bd43c3a9119f5b3ed1be">populateKrausSuperOperator2</a></div><div class="ttdeci">void populateKrausSuperOperator2(ComplexMatrix4 *superOp, ComplexMatrix2 *ops, int numOps)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00607">QuEST_common.c:607</a></div></div>
<div class="ttc" id="agroup__type_html_ggaa7d869b117ba5024d6b84938e8cdfc65a31a29a30f92fe34c35098f941adbbc93"><div class="ttname"><a href="group__type.html#ggaa7d869b117ba5024d6b84938e8cdfc65a31a29a30f92fe34c35098f941adbbc93">SCALED_PRODUCT</a></div><div class="ttdeci">@ SCALED_PRODUCT</div><div class="ttdef"><b>Definition:</b> <a href="QuEST_8h_source.html#l00233">QuEST.h:233</a></div></div>
<div class="ttc" id="aQuEST__internal_8h_html_aa4b77089145184d990848cd66993aa8b"><div class="ttname"><a href="QuEST__internal_8h.html#aa4b77089145184d990848cd66993aa8b">statevec_cloneQureg</a></div><div class="ttdeci">void statevec_cloneQureg(Qureg targetQureg, Qureg copyQureg)</div><div class="ttdoc">works for both statevectors and density matrices</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__cpu_8c_source.html#l01572">QuEST_cpu.c:1572</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_adf5fbea75826cde7d86fd7c31bcf3dec"><div class="ttname"><a href="QuEST__common_8c.html#adf5fbea75826cde7d86fd7c31bcf3dec">statevec_rotateAroundAxis</a></div><div class="ttdeci">void statevec_rotateAroundAxis(Qureg qureg, int rotQubit, qreal angle, Vector axis)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00314">QuEST_common.c:314</a></div></div>
<div class="ttc" id="astructComplexMatrix4_html_a1f6bb840b5a9e6252980cab7a14f1a35"><div class="ttname"><a href="structComplexMatrix4.html#a1f6bb840b5a9e6252980cab7a14f1a35">ComplexMatrix4::imag</a></div><div class="ttdeci">qreal imag[4][4]</div><div class="ttdef"><b>Definition:</b> <a href="QuEST_8h_source.html#l00178">QuEST.h:178</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_a9fd917537204e074119d4f283dea78e3"><div class="ttname"><a href="QuEST__common_8c.html#a9fd917537204e074119d4f283dea78e3">setConjugateMatrixN</a></div><div class="ttdeci">void setConjugateMatrixN(ComplexMatrixN m)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00115">QuEST_common.c:115</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_ab004038b5f906ea897621c88bf213129"><div class="ttname"><a href="QuEST__common_8c.html#ab004038b5f906ea897621c88bf213129">densmatr_mixTwoQubitKrausMap</a></div><div class="ttdeci">void densmatr_mixTwoQubitKrausMap(Qureg qureg, int target1, int target2, ComplexMatrix4 *ops, int numOps)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00682">QuEST_common.c:682</a></div></div>
<div class="ttc" id="aQuEST__internal_8h_html_a3eadcb581f2623e21e4967f791d0439c"><div class="ttname"><a href="QuEST__internal_8h.html#a3eadcb581f2623e21e4967f791d0439c">statevec_multiControlledTwoQubitUnitary</a></div><div class="ttdeci">void statevec_multiControlledTwoQubitUnitary(Qureg qureg, long long int ctrlMask, int targetQubit1, int targetQubit2, ComplexMatrix4 u)</div><div class="ttdoc">This calls swapQubitAmps only when it would involve a distributed communication; if the qubit chunks ...</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__cpu__distributed_8c_source.html#l01458">QuEST_cpu_distributed.c:1458</a></div></div>
<div class="ttc" id="astructPauliHamil_html_a737e46b40ef8990cf004d22eb9dbe2e4"><div class="ttname"><a href="structPauliHamil.html#a737e46b40ef8990cf004d22eb9dbe2e4">PauliHamil::numSumTerms</a></div><div class="ttdeci">int numSumTerms</div><div class="ttdoc">The number of terms in the weighted sum, or the number of Pauli products.</div><div class="ttdef"><b>Definition:</b> <a href="QuEST_8h_source.html#l00285">QuEST.h:285</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_acc45bb8300922ea5933e294921e25ea2"><div class="ttname"><a href="QuEST__common_8c.html#acc45bb8300922ea5933e294921e25ea2">getQubitBitMask</a></div><div class="ttdeci">long long int getQubitBitMask(int *qubits, int numQubits)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00050">QuEST_common.c:50</a></div></div>
<div class="ttc" id="aQuEST__internal_8h_html_aec63559e8b9ab0e17efe18d64d2ceca6"><div class="ttname"><a href="QuEST__internal_8h.html#aec63559e8b9ab0e17efe18d64d2ceca6">statevec_compactUnitary</a></div><div class="ttdeci">void statevec_compactUnitary(Qureg qureg, int targetQubit, Complex alpha, Complex beta)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__cpu__distributed_8c_source.html#l00858">QuEST_cpu_distributed.c:858</a></div></div>
<div class="ttc" id="agroup__type_html_gga1c703cf89629e4e9c7023cd402d67028a2c64a9aaee27f6bcfef894d159e759bd"><div class="ttname"><a href="group__type.html#gga1c703cf89629e4e9c7023cd402d67028a2c64a9aaee27f6bcfef894d159e759bd">PAULI_Y</a></div><div class="ttdeci">@ PAULI_Y</div><div class="ttdef"><b>Definition:</b> <a href="QuEST_8h_source.html#l00096">QuEST.h:96</a></div></div>
<div class="ttc" id="astructPauliHamil_html"><div class="ttname"><a href="structPauliHamil.html">PauliHamil</a></div><div class="ttdoc">A Pauli Hamiltonian, expressed as a real-weighted sum of pauli products, and which can hence represen...</div><div class="ttdef"><b>Definition:</b> <a href="QuEST_8h_source.html#l00277">QuEST.h:277</a></div></div>
<div class="ttc" id="amt19937ar_8c_html_ac94ab75771800274ed1a2bedeca86f04"><div class="ttname"><a href="mt19937ar_8c.html#ac94ab75771800274ed1a2bedeca86f04">genrand_real1</a></div><div class="ttdeci">double genrand_real1(void)</div><div class="ttdef"><b>Definition:</b> <a href="mt19937ar_8c_source.html#l00150">mt19937ar.c:150</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_a2153d8611247dd3a2ef483355265cc99"><div class="ttname"><a href="QuEST__common_8c.html#a2153d8611247dd3a2ef483355265cc99">statevec_multiRotatePauli</a></div><div class="ttdeci">void statevec_multiRotatePauli(Qureg qureg, int *targetQubits, enum pauliOpType *targetPaulis, int numTargets, qreal angle, int applyConj)</div><div class="ttdoc">applyConj=1 will apply conjugate operation, else applyConj=0</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00414">QuEST_common.c:414</a></div></div>
<div class="ttc" id="aQuEST__qasm_8c_html_a9a385fde708958617bbeb599064b8dcb"><div class="ttname"><a href="QuEST__qasm_8c.html#a9a385fde708958617bbeb599064b8dcb">qasm_recordComment</a></div><div class="ttdeci">void qasm_recordComment(Qureg qureg, char *comment,...)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__qasm_8c_source.html#l00121">QuEST_qasm.c:121</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_a7a76950824f17e821415070a4bfdad5b"><div class="ttname"><a href="QuEST__common_8c.html#a7a76950824f17e821415070a4bfdad5b">statevec_multiQubitUnitary</a></div><div class="ttdeci">void statevec_multiQubitUnitary(Qureg qureg, int *targets, int numTargets, ComplexMatrixN u)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00573">QuEST_common.c:573</a></div></div>
<div class="ttc" id="astructComplexMatrixN_html_a446e8a077597e4f76049d547df0c9738"><div class="ttname"><a href="structComplexMatrixN.html#a446e8a077597e4f76049d547df0c9738">ComplexMatrixN::real</a></div><div class="ttdeci">qreal ** real</div><div class="ttdef"><b>Definition:</b> <a href="QuEST_8h_source.html#l00189">QuEST.h:189</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_a336c3286d2434fc892ebe20e9fa22a24"><div class="ttname"><a href="QuEST__common_8c.html#a336c3286d2434fc892ebe20e9fa22a24">statevec_controlledRotateAroundAxis</a></div><div class="ttdeci">void statevec_controlledRotateAroundAxis(Qureg qureg, int controlQubit, int targetQubit, qreal angle, Vector axis)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00330">QuEST_common.c:330</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_af71c2470694795a8d639be8ff9c34ed5"><div class="ttname"><a href="QuEST__common_8c.html#af71c2470694795a8d639be8ff9c34ed5">statevec_applyPauliSum</a></div><div class="ttdeci">void statevec_applyPauliSum(Qureg inQureg, enum pauliOpType *allCodes, qreal *termCoeffs, int numSumTerms, Qureg outQureg)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00538">QuEST_common.c:538</a></div></div>
<div class="ttc" id="aQuEST__internal_8h_html_a7ebd3198a198f4cd20840f64fd8b84d0"><div class="ttname"><a href="QuEST__internal_8h.html#a7ebd3198a198f4cd20840f64fd8b84d0">statevec_calcInnerProduct</a></div><div class="ttdeci">Complex statevec_calcInnerProduct(Qureg bra, Qureg ket)</div><div class="ttdoc">Terrible code which unnecessarily individually computes and sums the real and imaginary components of...</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__cpu__distributed_8c_source.html#l00035">QuEST_cpu_distributed.c:35</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_a2cbfe9f416d4d8acb56213841f62cb53"><div class="ttname"><a href="QuEST__common_8c.html#a2cbfe9f416d4d8acb56213841f62cb53">statevec_rotateAroundAxisConj</a></div><div class="ttdeci">void statevec_rotateAroundAxisConj(Qureg qureg, int rotQubit, qreal angle, Vector axis)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00321">QuEST_common.c:321</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_a142b424d903dfaca7a7ac166076f543a"><div class="ttname"><a href="QuEST__common_8c.html#a142b424d903dfaca7a7ac166076f543a">statevec_sGateConj</a></div><div class="ttdeci">void statevec_sGateConj(Qureg qureg, int targetQubit)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00282">QuEST_common.c:282</a></div></div>
<div class="ttc" id="aQuEST__internal_8h_html_ad64a23933361fafc1b6b1787759c2663"><div class="ttname"><a href="QuEST__internal_8h.html#ad64a23933361fafc1b6b1787759c2663">statevec_setWeightedQureg</a></div><div class="ttdeci">void statevec_setWeightedQureg(Complex fac1, Qureg qureg1, Complex fac2, Qureg qureg2, Complex facOut, Qureg out)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__cpu_8c_source.html#l04005">QuEST_cpu.c:4005</a></div></div>
<div class="ttc" id="aQuEST__internal_8h_html_a008becc4a18f868340836e0ab9fd6df6"><div class="ttname"><a href="QuEST__internal_8h.html#a008becc4a18f868340836e0ab9fd6df6">statevec_collapseToKnownProbOutcome</a></div><div class="ttdeci">void statevec_collapseToKnownProbOutcome(Qureg qureg, int measureQubit, int outcome, qreal outcomeProb)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__cpu__distributed_8c_source.html#l01368">QuEST_cpu_distributed.c:1368</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_aecc3c36b275c53b321de1611c2f06f78"><div class="ttname"><a href="QuEST__common_8c.html#aecc3c36b275c53b321de1611c2f06f78">shiftSubregIndices</a></div><div class="ttdeci">void shiftSubregIndices(int *allInds, int *numIndsPerReg, int numRegs, int shift)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00161">QuEST_common.c:161</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_a6fa6f1bc617e9fef3cd87ecc3677a02b"><div class="ttname"><a href="QuEST__common_8c.html#a6fa6f1bc617e9fef3cd87ecc3677a02b">getControlFlipMask</a></div><div class="ttdeci">long long int getControlFlipMask(int *controlQubits, int *controlState, int numControlQubits)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00060">QuEST_common.c:60</a></div></div>
<div class="ttc" id="aQuEST__internal_8h_html_abd509244d57657e148e4084c5ab5d28f"><div class="ttname"><a href="QuEST__internal_8h.html#abd509244d57657e148e4084c5ab5d28f">statevec_getImagAmp</a></div><div class="ttdeci">qreal statevec_getImagAmp(Qureg qureg, long long int index)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__cpu__distributed_8c_source.html#l00222">QuEST_cpu_distributed.c:222</a></div></div>
<div class="ttc" id="astructQureg_html"><div class="ttname"><a href="structQureg.html">Qureg</a></div><div class="ttdoc">Represents a system of qubits.</div><div class="ttdef"><b>Definition:</b> <a href="QuEST_8h_source.html#l00322">QuEST.h:322</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_ac2eb81c71d0e0fbec7be5b45b1db1448"><div class="ttname"><a href="QuEST__common_8c.html#ac2eb81c71d0e0fbec7be5b45b1db1448">statevec_measureWithStats</a></div><div class="ttdeci">int statevec_measureWithStats(Qureg qureg, int measureQubit, qreal *outcomeProb)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00364">QuEST_common.c:364</a></div></div>
<div class="ttc" id="astructComplexMatrixN_html_afde13e469a3b4f179f45e53d1ec0a7cf"><div class="ttname"><a href="structComplexMatrixN.html#afde13e469a3b4f179f45e53d1ec0a7cf">ComplexMatrixN::imag</a></div><div class="ttdeci">qreal ** imag</div><div class="ttdef"><b>Definition:</b> <a href="QuEST_8h_source.html#l00190">QuEST.h:190</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_a404d177b894ce008289e44a599ad77d6"><div class="ttname"><a href="QuEST__common_8c.html#a404d177b894ce008289e44a599ad77d6">getZYZRotAnglesFromComplexPair</a></div><div class="ttdeci">void getZYZRotAnglesFromComplexPair(Complex alpha, Complex beta, qreal *rz2, qreal *ry, qreal *rz1)</div><div class="ttdoc">maps U(alpha, beta) to Rz(rz2) Ry(ry) Rz(rz1)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00130">QuEST_common.c:130</a></div></div>
<div class="ttc" id="aQuEST__precision_8h_html"><div class="ttname"><a href="QuEST__precision_8h.html">QuEST_precision.h</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_aa32eed5bc38e85193e1aa23b74e8541b"><div class="ttname"><a href="QuEST__common_8c.html#aa32eed5bc38e85193e1aa23b74e8541b">statevec_controlledRotateX</a></div><div class="ttdeci">void statevec_controlledRotateX(Qureg qureg, int controlQubit, int targetQubit, qreal angle)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00346">QuEST_common.c:346</a></div></div>
<div class="ttc" id="astructQureg_html_a441e4cacef6bd17adb9813c7442d42fe"><div class="ttname"><a href="structQureg.html#a441e4cacef6bd17adb9813c7442d42fe">Qureg::stateVec</a></div><div class="ttdeci">ComplexArray stateVec</div><div class="ttdoc">Computational state amplitudes - a subset thereof in the MPI version.</div><div class="ttdef"><b>Definition:</b> <a href="QuEST_8h_source.html#l00341">QuEST.h:341</a></div></div>
<div class="ttc" id="astructComplexMatrix2_html_a3ea98d545b67a190ee6abe936c72076f"><div class="ttname"><a href="structComplexMatrix2.html#a3ea98d545b67a190ee6abe936c72076f">ComplexMatrix2::real</a></div><div class="ttdeci">qreal real[2][2]</div><div class="ttdef"><b>Definition:</b> <a href="QuEST_8h_source.html#l00139">QuEST.h:139</a></div></div>
<div class="ttc" id="astructQureg_html_acf78445e9435d09f44f0cc832c6aee79"><div class="ttname"><a href="structQureg.html#acf78445e9435d09f44f0cc832c6aee79">Qureg::isDensityMatrix</a></div><div class="ttdeci">int isDensityMatrix</div><div class="ttdoc">Whether this instance is a density-state representation.</div><div class="ttdef"><b>Definition:</b> <a href="QuEST_8h_source.html#l00325">QuEST.h:325</a></div></div>
<div class="ttc" id="aQuEST__internal_8h_html_af3fc7bd149c70f083b149050242cd55e"><div class="ttname"><a href="QuEST__internal_8h.html#af3fc7bd149c70f083b149050242cd55e">statevec_hadamard</a></div><div class="ttdeci">void statevec_hadamard(Qureg qureg, int targetQubit)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__cpu__distributed_8c_source.html#l01258">QuEST_cpu_distributed.c:1258</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_a1f0db2800d13b6c2b120e05129cd13cb"><div class="ttname"><a href="QuEST__common_8c.html#a1f0db2800d13b6c2b120e05129cd13cb">statevec_rotateY</a></div><div class="ttdeci">void statevec_rotateY(Qureg qureg, int rotQubit, qreal angle)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00302">QuEST_common.c:302</a></div></div>
<div class="ttc" id="astructComplexMatrixN_html_a606f7afbc3f65a596bdca21cb43be6a8"><div class="ttname"><a href="structComplexMatrixN.html#a606f7afbc3f65a596bdca21cb43be6a8">ComplexMatrixN::numQubits</a></div><div class="ttdeci">int numQubits</div><div class="ttdef"><b>Definition:</b> <a href="QuEST_8h_source.html#l00188">QuEST.h:188</a></div></div>
<div class="ttc" id="aQuEST__internal_8h_html_a600bda10463e2078ce746ec28f2d39c9"><div class="ttname"><a href="QuEST__internal_8h.html#a600bda10463e2078ce746ec28f2d39c9">statevec_multiControlledUnitary</a></div><div class="ttdeci">void statevec_multiControlledUnitary(Qureg qureg, long long int ctrlQubitsMask, long long int ctrlFlipMask, int targetQubit, ComplexMatrix2 u)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__cpu__distributed_8c_source.html#l01011">QuEST_cpu_distributed.c:1011</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_a3cf54148e2162b04aeb634badb90826c"><div class="ttname"><a href="QuEST__common_8c.html#a3cf54148e2162b04aeb634badb90826c">getVectorMagnitude</a></div><div class="ttdeci">qreal getVectorMagnitude(Vector vec)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00079">QuEST_common.c:79</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_a3d39314af5d88822d4aefb5e9cc4e2f6"><div class="ttname"><a href="QuEST__common_8c.html#a3d39314af5d88822d4aefb5e9cc4e2f6">densmatr_applyTwoQubitKrausSuperoperator</a></div><div class="ttdeci">void densmatr_applyTwoQubitKrausSuperoperator(Qureg qureg, int target1, int target2, ComplexMatrixN superOp)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00626">QuEST_common.c:626</a></div></div>
<div class="ttc" id="aQuEST__internal_8h_html_a535cc3c9962d511078691e4d64b4e272"><div class="ttname"><a href="QuEST__internal_8h.html#a535cc3c9962d511078691e4d64b4e272">statevec_phaseShiftByTerm</a></div><div class="ttdeci">void statevec_phaseShiftByTerm(Qureg qureg, int targetQubit, Complex term)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__cpu_8c_source.html#l03185">QuEST_cpu.c:3185</a></div></div>
<div class="ttc" id="astructPauliHamil_html_acb49fff0e3160ea1904c9ca34408b4e1"><div class="ttname"><a href="structPauliHamil.html#acb49fff0e3160ea1904c9ca34408b4e1">PauliHamil::numQubits</a></div><div class="ttdeci">int numQubits</div><div class="ttdoc">The number of qubits informing the Hilbert dimension of the Hamiltonian.</div><div class="ttdef"><b>Definition:</b> <a href="QuEST_8h_source.html#l00287">QuEST.h:287</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_a33cde9b6518d55a544eff20301a5d844"><div class="ttname"><a href="QuEST__common_8c.html#a33cde9b6518d55a544eff20301a5d844">statevec_sqrtSwapGateConj</a></div><div class="ttdeci">void statevec_sqrtSwapGateConj(Qureg qureg, int qb1, int qb2)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00400">QuEST_common.c:400</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_a3e11a9b9fa5ee16d4738a97f4590ee61"><div class="ttname"><a href="QuEST__common_8c.html#a3e11a9b9fa5ee16d4738a97f4590ee61">statevec_controlledRotateZ</a></div><div class="ttdeci">void statevec_controlledRotateZ(Qureg qureg, int controlQubit, int targetQubit, qreal angle)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00358">QuEST_common.c:358</a></div></div>
<div class="ttc" id="astructQureg_html_ad08dff5316b8937f4b2a1417591543dc"><div class="ttname"><a href="structQureg.html#ad08dff5316b8937f4b2a1417591543dc">Qureg::numQubitsRepresented</a></div><div class="ttdeci">int numQubitsRepresented</div><div class="ttdoc">The number of qubits represented in either the state-vector or density matrix.</div><div class="ttdef"><b>Definition:</b> <a href="QuEST_8h_source.html#l00327">QuEST.h:327</a></div></div>
<div class="ttc" id="aQuEST__validation_8h_html"><div class="ttname"><a href="QuEST__validation_8h.html">QuEST_validation.h</a></div></div>
<div class="ttc" id="aQuEST__internal_8h_html_a1c82b5e914a8a39c5c7b782dbc16fa76"><div class="ttname"><a href="QuEST__internal_8h.html#a1c82b5e914a8a39c5c7b782dbc16fa76">statevec_pauliX</a></div><div class="ttdeci">void statevec_pauliX(Qureg qureg, int targetQubit)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__cpu__distributed_8c_source.html#l01048">QuEST_cpu_distributed.c:1048</a></div></div>
<div class="ttc" id="aQuEST__qasm_8h_html_a32bbd0276b80668812ac248cd81c725fab12d052de2b7c59c6ab87710fd7d3c78"><div class="ttname"><a href="QuEST__qasm_8h.html#a32bbd0276b80668812ac248cd81c725fab12d052de2b7c59c6ab87710fd7d3c78">GATE_SWAP</a></div><div class="ttdeci">@ GATE_SWAP</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__qasm_8h_source.html#l00033">QuEST_qasm.h:33</a></div></div>
<div class="ttc" id="astructComplex_html_ab5b2e2eca02131fc74b289a83636cbe3"><div class="ttname"><a href="structComplex.html#ab5b2e2eca02131fc74b289a83636cbe3">Complex::real</a></div><div class="ttdeci">qreal real</div><div class="ttdef"><b>Definition:</b> <a href="QuEST_8h_source.html#l00105">QuEST.h:105</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_a514ca44595676abf4bc2e10fc48ef970"><div class="ttname"><a href="QuEST__common_8c.html#a514ca44595676abf4bc2e10fc48ef970">densmatr_applyMultiQubitKrausSuperoperator</a></div><div class="ttdeci">void densmatr_applyMultiQubitKrausSuperoperator(Qureg qureg, int *targets, int numTargets, ComplexMatrixN superOp)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00634">QuEST_common.c:634</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_a63bf6575f12d7113fb854a2b6cc1a756"><div class="ttname"><a href="QuEST__common_8c.html#a63bf6575f12d7113fb854a2b6cc1a756">densmatr_mixPauli</a></div><div class="ttdeci">void densmatr_mixPauli(Qureg qureg, int qubit, qreal probX, qreal probY, qreal probZ)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00743">QuEST_common.c:743</a></div></div>
<div class="ttc" id="astructComplex_html_a84f5439aad0ef495efdd3b4c1c02d27e"><div class="ttname"><a href="structComplex.html#a84f5439aad0ef495efdd3b4c1c02d27e">Complex::imag</a></div><div class="ttdeci">qreal imag</div><div class="ttdef"><b>Definition:</b> <a href="QuEST_8h_source.html#l00106">QuEST.h:106</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_ac8ea84dd85eca525daeeb30b3c96b673"><div class="ttname"><a href="QuEST__common_8c.html#ac8ea84dd85eca525daeeb30b3c96b673">ensureIndsIncrease</a></div><div class="ttdeci">void ensureIndsIncrease(int *ind1, int *ind2)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00070">QuEST_common.c:70</a></div></div>
<div class="ttc" id="astructComplex_html"><div class="ttname"><a href="structComplex.html">Complex</a></div><div class="ttdoc">Represents one complex number.</div><div class="ttdef"><b>Definition:</b> <a href="QuEST_8h_source.html#l00103">QuEST.h:103</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_ade67874c89179541f616b3685595d3ae"><div class="ttname"><a href="QuEST__common_8c.html#ade67874c89179541f616b3685595d3ae">statevec_controlledRotateY</a></div><div class="ttdeci">void statevec_controlledRotateY(Qureg qureg, int controlQubit, int targetQubit, qreal angle)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00352">QuEST_common.c:352</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_a5ea2384a9b33c94eac891aacbfa5da6a"><div class="ttname"><a href="QuEST__common_8c.html#a5ea2384a9b33c94eac891aacbfa5da6a">statevec_controlledRotateAroundAxisConj</a></div><div class="ttdeci">void statevec_controlledRotateAroundAxisConj(Qureg qureg, int controlQubit, int targetQubit, qreal angle, Vector axis)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00337">QuEST_common.c:337</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_ab9babafb4a8cc2af97e3d815eadec9db"><div class="ttname"><a href="QuEST__common_8c.html#ab9babafb4a8cc2af97e3d815eadec9db">macro_populateKrausOperator</a></div><div class="ttdeci">#define macro_populateKrausOperator(superOp, ops, numOps, opDim)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00585">QuEST_common.c:585</a></div></div>
<div class="ttc" id="aQuEST__qasm_8c_html_a81749d23b8c3c1acf87c9db675872c24"><div class="ttname"><a href="QuEST__qasm_8c.html#a81749d23b8c3c1acf87c9db675872c24">qasm_recordGate</a></div><div class="ttdeci">void qasm_recordGate(Qureg qureg, TargetGate gate, int targetQubit)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__qasm_8c_source.html#l00179">QuEST_qasm.c:179</a></div></div>
<div class="ttc" id="aQuEST__internal_8h_html_ae37b3eff5d3affbc1edfe216cabc1c53"><div class="ttname"><a href="QuEST__internal_8h.html#ae37b3eff5d3affbc1edfe216cabc1c53">densmatr_calcProbOfOutcome</a></div><div class="ttdeci">qreal densmatr_calcProbOfOutcome(Qureg qureg, int measureQubit, int outcome)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__cpu__distributed_8c_source.html#l01328">QuEST_cpu_distributed.c:1328</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_a2b2b8c69566a49cdf4a3fa2cb4f5e03c"><div class="ttname"><a href="QuEST__common_8c.html#a2b2b8c69566a49cdf4a3fa2cb4f5e03c">statevec_rotateZ</a></div><div class="ttdeci">void statevec_rotateZ(Qureg qureg, int rotQubit, qreal angle)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00308">QuEST_common.c:308</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_ae71449b1cc6e6250b91f539153a7a0d3"><div class="ttname"><a href="QuEST__common_8c.html#ae71449b1cc6e6250b91f539153a7a0d3">M_PI</a></div><div class="ttdeci">#define M_PI</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00041">QuEST_common.c:41</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_a4f419c557e4505da481879407a920541"><div class="ttname"><a href="QuEST__common_8c.html#a4f419c557e4505da481879407a920541">applySymmetrizedTrotterCircuit</a></div><div class="ttdeci">void applySymmetrizedTrotterCircuit(Qureg qureg, PauliHamil hamil, qreal time, int order)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00820">QuEST_common.c:820</a></div></div>
<div class="ttc" id="aQuEST__qasm_8c_html_ab60a4a8af0530cbfa8cba1b02692d448"><div class="ttname"><a href="QuEST__qasm_8c.html#ab60a4a8af0530cbfa8cba1b02692d448">qasm_recordNamedPhaseFunc</a></div><div class="ttdeci">void qasm_recordNamedPhaseFunc(Qureg qureg, int *qubits, int *numQubitsPerReg, int numRegs, enum bitEncoding encoding, enum phaseFunc funcName, qreal *params, int numParams, long long int *overrideInds, qreal *overridePhases, int numOverrides)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__qasm_8c_source.html#l00726">QuEST_qasm.c:726</a></div></div>
<div class="ttc" id="aQuEST__internal_8h_html_a1275d604674224d87a173fb5bac78835"><div class="ttname"><a href="QuEST__internal_8h.html#a1275d604674224d87a173fb5bac78835">statevec_controlledCompactUnitary</a></div><div class="ttdeci">void statevec_controlledCompactUnitary(Qureg qureg, int controlQubit, int targetQubit, Complex alpha, Complex beta)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__cpu__distributed_8c_source.html#l00934">QuEST_cpu_distributed.c:934</a></div></div>
<div class="ttc" id="aQuEST__internal_8h_html_abc9a9ef4344c7faaaf28ac25c76649b9"><div class="ttname"><a href="QuEST__internal_8h.html#abc9a9ef4344c7faaaf28ac25c76649b9">statevec_getRealAmp</a></div><div class="ttdeci">qreal statevec_getRealAmp(Qureg qureg, long long int index)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__cpu__distributed_8c_source.html#l00212">QuEST_cpu_distributed.c:212</a></div></div>
<div class="ttc" id="aQuEST__qasm_8h_html"><div class="ttname"><a href="QuEST__qasm_8h.html">QuEST_qasm.h</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_ae4ebc1f4dd48922bfb30ad1692310269"><div class="ttname"><a href="QuEST__common_8c.html#ae4ebc1f4dd48922bfb30ad1692310269">densmatr_measureWithStats</a></div><div class="ttdeci">int densmatr_measureWithStats(Qureg qureg, int measureQubit, qreal *outcomeProb)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00372">QuEST_common.c:372</a></div></div>
<div class="ttc" id="astructComplexMatrix2_html"><div class="ttname"><a href="structComplexMatrix2.html">ComplexMatrix2</a></div><div class="ttdoc">Represents a 2x2 matrix of complex numbers.</div><div class="ttdef"><b>Definition:</b> <a href="QuEST_8h_source.html#l00137">QuEST.h:137</a></div></div>
<div class="ttc" id="aQuEST__common_8c_html_a73300e3cdca56a7615980274905ba097"><div class="ttname"><a href="QuEST__common_8c.html#a73300e3cdca56a7615980274905ba097">getComplexPairAndPhaseFromUnitary</a></div><div class="ttdeci">void getComplexPairAndPhaseFromUnitary(ComplexMatrix2 u, Complex *alpha, Complex *beta, qreal *globalPhase)</div><div class="ttdoc">maps U(r0c0, r0c1, r1c0, r1c1) to exp(i globalPhase) U(alpha, beta)</div><div class="ttdef"><b>Definition:</b> <a href="QuEST__common_8c_source.html#l00142">QuEST_common.c:142</a></div></div>
<!-- HTML footer for doxygen 1.8.8-->
<!-- start footer part -->
</div>
</div>
</div>
</div>
</div>
</body>
<script type="text/javascript" src="boot.js"></script>
</html>
