moflib-2.0/doc/html/basic__matrix_8hpp_source.html

   1 <!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd">
   2 <html xmlns="http://www.w3.org/1999/xhtml">
   3 <head>
   4 <meta http-equiv="Content-Type" content="text/xhtml;charset=UTF-8"/>
   5 <title>MotifOfFunctionsLIBrary: /home/yasutomo/devel/moflib/src/mof/math/basic_matrix.hpp ソースファイル</title>
   6 <link href="tabs.css" rel="stylesheet" type="text/css"/>
   7 <link href="doxygen.css" rel="stylesheet" type="text/css"/>
   8 </head>
   9 <body>
  10 <!-- 作成： Doxygen 1.6.3 -->
  11 <div class="navigation" id="top">
  12   <div class="tabs">
  13     <ul>
  14       <li><a href="index.html"><span>メインページ</span></a></li>
  15       <li><a href="namespaces.html"><span>ネームスペース</span></a></li>
  16       <li><a href="annotated.html"><span>クラス</span></a></li>
  17       <li class="current"><a href="files.html"><span>ファイル</span></a></li>
  18     </ul>
  19   </div>
  20   <div class="tabs">
  21     <ul>
  22       <li><a href="files.html"><span>ファイル一覧</span></a></li>
  23       <li><a href="globals.html"><span>ファイルメンバ</span></a></li>
  24     </ul>
  25   </div>
  26 <h1>/home/yasutomo/devel/moflib/src/mof/math/basic_matrix.hpp</h1><a href="basic__matrix_8hpp.html">説明を見る。</a><div class="fragment"><pre class="fragment"><a name="l00001"></a>00001 <span class="preprocessor">#pragma once</span>
  27 <a name="l00002"></a>00002 <span class="preprocessor"></span><span class="preprocessor">#include &lt;<a class="code" href="mofdef_8hpp.html">mof/base/mofdef.hpp</a>&gt;</span>
  28 <a name="l00003"></a>00003 <span class="preprocessor">#include &lt;<a class="code" href="threshold_8hpp.html">mof/math/threshold.hpp</a>&gt;</span>
  29 <a name="l00004"></a>00004 <span class="preprocessor">#include &lt;boost/operators.hpp&gt;</span>
  30 <a name="l00005"></a>00005 <span class="preprocessor">#include &lt;ostream&gt;</span>
  31 <a name="l00006"></a>00006 <span class="preprocessor">#include &lt;iomanip&gt;</span>
  32 <a name="l00007"></a>00007 <span class="preprocessor">#include &lt;cmath&gt;</span>
  33 <a name="l00008"></a>00008
  34 <a name="l00009"></a>00009 <span class="keyword">namespace </span>mof
  35 <a name="l00010"></a>00010 {
  36 <a name="l00011"></a><a class="code" href="namespacemof_1_1math.html">00011</a> <span class="keyword">namespace </span>math
  37 <a name="l00012"></a>00012 {
  38 <a name="l00013"></a>00013
  39 <a name="l00014"></a>00014     <span class="keyword">template</span> &lt;<span class="keywordtype">size_t</span> Dim&gt; <span class="keyword">struct </span><a class="code" href="structmof_1_1math_1_1row__of__matrix.html" title="M[i][j]のように行列の要素を参照できるようにするための補助クラス...">row_of_matrix</a>;
  40 <a name="l00015"></a>00015
  41 <a name="l00025"></a>00025     <span class="keyword">template</span> &lt;<span class="keywordtype">size_t</span> Dim, <span class="keyword">typename</span> Derived, <span class="keyword">typename</span> Coordinate&gt;
  42 <a name="l00026"></a><a class="code" href="classmof_1_1math_1_1basic__matrix.html">00026</a>     <span class="keyword">class </span><a class="code" href="classmof_1_1math_1_1basic__matrix.html" title="同次座標変換行列テンプレートクラス">basic_matrix</a>
  43 <a name="l00027"></a>00027         : boost::<a class="code" href="classboost_1_1addable.html">addable</a>&lt; Derived
  44 <a name="l00028"></a>00028         , boost::addable2&lt; Derived, float
  45 <a name="l00029"></a>00029         , boost::subtractable&lt; Derived
  46 <a name="l00030"></a>00030         , boost::subtractable2&lt; Derived, float
  47 <a name="l00031"></a>00031         , boost::multipliable&lt; Derived
  48 <a name="l00032"></a>00032         , boost::multipliable2&lt; Derived, float
  49 <a name="l00033"></a>00033         , boost::dividable2&lt; Derived, float
  50 <a name="l00034"></a>00034         , boost::equality_comparable&lt; Derived
  51 <a name="l00035"></a>00035         &gt; &gt; &gt; &gt; &gt; &gt; &gt; &gt;
  52 <a name="l00036"></a>00036     {
  53 <a name="l00037"></a>00037     <span class="keyword">protected</span>:
  54 <a name="l00038"></a>00038 <span class="comment">//{{{ last_index</span>
  55 <a name="l00042"></a><a class="code" href="classmof_1_1math_1_1basic__matrix.html#a31d7ff86afe909d714d66c8a697d4486">00042</a> <span class="comment"></span>        <span class="keywordtype">size_t</span> <a class="code" href="classmof_1_1math_1_1basic__matrix.html#a31d7ff86afe909d714d66c8a697d4486" title="elements_の最後の添字を得る">last_index</a>()<span class="keyword"> const</span>
  56 <a name="l00043"></a>00043 <span class="keyword">        </span>{
  57 <a name="l00044"></a>00044             <span class="keywordflow">return</span> (Dim + 1) * (Dim + 1) - 1;
  58 <a name="l00045"></a>00045         }
  59 <a name="l00046"></a>00046 <span class="comment">//}}}</span>
  60 <a name="l00047"></a><a class="code" href="classmof_1_1math_1_1basic__matrix.html#a31738810fd6ff013ddcb0dd3d78353f0">00047</a>         <span class="keywordtype">float</span> <a class="code" href="classmof_1_1math_1_1basic__matrix.html#a31738810fd6ff013ddcb0dd3d78353f0" title="要素の配列">elements_</a>[(Dim + 1) * (Dim + 1)];
  61 <a name="l00048"></a>00048     <span class="keyword">public</span>:
  62 <a name="l00049"></a>00049         <span class="comment">// コンストラクタ，デストラクタはデフォルトのものを使う</span>
  63 <a name="l00050"></a>00050         <span class="comment">// 代入演算子，コピーコンストラクタはデフォルトのものを使う</span>
  64 <a name="l00051"></a>00051 <span class="comment">//{{{ operator +=</span>
  65 <a name="l00052"></a><a class="code" href="classmof_1_1math_1_1basic__matrix.html#ade7472764844aa239cd81bc6576d3d01">00052</a>         Derived&amp; <a class="code" href="classmof_1_1math_1_1basic__matrix.html#ade7472764844aa239cd81bc6576d3d01">operator+=</a>(<span class="keyword">const</span> Derived&amp; rhs)
  66 <a name="l00053"></a>00053         {
  67 <a name="l00054"></a>00054             <span class="keyword">const</span> <span class="keywordtype">size_t</span> END = <a class="code" href="classmof_1_1math_1_1basic__matrix.html#a31d7ff86afe909d714d66c8a697d4486" title="elements_の最後の添字を得る">last_index</a>() - 1;<span class="comment">// 最後の要素は1で保存</span>
  68 <a name="l00055"></a>00055             <span class="keywordflow">for</span> (<span class="keywordtype">size_t</span> i = 0; i &lt;= END; ++i) {
  69 <a name="l00056"></a>00056                 <a class="code" href="classmof_1_1math_1_1basic__matrix.html#a31738810fd6ff013ddcb0dd3d78353f0" title="要素の配列">elements_</a>[i] += rhs.elements_[i];
  70 <a name="l00057"></a>00057             }
  71 <a name="l00058"></a>00058             <span class="keywordflow">return</span> *<span class="keyword">reinterpret_cast&lt;</span>Derived*<span class="keyword">&gt;</span>(<span class="keyword">this</span>);<span class="comment">//thisがDerived型であることは保証されている．</span>
  72 <a name="l00059"></a>00059         }
  73 <a name="l00060"></a>00060
  74 <a name="l00061"></a><a class="code" href="classmof_1_1math_1_1basic__matrix.html#a57137691ca3efba9e3c59be45ec112a6">00061</a>         Derived&amp; <a class="code" href="classmof_1_1math_1_1basic__matrix.html#ade7472764844aa239cd81bc6576d3d01">operator+=</a>(<span class="keywordtype">float</span> rhs)
  75 <a name="l00062"></a>00062         {
  76 <a name="l00063"></a>00063             <span class="keyword">const</span> <span class="keywordtype">size_t</span> END = <a class="code" href="classmof_1_1math_1_1basic__matrix.html#a31d7ff86afe909d714d66c8a697d4486" title="elements_の最後の添字を得る">last_index</a>() - 1;<span class="comment">// 最後の要素は1で保存</span>
  77 <a name="l00064"></a>00064             <span class="keywordflow">for</span> (<span class="keywordtype">size_t</span> i = 0; i &lt;= END; ++i) {
  78 <a name="l00065"></a>00065                 <a class="code" href="classmof_1_1math_1_1basic__matrix.html#a31738810fd6ff013ddcb0dd3d78353f0" title="要素の配列">elements_</a>[i] += rhs;
  79 <a name="l00066"></a>00066             }
  80 <a name="l00067"></a>00067             <span class="keywordflow">return</span> *<span class="keyword">reinterpret_cast&lt;</span>Derived*<span class="keyword">&gt;</span>(<span class="keyword">this</span>);<span class="comment">//thisがDerived型であることは保証されている．</span>
  81 <a name="l00068"></a>00068         }
  82 <a name="l00069"></a>00069
  83 <a name="l00070"></a><a class="code" href="classmof_1_1math_1_1basic__matrix.html#ab51200fa1e6bb7c9a1cfad79954049f6">00070</a>         <span class="keyword">friend</span> Derived <a class="code" href="classmof_1_1math_1_1basic__matrix.html#ab51200fa1e6bb7c9a1cfad79954049f6">operator+</a>(<span class="keywordtype">float</span> rhs1, Derived&amp; rhs2)
  84 <a name="l00071"></a>00071         {
  85 <a name="l00072"></a>00072             Derived retval;
  86 <a name="l00073"></a>00073             <span class="keyword">const</span> <span class="keywordtype">size_t</span> END = <a class="code" href="classmof_1_1math_1_1basic__matrix.html#a31d7ff86afe909d714d66c8a697d4486" title="elements_の最後の添字を得る">last_index</a>() - 1;<span class="comment">// 最後の要素は1で保存</span>
  87 <a name="l00074"></a>00074             <span class="keywordflow">for</span> (<span class="keywordtype">size_t</span> i = 0; i &lt;= END; ++i) {
  88 <a name="l00075"></a>00075                 retval.elements_[i] = rhs1 + rhs2.elements_[i];
  89 <a name="l00076"></a>00076             }
  90 <a name="l00077"></a>00077             <span class="keywordflow">return</span> retval;
  91 <a name="l00078"></a>00078         }
  92 <a name="l00079"></a>00079 <span class="comment">//}}}</span>
  93 <a name="l00080"></a>00080 <span class="comment">//{{{ operator -=</span>
  94 <a name="l00081"></a><a class="code" href="classmof_1_1math_1_1basic__matrix.html#a55c53c6c312208a6f7c716542b2dcf7e">00081</a>         Derived&amp; <a class="code" href="classmof_1_1math_1_1basic__matrix.html#a55c53c6c312208a6f7c716542b2dcf7e">operator-=</a>(<span class="keyword">const</span> Derived&amp; rhs)
  95 <a name="l00082"></a>00082         {
  96 <a name="l00083"></a>00083             <span class="keyword">const</span> <span class="keywordtype">size_t</span> END = <a class="code" href="classmof_1_1math_1_1basic__matrix.html#a31d7ff86afe909d714d66c8a697d4486" title="elements_の最後の添字を得る">last_index</a>() - 1;<span class="comment">// 最後の要素は1で保存</span>
  97 <a name="l00084"></a>00084             <span class="keywordflow">for</span> (<span class="keywordtype">size_t</span> i = 0; i &lt;= END; ++i) {
  98 <a name="l00085"></a>00085                 <a class="code" href="classmof_1_1math_1_1basic__matrix.html#a31738810fd6ff013ddcb0dd3d78353f0" title="要素の配列">elements_</a>[i] -= rhs.elements_[i];
  99 <a name="l00086"></a>00086             }
 100 <a name="l00087"></a>00087             <span class="keywordflow">return</span> *<span class="keyword">reinterpret_cast&lt;</span>Derived*<span class="keyword">&gt;</span>(<span class="keyword">this</span>);<span class="comment">//thisがDerived型であることは保証されている．</span>
 101 <a name="l00088"></a>00088         }
 102 <a name="l00089"></a>00089
 103 <a name="l00090"></a><a class="code" href="classmof_1_1math_1_1basic__matrix.html#a13104fef1f1ba8f89b885d32110ae946">00090</a>         Derived&amp; <a class="code" href="classmof_1_1math_1_1basic__matrix.html#a55c53c6c312208a6f7c716542b2dcf7e">operator-=</a>(<span class="keywordtype">float</span> rhs)
 104 <a name="l00091"></a>00091         {
 105 <a name="l00092"></a>00092             <span class="keyword">const</span> <span class="keywordtype">size_t</span> END = <a class="code" href="classmof_1_1math_1_1basic__matrix.html#a31d7ff86afe909d714d66c8a697d4486" title="elements_の最後の添字を得る">last_index</a>() - 1;<span class="comment">// 最後の要素は1で保存</span>
 106 <a name="l00093"></a>00093             <span class="keywordflow">for</span> (<span class="keywordtype">size_t</span> i = 0; i &lt;= END; ++i) {
 107 <a name="l00094"></a>00094                 <a class="code" href="classmof_1_1math_1_1basic__matrix.html#a31738810fd6ff013ddcb0dd3d78353f0" title="要素の配列">elements_</a>[i] -= rhs;
 108 <a name="l00095"></a>00095             }
 109 <a name="l00096"></a>00096             <span class="keywordflow">return</span> *<span class="keyword">reinterpret_cast&lt;</span>Derived*<span class="keyword">&gt;</span>(<span class="keyword">this</span>);<span class="comment">//thisがDerived型であることは保証されている．</span>
 110 <a name="l00097"></a>00097         }
 111 <a name="l00098"></a>00098
 112 <a name="l00099"></a><a class="code" href="classmof_1_1math_1_1basic__matrix.html#a87b3068df3cfd1cb062ab11d6700178a">00099</a>         <span class="keyword">friend</span> Derived <a class="code" href="classmof_1_1math_1_1basic__matrix.html#a87b3068df3cfd1cb062ab11d6700178a">operator-</a>(<span class="keywordtype">float</span> rhs1, Derived&amp; rhs2)
 113 <a name="l00100"></a>00100         {
 114 <a name="l00101"></a>00101             Derived retval;
 115 <a name="l00102"></a>00102             <span class="keyword">const</span> <span class="keywordtype">size_t</span> END = <a class="code" href="classmof_1_1math_1_1basic__matrix.html#a31d7ff86afe909d714d66c8a697d4486" title="elements_の最後の添字を得る">last_index</a>() - 1;<span class="comment">// 最後の要素は1で保存</span>
 116 <a name="l00103"></a>00103             <span class="keywordflow">for</span> (<span class="keywordtype">size_t</span> i = 0; i &lt;= END; ++i) {
 117 <a name="l00104"></a>00104                 retval.elements_[i] = rhs1 - rhs2.elements_[i];
 118 <a name="l00105"></a>00105             }
 119 <a name="l00106"></a>00106             <span class="keywordflow">return</span> retval;
 120 <a name="l00107"></a>00107         }
 121 <a name="l00108"></a>00108 <span class="comment">//}}}</span>
 122 <a name="l00109"></a>00109 <span class="comment">//{{{ operator *=   </span>
 123 <a name="l00113"></a><a class="code" href="classmof_1_1math_1_1basic__matrix.html#ab9a8a7be1b789c318ea52eeaa1d2a55c">00113</a> <span class="comment"></span>        Derived&amp; <a class="code" href="classmof_1_1math_1_1basic__matrix.html#ab9a8a7be1b789c318ea52eeaa1d2a55c" title="行列の積を計算し，最後の要素が1になるように定数倍する．...">operator*=</a>(<span class="keyword">const</span> Derived&amp; rhs)
 124 <a name="l00114"></a>00114         {
 125 <a name="l00115"></a>00115             Derived retval;
 126 <a name="l00116"></a>00116             <span class="keyword">const</span> <span class="keywordtype">int</span> SIZE = Dim + 1;
 127 <a name="l00117"></a>00117
 128 <a name="l00118"></a>00118             <span class="comment">// calculate the last element previously</span>
 129 <a name="l00119"></a>00119             <span class="keywordtype">int</span> b = <a class="code" href="classmof_1_1math_1_1basic__matrix.html#a31d7ff86afe909d714d66c8a697d4486" title="elements_の最後の添字を得る">last_index</a>() - Dim;
 130 <a name="l00120"></a>00120             <span class="keywordtype">int</span> c = Dim;
 131 <a name="l00121"></a>00121             <span class="keywordtype">float</span> last_sum = 0;
 132 <a name="l00122"></a>00122             <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i = 0; i &lt; SIZE; ++i) {
 133 <a name="l00123"></a>00123                 last_sum += <a class="code" href="classmof_1_1math_1_1basic__matrix.html#a31738810fd6ff013ddcb0dd3d78353f0" title="要素の配列">elements_</a>[b + i] * rhs.elements_[c + i * SIZE];
 134 <a name="l00124"></a>00124             }
 135 <a name="l00125"></a>00125             retval.elements_[<a class="code" href="classmof_1_1math_1_1basic__matrix.html#a31d7ff86afe909d714d66c8a697d4486" title="elements_の最後の添字を得る">last_index</a>()] = 1;
 136 <a name="l00126"></a>00126
 137 <a name="l00127"></a>00127             <span class="keywordflow">for</span> (<span class="keywordtype">int</span> a = <a class="code" href="classmof_1_1math_1_1basic__matrix.html#a31d7ff86afe909d714d66c8a697d4486" title="elements_の最後の添字を得る">last_index</a>() - 1; a &gt;= 0; --a) {
 138 <a name="l00128"></a>00128                 <span class="keywordtype">int</span> b = a / SIZE * SIZE;
 139 <a name="l00129"></a>00129                 <span class="keywordtype">int</span> c = a % SIZE;
 140 <a name="l00130"></a>00130                 <span class="keywordtype">float</span> sum = 0;
 141 <a name="l00131"></a>00131                 <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i = 0; i &lt; SIZE; ++i) {
 142 <a name="l00132"></a>00132                     sum += <a class="code" href="classmof_1_1math_1_1basic__matrix.html#a31738810fd6ff013ddcb0dd3d78353f0" title="要素の配列">elements_</a>[b + i] * rhs.elements_[c + i * SIZE];
 143 <a name="l00133"></a>00133                 }
 144 <a name="l00134"></a>00134                 retval.elements_[a] = sum / last_sum;
 145 <a name="l00135"></a>00135             }
 146 <a name="l00136"></a>00136             *<span class="keyword">this</span> = retval;
 147 <a name="l00137"></a>00137             <span class="keywordflow">return</span> *<span class="keyword">reinterpret_cast&lt;</span>Derived*<span class="keyword">&gt;</span>(<span class="keyword">this</span>);<span class="comment">//thisがDerived型であることは保証されている．</span>
 148 <a name="l00138"></a>00138         }
 149 <a name="l00139"></a>00139
 150 <a name="l00140"></a><a class="code" href="classmof_1_1math_1_1basic__matrix.html#ab3aaa6779c7a24ea87a83be329969da6">00140</a>         Derived&amp; <a class="code" href="classmof_1_1math_1_1basic__matrix.html#ab9a8a7be1b789c318ea52eeaa1d2a55c" title="行列の積を計算し，最後の要素が1になるように定数倍する．...">operator*=</a>(<span class="keywordtype">float</span> rhs)
 151 <a name="l00141"></a>00141         {
 152 <a name="l00142"></a>00142             <span class="keyword">const</span> <span class="keywordtype">size_t</span> END = <a class="code" href="classmof_1_1math_1_1basic__matrix.html#a31d7ff86afe909d714d66c8a697d4486" title="elements_の最後の添字を得る">last_index</a>() - 1;<span class="comment">// 最後の要素は1で保存</span>
 153 <a name="l00143"></a>00143             <span class="keywordflow">for</span> (<span class="keywordtype">size_t</span> i = 0; i &lt;= END; ++i) {
 154 <a name="l00144"></a>00144                 <a class="code" href="classmof_1_1math_1_1basic__matrix.html#a31738810fd6ff013ddcb0dd3d78353f0" title="要素の配列">elements_</a>[i] *= rhs;
 155 <a name="l00145"></a>00145             }
 156 <a name="l00146"></a>00146             <span class="keywordflow">return</span> *<span class="keyword">reinterpret_cast&lt;</span>Derived*<span class="keyword">&gt;</span>(<span class="keyword">this</span>);<span class="comment">//thisがDerived型であることは保証されている．</span>
 157 <a name="l00147"></a>00147         }
 158 <a name="l00148"></a>00148
 159 <a name="l00152"></a><a class="code" href="classmof_1_1math_1_1basic__matrix.html#a66e89f8a7eb8da98314b4218278bfcbf">00152</a>         Coordinate <a class="code" href="classmof_1_1math_1_1basic__matrix.html#a66e89f8a7eb8da98314b4218278bfcbf" title="ベクトルを変換し，最後の要素が1になるように定数倍する．...">operator*</a>(<span class="keyword">const</span> Coordinate&amp; rhs)
 160 <a name="l00153"></a>00153         {
 161 <a name="l00154"></a>00154             <span class="keyword">const</span> <span class="keywordtype">size_t</span> SIZE = Dim + 1;
 162 <a name="l00155"></a>00155
 163 <a name="l00156"></a>00156             <span class="keywordtype">float</span> last_sum = 0;
 164 <a name="l00157"></a>00157             <span class="keywordflow">for</span> (<span class="keywordtype">size_t</span> i = 0; i &lt; SIZE; ++i) {
 165 <a name="l00158"></a>00158                 last_sum += <a class="code" href="classmof_1_1math_1_1basic__matrix.html#a31738810fd6ff013ddcb0dd3d78353f0" title="要素の配列">elements_</a>[Dim * SIZE + i] * rhs[i];
 166 <a name="l00159"></a>00159             }
 167 <a name="l00160"></a>00160
 168 <a name="l00161"></a>00161             <span class="keywordtype">float</span> arr[Dim];
 169 <a name="l00162"></a>00162             <span class="keywordflow">for</span> (<span class="keywordtype">size_t</span> a = 0; a &lt; Dim; ++a) {
 170 <a name="l00163"></a>00163                 <span class="keywordtype">float</span> sum = 0;
 171 <a name="l00164"></a>00164                 <span class="keywordflow">for</span> (<span class="keywordtype">size_t</span> i = 0; i &lt; SIZE; ++i) {
 172 <a name="l00165"></a>00165                     sum += <a class="code" href="classmof_1_1math_1_1basic__matrix.html#a31738810fd6ff013ddcb0dd3d78353f0" title="要素の配列">elements_</a>[a * SIZE + i] * rhs[i];
 173 <a name="l00166"></a>00166                 }
 174 <a name="l00167"></a>00167                 arr[a] = sum / last_sum;
 175 <a name="l00168"></a>00168             }
 176 <a name="l00169"></a>00169             <span class="keywordflow">return</span> Coordinate(arr);
 177 <a name="l00170"></a>00170         }
 178 <a name="l00171"></a>00171
 179 <a name="l00172"></a>00172
 180 <a name="l00173"></a><a class="code" href="classmof_1_1math_1_1basic__matrix.html#a83780bceadc61ec021a6d0ceae3da02a">00173</a>         <span class="keyword">friend</span> Derived <a class="code" href="classmof_1_1math_1_1basic__matrix.html#a66e89f8a7eb8da98314b4218278bfcbf" title="ベクトルを変換し，最後の要素が1になるように定数倍する．...">operator*</a>(<span class="keywordtype">float</span> rhs1, Derived&amp; rhs2)
 181 <a name="l00174"></a>00174         {
 182 <a name="l00175"></a>00175             Derived retval;
 183 <a name="l00176"></a>00176             <span class="keyword">const</span> <span class="keywordtype">size_t</span> END = rhs2.last_index() - 1;<span class="comment">// 最後の要素は1で保存</span>
 184 <a name="l00177"></a>00177             <span class="keywordflow">for</span> (<span class="keywordtype">size_t</span> i = 0; i &lt;= END; ++i) {
 185 <a name="l00178"></a>00178                 retval.elements_[i] = rhs1 * rhs2.elements_[i];
 186 <a name="l00179"></a>00179             }
 187 <a name="l00180"></a>00180             <span class="keywordflow">return</span> retval;
 188 <a name="l00181"></a>00181         }
 189 <a name="l00182"></a>00182 <span class="comment">//}}}</span>
 190 <a name="l00183"></a>00183 <span class="comment">//{{{ operator /=</span>
 191 <a name="l00184"></a><a class="code" href="classmof_1_1math_1_1basic__matrix.html#adcd8c94f420be936816d06ba99b6111c">00184</a>         Derived&amp; <a class="code" href="classmof_1_1math_1_1basic__matrix.html#adcd8c94f420be936816d06ba99b6111c">operator/=</a>(<span class="keywordtype">float</span> rhs)
 192 <a name="l00185"></a>00185         {
 193 <a name="l00186"></a>00186             <span class="keyword">const</span> <span class="keywordtype">size_t</span> END = <a class="code" href="classmof_1_1math_1_1basic__matrix.html#a31d7ff86afe909d714d66c8a697d4486" title="elements_の最後の添字を得る">last_index</a>() - 1;<span class="comment">// 最後の要素は1で保存</span>
 194 <a name="l00187"></a>00187             <span class="keywordflow">for</span> (<span class="keywordtype">size_t</span> i = 0; i &lt;= END; ++i) {
 195 <a name="l00188"></a>00188                 <a class="code" href="classmof_1_1math_1_1basic__matrix.html#a31738810fd6ff013ddcb0dd3d78353f0" title="要素の配列">elements_</a>[i] /= rhs;
 196 <a name="l00189"></a>00189             }
 197 <a name="l00190"></a>00190             <span class="keywordflow">return</span> *<span class="keyword">reinterpret_cast&lt;</span>Derived*<span class="keyword">&gt;</span>(<span class="keyword">this</span>);<span class="comment">//thisがDerived型であることは保証されている．</span>
 198 <a name="l00191"></a>00191         }
 199 <a name="l00192"></a>00192 <span class="comment">//}}}</span>
 200 <a name="l00193"></a>00193 <span class="comment">//{{{ operator ==</span>
 201 <a name="l00194"></a><a class="code" href="classmof_1_1math_1_1basic__matrix.html#a25ee023b18ea540e2d1227189d1bb898">00194</a>         <span class="keywordtype">bool</span> <a class="code" href="classmof_1_1math_1_1basic__matrix.html#a25ee023b18ea540e2d1227189d1bb898">operator==</a>(<span class="keyword">const</span> Derived&amp; rhs)<span class="keyword"> const</span>
 202 <a name="l00195"></a>00195 <span class="keyword">        </span>{
 203 <a name="l00196"></a>00196             <span class="keyword">const</span> <span class="keywordtype">size_t</span> END = <a class="code" href="classmof_1_1math_1_1basic__matrix.html#a31d7ff86afe909d714d66c8a697d4486" title="elements_の最後の添字を得る">last_index</a>() - 1;<span class="comment">// 最後の要素は1で保存されている</span>
 204 <a name="l00197"></a>00197             <span class="keywordflow">for</span> (<span class="keywordtype">size_t</span> i = 0; i &lt;= END; ++i) {
 205 <a name="l00198"></a>00198                 <span class="keywordflow">if</span> (std::abs(<a class="code" href="classmof_1_1math_1_1basic__matrix.html#a31738810fd6ff013ddcb0dd3d78353f0" title="要素の配列">elements_</a>[i] - rhs.elements_[i]) &gt; <a class="code" href="threshold_8hpp.html#acefd3ede4bb1937da93c49823d4d008e">MOF_ERROR_THRESHOLD</a>) <span class="keywordflow">return</span> <span class="keyword">false</span>;
 206 <a name="l00199"></a>00199             }
 207 <a name="l00200"></a>00200             <span class="keywordflow">return</span> <span class="keyword">true</span>;
 208 <a name="l00201"></a>00201         }
 209 <a name="l00202"></a>00202 <span class="comment">//}}}</span>
 210 <a name="l00203"></a>00203 <span class="comment">//{{{ operator []</span>
 211 <a name="l00208"></a><a class="code" href="classmof_1_1math_1_1basic__matrix.html#adf2ab1f146d9ff006522a437ab304617">00208</a> <span class="comment"></span>        <a class="code" href="structmof_1_1math_1_1row__of__matrix.html" title="M[i][j]のように行列の要素を参照できるようにするための補助クラス...">row_of_matrix&lt;Dim&gt;</a> <a class="code" href="classmof_1_1math_1_1basic__matrix.html#adf2ab1f146d9ff006522a437ab304617">operator []</a>(<span class="keywordtype">size_t</span> i)<span class="keyword"> const</span>
 212 <a name="l00209"></a>00209 <span class="keyword">        </span>{
 213 <a name="l00210"></a>00210             <span class="keywordflow">return</span> <a class="code" href="structmof_1_1math_1_1row__of__matrix.html" title="M[i][j]のように行列の要素を参照できるようにするための補助クラス...">row_of_matrix&lt;Dim&gt;</a>(<a class="code" href="classmof_1_1math_1_1basic__matrix.html#a31738810fd6ff013ddcb0dd3d78353f0" title="要素の配列">elements_</a>, i);
 214 <a name="l00211"></a>00211         }
 215 <a name="l00212"></a>00212 <span class="comment">//}}}</span>
 216 <a name="l00213"></a>00213 <span class="comment">//{{{ operator &lt;&lt;</span>
 217 <a name="l00220"></a>00220 <span class="comment"></span>        <span class="keyword">friend</span> std::ostream&amp; <span class="keyword">operator</span>&lt;&lt;
 218 <a name="l00221"></a><a class="code" href="classmof_1_1math_1_1basic__matrix.html#a46d2540a8abf6ca59db3ee1551769ab9">00221</a>         (
 219 <a name="l00222"></a>00222             std::ostream&amp; stream,
 220 <a name="l00223"></a>00223             <span class="keyword">const</span> Derived&amp; rhs
 221 <a name="l00224"></a>00224         )
 222 <a name="l00225"></a>00225         {
 223 <a name="l00226"></a>00226             <span class="keywordflow">for</span> (<span class="keywordtype">size_t</span> i = 0; i &lt; Dim + 1; ++i) {
 224 <a name="l00227"></a>00227                 <span class="keywordflow">for</span> (<span class="keywordtype">size_t</span> j = 0; j &lt; Dim + 1; ++j) {
 225 <a name="l00228"></a>00228                     <span class="keywordflow">if</span> (j != 0) stream &lt;&lt; <span class="stringliteral">&quot;, &quot;</span>;
 226 <a name="l00229"></a>00229                     <span class="keywordflow">else</span> <span class="keywordflow">if</span> (j != Dim) stream &lt;&lt; <span class="stringliteral">&quot;\n&quot;</span>;
 227 <a name="l00230"></a>00230                     stream &lt;&lt; std::setw(5) &lt;&lt; rhs.elements_[i * (Dim + 1) + j];
 228 <a name="l00231"></a>00231                 }
 229 <a name="l00232"></a>00232             }
 230 <a name="l00233"></a>00233             <span class="keywordflow">return</span> stream;
 231 <a name="l00234"></a>00234         }
 232 <a name="l00235"></a>00235 <span class="comment">//}}}</span>
 233 <a name="l00236"></a>00236 <span class="comment">//{{{ at</span>
 234 <a name="l00243"></a><a class="code" href="classmof_1_1math_1_1basic__matrix.html#a02dbe242ff1395bc2577c67573433efd">00243</a> <span class="comment"></span>        <span class="keywordtype">float</span> <a class="code" href="classmof_1_1math_1_1basic__matrix.html#a02dbe242ff1395bc2577c67573433efd" title="行列の要素にアクセス">at</a>(<span class="keywordtype">size_t</span> i, <span class="keywordtype">size_t</span> j)<span class="keyword"> const</span>
 235 <a name="l00244"></a>00244 <span class="keyword">        </span>{
 236 <a name="l00245"></a>00245             <span class="keywordflow">return</span> <a class="code" href="classmof_1_1math_1_1basic__matrix.html#a31738810fd6ff013ddcb0dd3d78353f0" title="要素の配列">elements_</a>[i * (Dim + 1) + j];
 237 <a name="l00246"></a>00246         }
 238 <a name="l00247"></a>00247 <span class="comment">//}}}</span>
 239 <a name="l00248"></a>00248     };
 240 <a name="l00249"></a>00249 <span class="comment">//{{{ row_of_matrix</span>
 241 <a name="l00253"></a>00253 <span class="comment"></span>    <span class="keyword">template</span> &lt;<span class="keywordtype">size_t</span> Dim&gt;
 242 <a name="l00254"></a><a class="code" href="structmof_1_1math_1_1row__of__matrix.html">00254</a>     <span class="keyword">struct </span><a class="code" href="structmof_1_1math_1_1row__of__matrix.html" title="M[i][j]のように行列の要素を参照できるようにするための補助クラス...">row_of_matrix</a>
 243 <a name="l00255"></a>00255     {
 244 <a name="l00256"></a><a class="code" href="structmof_1_1math_1_1row__of__matrix.html#a0bdc04595281825f274b1016728f6a72">00256</a>         <span class="keyword">const</span> <span class="keywordtype">float</span>* <a class="code" href="structmof_1_1math_1_1row__of__matrix.html#a0bdc04595281825f274b1016728f6a72">elements_</a>;
 245 <a name="l00257"></a><a class="code" href="structmof_1_1math_1_1row__of__matrix.html#af1d3a14a09bb4bbe6ddb5d8e37730b9e">00257</a>         <span class="keywordtype">size_t</span> <a class="code" href="structmof_1_1math_1_1row__of__matrix.html#af1d3a14a09bb4bbe6ddb5d8e37730b9e">row_index_</a>;
 246 <a name="l00258"></a>00258
 247 <a name="l00259"></a><a class="code" href="structmof_1_1math_1_1row__of__matrix.html#add819c408998bb73fd68da047504684c">00259</a>         <a class="code" href="structmof_1_1math_1_1row__of__matrix.html#add819c408998bb73fd68da047504684c">row_of_matrix</a>(<span class="keyword">const</span> <span class="keywordtype">float</span>* elements, <span class="keywordtype">size_t</span> row_index)
 248 <a name="l00260"></a>00260             : <a class="code" href="structmof_1_1math_1_1row__of__matrix.html#a0bdc04595281825f274b1016728f6a72">elements_</a>(elements), <a class="code" href="structmof_1_1math_1_1row__of__matrix.html#af1d3a14a09bb4bbe6ddb5d8e37730b9e">row_index_</a>(row_index)
 249 <a name="l00261"></a>00261         {
 250 <a name="l00262"></a>00262         }
 251 <a name="l00263"></a>00263
 252 <a name="l00264"></a><a class="code" href="structmof_1_1math_1_1row__of__matrix.html#ad26e390fd29458b4ccb1a6efc2b9fc2d">00264</a>         <span class="keywordtype">float</span> <a class="code" href="structmof_1_1math_1_1row__of__matrix.html#ad26e390fd29458b4ccb1a6efc2b9fc2d">operator []</a>(<span class="keywordtype">size_t</span> index)<span class="keyword"> const </span>{<span class="keywordflow">return</span> <a class="code" href="structmof_1_1math_1_1row__of__matrix.html#a0bdc04595281825f274b1016728f6a72">elements_</a>[<a class="code" href="structmof_1_1math_1_1row__of__matrix.html#af1d3a14a09bb4bbe6ddb5d8e37730b9e">row_index_</a> * (Dim + 1) + index];}
 253 <a name="l00265"></a>00265     };
 254 <a name="l00266"></a>00266 <span class="comment">//}}}</span>
 255 <a name="l00267"></a>00267
 256 <a name="l00268"></a>00268 }<span class="comment">// namespace math</span>
 257 <a name="l00269"></a>00269 }<span class="comment">// namespace mof</span>
 258 </pre></div></div>
 259 <hr class="footer"/><address style="text-align: right;"><small>MotifOfFunctionsLIBraryに対してThu Aug 19 09:37:38 2010に生成されました。&nbsp;
 260 <a href="http://www.doxygen.org/index.html">
 261 <img class="footer" src="doxygen.png" alt="doxygen"/></a> 1.6.3 </small></address>
 262 </body>
 263 </html>