Google Git
Sign in
apache / chemistry-site / 19af7d0fa4e1872c926a0e153c3f24d6bc996bae / . / content / java / 0.13.0 / maven / chemistry-opencmis-test / chemistry-opencmis-test-util / xref / org / apache / chemistry / opencmis / util / content / fractal / FractalCalculator.html
blob: 29352f11316e7d5e027f1cf53316ee9c63d0450a [file] [log] [blame]
<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd">
<html xmlns="http://www.w3.org/1999/xhtml" xml:lang="en" lang="en">
<head>
<meta http-equiv="content-type" content="text/html; charset=UTF-8" />
<title>FractalCalculator xref</title>
<link type="text/css" rel="stylesheet" href="../../../../../../../stylesheet.css" />
</head>
<body>
<div id="overview"><a href="../../../../../../../../apidocs/org/apache/chemistry/opencmis/util/content/fractal/FractalCalculator.html">View Javadoc</a></div><pre>
<a class="jxr_linenumber" name="1" href="#1">1</a> <em class="jxr_comment">////////////////////////////////////////////////////////////////////////////////</em>
<a class="jxr_linenumber" name="2" href="#2">2</a> <em class="jxr_comment">////////////////////////////////////////////////////////////////////////////////</em>
<a class="jxr_linenumber" name="3" href="#3">3</a> <em class="jxr_comment">/*</em>
<a class="jxr_linenumber" name="4" href="#4">4</a> <em class="jxr_comment"> * Licensed to the Apache Software Foundation (ASF) under one</em>
<a class="jxr_linenumber" name="5" href="#5">5</a> <em class="jxr_comment"> * or more contributor license agreements. See the NOTICE file</em>
<a class="jxr_linenumber" name="6" href="#6">6</a> <em class="jxr_comment"> * distributed with this work for additional information</em>
<a class="jxr_linenumber" name="7" href="#7">7</a> <em class="jxr_comment"> * regarding copyright ownership. The ASF licenses this file</em>
<a class="jxr_linenumber" name="8" href="#8">8</a> <em class="jxr_comment"> * to you under the Apache License, Version 2.0 (the</em>
<a class="jxr_linenumber" name="9" href="#9">9</a> <em class="jxr_comment"> * "License"); you may not use this file except in compliance</em>
<a class="jxr_linenumber" name="10" href="#10">10</a> <em class="jxr_comment"> * with the License. You may obtain a copy of the License at</em>
<a class="jxr_linenumber" name="11" href="#11">11</a> <em class="jxr_comment"> *</em>
<a class="jxr_linenumber" name="12" href="#12">12</a> <em class="jxr_comment"> * <a href="http://www.apache.org/licenses/LICENSE-2.0" target="alexandria_uri">http://www.apache.org/licenses/LICENSE-2.0</a></em>
<a class="jxr_linenumber" name="13" href="#13">13</a> <em class="jxr_comment"> *</em>
<a class="jxr_linenumber" name="14" href="#14">14</a> <em class="jxr_comment"> * Unless required by applicable law or agreed to in writing,</em>
<a class="jxr_linenumber" name="15" href="#15">15</a> <em class="jxr_comment"> * software distributed under the License is distributed on an</em>
<a class="jxr_linenumber" name="16" href="#16">16</a> <em class="jxr_comment"> * "AS IS" BASIS, WITHOUT WARRANTIES OR CONDITIONS OF ANY</em>
<a class="jxr_linenumber" name="17" href="#17">17</a> <em class="jxr_comment"> * KIND, either express or implied. See the License for the</em>
<a class="jxr_linenumber" name="18" href="#18">18</a> <em class="jxr_comment"> * specific language governing permissions and limitations</em>
<a class="jxr_linenumber" name="19" href="#19">19</a> <em class="jxr_comment"> * under the License.</em>
<a class="jxr_linenumber" name="20" href="#20">20</a> <em class="jxr_comment"> */</em>
<a class="jxr_linenumber" name="21" href="#21">21</a>
<a class="jxr_linenumber" name="22" href="#22">22</a> <strong class="jxr_keyword">package</strong> org.apache.chemistry.opencmis.util.content.fractal;
<a class="jxr_linenumber" name="23" href="#23">23</a>
<a class="jxr_linenumber" name="24" href="#24">24</a> <strong class="jxr_keyword">import</strong> java.awt.Color;
<a class="jxr_linenumber" name="25" href="#25">25</a> <strong class="jxr_keyword">import</strong> java.awt.image.BufferedImage;
<a class="jxr_linenumber" name="26" href="#26">26</a> <strong class="jxr_keyword">import</strong> java.util.Arrays;
<a class="jxr_linenumber" name="27" href="#27">27</a>
<a class="jxr_linenumber" name="28" href="#28">28</a> <strong class="jxr_keyword">final</strong> <strong class="jxr_keyword">class</strong> <a href="../../../../../../../org/apache/chemistry/opencmis/util/content/fractal/FractalCalculator.html">FractalCalculator</a> {
<a class="jxr_linenumber" name="29" href="#29">29</a> <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">int</strong>[] colorMap;
<a class="jxr_linenumber" name="30" href="#30">30</a> <strong class="jxr_keyword">protected</strong> <strong class="jxr_keyword">int</strong>[][] noIterations;
<a class="jxr_linenumber" name="31" href="#31">31</a> <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">double</strong> delta;
<a class="jxr_linenumber" name="32" href="#32">32</a> <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">double</strong> iRangeMax;
<a class="jxr_linenumber" name="33" href="#33">33</a> <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">double</strong> iRangeMin;
<a class="jxr_linenumber" name="34" href="#34">34</a> <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">int</strong> maxIterations;
<a class="jxr_linenumber" name="35" href="#35">35</a> <strong class="jxr_keyword">private</strong> <a href="../../../../../../../org/apache/chemistry/opencmis/util/content/fractal/ComplexRectangle.html">ComplexRectangle</a> newRect;
<a class="jxr_linenumber" name="36" href="#36">36</a> <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">int</strong> numColors;
<a class="jxr_linenumber" name="37" href="#37">37</a> <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">int</strong> imageHeight;
<a class="jxr_linenumber" name="38" href="#38">38</a> <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">int</strong> imageWidth;
<a class="jxr_linenumber" name="39" href="#39">39</a> <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">double</strong> rRangeMax;
<a class="jxr_linenumber" name="40" href="#40">40</a> <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">double</strong> rRangeMin;
<a class="jxr_linenumber" name="41" href="#41">41</a> <em class="jxr_comment">// For Julia set:</em>
<a class="jxr_linenumber" name="42" href="#42">42</a> <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">double</strong> cJuliaPointR = 0.0; <em class="jxr_comment">// Real</em>
<a class="jxr_linenumber" name="43" href="#43">43</a> <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">double</strong> cJuliaPointI = 0.0; <em class="jxr_comment">// Imaginary</em>
<a class="jxr_linenumber" name="44" href="#44">44</a> <strong class="jxr_keyword">boolean</strong> useJulia = false;
<a class="jxr_linenumber" name="45" href="#45">45</a>
<a class="jxr_linenumber" name="46" href="#46">46</a> <strong class="jxr_keyword">public</strong> <a href="../../../../../../../org/apache/chemistry/opencmis/util/content/fractal/FractalCalculator.html">FractalCalculator</a>(<a href="../../../../../../../org/apache/chemistry/opencmis/util/content/fractal/ComplexRectangle.html">ComplexRectangle</a> complRect, <strong class="jxr_keyword">int</strong> maxIters, <strong class="jxr_keyword">int</strong> imgWidth, <strong class="jxr_keyword">int</strong> imgHeight, <strong class="jxr_keyword">int</strong>[] colMap,
<a class="jxr_linenumber" name="47" href="#47">47</a> <a href="../../../../../../../org/apache/chemistry/opencmis/util/content/fractal/ComplexPoint.html">ComplexPoint</a> juliaPoint) {
<a class="jxr_linenumber" name="48" href="#48">48</a> maxIterations = maxIters;
<a class="jxr_linenumber" name="49" href="#49">49</a> newRect = complRect;
<a class="jxr_linenumber" name="50" href="#50">50</a> imageWidth = imgWidth;
<a class="jxr_linenumber" name="51" href="#51">51</a> imageHeight = imgHeight;
<a class="jxr_linenumber" name="52" href="#52">52</a> colorMap = Arrays.copyOf(colMap, colMap.length);
<a class="jxr_linenumber" name="53" href="#53">53</a> numColors = colorMap.length;
<a class="jxr_linenumber" name="54" href="#54">54</a> rRangeMin = newRect.getRMin();
<a class="jxr_linenumber" name="55" href="#55">55</a> rRangeMax = newRect.getRMax();
<a class="jxr_linenumber" name="56" href="#56">56</a> iRangeMin = newRect.getIMin();
<a class="jxr_linenumber" name="57" href="#57">57</a> iRangeMax = newRect.getIMax();
<a class="jxr_linenumber" name="58" href="#58">58</a> delta = (rRangeMax - rRangeMin) / imageWidth;
<a class="jxr_linenumber" name="59" href="#59">59</a> <strong class="jxr_keyword">if</strong> (<strong class="jxr_keyword">null</strong> != juliaPoint) {
<a class="jxr_linenumber" name="60" href="#60">60</a> cJuliaPointR = juliaPoint.getReal();
<a class="jxr_linenumber" name="61" href="#61">61</a> cJuliaPointI = juliaPoint.getImaginary();
<a class="jxr_linenumber" name="62" href="#62">62</a> useJulia = <strong class="jxr_keyword">true</strong>;
<a class="jxr_linenumber" name="63" href="#63">63</a> }
<a class="jxr_linenumber" name="64" href="#64">64</a> }
<a class="jxr_linenumber" name="65" href="#65">65</a>
<a class="jxr_linenumber" name="66" href="#66">66</a> <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">int</strong>[][] calcFractal() {
<a class="jxr_linenumber" name="67" href="#67">67</a> noIterations = <strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">int</strong>[ imageWidth ][ imageHeight ];
<a class="jxr_linenumber" name="68" href="#68">68</a>
<a class="jxr_linenumber" name="69" href="#69">69</a> <em class="jxr_comment">// For each pixel...</em>
<a class="jxr_linenumber" name="70" href="#70">70</a> <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">int</strong> x = 0; x &lt; imageWidth; x++) {
<a class="jxr_linenumber" name="71" href="#71">71</a> <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">int</strong> y = 0; y &lt; imageHeight; y++) {
<a class="jxr_linenumber" name="72" href="#72">72</a> <strong class="jxr_keyword">double</strong> zR = rRangeMin + x * delta;
<a class="jxr_linenumber" name="73" href="#73">73</a> <strong class="jxr_keyword">double</strong> zI = iRangeMin + (imageHeight - y) * delta;
<a class="jxr_linenumber" name="74" href="#74">74</a>
<a class="jxr_linenumber" name="75" href="#75">75</a> <em class="jxr_comment">// Is the point inside the set?</em>
<a class="jxr_linenumber" name="76" href="#76">76</a> <strong class="jxr_keyword">if</strong> (useJulia) {
<a class="jxr_linenumber" name="77" href="#77">77</a> noIterations[x][y] = testPointJuliaSet(zR, zI, maxIterations);
<a class="jxr_linenumber" name="78" href="#78">78</a> } <strong class="jxr_keyword">else</strong> {
<a class="jxr_linenumber" name="79" href="#79">79</a> noIterations[x][y] = testPointMandelbrot(zR, zI, maxIterations);
<a class="jxr_linenumber" name="80" href="#80">80</a> }
<a class="jxr_linenumber" name="81" href="#81">81</a> }
<a class="jxr_linenumber" name="82" href="#82">82</a> }
<a class="jxr_linenumber" name="83" href="#83">83</a> <strong class="jxr_keyword">return</strong> noIterations;
<a class="jxr_linenumber" name="84" href="#84">84</a> }
<a class="jxr_linenumber" name="85" href="#85">85</a>
<a class="jxr_linenumber" name="86" href="#86">86</a> <strong class="jxr_keyword">public</strong> BufferedImage mapItersToColors(<strong class="jxr_keyword">int</strong>[][] iterations) {
<a class="jxr_linenumber" name="87" href="#87">87</a>
<a class="jxr_linenumber" name="88" href="#88">88</a> <em class="jxr_comment">// Assign a color to every pixel ( x , y ) in the Image, corresponding</em>
<a class="jxr_linenumber" name="89" href="#89">89</a> <em class="jxr_comment">// to</em>
<a class="jxr_linenumber" name="90" href="#90">90</a> <em class="jxr_comment">// one point, z, in the imaginary plane ( zr, zi ).</em>
<a class="jxr_linenumber" name="91" href="#91">91</a> BufferedImage image = <strong class="jxr_keyword">new</strong> BufferedImage(imageWidth, imageHeight, BufferedImage.TYPE_3BYTE_BGR );
<a class="jxr_linenumber" name="92" href="#92">92</a>
<a class="jxr_linenumber" name="93" href="#93">93</a> <em class="jxr_comment">// For each pixel...</em>
<a class="jxr_linenumber" name="94" href="#94">94</a> <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">int</strong> x = 0; x &lt; imageWidth; x++) {
<a class="jxr_linenumber" name="95" href="#95">95</a> <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">int</strong> y = 0; y &lt; imageHeight; y++) {
<a class="jxr_linenumber" name="96" href="#96">96</a> <strong class="jxr_keyword">int</strong> color = getColor(iterations[x][y]);
<a class="jxr_linenumber" name="97" href="#97">97</a> image.setRGB(x, y, color);
<a class="jxr_linenumber" name="98" href="#98">98</a> }
<a class="jxr_linenumber" name="99" href="#99">99</a> }
<a class="jxr_linenumber" name="100" href="#100">100</a> <strong class="jxr_keyword">return</strong> image;
<a class="jxr_linenumber" name="101" href="#101">101</a> }
<a class="jxr_linenumber" name="102" href="#102">102</a>
<a class="jxr_linenumber" name="103" href="#103">103</a> <strong class="jxr_keyword">protected</strong> <strong class="jxr_keyword">int</strong> getColor(<strong class="jxr_keyword">int</strong> numIterations) {
<a class="jxr_linenumber" name="104" href="#104">104</a> <strong class="jxr_keyword">int</strong> c = Color.black.getRGB();
<a class="jxr_linenumber" name="105" href="#105">105</a>
<a class="jxr_linenumber" name="106" href="#106">106</a> <strong class="jxr_keyword">if</strong> (numIterations != 0) {
<a class="jxr_linenumber" name="107" href="#107">107</a> <em class="jxr_comment">// The point is outside the set. It gets a color based on the number</em>
<a class="jxr_linenumber" name="108" href="#108">108</a> <em class="jxr_comment">// of iterations it took to know this.</em>
<a class="jxr_linenumber" name="109" href="#109">109</a> <strong class="jxr_keyword">int</strong> colorNum = (<strong class="jxr_keyword">int</strong>) (numColors * (1.0 - (<strong class="jxr_keyword">float</strong>) numIterations / (<strong class="jxr_keyword">float</strong>) maxIterations));
<a class="jxr_linenumber" name="110" href="#110">110</a> colorNum = (colorNum == numColors) ? 0 : colorNum;
<a class="jxr_linenumber" name="111" href="#111">111</a>
<a class="jxr_linenumber" name="112" href="#112">112</a> c = colorMap[colorNum];
<a class="jxr_linenumber" name="113" href="#113">113</a> }
<a class="jxr_linenumber" name="114" href="#114">114</a> <strong class="jxr_keyword">return</strong> c;
<a class="jxr_linenumber" name="115" href="#115">115</a> }
<a class="jxr_linenumber" name="116" href="#116">116</a>
<a class="jxr_linenumber" name="117" href="#117">117</a> <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">int</strong> testPointMandelbrot(<strong class="jxr_keyword">double</strong> cR, <strong class="jxr_keyword">double</strong> cI, <strong class="jxr_keyword">int</strong> maxIterations) {
<a class="jxr_linenumber" name="118" href="#118">118</a> <em class="jxr_comment">// Is the given complex point, (cR, cI), in the Mandelbrot set?</em>
<a class="jxr_linenumber" name="119" href="#119">119</a> <em class="jxr_comment">// Use the formula: z &lt;= z*z + c, where z is initially equal to c.</em>
<a class="jxr_linenumber" name="120" href="#120">120</a> <em class="jxr_comment">// If |z| &gt;= 2, then the point is not in the set.</em>
<a class="jxr_linenumber" name="121" href="#121">121</a> <em class="jxr_comment">// Return 0 if the point is in the set; else return the number of</em>
<a class="jxr_linenumber" name="122" href="#122">122</a> <em class="jxr_comment">// iterations it took to decide that the point is not in the set.</em>
<a class="jxr_linenumber" name="123" href="#123">123</a> <strong class="jxr_keyword">double</strong> zR = cR;
<a class="jxr_linenumber" name="124" href="#124">124</a> <strong class="jxr_keyword">double</strong> zI = cI;
<a class="jxr_linenumber" name="125" href="#125">125</a>
<a class="jxr_linenumber" name="126" href="#126">126</a> <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">int</strong> i = 1; i &lt;= maxIterations; i++) {
<a class="jxr_linenumber" name="127" href="#127">127</a> <em class="jxr_comment">// To square a complex number: (a+bi)(a+bi) = a*a - b*b + 2abi</em>
<a class="jxr_linenumber" name="128" href="#128">128</a> <strong class="jxr_keyword">double</strong> zROld = zR;
<a class="jxr_linenumber" name="129" href="#129">129</a> zR = zR * zR - zI * zI + cR;
<a class="jxr_linenumber" name="130" href="#130">130</a> zI = 2 * zROld * zI + cI;
<a class="jxr_linenumber" name="131" href="#131">131</a>
<a class="jxr_linenumber" name="132" href="#132">132</a> <em class="jxr_comment">// We know that if the distance from z to the origin is &gt;= 2</em>
<a class="jxr_linenumber" name="133" href="#133">133</a> <em class="jxr_comment">// then the point is out of the set. To avoid a square root,</em>
<a class="jxr_linenumber" name="134" href="#134">134</a> <em class="jxr_comment">// we'll instead check if the distance squared &gt;= 4.</em>
<a class="jxr_linenumber" name="135" href="#135">135</a> <strong class="jxr_keyword">double</strong> distSquared = zR * zR + zI * zI;
<a class="jxr_linenumber" name="136" href="#136">136</a> <strong class="jxr_keyword">if</strong> (distSquared &gt;= 4) {
<a class="jxr_linenumber" name="137" href="#137">137</a> <strong class="jxr_keyword">return</strong> i;
<a class="jxr_linenumber" name="138" href="#138">138</a> }
<a class="jxr_linenumber" name="139" href="#139">139</a> }
<a class="jxr_linenumber" name="140" href="#140">140</a> <strong class="jxr_keyword">return</strong> 0;
<a class="jxr_linenumber" name="141" href="#141">141</a> }
<a class="jxr_linenumber" name="142" href="#142">142</a>
<a class="jxr_linenumber" name="143" href="#143">143</a> <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">int</strong> testPointJuliaSet(<strong class="jxr_keyword">double</strong> zR, <strong class="jxr_keyword">double</strong> zI, <strong class="jxr_keyword">int</strong> maxIterations) {
<a class="jxr_linenumber" name="144" href="#144">144</a> <em class="jxr_comment">// Is the given complex point, (zR, zI), in the Julia set?</em>
<a class="jxr_linenumber" name="145" href="#145">145</a> <em class="jxr_comment">// Use the formula: z &lt;= z*z + c, where z is the point being tested,</em>
<a class="jxr_linenumber" name="146" href="#146">146</a> <em class="jxr_comment">// and c is the Julia Set constant.</em>
<a class="jxr_linenumber" name="147" href="#147">147</a> <em class="jxr_comment">// If |z| &gt;= 2, then the point is not in the set.</em>
<a class="jxr_linenumber" name="148" href="#148">148</a> <em class="jxr_comment">// Return 0 if the point is in the set; else return the number of</em>
<a class="jxr_linenumber" name="149" href="#149">149</a> <em class="jxr_comment">// iterations it took to decide that the point is not in the set.</em>
<a class="jxr_linenumber" name="150" href="#150">150</a> <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">int</strong> i = 1; i &lt;= maxIterations; i++) {
<a class="jxr_linenumber" name="151" href="#151">151</a> <strong class="jxr_keyword">double</strong> zROld = zR;
<a class="jxr_linenumber" name="152" href="#152">152</a> <em class="jxr_comment">// To square a complex number: (a+bi)(a+bi) = a*a - b*b + 2abi</em>
<a class="jxr_linenumber" name="153" href="#153">153</a> zR = zR * zR - zI * zI + cJuliaPointR;
<a class="jxr_linenumber" name="154" href="#154">154</a> zI = 2 * zROld * zI + cJuliaPointI;
<a class="jxr_linenumber" name="155" href="#155">155</a> <em class="jxr_comment">// We know that if the distance from z to the origin is &gt;= 2</em>
<a class="jxr_linenumber" name="156" href="#156">156</a> <em class="jxr_comment">// then the point is out of the set. To avoid a square root,</em>
<a class="jxr_linenumber" name="157" href="#157">157</a> <em class="jxr_comment">// we'll instead check if the distance squared &gt;= 4.</em>
<a class="jxr_linenumber" name="158" href="#158">158</a> <strong class="jxr_keyword">double</strong> distSquared = zR * zR + zI * zI;
<a class="jxr_linenumber" name="159" href="#159">159</a> <strong class="jxr_keyword">if</strong> (distSquared &gt;= 4) {
<a class="jxr_linenumber" name="160" href="#160">160</a> <strong class="jxr_keyword">return</strong> i;
<a class="jxr_linenumber" name="161" href="#161">161</a> }
<a class="jxr_linenumber" name="162" href="#162">162</a> }
<a class="jxr_linenumber" name="163" href="#163">163</a> <strong class="jxr_keyword">return</strong> 0;
<a class="jxr_linenumber" name="164" href="#164">164</a> }
<a class="jxr_linenumber" name="165" href="#165">165</a> }
</pre>
<hr/><div id="footer">This page was automatically generated by <a href="http://maven.apache.org/">Maven</a></div></body>
</html>