697 lines
68 KiB
HTML
697 lines
68 KiB
HTML
<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
|
|
|
|
<html lang="en">
|
|
|
|
<head>
|
|
<meta http-equiv="Content-Type" content="text/html; charset=UTF-8">
|
|
<title>LCOV - coverage.info - libtorch/include/c10/util/complex.h</title>
|
|
<link rel="stylesheet" type="text/css" href="../../../../gcov.css">
|
|
</head>
|
|
|
|
<body>
|
|
|
|
<table width="100%" border=0 cellspacing=0 cellpadding=0>
|
|
<tr><td class="title">LCOV - code coverage report</td></tr>
|
|
<tr><td class="ruler"><img src="../../../../glass.png" width=3 height=3 alt=""></td></tr>
|
|
|
|
<tr>
|
|
<td width="100%">
|
|
<table cellpadding=1 border=0 width="100%">
|
|
<tr>
|
|
<td width="10%" class="headerItem">Current view:</td>
|
|
<td width="10%" class="headerValue"><a href="../../../../index.html">top level</a> - <a href="index.html">libtorch/include/c10/util</a> - complex.h<span style="font-size: 80%;"> (source / <a href="complex.h.func-c.html">functions</a>)</span></td>
|
|
<td width="5%"></td>
|
|
<td width="5%"></td>
|
|
<td width="5%" class="headerCovTableHead">Coverage</td>
|
|
<td width="5%" class="headerCovTableHead" title="Covered + Uncovered code">Total</td>
|
|
<td width="5%" class="headerCovTableHead" title="Exercised code only">Hit</td>
|
|
</tr>
|
|
<tr>
|
|
<td class="headerItem">Test:</td>
|
|
<td class="headerValue">coverage.info</td>
|
|
<td></td>
|
|
<td class="headerItem">Lines:</td>
|
|
<td class="headerCovTableEntryLo">0.0 %</td>
|
|
<td class="headerCovTableEntry">4</td>
|
|
<td class="headerCovTableEntry">0</td>
|
|
</tr>
|
|
<tr>
|
|
<td class="headerItem">Test Date:</td>
|
|
<td class="headerValue">2024-04-30 13:17:26</td>
|
|
<td></td>
|
|
<td class="headerItem">Functions:</td>
|
|
<td class="headerCovTableEntryLo">0.0 %</td>
|
|
<td class="headerCovTableEntry">2</td>
|
|
<td class="headerCovTableEntry">0</td>
|
|
</tr>
|
|
<tr><td><img src="../../../../glass.png" width=3 height=3 alt=""></td></tr>
|
|
</table>
|
|
</td>
|
|
</tr>
|
|
|
|
<tr><td class="ruler"><img src="../../../../glass.png" width=3 height=3 alt=""></td></tr>
|
|
</table>
|
|
|
|
<table cellpadding=0 cellspacing=0 border=0>
|
|
<tr>
|
|
<td><br></td>
|
|
</tr>
|
|
<tr>
|
|
<td>
|
|
<pre class="sourceHeading"> Line data Source code</pre>
|
|
<pre class="source">
|
|
<span id="L1"><span class="lineNum"> 1</span> : #pragma once</span>
|
|
<span id="L2"><span class="lineNum"> 2</span> : </span>
|
|
<span id="L3"><span class="lineNum"> 3</span> : #include <complex></span>
|
|
<span id="L4"><span class="lineNum"> 4</span> : </span>
|
|
<span id="L5"><span class="lineNum"> 5</span> : #include <c10/macros/Macros.h></span>
|
|
<span id="L6"><span class="lineNum"> 6</span> : </span>
|
|
<span id="L7"><span class="lineNum"> 7</span> : #if defined(__CUDACC__) || defined(__HIPCC__)</span>
|
|
<span id="L8"><span class="lineNum"> 8</span> : #include <thrust/complex.h></span>
|
|
<span id="L9"><span class="lineNum"> 9</span> : #endif</span>
|
|
<span id="L10"><span class="lineNum"> 10</span> : </span>
|
|
<span id="L11"><span class="lineNum"> 11</span> : C10_CLANG_DIAGNOSTIC_PUSH()</span>
|
|
<span id="L12"><span class="lineNum"> 12</span> : #if C10_CLANG_HAS_WARNING("-Wimplicit-float-conversion")</span>
|
|
<span id="L13"><span class="lineNum"> 13</span> : C10_CLANG_DIAGNOSTIC_IGNORE("-Wimplicit-float-conversion")</span>
|
|
<span id="L14"><span class="lineNum"> 14</span> : #endif</span>
|
|
<span id="L15"><span class="lineNum"> 15</span> : #if C10_CLANG_HAS_WARNING("-Wfloat-conversion")</span>
|
|
<span id="L16"><span class="lineNum"> 16</span> : C10_CLANG_DIAGNOSTIC_IGNORE("-Wfloat-conversion")</span>
|
|
<span id="L17"><span class="lineNum"> 17</span> : #endif</span>
|
|
<span id="L18"><span class="lineNum"> 18</span> : </span>
|
|
<span id="L19"><span class="lineNum"> 19</span> : namespace c10 {</span>
|
|
<span id="L20"><span class="lineNum"> 20</span> : </span>
|
|
<span id="L21"><span class="lineNum"> 21</span> : // c10::complex is an implementation of complex numbers that aims</span>
|
|
<span id="L22"><span class="lineNum"> 22</span> : // to work on all devices supported by PyTorch</span>
|
|
<span id="L23"><span class="lineNum"> 23</span> : //</span>
|
|
<span id="L24"><span class="lineNum"> 24</span> : // Most of the APIs duplicates std::complex</span>
|
|
<span id="L25"><span class="lineNum"> 25</span> : // Reference: https://en.cppreference.com/w/cpp/numeric/complex</span>
|
|
<span id="L26"><span class="lineNum"> 26</span> : //</span>
|
|
<span id="L27"><span class="lineNum"> 27</span> : // [NOTE: Complex Operator Unification]</span>
|
|
<span id="L28"><span class="lineNum"> 28</span> : // Operators currently use a mix of std::complex, thrust::complex, and</span>
|
|
<span id="L29"><span class="lineNum"> 29</span> : // c10::complex internally. The end state is that all operators will use</span>
|
|
<span id="L30"><span class="lineNum"> 30</span> : // c10::complex internally. Until then, there may be some hacks to support all</span>
|
|
<span id="L31"><span class="lineNum"> 31</span> : // variants.</span>
|
|
<span id="L32"><span class="lineNum"> 32</span> : //</span>
|
|
<span id="L33"><span class="lineNum"> 33</span> : //</span>
|
|
<span id="L34"><span class="lineNum"> 34</span> : // [Note on Constructors]</span>
|
|
<span id="L35"><span class="lineNum"> 35</span> : //</span>
|
|
<span id="L36"><span class="lineNum"> 36</span> : // The APIs of constructors are mostly copied from C++ standard:</span>
|
|
<span id="L37"><span class="lineNum"> 37</span> : // https://en.cppreference.com/w/cpp/numeric/complex/complex</span>
|
|
<span id="L38"><span class="lineNum"> 38</span> : //</span>
|
|
<span id="L39"><span class="lineNum"> 39</span> : // Since C++14, all constructors are constexpr in std::complex</span>
|
|
<span id="L40"><span class="lineNum"> 40</span> : //</span>
|
|
<span id="L41"><span class="lineNum"> 41</span> : // There are three types of constructors:</span>
|
|
<span id="L42"><span class="lineNum"> 42</span> : // - initializing from real and imag:</span>
|
|
<span id="L43"><span class="lineNum"> 43</span> : // `constexpr complex( const T& re = T(), const T& im = T() );`</span>
|
|
<span id="L44"><span class="lineNum"> 44</span> : // - implicitly-declared copy constructor</span>
|
|
<span id="L45"><span class="lineNum"> 45</span> : // - converting constructors</span>
|
|
<span id="L46"><span class="lineNum"> 46</span> : //</span>
|
|
<span id="L47"><span class="lineNum"> 47</span> : // Converting constructors:</span>
|
|
<span id="L48"><span class="lineNum"> 48</span> : // - std::complex defines converting constructor between float/double/long</span>
|
|
<span id="L49"><span class="lineNum"> 49</span> : // double,</span>
|
|
<span id="L50"><span class="lineNum"> 50</span> : // while we define converting constructor between float/double.</span>
|
|
<span id="L51"><span class="lineNum"> 51</span> : // - For these converting constructors, upcasting is implicit, downcasting is</span>
|
|
<span id="L52"><span class="lineNum"> 52</span> : // explicit.</span>
|
|
<span id="L53"><span class="lineNum"> 53</span> : // - We also define explicit casting from std::complex/thrust::complex</span>
|
|
<span id="L54"><span class="lineNum"> 54</span> : // - Note that the conversion from thrust is not constexpr, because</span>
|
|
<span id="L55"><span class="lineNum"> 55</span> : // thrust does not define them as constexpr ????</span>
|
|
<span id="L56"><span class="lineNum"> 56</span> : //</span>
|
|
<span id="L57"><span class="lineNum"> 57</span> : //</span>
|
|
<span id="L58"><span class="lineNum"> 58</span> : // [Operator =]</span>
|
|
<span id="L59"><span class="lineNum"> 59</span> : //</span>
|
|
<span id="L60"><span class="lineNum"> 60</span> : // The APIs of operator = are mostly copied from C++ standard:</span>
|
|
<span id="L61"><span class="lineNum"> 61</span> : // https://en.cppreference.com/w/cpp/numeric/complex/operator%3D</span>
|
|
<span id="L62"><span class="lineNum"> 62</span> : //</span>
|
|
<span id="L63"><span class="lineNum"> 63</span> : // Since C++20, all operator= are constexpr. Although we are not building with</span>
|
|
<span id="L64"><span class="lineNum"> 64</span> : // C++20, we also obey this behavior.</span>
|
|
<span id="L65"><span class="lineNum"> 65</span> : //</span>
|
|
<span id="L66"><span class="lineNum"> 66</span> : // There are three types of assign operator:</span>
|
|
<span id="L67"><span class="lineNum"> 67</span> : // - Assign a real value from the same scalar type</span>
|
|
<span id="L68"><span class="lineNum"> 68</span> : // - In std, this is templated as complex& operator=(const T& x)</span>
|
|
<span id="L69"><span class="lineNum"> 69</span> : // with specialization `complex& operator=(T x)` for float/double/long</span>
|
|
<span id="L70"><span class="lineNum"> 70</span> : // double Since we only support float and double, on will use `complex&</span>
|
|
<span id="L71"><span class="lineNum"> 71</span> : // operator=(T x)`</span>
|
|
<span id="L72"><span class="lineNum"> 72</span> : // - Copy assignment operator and converting assignment operator</span>
|
|
<span id="L73"><span class="lineNum"> 73</span> : // - There is no specialization of converting assignment operators, which type</span>
|
|
<span id="L74"><span class="lineNum"> 74</span> : // is</span>
|
|
<span id="L75"><span class="lineNum"> 75</span> : // convertible is solely dependent on whether the scalar type is convertible</span>
|
|
<span id="L76"><span class="lineNum"> 76</span> : //</span>
|
|
<span id="L77"><span class="lineNum"> 77</span> : // In addition to the standard assignment, we also provide assignment operators</span>
|
|
<span id="L78"><span class="lineNum"> 78</span> : // with std and thrust</span>
|
|
<span id="L79"><span class="lineNum"> 79</span> : //</span>
|
|
<span id="L80"><span class="lineNum"> 80</span> : //</span>
|
|
<span id="L81"><span class="lineNum"> 81</span> : // [Casting operators]</span>
|
|
<span id="L82"><span class="lineNum"> 82</span> : //</span>
|
|
<span id="L83"><span class="lineNum"> 83</span> : // std::complex does not have casting operators. We define casting operators</span>
|
|
<span id="L84"><span class="lineNum"> 84</span> : // casting to std::complex and thrust::complex</span>
|
|
<span id="L85"><span class="lineNum"> 85</span> : //</span>
|
|
<span id="L86"><span class="lineNum"> 86</span> : //</span>
|
|
<span id="L87"><span class="lineNum"> 87</span> : // [Operator ""]</span>
|
|
<span id="L88"><span class="lineNum"> 88</span> : //</span>
|
|
<span id="L89"><span class="lineNum"> 89</span> : // std::complex has custom literals `i`, `if` and `il` defined in namespace</span>
|
|
<span id="L90"><span class="lineNum"> 90</span> : // `std::literals::complex_literals`. We define our own custom literals in the</span>
|
|
<span id="L91"><span class="lineNum"> 91</span> : // namespace `c10::complex_literals`. Our custom literals does not follow the</span>
|
|
<span id="L92"><span class="lineNum"> 92</span> : // same behavior as in std::complex, instead, we define _if, _id to construct</span>
|
|
<span id="L93"><span class="lineNum"> 93</span> : // float/double complex literals.</span>
|
|
<span id="L94"><span class="lineNum"> 94</span> : //</span>
|
|
<span id="L95"><span class="lineNum"> 95</span> : //</span>
|
|
<span id="L96"><span class="lineNum"> 96</span> : // [real() and imag()]</span>
|
|
<span id="L97"><span class="lineNum"> 97</span> : //</span>
|
|
<span id="L98"><span class="lineNum"> 98</span> : // In C++20, there are two overload of these functions, one it to return the</span>
|
|
<span id="L99"><span class="lineNum"> 99</span> : // real/imag, another is to set real/imag, they are both constexpr. We follow</span>
|
|
<span id="L100"><span class="lineNum"> 100</span> : // this design.</span>
|
|
<span id="L101"><span class="lineNum"> 101</span> : //</span>
|
|
<span id="L102"><span class="lineNum"> 102</span> : //</span>
|
|
<span id="L103"><span class="lineNum"> 103</span> : // [Operator +=,-=,*=,/=]</span>
|
|
<span id="L104"><span class="lineNum"> 104</span> : //</span>
|
|
<span id="L105"><span class="lineNum"> 105</span> : // Since C++20, these operators become constexpr. In our implementation, they</span>
|
|
<span id="L106"><span class="lineNum"> 106</span> : // are also constexpr.</span>
|
|
<span id="L107"><span class="lineNum"> 107</span> : //</span>
|
|
<span id="L108"><span class="lineNum"> 108</span> : // There are two types of such operators: operating with a real number, or</span>
|
|
<span id="L109"><span class="lineNum"> 109</span> : // operating with another complex number. For the operating with a real number,</span>
|
|
<span id="L110"><span class="lineNum"> 110</span> : // the generic template form has argument type `const T &`, while the overload</span>
|
|
<span id="L111"><span class="lineNum"> 111</span> : // for float/double/long double has `T`. We will follow the same type as</span>
|
|
<span id="L112"><span class="lineNum"> 112</span> : // float/double/long double in std.</span>
|
|
<span id="L113"><span class="lineNum"> 113</span> : //</span>
|
|
<span id="L114"><span class="lineNum"> 114</span> : // [Unary operator +-]</span>
|
|
<span id="L115"><span class="lineNum"> 115</span> : //</span>
|
|
<span id="L116"><span class="lineNum"> 116</span> : // Since C++20, they are constexpr. We also make them expr</span>
|
|
<span id="L117"><span class="lineNum"> 117</span> : //</span>
|
|
<span id="L118"><span class="lineNum"> 118</span> : // [Binary operators +-*/]</span>
|
|
<span id="L119"><span class="lineNum"> 119</span> : //</span>
|
|
<span id="L120"><span class="lineNum"> 120</span> : // Each operator has three versions (taking + as example):</span>
|
|
<span id="L121"><span class="lineNum"> 121</span> : // - complex + complex</span>
|
|
<span id="L122"><span class="lineNum"> 122</span> : // - complex + real</span>
|
|
<span id="L123"><span class="lineNum"> 123</span> : // - real + complex</span>
|
|
<span id="L124"><span class="lineNum"> 124</span> : //</span>
|
|
<span id="L125"><span class="lineNum"> 125</span> : // [Operator ==, !=]</span>
|
|
<span id="L126"><span class="lineNum"> 126</span> : //</span>
|
|
<span id="L127"><span class="lineNum"> 127</span> : // Each operator has three versions (taking == as example):</span>
|
|
<span id="L128"><span class="lineNum"> 128</span> : // - complex == complex</span>
|
|
<span id="L129"><span class="lineNum"> 129</span> : // - complex == real</span>
|
|
<span id="L130"><span class="lineNum"> 130</span> : // - real == complex</span>
|
|
<span id="L131"><span class="lineNum"> 131</span> : //</span>
|
|
<span id="L132"><span class="lineNum"> 132</span> : // Some of them are removed on C++20, but we decide to keep them</span>
|
|
<span id="L133"><span class="lineNum"> 133</span> : //</span>
|
|
<span id="L134"><span class="lineNum"> 134</span> : // [Operator <<, >>]</span>
|
|
<span id="L135"><span class="lineNum"> 135</span> : //</span>
|
|
<span id="L136"><span class="lineNum"> 136</span> : // These are implemented by casting to std::complex</span>
|
|
<span id="L137"><span class="lineNum"> 137</span> : //</span>
|
|
<span id="L138"><span class="lineNum"> 138</span> : //</span>
|
|
<span id="L139"><span class="lineNum"> 139</span> : //</span>
|
|
<span id="L140"><span class="lineNum"> 140</span> : // TODO(@zasdfgbnm): c10::complex<c10::Half> is not currently supported,</span>
|
|
<span id="L141"><span class="lineNum"> 141</span> : // because:</span>
|
|
<span id="L142"><span class="lineNum"> 142</span> : // - lots of members and functions of c10::Half are not constexpr</span>
|
|
<span id="L143"><span class="lineNum"> 143</span> : // - thrust::complex only support float and double</span>
|
|
<span id="L144"><span class="lineNum"> 144</span> : </span>
|
|
<span id="L145"><span class="lineNum"> 145</span> : template <typename T></span>
|
|
<span id="L146"><span class="lineNum"> 146</span> : struct alignas(sizeof(T) * 2) complex {</span>
|
|
<span id="L147"><span class="lineNum"> 147</span> : using value_type = T;</span>
|
|
<span id="L148"><span class="lineNum"> 148</span> : </span>
|
|
<span id="L149"><span class="lineNum"> 149</span> : T real_ = T(0);</span>
|
|
<span id="L150"><span class="lineNum"> 150</span> : T imag_ = T(0);</span>
|
|
<span id="L151"><span class="lineNum"> 151</span> : </span>
|
|
<span id="L152"><span class="lineNum"> 152</span> : constexpr complex() = default;</span>
|
|
<span id="L153"><span class="lineNum"> 153</span> : C10_HOST_DEVICE constexpr complex(const T& re, const T& im = T())</span>
|
|
<span id="L154"><span class="lineNum"> 154</span> : : real_(re), imag_(im) {}</span>
|
|
<span id="L155"><span class="lineNum"> 155</span> : template <typename U></span>
|
|
<span id="L156"><span class="lineNum"> 156</span> : explicit constexpr complex(const std::complex<U>& other)</span>
|
|
<span id="L157"><span class="lineNum"> 157</span> : : complex(other.real(), other.imag()) {}</span>
|
|
<span id="L158"><span class="lineNum"> 158</span> : #if defined(__CUDACC__) || defined(__HIPCC__)</span>
|
|
<span id="L159"><span class="lineNum"> 159</span> : template <typename U></span>
|
|
<span id="L160"><span class="lineNum"> 160</span> : explicit C10_HOST_DEVICE complex(const thrust::complex<U>& other)</span>
|
|
<span id="L161"><span class="lineNum"> 161</span> : : real_(other.real()), imag_(other.imag()) {}</span>
|
|
<span id="L162"><span class="lineNum"> 162</span> : // NOTE can not be implemented as follow due to ROCm bug:</span>
|
|
<span id="L163"><span class="lineNum"> 163</span> : // explicit C10_HOST_DEVICE complex(const thrust::complex<U> &other):</span>
|
|
<span id="L164"><span class="lineNum"> 164</span> : // complex(other.real(), other.imag()) {}</span>
|
|
<span id="L165"><span class="lineNum"> 165</span> : #endif</span>
|
|
<span id="L166"><span class="lineNum"> 166</span> : </span>
|
|
<span id="L167"><span class="lineNum"> 167</span> : // Use SFINAE to specialize casting constructor for c10::complex<float> and</span>
|
|
<span id="L168"><span class="lineNum"> 168</span> : // c10::complex<double></span>
|
|
<span id="L169"><span class="lineNum"> 169</span> : template <typename U = T></span>
|
|
<span id="L170"><span class="lineNum"> 170</span> : C10_HOST_DEVICE explicit constexpr complex(</span>
|
|
<span id="L171"><span class="lineNum"> 171</span> : const std::enable_if_t<std::is_same<U, float>::value, complex<double>>&</span>
|
|
<span id="L172"><span class="lineNum"> 172</span> : other)</span>
|
|
<span id="L173"><span class="lineNum"> 173</span> : : real_(other.real_), imag_(other.imag_) {}</span>
|
|
<span id="L174"><span class="lineNum"> 174</span> : template <typename U = T></span>
|
|
<span id="L175"><span class="lineNum"> 175</span> : C10_HOST_DEVICE constexpr complex(</span>
|
|
<span id="L176"><span class="lineNum"> 176</span> : const std::enable_if_t<std::is_same<U, double>::value, complex<float>>&</span>
|
|
<span id="L177"><span class="lineNum"> 177</span> : other)</span>
|
|
<span id="L178"><span class="lineNum"> 178</span> : : real_(other.real_), imag_(other.imag_) {}</span>
|
|
<span id="L179"><span class="lineNum"> 179</span> : </span>
|
|
<span id="L180"><span class="lineNum"> 180</span> : constexpr complex<T>& operator=(T re) {</span>
|
|
<span id="L181"><span class="lineNum"> 181</span> : real_ = re;</span>
|
|
<span id="L182"><span class="lineNum"> 182</span> : imag_ = 0;</span>
|
|
<span id="L183"><span class="lineNum"> 183</span> : return *this;</span>
|
|
<span id="L184"><span class="lineNum"> 184</span> : }</span>
|
|
<span id="L185"><span class="lineNum"> 185</span> : </span>
|
|
<span id="L186"><span class="lineNum"> 186</span> : constexpr complex<T>& operator+=(T re) {</span>
|
|
<span id="L187"><span class="lineNum"> 187</span> : real_ += re;</span>
|
|
<span id="L188"><span class="lineNum"> 188</span> : return *this;</span>
|
|
<span id="L189"><span class="lineNum"> 189</span> : }</span>
|
|
<span id="L190"><span class="lineNum"> 190</span> : </span>
|
|
<span id="L191"><span class="lineNum"> 191</span> : constexpr complex<T>& operator-=(T re) {</span>
|
|
<span id="L192"><span class="lineNum"> 192</span> : real_ -= re;</span>
|
|
<span id="L193"><span class="lineNum"> 193</span> : return *this;</span>
|
|
<span id="L194"><span class="lineNum"> 194</span> : }</span>
|
|
<span id="L195"><span class="lineNum"> 195</span> : </span>
|
|
<span id="L196"><span class="lineNum"> 196</span> : constexpr complex<T>& operator*=(T re) {</span>
|
|
<span id="L197"><span class="lineNum"> 197</span> : real_ *= re;</span>
|
|
<span id="L198"><span class="lineNum"> 198</span> : imag_ *= re;</span>
|
|
<span id="L199"><span class="lineNum"> 199</span> : return *this;</span>
|
|
<span id="L200"><span class="lineNum"> 200</span> : }</span>
|
|
<span id="L201"><span class="lineNum"> 201</span> : </span>
|
|
<span id="L202"><span class="lineNum"> 202</span> : constexpr complex<T>& operator/=(T re) {</span>
|
|
<span id="L203"><span class="lineNum"> 203</span> : real_ /= re;</span>
|
|
<span id="L204"><span class="lineNum"> 204</span> : imag_ /= re;</span>
|
|
<span id="L205"><span class="lineNum"> 205</span> : return *this;</span>
|
|
<span id="L206"><span class="lineNum"> 206</span> : }</span>
|
|
<span id="L207"><span class="lineNum"> 207</span> : </span>
|
|
<span id="L208"><span class="lineNum"> 208</span> : template <typename U></span>
|
|
<span id="L209"><span class="lineNum"> 209</span> : constexpr complex<T>& operator=(const complex<U>& rhs) {</span>
|
|
<span id="L210"><span class="lineNum"> 210</span> : real_ = rhs.real();</span>
|
|
<span id="L211"><span class="lineNum"> 211</span> : imag_ = rhs.imag();</span>
|
|
<span id="L212"><span class="lineNum"> 212</span> : return *this;</span>
|
|
<span id="L213"><span class="lineNum"> 213</span> : }</span>
|
|
<span id="L214"><span class="lineNum"> 214</span> : </span>
|
|
<span id="L215"><span class="lineNum"> 215</span> : template <typename U></span>
|
|
<span id="L216"><span class="lineNum"> 216</span> : constexpr complex<T>& operator+=(const complex<U>& rhs) {</span>
|
|
<span id="L217"><span class="lineNum"> 217</span> : real_ += rhs.real();</span>
|
|
<span id="L218"><span class="lineNum"> 218</span> : imag_ += rhs.imag();</span>
|
|
<span id="L219"><span class="lineNum"> 219</span> : return *this;</span>
|
|
<span id="L220"><span class="lineNum"> 220</span> : }</span>
|
|
<span id="L221"><span class="lineNum"> 221</span> : </span>
|
|
<span id="L222"><span class="lineNum"> 222</span> : template <typename U></span>
|
|
<span id="L223"><span class="lineNum"> 223</span> : constexpr complex<T>& operator-=(const complex<U>& rhs) {</span>
|
|
<span id="L224"><span class="lineNum"> 224</span> : real_ -= rhs.real();</span>
|
|
<span id="L225"><span class="lineNum"> 225</span> : imag_ -= rhs.imag();</span>
|
|
<span id="L226"><span class="lineNum"> 226</span> : return *this;</span>
|
|
<span id="L227"><span class="lineNum"> 227</span> : }</span>
|
|
<span id="L228"><span class="lineNum"> 228</span> : </span>
|
|
<span id="L229"><span class="lineNum"> 229</span> : template <typename U></span>
|
|
<span id="L230"><span class="lineNum"> 230</span> : constexpr complex<T>& operator*=(const complex<U>& rhs) {</span>
|
|
<span id="L231"><span class="lineNum"> 231</span> : // (a + bi) * (c + di) = (a*c - b*d) + (a * d + b * c) i</span>
|
|
<span id="L232"><span class="lineNum"> 232</span> : T a = real_;</span>
|
|
<span id="L233"><span class="lineNum"> 233</span> : T b = imag_;</span>
|
|
<span id="L234"><span class="lineNum"> 234</span> : U c = rhs.real();</span>
|
|
<span id="L235"><span class="lineNum"> 235</span> : U d = rhs.imag();</span>
|
|
<span id="L236"><span class="lineNum"> 236</span> : real_ = a * c - b * d;</span>
|
|
<span id="L237"><span class="lineNum"> 237</span> : imag_ = a * d + b * c;</span>
|
|
<span id="L238"><span class="lineNum"> 238</span> : return *this;</span>
|
|
<span id="L239"><span class="lineNum"> 239</span> : }</span>
|
|
<span id="L240"><span class="lineNum"> 240</span> : </span>
|
|
<span id="L241"><span class="lineNum"> 241</span> : #ifdef __APPLE__</span>
|
|
<span id="L242"><span class="lineNum"> 242</span> : #define FORCE_INLINE_APPLE __attribute__((always_inline))</span>
|
|
<span id="L243"><span class="lineNum"> 243</span> : #else</span>
|
|
<span id="L244"><span class="lineNum"> 244</span> : #define FORCE_INLINE_APPLE</span>
|
|
<span id="L245"><span class="lineNum"> 245</span> : #endif</span>
|
|
<span id="L246"><span class="lineNum"> 246</span> : template <typename U></span>
|
|
<span id="L247"><span class="lineNum"> 247</span> : constexpr FORCE_INLINE_APPLE complex<T>& operator/=(const complex<U>& rhs)</span>
|
|
<span id="L248"><span class="lineNum"> 248</span> : __ubsan_ignore_float_divide_by_zero__ {</span>
|
|
<span id="L249"><span class="lineNum"> 249</span> : // (a + bi) / (c + di) = (ac + bd)/(c^2 + d^2) + (bc - ad)/(c^2 + d^2) i</span>
|
|
<span id="L250"><span class="lineNum"> 250</span> : // the calculation below follows numpy's complex division</span>
|
|
<span id="L251"><span class="lineNum"> 251</span> : T a = real_;</span>
|
|
<span id="L252"><span class="lineNum"> 252</span> : T b = imag_;</span>
|
|
<span id="L253"><span class="lineNum"> 253</span> : U c = rhs.real();</span>
|
|
<span id="L254"><span class="lineNum"> 254</span> : U d = rhs.imag();</span>
|
|
<span id="L255"><span class="lineNum"> 255</span> : </span>
|
|
<span id="L256"><span class="lineNum"> 256</span> : #if defined(__GNUC__) && !defined(__clang__)</span>
|
|
<span id="L257"><span class="lineNum"> 257</span> : // std::abs is already constexpr by gcc</span>
|
|
<span id="L258"><span class="lineNum"> 258</span> : auto abs_c = std::abs(c);</span>
|
|
<span id="L259"><span class="lineNum"> 259</span> : auto abs_d = std::abs(d);</span>
|
|
<span id="L260"><span class="lineNum"> 260</span> : #else</span>
|
|
<span id="L261"><span class="lineNum"> 261</span> : auto abs_c = c < 0 ? -c : c;</span>
|
|
<span id="L262"><span class="lineNum"> 262</span> : auto abs_d = d < 0 ? -d : d;</span>
|
|
<span id="L263"><span class="lineNum"> 263</span> : #endif</span>
|
|
<span id="L264"><span class="lineNum"> 264</span> : </span>
|
|
<span id="L265"><span class="lineNum"> 265</span> : if (abs_c >= abs_d) {</span>
|
|
<span id="L266"><span class="lineNum"> 266</span> : if (abs_c == 0 && abs_d == 0) {</span>
|
|
<span id="L267"><span class="lineNum"> 267</span> : /* divide by zeros should yield a complex inf or nan */</span>
|
|
<span id="L268"><span class="lineNum"> 268</span> : real_ = a / abs_c;</span>
|
|
<span id="L269"><span class="lineNum"> 269</span> : imag_ = b / abs_d;</span>
|
|
<span id="L270"><span class="lineNum"> 270</span> : } else {</span>
|
|
<span id="L271"><span class="lineNum"> 271</span> : auto rat = d / c;</span>
|
|
<span id="L272"><span class="lineNum"> 272</span> : auto scl = 1.0 / (c + d * rat);</span>
|
|
<span id="L273"><span class="lineNum"> 273</span> : real_ = (a + b * rat) * scl;</span>
|
|
<span id="L274"><span class="lineNum"> 274</span> : imag_ = (b - a * rat) * scl;</span>
|
|
<span id="L275"><span class="lineNum"> 275</span> : }</span>
|
|
<span id="L276"><span class="lineNum"> 276</span> : } else {</span>
|
|
<span id="L277"><span class="lineNum"> 277</span> : auto rat = c / d;</span>
|
|
<span id="L278"><span class="lineNum"> 278</span> : auto scl = 1.0 / (d + c * rat);</span>
|
|
<span id="L279"><span class="lineNum"> 279</span> : real_ = (a * rat + b) * scl;</span>
|
|
<span id="L280"><span class="lineNum"> 280</span> : imag_ = (b * rat - a) * scl;</span>
|
|
<span id="L281"><span class="lineNum"> 281</span> : }</span>
|
|
<span id="L282"><span class="lineNum"> 282</span> : return *this;</span>
|
|
<span id="L283"><span class="lineNum"> 283</span> : }</span>
|
|
<span id="L284"><span class="lineNum"> 284</span> : #undef FORCE_INLINE_APPLE</span>
|
|
<span id="L285"><span class="lineNum"> 285</span> : </span>
|
|
<span id="L286"><span class="lineNum"> 286</span> : template <typename U></span>
|
|
<span id="L287"><span class="lineNum"> 287</span> : constexpr complex<T>& operator=(const std::complex<U>& rhs) {</span>
|
|
<span id="L288"><span class="lineNum"> 288</span> : real_ = rhs.real();</span>
|
|
<span id="L289"><span class="lineNum"> 289</span> : imag_ = rhs.imag();</span>
|
|
<span id="L290"><span class="lineNum"> 290</span> : return *this;</span>
|
|
<span id="L291"><span class="lineNum"> 291</span> : }</span>
|
|
<span id="L292"><span class="lineNum"> 292</span> : </span>
|
|
<span id="L293"><span class="lineNum"> 293</span> : #if defined(__CUDACC__) || defined(__HIPCC__)</span>
|
|
<span id="L294"><span class="lineNum"> 294</span> : template <typename U></span>
|
|
<span id="L295"><span class="lineNum"> 295</span> : C10_HOST_DEVICE complex<T>& operator=(const thrust::complex<U>& rhs) {</span>
|
|
<span id="L296"><span class="lineNum"> 296</span> : real_ = rhs.real();</span>
|
|
<span id="L297"><span class="lineNum"> 297</span> : imag_ = rhs.imag();</span>
|
|
<span id="L298"><span class="lineNum"> 298</span> : return *this;</span>
|
|
<span id="L299"><span class="lineNum"> 299</span> : }</span>
|
|
<span id="L300"><span class="lineNum"> 300</span> : #endif</span>
|
|
<span id="L301"><span class="lineNum"> 301</span> : </span>
|
|
<span id="L302"><span class="lineNum"> 302</span> : template <typename U></span>
|
|
<span id="L303"><span class="lineNum"> 303</span> : explicit constexpr operator std::complex<U>() const {</span>
|
|
<span id="L304"><span class="lineNum"> 304</span> : return std::complex<U>(std::complex<T>(real(), imag()));</span>
|
|
<span id="L305"><span class="lineNum"> 305</span> : }</span>
|
|
<span id="L306"><span class="lineNum"> 306</span> : </span>
|
|
<span id="L307"><span class="lineNum"> 307</span> : #if defined(__CUDACC__) || defined(__HIPCC__)</span>
|
|
<span id="L308"><span class="lineNum"> 308</span> : template <typename U></span>
|
|
<span id="L309"><span class="lineNum"> 309</span> : C10_HOST_DEVICE explicit operator thrust::complex<U>() const {</span>
|
|
<span id="L310"><span class="lineNum"> 310</span> : return static_cast<thrust::complex<U>>(thrust::complex<T>(real(), imag()));</span>
|
|
<span id="L311"><span class="lineNum"> 311</span> : }</span>
|
|
<span id="L312"><span class="lineNum"> 312</span> : #endif</span>
|
|
<span id="L313"><span class="lineNum"> 313</span> : </span>
|
|
<span id="L314"><span class="lineNum"> 314</span> : // consistent with NumPy behavior</span>
|
|
<span id="L315"><span class="lineNum"> 315</span> : explicit constexpr operator bool() const {</span>
|
|
<span id="L316"><span class="lineNum"> 316</span> : return real() || imag();</span>
|
|
<span id="L317"><span class="lineNum"> 317</span> : }</span>
|
|
<span id="L318"><span class="lineNum"> 318</span> : </span>
|
|
<span id="L319"><span class="lineNum"> 319</span> <span class="tlaUNC tlaBgUNC"> 0 : C10_HOST_DEVICE constexpr T real() const {</span></span>
|
|
<span id="L320"><span class="lineNum"> 320</span> <span class="tlaUNC"> 0 : return real_;</span></span>
|
|
<span id="L321"><span class="lineNum"> 321</span> : }</span>
|
|
<span id="L322"><span class="lineNum"> 322</span> : constexpr void real(T value) {</span>
|
|
<span id="L323"><span class="lineNum"> 323</span> : real_ = value;</span>
|
|
<span id="L324"><span class="lineNum"> 324</span> : }</span>
|
|
<span id="L325"><span class="lineNum"> 325</span> <span class="tlaUNC"> 0 : constexpr T imag() const {</span></span>
|
|
<span id="L326"><span class="lineNum"> 326</span> <span class="tlaUNC"> 0 : return imag_;</span></span>
|
|
<span id="L327"><span class="lineNum"> 327</span> : }</span>
|
|
<span id="L328"><span class="lineNum"> 328</span> : constexpr void imag(T value) {</span>
|
|
<span id="L329"><span class="lineNum"> 329</span> : imag_ = value;</span>
|
|
<span id="L330"><span class="lineNum"> 330</span> : }</span>
|
|
<span id="L331"><span class="lineNum"> 331</span> : };</span>
|
|
<span id="L332"><span class="lineNum"> 332</span> : </span>
|
|
<span id="L333"><span class="lineNum"> 333</span> : namespace complex_literals {</span>
|
|
<span id="L334"><span class="lineNum"> 334</span> : </span>
|
|
<span id="L335"><span class="lineNum"> 335</span> : constexpr complex<float> operator"" _if(long double imag) {</span>
|
|
<span id="L336"><span class="lineNum"> 336</span> : return complex<float>(0.0f, static_cast<float>(imag));</span>
|
|
<span id="L337"><span class="lineNum"> 337</span> : }</span>
|
|
<span id="L338"><span class="lineNum"> 338</span> : </span>
|
|
<span id="L339"><span class="lineNum"> 339</span> : constexpr complex<double> operator"" _id(long double imag) {</span>
|
|
<span id="L340"><span class="lineNum"> 340</span> : return complex<double>(0.0, static_cast<double>(imag));</span>
|
|
<span id="L341"><span class="lineNum"> 341</span> : }</span>
|
|
<span id="L342"><span class="lineNum"> 342</span> : </span>
|
|
<span id="L343"><span class="lineNum"> 343</span> : constexpr complex<float> operator"" _if(unsigned long long imag) {</span>
|
|
<span id="L344"><span class="lineNum"> 344</span> : return complex<float>(0.0f, static_cast<float>(imag));</span>
|
|
<span id="L345"><span class="lineNum"> 345</span> : }</span>
|
|
<span id="L346"><span class="lineNum"> 346</span> : </span>
|
|
<span id="L347"><span class="lineNum"> 347</span> : constexpr complex<double> operator"" _id(unsigned long long imag) {</span>
|
|
<span id="L348"><span class="lineNum"> 348</span> : return complex<double>(0.0, static_cast<double>(imag));</span>
|
|
<span id="L349"><span class="lineNum"> 349</span> : }</span>
|
|
<span id="L350"><span class="lineNum"> 350</span> : </span>
|
|
<span id="L351"><span class="lineNum"> 351</span> : } // namespace complex_literals</span>
|
|
<span id="L352"><span class="lineNum"> 352</span> : </span>
|
|
<span id="L353"><span class="lineNum"> 353</span> : template <typename T></span>
|
|
<span id="L354"><span class="lineNum"> 354</span> : constexpr complex<T> operator+(const complex<T>& val) {</span>
|
|
<span id="L355"><span class="lineNum"> 355</span> : return val;</span>
|
|
<span id="L356"><span class="lineNum"> 356</span> : }</span>
|
|
<span id="L357"><span class="lineNum"> 357</span> : </span>
|
|
<span id="L358"><span class="lineNum"> 358</span> : template <typename T></span>
|
|
<span id="L359"><span class="lineNum"> 359</span> : constexpr complex<T> operator-(const complex<T>& val) {</span>
|
|
<span id="L360"><span class="lineNum"> 360</span> : return complex<T>(-val.real(), -val.imag());</span>
|
|
<span id="L361"><span class="lineNum"> 361</span> : }</span>
|
|
<span id="L362"><span class="lineNum"> 362</span> : </span>
|
|
<span id="L363"><span class="lineNum"> 363</span> : template <typename T></span>
|
|
<span id="L364"><span class="lineNum"> 364</span> : constexpr complex<T> operator+(const complex<T>& lhs, const complex<T>& rhs) {</span>
|
|
<span id="L365"><span class="lineNum"> 365</span> : complex<T> result = lhs;</span>
|
|
<span id="L366"><span class="lineNum"> 366</span> : return result += rhs;</span>
|
|
<span id="L367"><span class="lineNum"> 367</span> : }</span>
|
|
<span id="L368"><span class="lineNum"> 368</span> : </span>
|
|
<span id="L369"><span class="lineNum"> 369</span> : template <typename T></span>
|
|
<span id="L370"><span class="lineNum"> 370</span> : constexpr complex<T> operator+(const complex<T>& lhs, const T& rhs) {</span>
|
|
<span id="L371"><span class="lineNum"> 371</span> : complex<T> result = lhs;</span>
|
|
<span id="L372"><span class="lineNum"> 372</span> : return result += rhs;</span>
|
|
<span id="L373"><span class="lineNum"> 373</span> : }</span>
|
|
<span id="L374"><span class="lineNum"> 374</span> : </span>
|
|
<span id="L375"><span class="lineNum"> 375</span> : template <typename T></span>
|
|
<span id="L376"><span class="lineNum"> 376</span> : constexpr complex<T> operator+(const T& lhs, const complex<T>& rhs) {</span>
|
|
<span id="L377"><span class="lineNum"> 377</span> : return complex<T>(lhs + rhs.real(), rhs.imag());</span>
|
|
<span id="L378"><span class="lineNum"> 378</span> : }</span>
|
|
<span id="L379"><span class="lineNum"> 379</span> : </span>
|
|
<span id="L380"><span class="lineNum"> 380</span> : template <typename T></span>
|
|
<span id="L381"><span class="lineNum"> 381</span> : constexpr complex<T> operator-(const complex<T>& lhs, const complex<T>& rhs) {</span>
|
|
<span id="L382"><span class="lineNum"> 382</span> : complex<T> result = lhs;</span>
|
|
<span id="L383"><span class="lineNum"> 383</span> : return result -= rhs;</span>
|
|
<span id="L384"><span class="lineNum"> 384</span> : }</span>
|
|
<span id="L385"><span class="lineNum"> 385</span> : </span>
|
|
<span id="L386"><span class="lineNum"> 386</span> : template <typename T></span>
|
|
<span id="L387"><span class="lineNum"> 387</span> : constexpr complex<T> operator-(const complex<T>& lhs, const T& rhs) {</span>
|
|
<span id="L388"><span class="lineNum"> 388</span> : complex<T> result = lhs;</span>
|
|
<span id="L389"><span class="lineNum"> 389</span> : return result -= rhs;</span>
|
|
<span id="L390"><span class="lineNum"> 390</span> : }</span>
|
|
<span id="L391"><span class="lineNum"> 391</span> : </span>
|
|
<span id="L392"><span class="lineNum"> 392</span> : template <typename T></span>
|
|
<span id="L393"><span class="lineNum"> 393</span> : constexpr complex<T> operator-(const T& lhs, const complex<T>& rhs) {</span>
|
|
<span id="L394"><span class="lineNum"> 394</span> : complex<T> result = -rhs;</span>
|
|
<span id="L395"><span class="lineNum"> 395</span> : return result += lhs;</span>
|
|
<span id="L396"><span class="lineNum"> 396</span> : }</span>
|
|
<span id="L397"><span class="lineNum"> 397</span> : </span>
|
|
<span id="L398"><span class="lineNum"> 398</span> : template <typename T></span>
|
|
<span id="L399"><span class="lineNum"> 399</span> : constexpr complex<T> operator*(const complex<T>& lhs, const complex<T>& rhs) {</span>
|
|
<span id="L400"><span class="lineNum"> 400</span> : complex<T> result = lhs;</span>
|
|
<span id="L401"><span class="lineNum"> 401</span> : return result *= rhs;</span>
|
|
<span id="L402"><span class="lineNum"> 402</span> : }</span>
|
|
<span id="L403"><span class="lineNum"> 403</span> : </span>
|
|
<span id="L404"><span class="lineNum"> 404</span> : template <typename T></span>
|
|
<span id="L405"><span class="lineNum"> 405</span> : constexpr complex<T> operator*(const complex<T>& lhs, const T& rhs) {</span>
|
|
<span id="L406"><span class="lineNum"> 406</span> : complex<T> result = lhs;</span>
|
|
<span id="L407"><span class="lineNum"> 407</span> : return result *= rhs;</span>
|
|
<span id="L408"><span class="lineNum"> 408</span> : }</span>
|
|
<span id="L409"><span class="lineNum"> 409</span> : </span>
|
|
<span id="L410"><span class="lineNum"> 410</span> : template <typename T></span>
|
|
<span id="L411"><span class="lineNum"> 411</span> : constexpr complex<T> operator*(const T& lhs, const complex<T>& rhs) {</span>
|
|
<span id="L412"><span class="lineNum"> 412</span> : complex<T> result = rhs;</span>
|
|
<span id="L413"><span class="lineNum"> 413</span> : return result *= lhs;</span>
|
|
<span id="L414"><span class="lineNum"> 414</span> : }</span>
|
|
<span id="L415"><span class="lineNum"> 415</span> : </span>
|
|
<span id="L416"><span class="lineNum"> 416</span> : template <typename T></span>
|
|
<span id="L417"><span class="lineNum"> 417</span> : constexpr complex<T> operator/(const complex<T>& lhs, const complex<T>& rhs) {</span>
|
|
<span id="L418"><span class="lineNum"> 418</span> : complex<T> result = lhs;</span>
|
|
<span id="L419"><span class="lineNum"> 419</span> : return result /= rhs;</span>
|
|
<span id="L420"><span class="lineNum"> 420</span> : }</span>
|
|
<span id="L421"><span class="lineNum"> 421</span> : </span>
|
|
<span id="L422"><span class="lineNum"> 422</span> : template <typename T></span>
|
|
<span id="L423"><span class="lineNum"> 423</span> : constexpr complex<T> operator/(const complex<T>& lhs, const T& rhs) {</span>
|
|
<span id="L424"><span class="lineNum"> 424</span> : complex<T> result = lhs;</span>
|
|
<span id="L425"><span class="lineNum"> 425</span> : return result /= rhs;</span>
|
|
<span id="L426"><span class="lineNum"> 426</span> : }</span>
|
|
<span id="L427"><span class="lineNum"> 427</span> : </span>
|
|
<span id="L428"><span class="lineNum"> 428</span> : template <typename T></span>
|
|
<span id="L429"><span class="lineNum"> 429</span> : constexpr complex<T> operator/(const T& lhs, const complex<T>& rhs) {</span>
|
|
<span id="L430"><span class="lineNum"> 430</span> : complex<T> result(lhs, T());</span>
|
|
<span id="L431"><span class="lineNum"> 431</span> : return result /= rhs;</span>
|
|
<span id="L432"><span class="lineNum"> 432</span> : }</span>
|
|
<span id="L433"><span class="lineNum"> 433</span> : </span>
|
|
<span id="L434"><span class="lineNum"> 434</span> : // Define operators between integral scalars and c10::complex. std::complex does</span>
|
|
<span id="L435"><span class="lineNum"> 435</span> : // not support this when T is a floating-point number. This is useful because it</span>
|
|
<span id="L436"><span class="lineNum"> 436</span> : // saves a lot of "static_cast" when operate a complex and an integer. This</span>
|
|
<span id="L437"><span class="lineNum"> 437</span> : // makes the code both less verbose and potentially more efficient.</span>
|
|
<span id="L438"><span class="lineNum"> 438</span> : #define COMPLEX_INTEGER_OP_TEMPLATE_CONDITION \</span>
|
|
<span id="L439"><span class="lineNum"> 439</span> : typename std::enable_if_t< \</span>
|
|
<span id="L440"><span class="lineNum"> 440</span> : std::is_floating_point<fT>::value && std::is_integral<iT>::value, \</span>
|
|
<span id="L441"><span class="lineNum"> 441</span> : int> = 0</span>
|
|
<span id="L442"><span class="lineNum"> 442</span> : </span>
|
|
<span id="L443"><span class="lineNum"> 443</span> : template <typename fT, typename iT, COMPLEX_INTEGER_OP_TEMPLATE_CONDITION></span>
|
|
<span id="L444"><span class="lineNum"> 444</span> : constexpr c10::complex<fT> operator+(const c10::complex<fT>& a, const iT& b) {</span>
|
|
<span id="L445"><span class="lineNum"> 445</span> : return a + static_cast<fT>(b);</span>
|
|
<span id="L446"><span class="lineNum"> 446</span> : }</span>
|
|
<span id="L447"><span class="lineNum"> 447</span> : </span>
|
|
<span id="L448"><span class="lineNum"> 448</span> : template <typename fT, typename iT, COMPLEX_INTEGER_OP_TEMPLATE_CONDITION></span>
|
|
<span id="L449"><span class="lineNum"> 449</span> : constexpr c10::complex<fT> operator+(const iT& a, const c10::complex<fT>& b) {</span>
|
|
<span id="L450"><span class="lineNum"> 450</span> : return static_cast<fT>(a) + b;</span>
|
|
<span id="L451"><span class="lineNum"> 451</span> : }</span>
|
|
<span id="L452"><span class="lineNum"> 452</span> : </span>
|
|
<span id="L453"><span class="lineNum"> 453</span> : template <typename fT, typename iT, COMPLEX_INTEGER_OP_TEMPLATE_CONDITION></span>
|
|
<span id="L454"><span class="lineNum"> 454</span> : constexpr c10::complex<fT> operator-(const c10::complex<fT>& a, const iT& b) {</span>
|
|
<span id="L455"><span class="lineNum"> 455</span> : return a - static_cast<fT>(b);</span>
|
|
<span id="L456"><span class="lineNum"> 456</span> : }</span>
|
|
<span id="L457"><span class="lineNum"> 457</span> : </span>
|
|
<span id="L458"><span class="lineNum"> 458</span> : template <typename fT, typename iT, COMPLEX_INTEGER_OP_TEMPLATE_CONDITION></span>
|
|
<span id="L459"><span class="lineNum"> 459</span> : constexpr c10::complex<fT> operator-(const iT& a, const c10::complex<fT>& b) {</span>
|
|
<span id="L460"><span class="lineNum"> 460</span> : return static_cast<fT>(a) - b;</span>
|
|
<span id="L461"><span class="lineNum"> 461</span> : }</span>
|
|
<span id="L462"><span class="lineNum"> 462</span> : </span>
|
|
<span id="L463"><span class="lineNum"> 463</span> : template <typename fT, typename iT, COMPLEX_INTEGER_OP_TEMPLATE_CONDITION></span>
|
|
<span id="L464"><span class="lineNum"> 464</span> : constexpr c10::complex<fT> operator*(const c10::complex<fT>& a, const iT& b) {</span>
|
|
<span id="L465"><span class="lineNum"> 465</span> : return a * static_cast<fT>(b);</span>
|
|
<span id="L466"><span class="lineNum"> 466</span> : }</span>
|
|
<span id="L467"><span class="lineNum"> 467</span> : </span>
|
|
<span id="L468"><span class="lineNum"> 468</span> : template <typename fT, typename iT, COMPLEX_INTEGER_OP_TEMPLATE_CONDITION></span>
|
|
<span id="L469"><span class="lineNum"> 469</span> : constexpr c10::complex<fT> operator*(const iT& a, const c10::complex<fT>& b) {</span>
|
|
<span id="L470"><span class="lineNum"> 470</span> : return static_cast<fT>(a) * b;</span>
|
|
<span id="L471"><span class="lineNum"> 471</span> : }</span>
|
|
<span id="L472"><span class="lineNum"> 472</span> : </span>
|
|
<span id="L473"><span class="lineNum"> 473</span> : template <typename fT, typename iT, COMPLEX_INTEGER_OP_TEMPLATE_CONDITION></span>
|
|
<span id="L474"><span class="lineNum"> 474</span> : constexpr c10::complex<fT> operator/(const c10::complex<fT>& a, const iT& b) {</span>
|
|
<span id="L475"><span class="lineNum"> 475</span> : return a / static_cast<fT>(b);</span>
|
|
<span id="L476"><span class="lineNum"> 476</span> : }</span>
|
|
<span id="L477"><span class="lineNum"> 477</span> : </span>
|
|
<span id="L478"><span class="lineNum"> 478</span> : template <typename fT, typename iT, COMPLEX_INTEGER_OP_TEMPLATE_CONDITION></span>
|
|
<span id="L479"><span class="lineNum"> 479</span> : constexpr c10::complex<fT> operator/(const iT& a, const c10::complex<fT>& b) {</span>
|
|
<span id="L480"><span class="lineNum"> 480</span> : return static_cast<fT>(a) / b;</span>
|
|
<span id="L481"><span class="lineNum"> 481</span> : }</span>
|
|
<span id="L482"><span class="lineNum"> 482</span> : </span>
|
|
<span id="L483"><span class="lineNum"> 483</span> : #undef COMPLEX_INTEGER_OP_TEMPLATE_CONDITION</span>
|
|
<span id="L484"><span class="lineNum"> 484</span> : </span>
|
|
<span id="L485"><span class="lineNum"> 485</span> : template <typename T></span>
|
|
<span id="L486"><span class="lineNum"> 486</span> : constexpr bool operator==(const complex<T>& lhs, const complex<T>& rhs) {</span>
|
|
<span id="L487"><span class="lineNum"> 487</span> : return (lhs.real() == rhs.real()) && (lhs.imag() == rhs.imag());</span>
|
|
<span id="L488"><span class="lineNum"> 488</span> : }</span>
|
|
<span id="L489"><span class="lineNum"> 489</span> : </span>
|
|
<span id="L490"><span class="lineNum"> 490</span> : template <typename T></span>
|
|
<span id="L491"><span class="lineNum"> 491</span> : constexpr bool operator==(const complex<T>& lhs, const T& rhs) {</span>
|
|
<span id="L492"><span class="lineNum"> 492</span> : return (lhs.real() == rhs) && (lhs.imag() == T());</span>
|
|
<span id="L493"><span class="lineNum"> 493</span> : }</span>
|
|
<span id="L494"><span class="lineNum"> 494</span> : </span>
|
|
<span id="L495"><span class="lineNum"> 495</span> : template <typename T></span>
|
|
<span id="L496"><span class="lineNum"> 496</span> : constexpr bool operator==(const T& lhs, const complex<T>& rhs) {</span>
|
|
<span id="L497"><span class="lineNum"> 497</span> : return (lhs == rhs.real()) && (T() == rhs.imag());</span>
|
|
<span id="L498"><span class="lineNum"> 498</span> : }</span>
|
|
<span id="L499"><span class="lineNum"> 499</span> : </span>
|
|
<span id="L500"><span class="lineNum"> 500</span> : template <typename T></span>
|
|
<span id="L501"><span class="lineNum"> 501</span> : constexpr bool operator!=(const complex<T>& lhs, const complex<T>& rhs) {</span>
|
|
<span id="L502"><span class="lineNum"> 502</span> : return !(lhs == rhs);</span>
|
|
<span id="L503"><span class="lineNum"> 503</span> : }</span>
|
|
<span id="L504"><span class="lineNum"> 504</span> : </span>
|
|
<span id="L505"><span class="lineNum"> 505</span> : template <typename T></span>
|
|
<span id="L506"><span class="lineNum"> 506</span> : constexpr bool operator!=(const complex<T>& lhs, const T& rhs) {</span>
|
|
<span id="L507"><span class="lineNum"> 507</span> : return !(lhs == rhs);</span>
|
|
<span id="L508"><span class="lineNum"> 508</span> : }</span>
|
|
<span id="L509"><span class="lineNum"> 509</span> : </span>
|
|
<span id="L510"><span class="lineNum"> 510</span> : template <typename T></span>
|
|
<span id="L511"><span class="lineNum"> 511</span> : constexpr bool operator!=(const T& lhs, const complex<T>& rhs) {</span>
|
|
<span id="L512"><span class="lineNum"> 512</span> : return !(lhs == rhs);</span>
|
|
<span id="L513"><span class="lineNum"> 513</span> : }</span>
|
|
<span id="L514"><span class="lineNum"> 514</span> : </span>
|
|
<span id="L515"><span class="lineNum"> 515</span> : template <typename T, typename CharT, typename Traits></span>
|
|
<span id="L516"><span class="lineNum"> 516</span> : std::basic_ostream<CharT, Traits>& operator<<(</span>
|
|
<span id="L517"><span class="lineNum"> 517</span> : std::basic_ostream<CharT, Traits>& os,</span>
|
|
<span id="L518"><span class="lineNum"> 518</span> : const complex<T>& x) {</span>
|
|
<span id="L519"><span class="lineNum"> 519</span> : return (os << static_cast<std::complex<T>>(x));</span>
|
|
<span id="L520"><span class="lineNum"> 520</span> : }</span>
|
|
<span id="L521"><span class="lineNum"> 521</span> : </span>
|
|
<span id="L522"><span class="lineNum"> 522</span> : template <typename T, typename CharT, typename Traits></span>
|
|
<span id="L523"><span class="lineNum"> 523</span> : std::basic_istream<CharT, Traits>& operator>>(</span>
|
|
<span id="L524"><span class="lineNum"> 524</span> : std::basic_istream<CharT, Traits>& is,</span>
|
|
<span id="L525"><span class="lineNum"> 525</span> : complex<T>& x) {</span>
|
|
<span id="L526"><span class="lineNum"> 526</span> : std::complex<T> tmp;</span>
|
|
<span id="L527"><span class="lineNum"> 527</span> : is >> tmp;</span>
|
|
<span id="L528"><span class="lineNum"> 528</span> : x = tmp;</span>
|
|
<span id="L529"><span class="lineNum"> 529</span> : return is;</span>
|
|
<span id="L530"><span class="lineNum"> 530</span> : }</span>
|
|
<span id="L531"><span class="lineNum"> 531</span> : </span>
|
|
<span id="L532"><span class="lineNum"> 532</span> : } // namespace c10</span>
|
|
<span id="L533"><span class="lineNum"> 533</span> : </span>
|
|
<span id="L534"><span class="lineNum"> 534</span> : // std functions</span>
|
|
<span id="L535"><span class="lineNum"> 535</span> : //</span>
|
|
<span id="L536"><span class="lineNum"> 536</span> : // The implementation of these functions also follow the design of C++20</span>
|
|
<span id="L537"><span class="lineNum"> 537</span> : </span>
|
|
<span id="L538"><span class="lineNum"> 538</span> : namespace std {</span>
|
|
<span id="L539"><span class="lineNum"> 539</span> : </span>
|
|
<span id="L540"><span class="lineNum"> 540</span> : template <typename T></span>
|
|
<span id="L541"><span class="lineNum"> 541</span> : constexpr T real(const c10::complex<T>& z) {</span>
|
|
<span id="L542"><span class="lineNum"> 542</span> : return z.real();</span>
|
|
<span id="L543"><span class="lineNum"> 543</span> : }</span>
|
|
<span id="L544"><span class="lineNum"> 544</span> : </span>
|
|
<span id="L545"><span class="lineNum"> 545</span> : template <typename T></span>
|
|
<span id="L546"><span class="lineNum"> 546</span> : constexpr T imag(const c10::complex<T>& z) {</span>
|
|
<span id="L547"><span class="lineNum"> 547</span> : return z.imag();</span>
|
|
<span id="L548"><span class="lineNum"> 548</span> : }</span>
|
|
<span id="L549"><span class="lineNum"> 549</span> : </span>
|
|
<span id="L550"><span class="lineNum"> 550</span> : template <typename T></span>
|
|
<span id="L551"><span class="lineNum"> 551</span> : C10_HOST_DEVICE T abs(const c10::complex<T>& z) {</span>
|
|
<span id="L552"><span class="lineNum"> 552</span> : #if defined(__CUDACC__) || defined(__HIPCC__)</span>
|
|
<span id="L553"><span class="lineNum"> 553</span> : return thrust::abs(static_cast<thrust::complex<T>>(z));</span>
|
|
<span id="L554"><span class="lineNum"> 554</span> : #else</span>
|
|
<span id="L555"><span class="lineNum"> 555</span> : return std::abs(static_cast<std::complex<T>>(z));</span>
|
|
<span id="L556"><span class="lineNum"> 556</span> : #endif</span>
|
|
<span id="L557"><span class="lineNum"> 557</span> : }</span>
|
|
<span id="L558"><span class="lineNum"> 558</span> : </span>
|
|
<span id="L559"><span class="lineNum"> 559</span> : #if defined(USE_ROCM)</span>
|
|
<span id="L560"><span class="lineNum"> 560</span> : #define ROCm_Bug(x)</span>
|
|
<span id="L561"><span class="lineNum"> 561</span> : #else</span>
|
|
<span id="L562"><span class="lineNum"> 562</span> : #define ROCm_Bug(x) x</span>
|
|
<span id="L563"><span class="lineNum"> 563</span> : #endif</span>
|
|
<span id="L564"><span class="lineNum"> 564</span> : </span>
|
|
<span id="L565"><span class="lineNum"> 565</span> : template <typename T></span>
|
|
<span id="L566"><span class="lineNum"> 566</span> : C10_HOST_DEVICE T arg(const c10::complex<T>& z) {</span>
|
|
<span id="L567"><span class="lineNum"> 567</span> : return ROCm_Bug(std)::atan2(std::imag(z), std::real(z));</span>
|
|
<span id="L568"><span class="lineNum"> 568</span> : }</span>
|
|
<span id="L569"><span class="lineNum"> 569</span> : </span>
|
|
<span id="L570"><span class="lineNum"> 570</span> : #undef ROCm_Bug</span>
|
|
<span id="L571"><span class="lineNum"> 571</span> : </span>
|
|
<span id="L572"><span class="lineNum"> 572</span> : template <typename T></span>
|
|
<span id="L573"><span class="lineNum"> 573</span> : constexpr T norm(const c10::complex<T>& z) {</span>
|
|
<span id="L574"><span class="lineNum"> 574</span> : return z.real() * z.real() + z.imag() * z.imag();</span>
|
|
<span id="L575"><span class="lineNum"> 575</span> : }</span>
|
|
<span id="L576"><span class="lineNum"> 576</span> : </span>
|
|
<span id="L577"><span class="lineNum"> 577</span> : // For std::conj, there are other versions of it:</span>
|
|
<span id="L578"><span class="lineNum"> 578</span> : // constexpr std::complex<float> conj( float z );</span>
|
|
<span id="L579"><span class="lineNum"> 579</span> : // template< class DoubleOrInteger ></span>
|
|
<span id="L580"><span class="lineNum"> 580</span> : // constexpr std::complex<double> conj( DoubleOrInteger z );</span>
|
|
<span id="L581"><span class="lineNum"> 581</span> : // constexpr std::complex<long double> conj( long double z );</span>
|
|
<span id="L582"><span class="lineNum"> 582</span> : // These are not implemented</span>
|
|
<span id="L583"><span class="lineNum"> 583</span> : // TODO(@zasdfgbnm): implement them as c10::conj</span>
|
|
<span id="L584"><span class="lineNum"> 584</span> : template <typename T></span>
|
|
<span id="L585"><span class="lineNum"> 585</span> : constexpr c10::complex<T> conj(const c10::complex<T>& z) {</span>
|
|
<span id="L586"><span class="lineNum"> 586</span> : return c10::complex<T>(z.real(), -z.imag());</span>
|
|
<span id="L587"><span class="lineNum"> 587</span> : }</span>
|
|
<span id="L588"><span class="lineNum"> 588</span> : </span>
|
|
<span id="L589"><span class="lineNum"> 589</span> : // Thrust does not have complex --> complex version of thrust::proj,</span>
|
|
<span id="L590"><span class="lineNum"> 590</span> : // so this function is not implemented at c10 right now.</span>
|
|
<span id="L591"><span class="lineNum"> 591</span> : // TODO(@zasdfgbnm): implement it by ourselves</span>
|
|
<span id="L592"><span class="lineNum"> 592</span> : </span>
|
|
<span id="L593"><span class="lineNum"> 593</span> : // There is no c10 version of std::polar, because std::polar always</span>
|
|
<span id="L594"><span class="lineNum"> 594</span> : // returns std::complex. Use c10::polar instead;</span>
|
|
<span id="L595"><span class="lineNum"> 595</span> : </span>
|
|
<span id="L596"><span class="lineNum"> 596</span> : } // namespace std</span>
|
|
<span id="L597"><span class="lineNum"> 597</span> : </span>
|
|
<span id="L598"><span class="lineNum"> 598</span> : namespace c10 {</span>
|
|
<span id="L599"><span class="lineNum"> 599</span> : </span>
|
|
<span id="L600"><span class="lineNum"> 600</span> : template <typename T></span>
|
|
<span id="L601"><span class="lineNum"> 601</span> : C10_HOST_DEVICE complex<T> polar(const T& r, const T& theta = T()) {</span>
|
|
<span id="L602"><span class="lineNum"> 602</span> : #if defined(__CUDACC__) || defined(__HIPCC__)</span>
|
|
<span id="L603"><span class="lineNum"> 603</span> : return static_cast<complex<T>>(thrust::polar(r, theta));</span>
|
|
<span id="L604"><span class="lineNum"> 604</span> : #else</span>
|
|
<span id="L605"><span class="lineNum"> 605</span> : // std::polar() requires r >= 0, so spell out the explicit implementation to</span>
|
|
<span id="L606"><span class="lineNum"> 606</span> : // avoid a branch.</span>
|
|
<span id="L607"><span class="lineNum"> 607</span> : return complex<T>(r * std::cos(theta), r * std::sin(theta));</span>
|
|
<span id="L608"><span class="lineNum"> 608</span> : #endif</span>
|
|
<span id="L609"><span class="lineNum"> 609</span> : }</span>
|
|
<span id="L610"><span class="lineNum"> 610</span> : </span>
|
|
<span id="L611"><span class="lineNum"> 611</span> : } // namespace c10</span>
|
|
<span id="L612"><span class="lineNum"> 612</span> : </span>
|
|
<span id="L613"><span class="lineNum"> 613</span> : C10_CLANG_DIAGNOSTIC_POP()</span>
|
|
<span id="L614"><span class="lineNum"> 614</span> : </span>
|
|
<span id="L615"><span class="lineNum"> 615</span> : #define C10_INTERNAL_INCLUDE_COMPLEX_REMAINING_H</span>
|
|
<span id="L616"><span class="lineNum"> 616</span> : // math functions are included in a separate file</span>
|
|
<span id="L617"><span class="lineNum"> 617</span> : #include <c10/util/complex_math.h> // IWYU pragma: keep</span>
|
|
<span id="L618"><span class="lineNum"> 618</span> : // utilities for complex types</span>
|
|
<span id="L619"><span class="lineNum"> 619</span> : #include <c10/util/complex_utils.h> // IWYU pragma: keep</span>
|
|
<span id="L620"><span class="lineNum"> 620</span> : #undef C10_INTERNAL_INCLUDE_COMPLEX_REMAINING_H</span>
|
|
</pre>
|
|
</td>
|
|
</tr>
|
|
</table>
|
|
<br>
|
|
|
|
<table width="100%" border=0 cellspacing=0 cellpadding=0>
|
|
<tr><td class="ruler"><img src="../../../../glass.png" width=3 height=3 alt=""></td></tr>
|
|
<tr><td class="versionInfo">Generated by: <a href="https://github.com//linux-test-project/lcov" target="_parent">LCOV version 2.0-1</a></td></tr>
|
|
</table>
|
|
<br>
|
|
|
|
</body>
|
|
</html>
|