html/warpcoord_8h-source.html

   1 <!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
   2 <html><head><meta http-equiv="Content-Type" content="text/html;charset=iso-8859-1">
   3 <title>WarpTree: src/warpcoord.h Source File</title>
   4 <link href="doxygen.css" rel="stylesheet" type="text/css">
   5 <link href="tabs.css" rel="stylesheet" type="text/css">
   6 </head><body>
   7 <!-- Generated by Doxygen 1.5.0 -->
   8 <div class="tabs">
   9   <ul>
  10     <li><a href="index.html"><span>Main&nbsp;Page</span></a></li>
  11     <li><a href="annotated.html"><span>Classes</span></a></li>
  12     <li id="current"><a href="files.html"><span>Files</span></a></li>
  13   </ul></div>
  14 <h1>src/warpcoord.h</h1><div class="fragment"><pre class="fragment"><a name="l00001"></a>00001 <span class="comment">/* This file is part of WarpTree</span>
  15 <a name="l00002"></a>00002 <span class="comment"> *</span>
  16 <a name="l00003"></a>00003 <span class="comment"> * Copyright (C) 2007 Jos van den Oever &lt;jos@vandenoever.info&gt;</span>
  17 <a name="l00004"></a>00004 <span class="comment"> *</span>
  18 <a name="l00005"></a>00005 <span class="comment"> * This library is free software; you can redistribute it and/or</span>
  19 <a name="l00006"></a>00006 <span class="comment"> * modify it under the terms of the GNU Library General Public</span>
  20 <a name="l00007"></a>00007 <span class="comment"> * License as published by the Free Software Foundation; either</span>
  21 <a name="l00008"></a>00008 <span class="comment"> * version 2 of the License, or (at your option) any later version.</span>
  22 <a name="l00009"></a>00009 <span class="comment"> *</span>
  23 <a name="l00010"></a>00010 <span class="comment"> * This library is distributed in the hope that it will be useful,</span>
  24 <a name="l00011"></a>00011 <span class="comment"> * but WITHOUT ANY WARRANTY; without even the implied warranty of</span>
  25 <a name="l00012"></a>00012 <span class="comment"> * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU</span>
  26 <a name="l00013"></a>00013 <span class="comment"> * Library General Public License for more details.</span>
  27 <a name="l00014"></a>00014 <span class="comment"> *</span>
  28 <a name="l00015"></a>00015 <span class="comment"> * You should have received a copy of the GNU Library General Public License</span>
  29 <a name="l00016"></a>00016 <span class="comment"> * along with this library; see the file COPYING.LIB.  If not, write to</span>
  30 <a name="l00017"></a>00017 <span class="comment"> * the Free Software Foundation, Inc., 51 Franklin Street, Fifth Floor,</span>
  31 <a name="l00018"></a>00018 <span class="comment"> * Boston, MA 02110-1301, USA.</span>
  32 <a name="l00019"></a>00019 <span class="comment"> */</span>
  33 <a name="l00020"></a>00020
  34 <a name="l00021"></a>00021 <span class="comment">// code for transformations was adapted from HyperTree published under</span>
  35 <a name="l00022"></a>00022 <span class="comment">// MIT License 2001 www.bouthier.net</span>
  36 <a name="l00023"></a>00023
  37 <a name="l00024"></a>00024 <span class="preprocessor">#ifndef WARPTREE_H</span>
  38 <a name="l00025"></a>00025 <span class="preprocessor"></span><span class="preprocessor">#define WARPTREE_H</span>
  39 <a name="l00026"></a>00026 <span class="preprocessor"></span>
  40 <a name="l00027"></a>00027 <span class="preprocessor">#include &lt;iostream&gt;</span>
  41 <a name="l00028"></a>00028 <span class="preprocessor">#include &lt;cmath&gt;</span>
  42 <a name="l00029"></a>00029 <span class="preprocessor">#include &lt;cassert&gt;</span>
  43 <a name="l00030"></a>00030 <span class="preprocessor">#include &lt;cstdlib&gt;</span>
  44 <a name="l00031"></a>00031
  45 <a name="l00032"></a>00032 <span class="comment">/* This class provides the methods to handle the points in the graph.</span>
  46 <a name="l00033"></a>00033 <span class="comment"> * It allows points to be translated in the unit square.</span>
  47 <a name="l00034"></a>00034 <span class="comment"> */</span>
  48 <a name="l00035"></a>00035
  49 <a name="l00036"></a>00036 <span class="preprocessor">#define ZONE_LENGTH 4 // size of the zone</span>
  50 <a name="l00037"></a>00037 <span class="preprocessor"></span>
  51 <a name="l00038"></a>00038 <span class="comment">// if WARPCOORD is not defined, the standard Poincare disk projection is used</span>
  52 <a name="l00039"></a>00039 <span class="preprocessor">#define WARPCOORD</span>
  53 <a name="l00040"></a>00040 <span class="preprocessor"></span>
  54 <a name="l00041"></a>00041 <span class="keyword">class </span>WarpCoord;
  55 <a name="l00042"></a>00042 std::ostream&amp; operator&lt;&lt;(std::ostream&amp; o, <span class="keyword">const</span> WarpCoord&amp; e);
  56 <a name="l00043"></a>00043 <span class="keyword">class </span>WarpCoord {
  57 <a name="l00044"></a>00044 <span class="keyword">private</span>:
  58 <a name="l00045"></a>00045     <span class="keywordtype">float</span> mx;
  59 <a name="l00046"></a>00046     <span class="keywordtype">float</span> my;
  60 <a name="l00047"></a>00047     <span class="keywordtype">float</span> md;
  61 <a name="l00048"></a>00048
  62 <a name="l00049"></a>00049     <span class="comment">/* the value for d is cached in this class. This function updates the</span>
  63 <a name="l00050"></a>00050 <span class="comment">       value */</span>
  64 <a name="l00051"></a>00051     <span class="keyword">inline</span> <span class="keywordtype">void</span> updateD() {
  65 <a name="l00052"></a>00052         md = sqrt(d2());
  66 <a name="l00053"></a>00053     }
  67 <a name="l00054"></a>00054     <span class="comment">/* multiply this point with point z */</span>
  68 <a name="l00055"></a>00055     <span class="keywordtype">void</span>
  69 <a name="l00056"></a>00056     operator*=(<span class="keyword">const</span> WarpCoord&amp; z) {
  70 <a name="l00057"></a>00057         <span class="keywordtype">float</span> tx = mx;
  71 <a name="l00058"></a>00058         mx = tx*z.mx - my * z.my;
  72 <a name="l00059"></a>00059         my = tx*z.my + my * z.mx;
  73 <a name="l00060"></a>00060         updateD();
  74 <a name="l00061"></a>00061     }
  75 <a name="l00062"></a>00062     <span class="comment">/* divide this point by point z */</span>
  76 <a name="l00063"></a>00063     <span class="keywordtype">void</span>
  77 <a name="l00064"></a>00064     operator/=(<span class="keyword">const</span> WarpCoord&amp; z) {
  78 <a name="l00065"></a>00065         <span class="keywordtype">float</span> d2 = z.d2();
  79 <a name="l00066"></a>00066         <span class="keywordtype">float</span> tx = mx;
  80 <a name="l00067"></a>00067         mx = (tx*z.mx + my*z.my) / d2;
  81 <a name="l00068"></a>00068         my = (my*z.mx - tx*z.my) / d2;
  82 <a name="l00069"></a>00069         updateD();
  83 <a name="l00070"></a>00070     }
  84 <a name="l00071"></a>00071 <span class="keyword">public</span>:
  85 <a name="l00072"></a>00072     WarpCoord() :mx(0), my(0), md(0) {}
  86 <a name="l00073"></a>00073     WarpCoord(<span class="keywordtype">float</span> x, <span class="keywordtype">float</span> y) :mx(x), my(y) {
  87 <a name="l00074"></a>00074         updateD();
  88 <a name="l00075"></a>00075     }
  89 <a name="l00076"></a>00076     <span class="keyword">inline</span> <span class="keywordtype">float</span> x()<span class="keyword"> const </span>{ <span class="keywordflow">return</span> mx; }
  90 <a name="l00077"></a>00077     <span class="keyword">inline</span> <span class="keywordtype">float</span> y()<span class="keyword"> const </span>{ <span class="keywordflow">return</span> my; }
  91 <a name="l00078"></a>00078     <span class="keyword">inline</span> <span class="keywordtype">float</span> d()<span class="keyword"> const </span>{ <span class="keywordflow">return</span> md; }
  92 <a name="l00079"></a>00079     <span class="keyword">inline</span> <span class="keywordtype">float</span> d2()<span class="keyword"> const </span>{ <span class="keywordflow">return</span> mx*mx + my*my; }
  93 <a name="l00080"></a>00080     <span class="comment">/* calculate the arg of this point */</span>
  94 <a name="l00081"></a>00081     <span class="keywordtype">float</span> arg()<span class="keyword"> const </span>{
  95 <a name="l00082"></a>00082         <span class="keywordtype">float</span> a = atan(my / mx);
  96 <a name="l00083"></a>00083         <span class="keywordflow">if</span> (mx &lt; 0) {
  97 <a name="l00084"></a>00084             a += M_PI;
  98 <a name="l00085"></a>00085         } <span class="keywordflow">else</span> <span class="keywordflow">if</span> (my &lt; 0) {
  99 <a name="l00086"></a>00086             a += 2 * M_PI;
 100 <a name="l00087"></a>00087         }
 101 <a name="l00088"></a>00088         <span class="keywordflow">return</span> a;
 102 <a name="l00089"></a>00089     }
 103 <a name="l00090"></a>00090     <span class="comment">/* translate this point by adding the vector t */</span>
 104 <a name="l00091"></a>00091     <span class="keywordtype">void</span> translate(<span class="keyword">const</span> WarpCoord&amp; t) {
 105 <a name="l00092"></a>00092         <span class="keywordtype">float</span> denX = mx * t.mx + my * t.my + 1;
 106 <a name="l00093"></a>00093         <span class="keywordtype">float</span> denY = my * t.mx - mx * t.my;
 107 <a name="l00094"></a>00094         <span class="keywordtype">float</span> dd = denX * denX + denY * denY;
 108 <a name="l00095"></a>00095
 109 <a name="l00096"></a>00096         <span class="keywordtype">double</span> numX = mx + t.mx;
 110 <a name="l00097"></a>00097         <span class="keywordtype">double</span> numY = my + t.my;
 111 <a name="l00098"></a>00098
 112 <a name="l00099"></a>00099         mx = (numX * denX + numY * denY) / dd;
 113 <a name="l00100"></a>00100         my = (numY * denX - numX * denY) / dd;
 114 <a name="l00101"></a>00101         updateD();
 115 <a name="l00102"></a>00102     }
 116 <a name="l00103"></a>00103     <span class="comment">/* zoom this point. this means the point move closer to or away from the</span>
 117 <a name="l00104"></a>00104 <span class="comment">       center of the square */</span>
 118 <a name="l00105"></a>00105     <span class="keywordtype">void</span> zoom(<span class="keywordtype">float</span> f) {
 119 <a name="l00106"></a>00106         assert(f &gt; 0);
 120 <a name="l00107"></a>00107         <span class="keywordflow">if</span> (md &lt; 0.001) <span class="keywordflow">return</span>;
 121 <a name="l00108"></a>00108         <span class="keywordtype">float</span> r = md;
 122 <a name="l00109"></a>00109         <span class="keywordtype">float</span> c = pow((1+md)/(1-md),f);
 123 <a name="l00110"></a>00110         c = (c-1)/(c+1)/r;
 124 <a name="l00111"></a>00111         mx *= c;
 125 <a name="l00112"></a>00112         my *= c;
 126 <a name="l00113"></a>00113         md *= c;
 127 <a name="l00114"></a>00114     }
 128 <a name="l00115"></a>00115     <span class="comment">/* move this point from the disk positiong to the square position */</span>
 129 <a name="l00116"></a>00116     <span class="keywordtype">void</span> stretch() {
 130 <a name="l00117"></a>00117         <span class="keywordflow">if</span> (md &lt; 0.0001) <span class="keywordflow">return</span>;
 131 <a name="l00118"></a>00118 <span class="preprocessor">#ifdef WARPCOORD</span>
 132 <a name="l00119"></a>00119 <span class="preprocessor"></span>        <span class="keywordtype">float</span> fx = fabs(mx);
 133 <a name="l00120"></a>00120         <span class="keywordtype">float</span> fy = fabs(my);
 134 <a name="l00121"></a>00121         <span class="keywordflow">if</span> (fy &gt; 0 &amp;&amp; fy &gt; fx) {
 135 <a name="l00122"></a>00122             <span class="keywordtype">float</span> f = md/fy;
 136 <a name="l00123"></a>00123             my = (my&gt;0) ?md :-md;
 137 <a name="l00124"></a>00124             mx *= f;
 138 <a name="l00125"></a>00125             md *= f;
 139 <a name="l00126"></a>00126         } <span class="keywordflow">else</span> <span class="keywordflow">if</span> (fx &gt; 0) {
 140 <a name="l00127"></a>00127             <span class="keywordtype">float</span> f = md/fx;
 141 <a name="l00128"></a>00128             mx = (mx&gt;0) ?md :-md;
 142 <a name="l00129"></a>00129             my *= f;
 143 <a name="l00130"></a>00130             md *= f;
 144 <a name="l00131"></a>00131         }
 145 <a name="l00132"></a>00132 <span class="preprocessor">#endif</span>
 146 <a name="l00133"></a>00133 <span class="preprocessor"></span>    }
 147 <a name="l00134"></a>00134     <span class="comment">/* move this point from the square position to the disk position */</span>
 148 <a name="l00135"></a>00135     <span class="keywordtype">void</span> unstretch() {
 149 <a name="l00136"></a>00136         assert(md &gt;= 0);
 150 <a name="l00137"></a>00137         <span class="keywordflow">if</span> (md &lt; 0.0001) <span class="keywordflow">return</span>;
 151 <a name="l00138"></a>00138 <span class="preprocessor">#ifdef WARPCOORD</span>
 152 <a name="l00139"></a>00139 <span class="preprocessor"></span>        <span class="keywordtype">float</span> f = 1;
 153 <a name="l00140"></a>00140         <span class="keywordtype">float</span> fx = fabs(mx);
 154 <a name="l00141"></a>00141         <span class="keywordtype">float</span> fy = fabs(my);
 155 <a name="l00142"></a>00142         <span class="keywordflow">if</span> (fy &gt; 0 &amp;&amp; fy &gt; fx) {
 156 <a name="l00143"></a>00143             f = md/fy;
 157 <a name="l00144"></a>00144             md = (my&gt;0) ?my :-my;
 158 <a name="l00145"></a>00145         } <span class="keywordflow">else</span> <span class="keywordflow">if</span> (fx &gt; 0) {
 159 <a name="l00146"></a>00146             f = md/fx;
 160 <a name="l00147"></a>00147             md = (mx&gt;0) ?mx :-mx;
 161 <a name="l00148"></a>00148         }
 162 <a name="l00149"></a>00149         mx /= f;
 163 <a name="l00150"></a>00150         my /= f;
 164 <a name="l00151"></a>00151 <span class="preprocessor">#endif</span>
 165 <a name="l00152"></a>00152 <span class="preprocessor"></span>    }
 166 <a name="l00153"></a>00153     <span class="comment">/* is this point visible, in other words, is it not too close to the edge</span>
 167 <a name="l00154"></a>00154 <span class="comment">       of the unit square */</span>
 168 <a name="l00155"></a>00155     <span class="keyword">inline</span> <span class="keywordtype">bool</span> visible()<span class="keyword"> const </span>{
 169 <a name="l00156"></a>00156         <span class="keywordflow">return</span> mx &gt; -0.95 &amp;&amp; mx &lt; 0.95 &amp;&amp; my &gt; -0.95 &amp;&amp; my &lt; 0.95;
 170 <a name="l00157"></a>00157     }
 171 <a name="l00158"></a>00158     <span class="comment">/* is this point valid, in other words, does it lie inside the unit</span>
 172 <a name="l00159"></a>00159 <span class="comment">       square */</span>
 173 <a name="l00160"></a>00160     <span class="keyword">inline</span> <span class="keywordtype">bool</span> valid()<span class="keyword"> const </span>{
 174 <a name="l00161"></a>00161         <span class="keywordflow">return</span> mx &gt; -1 &amp;&amp; mx &lt; 1 &amp;&amp; my &gt; -1 &amp;&amp; my &lt; 1;
 175 <a name="l00162"></a>00162     }
 176 <a name="l00163"></a>00163     <span class="comment">/* return the negative value of this point */</span>
 177 <a name="l00164"></a>00164     <span class="keyword">inline</span> WarpCoord operator-()<span class="keyword"> const </span>{
 178 <a name="l00165"></a>00165         <span class="keywordflow">return</span> WarpCoord(-mx, -my);
 179 <a name="l00166"></a>00166     }
 180 <a name="l00167"></a>00167     <span class="comment">/* substract the value of point e from this point. Note that this is not</span>
 181 <a name="l00168"></a>00168 <span class="comment">       equal to translation. */</span>
 182 <a name="l00169"></a>00169     <span class="keyword">inline</span> <span class="keywordtype">void</span> operator-=(<span class="keyword">const</span> WarpCoord&amp; e) {
 183 <a name="l00170"></a>00170         mx -= e.mx;
 184 <a name="l00171"></a>00171         my -= e.my;
 185 <a name="l00172"></a>00172         updateD();
 186 <a name="l00173"></a>00173     }
 187 <a name="l00174"></a>00174     <span class="comment">/* apply the translation that is needed to get from point c1 to c2 to this</span>
 188 <a name="l00175"></a>00175 <span class="comment">       point */</span>
 189 <a name="l00176"></a>00176     <span class="keywordtype">void</span> transform(<span class="keyword">const</span> WarpCoord&amp; c1, <span class="keyword">const</span> WarpCoord&amp; c2) {
 190 <a name="l00177"></a>00177         WarpCoord p(c1.mx + c2.mx, c1.my + c2.my);
 191 <a name="l00178"></a>00178
 192 <a name="l00179"></a>00179         WarpCoord d(c2.mx, -c2.my);
 193 <a name="l00180"></a>00180         d *= c1;
 194 <a name="l00181"></a>00181         d.mx += 1;
 195 <a name="l00182"></a>00182         p /= d;
 196 <a name="l00183"></a>00183
 197 <a name="l00184"></a>00184         WarpCoord o(c1.mx, -c1.my);
 198 <a name="l00185"></a>00185         o *= c2;
 199 <a name="l00186"></a>00186         o.mx += 1;
 200 <a name="l00187"></a>00187         o /= d;
 201 <a name="l00188"></a>00188
 202 <a name="l00189"></a>00189         WarpCoord z(*<span class="keyword">this</span>);
 203 <a name="l00190"></a>00190         *<span class="keyword">this</span> *= o;
 204 <a name="l00191"></a>00191         mx += p.mx;
 205 <a name="l00192"></a>00192         my += p.my;
 206 <a name="l00193"></a>00193
 207 <a name="l00194"></a>00194         d.set(p.mx, -p.my);
 208 <a name="l00195"></a>00195         d *= z;
 209 <a name="l00196"></a>00196         d *= o;
 210 <a name="l00197"></a>00197         d.mx += 1;
 211 <a name="l00198"></a>00198
 212 <a name="l00199"></a>00199         *<span class="keyword">this</span> /= d;
 213 <a name="l00200"></a>00200         updateD();
 214 <a name="l00201"></a>00201     }
 215 <a name="l00202"></a>00202     <span class="comment">/* change the value of this point */</span>
 216 <a name="l00203"></a>00203     <span class="keyword">inline</span> <span class="keywordtype">void</span> set(<span class="keywordtype">float</span> x, <span class="keywordtype">float</span> y) {
 217 <a name="l00204"></a>00204         mx = x;
 218 <a name="l00205"></a>00205         my = y;
 219 <a name="l00206"></a>00206         updateD();
 220 <a name="l00207"></a>00207     }
 221 <a name="l00208"></a>00208 };
 222 <a name="l00209"></a>00209
 223 <a name="l00210"></a>00210 <span class="preprocessor">#endif </span>
 224 </pre></div><hr size="1"><address style="align: right;"><small>Generated on Sat Jun 2 19:18:06 2007 for WarpTree by&nbsp;
 225 <a href="http://www.doxygen.org/index.html">
 226 <img src="doxygen.png" alt="doxygen" align="middle" border="0"></a> 1.5.0 </small></address>
 227 </body>
 228 </html>