You are viewing a plain text version of this content. The canonical link for it is here.
Posted to commits@commons.apache.org by lu...@apache.org on 2013/04/07 09:52:13 UTC

svn commit: r857590 [13/48] - in /websites/production/commons/content/proper/commons-math/xref-test: ./ org/apache/commons/math3/ org/apache/commons/math3/analysis/ org/apache/commons/math3/analysis/differentiation/ org/apache/commons/math3/analysis/fu...

Modified: websites/production/commons/content/proper/commons-math/xref-test/org/apache/commons/math3/geometry/euclidean/threed/Vector3DTest.html
==============================================================================
--- websites/production/commons/content/proper/commons-math/xref-test/org/apache/commons/math3/geometry/euclidean/threed/Vector3DTest.html (original)
+++ websites/production/commons/content/proper/commons-math/xref-test/org/apache/commons/math3/geometry/euclidean/threed/Vector3DTest.html Sun Apr  7 07:52:05 2013
@@ -27,319 +27,399 @@
 <a class="jxr_linenumber" name="17" href="#17">17</a>  
 <a class="jxr_linenumber" name="18" href="#18">18</a>  <strong class="jxr_keyword">package</strong> org.apache.commons.math3.geometry.euclidean.threed;
 <a class="jxr_linenumber" name="19" href="#19">19</a>  
-<a class="jxr_linenumber" name="20" href="#20">20</a>  <strong class="jxr_keyword">import</strong> org.apache.commons.math3.exception.DimensionMismatchException;
-<a class="jxr_linenumber" name="21" href="#21">21</a>  <strong class="jxr_keyword">import</strong> org.apache.commons.math3.exception.MathArithmeticException;
-<a class="jxr_linenumber" name="22" href="#22">22</a>  <strong class="jxr_keyword">import</strong> org.apache.commons.math3.random.Well1024a;
-<a class="jxr_linenumber" name="23" href="#23">23</a>  <strong class="jxr_keyword">import</strong> org.apache.commons.math3.util.FastMath;
-<a class="jxr_linenumber" name="24" href="#24">24</a>  <strong class="jxr_keyword">import</strong> org.junit.Assert;
-<a class="jxr_linenumber" name="25" href="#25">25</a>  <strong class="jxr_keyword">import</strong> org.junit.Test;
-<a class="jxr_linenumber" name="26" href="#26">26</a>  
-<a class="jxr_linenumber" name="27" href="#27">27</a>  <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">class</strong> <a href="../../../../../../../org/apache/commons/math3/geometry/euclidean/threed/Vector3DTest.html">Vector3DTest</a> {
-<a class="jxr_linenumber" name="28" href="#28">28</a>      @Test
-<a class="jxr_linenumber" name="29" href="#29">29</a>      <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testConstructors() <strong class="jxr_keyword">throws</strong> DimensionMismatchException {
-<a class="jxr_linenumber" name="30" href="#30">30</a>          <strong class="jxr_keyword">double</strong> r = FastMath.sqrt(2) /2;
-<a class="jxr_linenumber" name="31" href="#31">31</a>          checkVector(<strong class="jxr_keyword">new</strong> Vector3D(2, <strong class="jxr_keyword">new</strong> Vector3D(FastMath.PI / 3, -FastMath.PI / 4)),
-<a class="jxr_linenumber" name="32" href="#32">32</a>                      r, r * FastMath.sqrt(3), -2 * r);
-<a class="jxr_linenumber" name="33" href="#33">33</a>          checkVector(<strong class="jxr_keyword">new</strong> Vector3D(2, Vector3D.PLUS_I,
-<a class="jxr_linenumber" name="34" href="#34">34</a>                                   -3, Vector3D.MINUS_K),
-<a class="jxr_linenumber" name="35" href="#35">35</a>                      2, 0, 3);
-<a class="jxr_linenumber" name="36" href="#36">36</a>          checkVector(<strong class="jxr_keyword">new</strong> Vector3D(2, Vector3D.PLUS_I,
-<a class="jxr_linenumber" name="37" href="#37">37</a>                                   5, Vector3D.PLUS_J,
-<a class="jxr_linenumber" name="38" href="#38">38</a>                                   -3, Vector3D.MINUS_K),
-<a class="jxr_linenumber" name="39" href="#39">39</a>                      2, 5, 3);
-<a class="jxr_linenumber" name="40" href="#40">40</a>          checkVector(<strong class="jxr_keyword">new</strong> Vector3D(2, Vector3D.PLUS_I,
-<a class="jxr_linenumber" name="41" href="#41">41</a>                                   5, Vector3D.PLUS_J,
-<a class="jxr_linenumber" name="42" href="#42">42</a>                                   5, Vector3D.MINUS_J,
-<a class="jxr_linenumber" name="43" href="#43">43</a>                                   -3, Vector3D.MINUS_K),
-<a class="jxr_linenumber" name="44" href="#44">44</a>                      2, 0, 3);
-<a class="jxr_linenumber" name="45" href="#45">45</a>          checkVector(<strong class="jxr_keyword">new</strong> Vector3D(<strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[] { 2,  5,  -3 }),
-<a class="jxr_linenumber" name="46" href="#46">46</a>                      2, 5, -3);
-<a class="jxr_linenumber" name="47" href="#47">47</a>      }
-<a class="jxr_linenumber" name="48" href="#48">48</a>  
-<a class="jxr_linenumber" name="49" href="#49">49</a>      @Test(expected=DimensionMismatchException.<strong class="jxr_keyword">class</strong>)
-<a class="jxr_linenumber" name="50" href="#50">50</a>      <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testWrongDimension() <strong class="jxr_keyword">throws</strong> DimensionMismatchException {
-<a class="jxr_linenumber" name="51" href="#51">51</a>          <strong class="jxr_keyword">new</strong> Vector3D(<strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[] { 2,  5 });
-<a class="jxr_linenumber" name="52" href="#52">52</a>      }
-<a class="jxr_linenumber" name="53" href="#53">53</a>  
-<a class="jxr_linenumber" name="54" href="#54">54</a>      @Test
-<a class="jxr_linenumber" name="55" href="#55">55</a>      <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testCoordinates() {
-<a class="jxr_linenumber" name="56" href="#56">56</a>          Vector3D v = <strong class="jxr_keyword">new</strong> Vector3D(1, 2, 3);
-<a class="jxr_linenumber" name="57" href="#57">57</a>          Assert.assertTrue(FastMath.abs(v.getX() - 1) &lt; 1.0e-12);
-<a class="jxr_linenumber" name="58" href="#58">58</a>          Assert.assertTrue(FastMath.abs(v.getY() - 2) &lt; 1.0e-12);
-<a class="jxr_linenumber" name="59" href="#59">59</a>          Assert.assertTrue(FastMath.abs(v.getZ() - 3) &lt; 1.0e-12);
-<a class="jxr_linenumber" name="60" href="#60">60</a>          <strong class="jxr_keyword">double</strong>[] coordinates = v.toArray();
-<a class="jxr_linenumber" name="61" href="#61">61</a>          Assert.assertTrue(FastMath.abs(coordinates[0] - 1) &lt; 1.0e-12);
-<a class="jxr_linenumber" name="62" href="#62">62</a>          Assert.assertTrue(FastMath.abs(coordinates[1] - 2) &lt; 1.0e-12);
-<a class="jxr_linenumber" name="63" href="#63">63</a>          Assert.assertTrue(FastMath.abs(coordinates[2] - 3) &lt; 1.0e-12);
+<a class="jxr_linenumber" name="20" href="#20">20</a>  <strong class="jxr_keyword">import</strong> java.text.DecimalFormat;
+<a class="jxr_linenumber" name="21" href="#21">21</a>  <strong class="jxr_keyword">import</strong> java.text.DecimalFormatSymbols;
+<a class="jxr_linenumber" name="22" href="#22">22</a>  <strong class="jxr_keyword">import</strong> java.text.NumberFormat;
+<a class="jxr_linenumber" name="23" href="#23">23</a>  <strong class="jxr_keyword">import</strong> java.util.Locale;
+<a class="jxr_linenumber" name="24" href="#24">24</a>  
+<a class="jxr_linenumber" name="25" href="#25">25</a>  <strong class="jxr_keyword">import</strong> org.apache.commons.math3.TestUtils;
+<a class="jxr_linenumber" name="26" href="#26">26</a>  <strong class="jxr_keyword">import</strong> org.apache.commons.math3.exception.DimensionMismatchException;
+<a class="jxr_linenumber" name="27" href="#27">27</a>  <strong class="jxr_keyword">import</strong> org.apache.commons.math3.exception.MathArithmeticException;
+<a class="jxr_linenumber" name="28" href="#28">28</a>  <strong class="jxr_keyword">import</strong> org.apache.commons.math3.geometry.Space;
+<a class="jxr_linenumber" name="29" href="#29">29</a>  <strong class="jxr_keyword">import</strong> org.apache.commons.math3.random.Well1024a;
+<a class="jxr_linenumber" name="30" href="#30">30</a>  <strong class="jxr_keyword">import</strong> org.apache.commons.math3.util.FastMath;
+<a class="jxr_linenumber" name="31" href="#31">31</a>  <strong class="jxr_keyword">import</strong> org.apache.commons.math3.util.Precision;
+<a class="jxr_linenumber" name="32" href="#32">32</a>  <strong class="jxr_keyword">import</strong> org.junit.Assert;
+<a class="jxr_linenumber" name="33" href="#33">33</a>  <strong class="jxr_keyword">import</strong> org.junit.Test;
+<a class="jxr_linenumber" name="34" href="#34">34</a>  
+<a class="jxr_linenumber" name="35" href="#35">35</a>  <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">class</strong> <a href="../../../../../../../org/apache/commons/math3/geometry/euclidean/threed/Vector3DTest.html">Vector3DTest</a> {
+<a class="jxr_linenumber" name="36" href="#36">36</a>      @Test
+<a class="jxr_linenumber" name="37" href="#37">37</a>      <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testConstructors() <strong class="jxr_keyword">throws</strong> DimensionMismatchException {
+<a class="jxr_linenumber" name="38" href="#38">38</a>          <strong class="jxr_keyword">double</strong> r = FastMath.sqrt(2) /2;
+<a class="jxr_linenumber" name="39" href="#39">39</a>          checkVector(<strong class="jxr_keyword">new</strong> Vector3D(2, <strong class="jxr_keyword">new</strong> Vector3D(FastMath.PI / 3, -FastMath.PI / 4)),
+<a class="jxr_linenumber" name="40" href="#40">40</a>                      r, r * FastMath.sqrt(3), -2 * r);
+<a class="jxr_linenumber" name="41" href="#41">41</a>          checkVector(<strong class="jxr_keyword">new</strong> Vector3D(2, Vector3D.PLUS_I,
+<a class="jxr_linenumber" name="42" href="#42">42</a>                                   -3, Vector3D.MINUS_K),
+<a class="jxr_linenumber" name="43" href="#43">43</a>                      2, 0, 3);
+<a class="jxr_linenumber" name="44" href="#44">44</a>          checkVector(<strong class="jxr_keyword">new</strong> Vector3D(2, Vector3D.PLUS_I,
+<a class="jxr_linenumber" name="45" href="#45">45</a>                                   5, Vector3D.PLUS_J,
+<a class="jxr_linenumber" name="46" href="#46">46</a>                                   -3, Vector3D.MINUS_K),
+<a class="jxr_linenumber" name="47" href="#47">47</a>                      2, 5, 3);
+<a class="jxr_linenumber" name="48" href="#48">48</a>          checkVector(<strong class="jxr_keyword">new</strong> Vector3D(2, Vector3D.PLUS_I,
+<a class="jxr_linenumber" name="49" href="#49">49</a>                                   5, Vector3D.PLUS_J,
+<a class="jxr_linenumber" name="50" href="#50">50</a>                                   5, Vector3D.MINUS_J,
+<a class="jxr_linenumber" name="51" href="#51">51</a>                                   -3, Vector3D.MINUS_K),
+<a class="jxr_linenumber" name="52" href="#52">52</a>                      2, 0, 3);
+<a class="jxr_linenumber" name="53" href="#53">53</a>          checkVector(<strong class="jxr_keyword">new</strong> Vector3D(<strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[] { 2,  5,  -3 }),
+<a class="jxr_linenumber" name="54" href="#54">54</a>                      2, 5, -3);
+<a class="jxr_linenumber" name="55" href="#55">55</a>      }
+<a class="jxr_linenumber" name="56" href="#56">56</a>  
+<a class="jxr_linenumber" name="57" href="#57">57</a>      @Test
+<a class="jxr_linenumber" name="58" href="#58">58</a>      <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testSpace() {
+<a class="jxr_linenumber" name="59" href="#59">59</a>          Space space = <strong class="jxr_keyword">new</strong> Vector3D(1, 2, 2).getSpace();
+<a class="jxr_linenumber" name="60" href="#60">60</a>          Assert.assertEquals(3, space.getDimension());
+<a class="jxr_linenumber" name="61" href="#61">61</a>          Assert.assertEquals(2, space.getSubSpace().getDimension());
+<a class="jxr_linenumber" name="62" href="#62">62</a>          Space deserialized = (Space) TestUtils.serializeAndRecover(space);
+<a class="jxr_linenumber" name="63" href="#63">63</a>          Assert.assertTrue(space == deserialized);
 <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>      @Test
-<a class="jxr_linenumber" name="67" href="#67">67</a>      <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testNorm1() {
-<a class="jxr_linenumber" name="68" href="#68">68</a>          Assert.assertEquals(0.0, Vector3D.ZERO.getNorm1(), 0);
-<a class="jxr_linenumber" name="69" href="#69">69</a>          Assert.assertEquals(6.0, <strong class="jxr_keyword">new</strong> Vector3D(1, -2, 3).getNorm1(), 0);
-<a class="jxr_linenumber" name="70" href="#70">70</a>      }
-<a class="jxr_linenumber" name="71" href="#71">71</a>  
-<a class="jxr_linenumber" name="72" href="#72">72</a>      @Test
-<a class="jxr_linenumber" name="73" href="#73">73</a>      <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testNorm() {
-<a class="jxr_linenumber" name="74" href="#74">74</a>          Assert.assertEquals(0.0, Vector3D.ZERO.getNorm(), 0);
-<a class="jxr_linenumber" name="75" href="#75">75</a>          Assert.assertEquals(FastMath.sqrt(14), <strong class="jxr_keyword">new</strong> Vector3D(1, 2, 3).getNorm(), 1.0e-12);
-<a class="jxr_linenumber" name="76" href="#76">76</a>      }
-<a class="jxr_linenumber" name="77" href="#77">77</a>  
-<a class="jxr_linenumber" name="78" href="#78">78</a>      @Test
-<a class="jxr_linenumber" name="79" href="#79">79</a>      <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testNormInf() {
-<a class="jxr_linenumber" name="80" href="#80">80</a>          Assert.assertEquals(0.0, Vector3D.ZERO.getNormInf(), 0);
-<a class="jxr_linenumber" name="81" href="#81">81</a>          Assert.assertEquals(3.0, <strong class="jxr_keyword">new</strong> Vector3D(1, -2, 3).getNormInf(), 0);
+<a class="jxr_linenumber" name="67" href="#67">67</a>      <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testZero() {
+<a class="jxr_linenumber" name="68" href="#68">68</a>          Assert.assertEquals(0, <strong class="jxr_keyword">new</strong> Vector3D(1, 2, 2).getZero().getNorm(), 1.0e-15);
+<a class="jxr_linenumber" name="69" href="#69">69</a>      }
+<a class="jxr_linenumber" name="70" href="#70">70</a>  
+<a class="jxr_linenumber" name="71" href="#71">71</a>      @Test
+<a class="jxr_linenumber" name="72" href="#72">72</a>      <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testEquals() {
+<a class="jxr_linenumber" name="73" href="#73">73</a>          Vector3D u1 = <strong class="jxr_keyword">new</strong> Vector3D(1, 2, 3);
+<a class="jxr_linenumber" name="74" href="#74">74</a>          Vector3D u2 = <strong class="jxr_keyword">new</strong> Vector3D(1, 2, 3);
+<a class="jxr_linenumber" name="75" href="#75">75</a>          Assert.assertTrue(u1.equals(u1));
+<a class="jxr_linenumber" name="76" href="#76">76</a>          Assert.assertTrue(u1.equals(u2));
+<a class="jxr_linenumber" name="77" href="#77">77</a>          Assert.assertFalse(u1.equals(<strong class="jxr_keyword">new</strong> Rotation(1, 0, 0, 0, false)));
+<a class="jxr_linenumber" name="78" href="#78">78</a>          Assert.assertFalse(u1.equals(<strong class="jxr_keyword">new</strong> Vector3D(1, 2, 3 + 10 * Precision.EPSILON)));
+<a class="jxr_linenumber" name="79" href="#79">79</a>          Assert.assertFalse(u1.equals(<strong class="jxr_keyword">new</strong> Vector3D(1, 2 + 10 * Precision.EPSILON, 3)));
+<a class="jxr_linenumber" name="80" href="#80">80</a>          Assert.assertFalse(u1.equals(<strong class="jxr_keyword">new</strong> Vector3D(1 + 10 * Precision.EPSILON, 2, 3)));
+<a class="jxr_linenumber" name="81" href="#81">81</a>          Assert.assertTrue(<strong class="jxr_keyword">new</strong> Vector3D(0, Double.NaN, 0).equals(<strong class="jxr_keyword">new</strong> Vector3D(0, 0, Double.NaN)));
 <a class="jxr_linenumber" name="82" href="#82">82</a>      }
 <a class="jxr_linenumber" name="83" href="#83">83</a>  
 <a class="jxr_linenumber" name="84" href="#84">84</a>      @Test
-<a class="jxr_linenumber" name="85" href="#85">85</a>      <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testDistance1() {
-<a class="jxr_linenumber" name="86" href="#86">86</a>          Vector3D v1 = <strong class="jxr_keyword">new</strong> Vector3D(1, -2, 3);
-<a class="jxr_linenumber" name="87" href="#87">87</a>          Vector3D v2 = <strong class="jxr_keyword">new</strong> Vector3D(-4, 2, 0);
-<a class="jxr_linenumber" name="88" href="#88">88</a>          Assert.assertEquals(0.0, Vector3D.distance1(Vector3D.MINUS_I, Vector3D.MINUS_I), 0);
-<a class="jxr_linenumber" name="89" href="#89">89</a>          Assert.assertEquals(12.0, Vector3D.distance1(v1, v2), 1.0e-12);
-<a class="jxr_linenumber" name="90" href="#90">90</a>          Assert.assertEquals(v1.subtract(v2).getNorm1(), Vector3D.distance1(v1, v2), 1.0e-12);
-<a class="jxr_linenumber" name="91" href="#91">91</a>      }
-<a class="jxr_linenumber" name="92" href="#92">92</a>  
-<a class="jxr_linenumber" name="93" href="#93">93</a>      @Test
-<a class="jxr_linenumber" name="94" href="#94">94</a>      <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testDistance() {
-<a class="jxr_linenumber" name="95" href="#95">95</a>          Vector3D v1 = <strong class="jxr_keyword">new</strong> Vector3D(1, -2, 3);
-<a class="jxr_linenumber" name="96" href="#96">96</a>          Vector3D v2 = <strong class="jxr_keyword">new</strong> Vector3D(-4, 2, 0);
-<a class="jxr_linenumber" name="97" href="#97">97</a>          Assert.assertEquals(0.0, Vector3D.distance(Vector3D.MINUS_I, Vector3D.MINUS_I), 0);
-<a class="jxr_linenumber" name="98" href="#98">98</a>          Assert.assertEquals(FastMath.sqrt(50), Vector3D.distance(v1, v2), 1.0e-12);
-<a class="jxr_linenumber" name="99" href="#99">99</a>          Assert.assertEquals(v1.subtract(v2).getNorm(), Vector3D.distance(v1, v2), 1.0e-12);
-<a class="jxr_linenumber" name="100" href="#100">100</a>     }
-<a class="jxr_linenumber" name="101" href="#101">101</a> 
-<a class="jxr_linenumber" name="102" href="#102">102</a>     @Test
-<a class="jxr_linenumber" name="103" href="#103">103</a>     <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testDistanceSq() {
-<a class="jxr_linenumber" name="104" href="#104">104</a>         Vector3D v1 = <strong class="jxr_keyword">new</strong> Vector3D(1, -2, 3);
-<a class="jxr_linenumber" name="105" href="#105">105</a>         Vector3D v2 = <strong class="jxr_keyword">new</strong> Vector3D(-4, 2, 0);
-<a class="jxr_linenumber" name="106" href="#106">106</a>         Assert.assertEquals(0.0, Vector3D.distanceSq(Vector3D.MINUS_I, Vector3D.MINUS_I), 0);
-<a class="jxr_linenumber" name="107" href="#107">107</a>         Assert.assertEquals(50.0, Vector3D.distanceSq(v1, v2), 1.0e-12);
-<a class="jxr_linenumber" name="108" href="#108">108</a>         Assert.assertEquals(Vector3D.distance(v1, v2) * Vector3D.distance(v1, v2),
-<a class="jxr_linenumber" name="109" href="#109">109</a>                             Vector3D.distanceSq(v1, v2), 1.0e-12);
-<a class="jxr_linenumber" name="110" href="#110">110</a>   }
-<a class="jxr_linenumber" name="111" href="#111">111</a> 
-<a class="jxr_linenumber" name="112" href="#112">112</a>     @Test
-<a class="jxr_linenumber" name="113" href="#113">113</a>     <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testDistanceInf() {
-<a class="jxr_linenumber" name="114" href="#114">114</a>         Vector3D v1 = <strong class="jxr_keyword">new</strong> Vector3D(1, -2, 3);
-<a class="jxr_linenumber" name="115" href="#115">115</a>         Vector3D v2 = <strong class="jxr_keyword">new</strong> Vector3D(-4, 2, 0);
-<a class="jxr_linenumber" name="116" href="#116">116</a>         Assert.assertEquals(0.0, Vector3D.distanceInf(Vector3D.MINUS_I, Vector3D.MINUS_I), 0);
-<a class="jxr_linenumber" name="117" href="#117">117</a>         Assert.assertEquals(5.0, Vector3D.distanceInf(v1, v2), 1.0e-12);
-<a class="jxr_linenumber" name="118" href="#118">118</a>         Assert.assertEquals(v1.subtract(v2).getNormInf(), Vector3D.distanceInf(v1, v2), 1.0e-12);
-<a class="jxr_linenumber" name="119" href="#119">119</a>     }
-<a class="jxr_linenumber" name="120" href="#120">120</a> 
-<a class="jxr_linenumber" name="121" href="#121">121</a>     @Test
-<a class="jxr_linenumber" name="122" href="#122">122</a>     <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testSubtract() {
-<a class="jxr_linenumber" name="123" href="#123">123</a>         Vector3D v1 = <strong class="jxr_keyword">new</strong> Vector3D(1, 2, 3);
-<a class="jxr_linenumber" name="124" href="#124">124</a>         Vector3D v2 = <strong class="jxr_keyword">new</strong> Vector3D(-3, -2, -1);
-<a class="jxr_linenumber" name="125" href="#125">125</a>         v1 = v1.subtract(v2);
-<a class="jxr_linenumber" name="126" href="#126">126</a>         checkVector(v1, 4, 4, 4);
-<a class="jxr_linenumber" name="127" href="#127">127</a> 
-<a class="jxr_linenumber" name="128" href="#128">128</a>         checkVector(v2.subtract(v1), -7, -6, -5);
-<a class="jxr_linenumber" name="129" href="#129">129</a>         checkVector(v2.subtract(3, v1), -15, -14, -13);
-<a class="jxr_linenumber" name="130" href="#130">130</a>     }
-<a class="jxr_linenumber" name="131" href="#131">131</a> 
-<a class="jxr_linenumber" name="132" href="#132">132</a>     @Test
-<a class="jxr_linenumber" name="133" href="#133">133</a>     <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testAdd() {
-<a class="jxr_linenumber" name="134" href="#134">134</a>         Vector3D v1 = <strong class="jxr_keyword">new</strong> Vector3D(1, 2, 3);
-<a class="jxr_linenumber" name="135" href="#135">135</a>         Vector3D v2 = <strong class="jxr_keyword">new</strong> Vector3D(-3, -2, -1);
-<a class="jxr_linenumber" name="136" href="#136">136</a>         v1 = v1.add(v2);
-<a class="jxr_linenumber" name="137" href="#137">137</a>         checkVector(v1, -2, 0, 2);
-<a class="jxr_linenumber" name="138" href="#138">138</a> 
-<a class="jxr_linenumber" name="139" href="#139">139</a>         checkVector(v2.add(v1), -5, -2, 1);
-<a class="jxr_linenumber" name="140" href="#140">140</a>         checkVector(v2.add(3, v1), -9, -2, 5);
-<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>     @Test
-<a class="jxr_linenumber" name="144" href="#144">144</a>     <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testScalarProduct() {
-<a class="jxr_linenumber" name="145" href="#145">145</a>         Vector3D v = <strong class="jxr_keyword">new</strong> Vector3D(1, 2, 3);
-<a class="jxr_linenumber" name="146" href="#146">146</a>         v = v.scalarMultiply(3);
-<a class="jxr_linenumber" name="147" href="#147">147</a>         checkVector(v, 3, 6, 9);
-<a class="jxr_linenumber" name="148" href="#148">148</a> 
-<a class="jxr_linenumber" name="149" href="#149">149</a>         checkVector(v.scalarMultiply(0.5), 1.5, 3, 4.5);
+<a class="jxr_linenumber" name="85" href="#85">85</a>      <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testHash() {
+<a class="jxr_linenumber" name="86" href="#86">86</a>          Assert.assertEquals(<strong class="jxr_keyword">new</strong> Vector3D(0, Double.NaN, 0).hashCode(), <strong class="jxr_keyword">new</strong> Vector3D(0, 0, Double.NaN).hashCode());
+<a class="jxr_linenumber" name="87" href="#87">87</a>          Vector3D u = <strong class="jxr_keyword">new</strong> Vector3D(1, 2, 3);
+<a class="jxr_linenumber" name="88" href="#88">88</a>          Vector3D v = <strong class="jxr_keyword">new</strong> Vector3D(1, 2, 3 + 10 * Precision.EPSILON);
+<a class="jxr_linenumber" name="89" href="#89">89</a>          Assert.assertTrue(u.hashCode() != v.hashCode());
+<a class="jxr_linenumber" name="90" href="#90">90</a>      }
+<a class="jxr_linenumber" name="91" href="#91">91</a>  
+<a class="jxr_linenumber" name="92" href="#92">92</a>      @Test
+<a class="jxr_linenumber" name="93" href="#93">93</a>      <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testInfinite() {
+<a class="jxr_linenumber" name="94" href="#94">94</a>          Assert.assertTrue(<strong class="jxr_keyword">new</strong> Vector3D(1, 1, Double.NEGATIVE_INFINITY).isInfinite());
+<a class="jxr_linenumber" name="95" href="#95">95</a>          Assert.assertTrue(<strong class="jxr_keyword">new</strong> Vector3D(1, Double.NEGATIVE_INFINITY, 1).isInfinite());
+<a class="jxr_linenumber" name="96" href="#96">96</a>          Assert.assertTrue(<strong class="jxr_keyword">new</strong> Vector3D(Double.NEGATIVE_INFINITY, 1, 1).isInfinite());
+<a class="jxr_linenumber" name="97" href="#97">97</a>          Assert.assertFalse(<strong class="jxr_keyword">new</strong> Vector3D(1, 1, 2).isInfinite());
+<a class="jxr_linenumber" name="98" href="#98">98</a>          Assert.assertFalse(<strong class="jxr_keyword">new</strong> Vector3D(1, Double.NaN, Double.NEGATIVE_INFINITY).isInfinite());
+<a class="jxr_linenumber" name="99" href="#99">99</a>      }
+<a class="jxr_linenumber" name="100" href="#100">100</a> 
+<a class="jxr_linenumber" name="101" href="#101">101</a>     @Test
+<a class="jxr_linenumber" name="102" href="#102">102</a>     <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testNaN() {
+<a class="jxr_linenumber" name="103" href="#103">103</a>         Assert.assertTrue(<strong class="jxr_keyword">new</strong> Vector3D(1, 1, Double.NaN).isNaN());
+<a class="jxr_linenumber" name="104" href="#104">104</a>         Assert.assertTrue(<strong class="jxr_keyword">new</strong> Vector3D(1, Double.NaN, 1).isNaN());
+<a class="jxr_linenumber" name="105" href="#105">105</a>         Assert.assertTrue(<strong class="jxr_keyword">new</strong> Vector3D(Double.NaN, 1, 1).isNaN());
+<a class="jxr_linenumber" name="106" href="#106">106</a>         Assert.assertFalse(<strong class="jxr_keyword">new</strong> Vector3D(1, 1, 2).isNaN());
+<a class="jxr_linenumber" name="107" href="#107">107</a>         Assert.assertFalse(<strong class="jxr_keyword">new</strong> Vector3D(1, 1, Double.NEGATIVE_INFINITY).isNaN());
+<a class="jxr_linenumber" name="108" href="#108">108</a>     }
+<a class="jxr_linenumber" name="109" href="#109">109</a> 
+<a class="jxr_linenumber" name="110" href="#110">110</a>     @Test
+<a class="jxr_linenumber" name="111" href="#111">111</a>     <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testToString() {
+<a class="jxr_linenumber" name="112" href="#112">112</a>         Assert.assertEquals(<span class="jxr_string">"{3; 2; 1}"</span>, <strong class="jxr_keyword">new</strong> Vector3D(3, 2, 1).toString());
+<a class="jxr_linenumber" name="113" href="#113">113</a>         NumberFormat format = <strong class="jxr_keyword">new</strong> DecimalFormat(<span class="jxr_string">"0.000"</span>, <strong class="jxr_keyword">new</strong> DecimalFormatSymbols(Locale.US));
+<a class="jxr_linenumber" name="114" href="#114">114</a>         Assert.assertEquals(<span class="jxr_string">"{3.000; 2.000; 1.000}"</span>, <strong class="jxr_keyword">new</strong> Vector3D(3, 2, 1).toString(format));
+<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>     @Test(expected=DimensionMismatchException.<strong class="jxr_keyword">class</strong>)
+<a class="jxr_linenumber" name="118" href="#118">118</a>     <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testWrongDimension() <strong class="jxr_keyword">throws</strong> DimensionMismatchException {
+<a class="jxr_linenumber" name="119" href="#119">119</a>         <strong class="jxr_keyword">new</strong> Vector3D(<strong class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[] { 2,  5 });
+<a class="jxr_linenumber" name="120" href="#120">120</a>     }
+<a class="jxr_linenumber" name="121" href="#121">121</a> 
+<a class="jxr_linenumber" name="122" href="#122">122</a>     @Test
+<a class="jxr_linenumber" name="123" href="#123">123</a>     <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testCoordinates() {
+<a class="jxr_linenumber" name="124" href="#124">124</a>         Vector3D v = <strong class="jxr_keyword">new</strong> Vector3D(1, 2, 3);
+<a class="jxr_linenumber" name="125" href="#125">125</a>         Assert.assertTrue(FastMath.abs(v.getX() - 1) &lt; 1.0e-12);
+<a class="jxr_linenumber" name="126" href="#126">126</a>         Assert.assertTrue(FastMath.abs(v.getY() - 2) &lt; 1.0e-12);
+<a class="jxr_linenumber" name="127" href="#127">127</a>         Assert.assertTrue(FastMath.abs(v.getZ() - 3) &lt; 1.0e-12);
+<a class="jxr_linenumber" name="128" href="#128">128</a>         <strong class="jxr_keyword">double</strong>[] coordinates = v.toArray();
+<a class="jxr_linenumber" name="129" href="#129">129</a>         Assert.assertTrue(FastMath.abs(coordinates[0] - 1) &lt; 1.0e-12);
+<a class="jxr_linenumber" name="130" href="#130">130</a>         Assert.assertTrue(FastMath.abs(coordinates[1] - 2) &lt; 1.0e-12);
+<a class="jxr_linenumber" name="131" href="#131">131</a>         Assert.assertTrue(FastMath.abs(coordinates[2] - 3) &lt; 1.0e-12);
+<a class="jxr_linenumber" name="132" href="#132">132</a>     }
+<a class="jxr_linenumber" name="133" href="#133">133</a> 
+<a class="jxr_linenumber" name="134" href="#134">134</a>     @Test
+<a class="jxr_linenumber" name="135" href="#135">135</a>     <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testNorm1() {
+<a class="jxr_linenumber" name="136" href="#136">136</a>         Assert.assertEquals(0.0, Vector3D.ZERO.getNorm1(), 0);
+<a class="jxr_linenumber" name="137" href="#137">137</a>         Assert.assertEquals(6.0, <strong class="jxr_keyword">new</strong> Vector3D(1, -2, 3).getNorm1(), 0);
+<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>     @Test
+<a class="jxr_linenumber" name="141" href="#141">141</a>     <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testNorm() {
+<a class="jxr_linenumber" name="142" href="#142">142</a>         Assert.assertEquals(0.0, Vector3D.ZERO.getNorm(), 0);
+<a class="jxr_linenumber" name="143" href="#143">143</a>         Assert.assertEquals(FastMath.sqrt(14), <strong class="jxr_keyword">new</strong> Vector3D(1, 2, 3).getNorm(), 1.0e-12);
+<a class="jxr_linenumber" name="144" href="#144">144</a>     }
+<a class="jxr_linenumber" name="145" href="#145">145</a> 
+<a class="jxr_linenumber" name="146" href="#146">146</a>     @Test
+<a class="jxr_linenumber" name="147" href="#147">147</a>     <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testNormSq() {
+<a class="jxr_linenumber" name="148" href="#148">148</a>         Assert.assertEquals(0.0, <strong class="jxr_keyword">new</strong> Vector3D(0, 0, 0).getNormSq(), 0);
+<a class="jxr_linenumber" name="149" href="#149">149</a>         Assert.assertEquals(14, <strong class="jxr_keyword">new</strong> Vector3D(1, 2, 3).getNormSq(), 1.0e-12);
 <a class="jxr_linenumber" name="150" href="#150">150</a>     }
 <a class="jxr_linenumber" name="151" href="#151">151</a> 
 <a class="jxr_linenumber" name="152" href="#152">152</a>     @Test
-<a class="jxr_linenumber" name="153" href="#153">153</a>     <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testVectorialProducts() {
-<a class="jxr_linenumber" name="154" href="#154">154</a>         Vector3D v1 = <strong class="jxr_keyword">new</strong> Vector3D(2, 1, -4);
-<a class="jxr_linenumber" name="155" href="#155">155</a>         Vector3D v2 = <strong class="jxr_keyword">new</strong> Vector3D(3, 1, -1);
-<a class="jxr_linenumber" name="156" href="#156">156</a> 
-<a class="jxr_linenumber" name="157" href="#157">157</a>         Assert.assertTrue(FastMath.abs(Vector3D.dotProduct(v1, v2) - 11) &lt; 1.0e-12);
-<a class="jxr_linenumber" name="158" href="#158">158</a> 
-<a class="jxr_linenumber" name="159" href="#159">159</a>         Vector3D v3 = Vector3D.crossProduct(v1, v2);
-<a class="jxr_linenumber" name="160" href="#160">160</a>         checkVector(v3, 3, -10, -1);
-<a class="jxr_linenumber" name="161" href="#161">161</a> 
-<a class="jxr_linenumber" name="162" href="#162">162</a>         Assert.assertTrue(FastMath.abs(Vector3D.dotProduct(v1, v3)) &lt; 1.0e-12);
-<a class="jxr_linenumber" name="163" href="#163">163</a>         Assert.assertTrue(FastMath.abs(Vector3D.dotProduct(v2, v3)) &lt; 1.0e-12);
-<a class="jxr_linenumber" name="164" href="#164">164</a>     }
-<a class="jxr_linenumber" name="165" href="#165">165</a> 
-<a class="jxr_linenumber" name="166" href="#166">166</a>     @Test
-<a class="jxr_linenumber" name="167" href="#167">167</a>     <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testCrossProductCancellation() {
-<a class="jxr_linenumber" name="168" href="#168">168</a>         Vector3D v1 = <strong class="jxr_keyword">new</strong> Vector3D(9070467121.0, 4535233560.0, 1);
-<a class="jxr_linenumber" name="169" href="#169">169</a>         Vector3D v2 = <strong class="jxr_keyword">new</strong> Vector3D(9070467123.0, 4535233561.0, 1);
-<a class="jxr_linenumber" name="170" href="#170">170</a>         checkVector(Vector3D.crossProduct(v1, v2), -1, 2, 1);
-<a class="jxr_linenumber" name="171" href="#171">171</a> 
-<a class="jxr_linenumber" name="172" href="#172">172</a>         <strong class="jxr_keyword">double</strong> scale    = FastMath.scalb(1.0, 100);
-<a class="jxr_linenumber" name="173" href="#173">173</a>         Vector3D big1   = <strong class="jxr_keyword">new</strong> Vector3D(scale, v1);
-<a class="jxr_linenumber" name="174" href="#174">174</a>         Vector3D small2 = <strong class="jxr_keyword">new</strong> Vector3D(1 / scale, v2);
-<a class="jxr_linenumber" name="175" href="#175">175</a>         checkVector(Vector3D.crossProduct(big1, small2), -1, 2, 1);
-<a class="jxr_linenumber" name="176" href="#176">176</a> 
-<a class="jxr_linenumber" name="177" href="#177">177</a>     }
-<a class="jxr_linenumber" name="178" href="#178">178</a> 
-<a class="jxr_linenumber" name="179" href="#179">179</a>     @Test
-<a class="jxr_linenumber" name="180" href="#180">180</a>     <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testAngular() {
-<a class="jxr_linenumber" name="181" href="#181">181</a>         Assert.assertEquals(0,           Vector3D.PLUS_I.getAlpha(), 1.0e-10);
-<a class="jxr_linenumber" name="182" href="#182">182</a>         Assert.assertEquals(0,           Vector3D.PLUS_I.getDelta(), 1.0e-10);
-<a class="jxr_linenumber" name="183" href="#183">183</a>         Assert.assertEquals(FastMath.PI / 2, Vector3D.PLUS_J.getAlpha(), 1.0e-10);
-<a class="jxr_linenumber" name="184" href="#184">184</a>         Assert.assertEquals(0,           Vector3D.PLUS_J.getDelta(), 1.0e-10);
-<a class="jxr_linenumber" name="185" href="#185">185</a>         Assert.assertEquals(0,           Vector3D.PLUS_K.getAlpha(), 1.0e-10);
-<a class="jxr_linenumber" name="186" href="#186">186</a>         Assert.assertEquals(FastMath.PI / 2, Vector3D.PLUS_K.getDelta(), 1.0e-10);
-<a class="jxr_linenumber" name="187" href="#187">187</a>       
-<a class="jxr_linenumber" name="188" href="#188">188</a>         Vector3D u = <strong class="jxr_keyword">new</strong> Vector3D(-1, 1, -1);
-<a class="jxr_linenumber" name="189" href="#189">189</a>         Assert.assertEquals(3 * FastMath.PI /4, u.getAlpha(), 1.0e-10);
-<a class="jxr_linenumber" name="190" href="#190">190</a>         Assert.assertEquals(-1.0 / FastMath.sqrt(3), FastMath.sin(u.getDelta()), 1.0e-10);
-<a class="jxr_linenumber" name="191" href="#191">191</a>     }
-<a class="jxr_linenumber" name="192" href="#192">192</a> 
-<a class="jxr_linenumber" name="193" href="#193">193</a>     @Test
-<a class="jxr_linenumber" name="194" href="#194">194</a>     <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testAngularSeparation() <strong class="jxr_keyword">throws</strong> MathArithmeticException {
-<a class="jxr_linenumber" name="195" href="#195">195</a>         Vector3D v1 = <strong class="jxr_keyword">new</strong> Vector3D(2, -1, 4);
-<a class="jxr_linenumber" name="196" href="#196">196</a> 
-<a class="jxr_linenumber" name="197" href="#197">197</a>         Vector3D  k = v1.normalize();
-<a class="jxr_linenumber" name="198" href="#198">198</a>         Vector3D  i = k.orthogonal();
-<a class="jxr_linenumber" name="199" href="#199">199</a>         Vector3D v2 = k.scalarMultiply(FastMath.cos(1.2)).add(i.scalarMultiply(FastMath.sin(1.2)));
-<a class="jxr_linenumber" name="200" href="#200">200</a> 
-<a class="jxr_linenumber" name="201" href="#201">201</a>         Assert.assertTrue(FastMath.abs(Vector3D.angle(v1, v2) - 1.2) &lt; 1.0e-12);
-<a class="jxr_linenumber" name="202" href="#202">202</a>   }
-<a class="jxr_linenumber" name="203" href="#203">203</a> 
-<a class="jxr_linenumber" name="204" href="#204">204</a>     @Test
-<a class="jxr_linenumber" name="205" href="#205">205</a>     <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testNormalize() <strong class="jxr_keyword">throws</strong> MathArithmeticException {
-<a class="jxr_linenumber" name="206" href="#206">206</a>         Assert.assertEquals(1.0, <strong class="jxr_keyword">new</strong> Vector3D(5, -4, 2).normalize().getNorm(), 1.0e-12);
-<a class="jxr_linenumber" name="207" href="#207">207</a>         <strong class="jxr_keyword">try</strong> {
-<a class="jxr_linenumber" name="208" href="#208">208</a>             Vector3D.ZERO.normalize();
-<a class="jxr_linenumber" name="209" href="#209">209</a>             Assert.fail(<span class="jxr_string">"an exception should have been thrown"</span>);
-<a class="jxr_linenumber" name="210" href="#210">210</a>         } <strong class="jxr_keyword">catch</strong> (MathArithmeticException ae) {
-<a class="jxr_linenumber" name="211" href="#211">211</a>             <em class="jxr_comment">// expected behavior</em>
-<a class="jxr_linenumber" name="212" href="#212">212</a>         }
-<a class="jxr_linenumber" name="213" href="#213">213</a>     }
-<a class="jxr_linenumber" name="214" href="#214">214</a> 
-<a class="jxr_linenumber" name="215" href="#215">215</a>     @Test
-<a class="jxr_linenumber" name="216" href="#216">216</a>     <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testOrthogonal() <strong class="jxr_keyword">throws</strong> MathArithmeticException {
-<a class="jxr_linenumber" name="217" href="#217">217</a>         Vector3D v1 = <strong class="jxr_keyword">new</strong> Vector3D(0.1, 2.5, 1.3);
-<a class="jxr_linenumber" name="218" href="#218">218</a>         Assert.assertEquals(0.0, Vector3D.dotProduct(v1, v1.orthogonal()), 1.0e-12);
-<a class="jxr_linenumber" name="219" href="#219">219</a>         Vector3D v2 = <strong class="jxr_keyword">new</strong> Vector3D(2.3, -0.003, 7.6);
-<a class="jxr_linenumber" name="220" href="#220">220</a>         Assert.assertEquals(0.0, Vector3D.dotProduct(v2, v2.orthogonal()), 1.0e-12);
-<a class="jxr_linenumber" name="221" href="#221">221</a>         Vector3D v3 = <strong class="jxr_keyword">new</strong> Vector3D(-1.7, 1.4, 0.2);
-<a class="jxr_linenumber" name="222" href="#222">222</a>         Assert.assertEquals(0.0, Vector3D.dotProduct(v3, v3.orthogonal()), 1.0e-12);
-<a class="jxr_linenumber" name="223" href="#223">223</a>         <strong class="jxr_keyword">try</strong> {
-<a class="jxr_linenumber" name="224" href="#224">224</a>             <strong class="jxr_keyword">new</strong> Vector3D(0, 0, 0).orthogonal();
-<a class="jxr_linenumber" name="225" href="#225">225</a>             Assert.fail(<span class="jxr_string">"an exception should have been thrown"</span>);
-<a class="jxr_linenumber" name="226" href="#226">226</a>         } <strong class="jxr_keyword">catch</strong> (MathArithmeticException ae) {
-<a class="jxr_linenumber" name="227" href="#227">227</a>             <em class="jxr_comment">// expected behavior</em>
-<a class="jxr_linenumber" name="228" href="#228">228</a>         }
-<a class="jxr_linenumber" name="229" href="#229">229</a>     }
+<a class="jxr_linenumber" name="153" href="#153">153</a>     <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testNormInf() {
+<a class="jxr_linenumber" name="154" href="#154">154</a>         Assert.assertEquals(0.0, Vector3D.ZERO.getNormInf(), 0);
+<a class="jxr_linenumber" name="155" href="#155">155</a>         Assert.assertEquals(3.0, <strong class="jxr_keyword">new</strong> Vector3D(1, -2, 3).getNormInf(), 0);
+<a class="jxr_linenumber" name="156" href="#156">156</a>     }
+<a class="jxr_linenumber" name="157" href="#157">157</a> 
+<a class="jxr_linenumber" name="158" href="#158">158</a>     @Test
+<a class="jxr_linenumber" name="159" href="#159">159</a>     <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testDistance1() {
+<a class="jxr_linenumber" name="160" href="#160">160</a>         Vector3D v1 = <strong class="jxr_keyword">new</strong> Vector3D(1, -2, 3);
+<a class="jxr_linenumber" name="161" href="#161">161</a>         Vector3D v2 = <strong class="jxr_keyword">new</strong> Vector3D(-4, 2, 0);
+<a class="jxr_linenumber" name="162" href="#162">162</a>         Assert.assertEquals(0.0, Vector3D.distance1(Vector3D.MINUS_I, Vector3D.MINUS_I), 0);
+<a class="jxr_linenumber" name="163" href="#163">163</a>         Assert.assertEquals(12.0, Vector3D.distance1(v1, v2), 1.0e-12);
+<a class="jxr_linenumber" name="164" href="#164">164</a>         Assert.assertEquals(v1.subtract(v2).getNorm1(), Vector3D.distance1(v1, v2), 1.0e-12);
+<a class="jxr_linenumber" name="165" href="#165">165</a>     }
+<a class="jxr_linenumber" name="166" href="#166">166</a> 
+<a class="jxr_linenumber" name="167" href="#167">167</a>     @Test
+<a class="jxr_linenumber" name="168" href="#168">168</a>     <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testDistance() {
+<a class="jxr_linenumber" name="169" href="#169">169</a>         Vector3D v1 = <strong class="jxr_keyword">new</strong> Vector3D(1, -2, 3);
+<a class="jxr_linenumber" name="170" href="#170">170</a>         Vector3D v2 = <strong class="jxr_keyword">new</strong> Vector3D(-4, 2, 0);
+<a class="jxr_linenumber" name="171" href="#171">171</a>         Assert.assertEquals(0.0, Vector3D.distance(Vector3D.MINUS_I, Vector3D.MINUS_I), 0);
+<a class="jxr_linenumber" name="172" href="#172">172</a>         Assert.assertEquals(FastMath.sqrt(50), Vector3D.distance(v1, v2), 1.0e-12);
+<a class="jxr_linenumber" name="173" href="#173">173</a>         Assert.assertEquals(v1.subtract(v2).getNorm(), Vector3D.distance(v1, v2), 1.0e-12);
+<a class="jxr_linenumber" name="174" href="#174">174</a>     }
+<a class="jxr_linenumber" name="175" href="#175">175</a> 
+<a class="jxr_linenumber" name="176" href="#176">176</a>     @Test
+<a class="jxr_linenumber" name="177" href="#177">177</a>     <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testDistanceSq() {
+<a class="jxr_linenumber" name="178" href="#178">178</a>         Vector3D v1 = <strong class="jxr_keyword">new</strong> Vector3D(1, -2, 3);
+<a class="jxr_linenumber" name="179" href="#179">179</a>         Vector3D v2 = <strong class="jxr_keyword">new</strong> Vector3D(-4, 2, 0);
+<a class="jxr_linenumber" name="180" href="#180">180</a>         Assert.assertEquals(0.0, Vector3D.distanceSq(Vector3D.MINUS_I, Vector3D.MINUS_I), 0);
+<a class="jxr_linenumber" name="181" href="#181">181</a>         Assert.assertEquals(50.0, Vector3D.distanceSq(v1, v2), 1.0e-12);
+<a class="jxr_linenumber" name="182" href="#182">182</a>         Assert.assertEquals(Vector3D.distance(v1, v2) * Vector3D.distance(v1, v2),
+<a class="jxr_linenumber" name="183" href="#183">183</a>                             Vector3D.distanceSq(v1, v2), 1.0e-12);
+<a class="jxr_linenumber" name="184" href="#184">184</a>   }
+<a class="jxr_linenumber" name="185" href="#185">185</a> 
+<a class="jxr_linenumber" name="186" href="#186">186</a>     @Test
+<a class="jxr_linenumber" name="187" href="#187">187</a>     <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testDistanceInf() {
+<a class="jxr_linenumber" name="188" href="#188">188</a>         Vector3D v1 = <strong class="jxr_keyword">new</strong> Vector3D(1, -2, 3);
+<a class="jxr_linenumber" name="189" href="#189">189</a>         Vector3D v2 = <strong class="jxr_keyword">new</strong> Vector3D(-4, 2, 0);
+<a class="jxr_linenumber" name="190" href="#190">190</a>         Assert.assertEquals(0.0, Vector3D.distanceInf(Vector3D.MINUS_I, Vector3D.MINUS_I), 0);
+<a class="jxr_linenumber" name="191" href="#191">191</a>         Assert.assertEquals(5.0, Vector3D.distanceInf(v1, v2), 1.0e-12);
+<a class="jxr_linenumber" name="192" href="#192">192</a>         Assert.assertEquals(v1.subtract(v2).getNormInf(), Vector3D.distanceInf(v1, v2), 1.0e-12);
+<a class="jxr_linenumber" name="193" href="#193">193</a>     }
+<a class="jxr_linenumber" name="194" href="#194">194</a> 
+<a class="jxr_linenumber" name="195" href="#195">195</a>     @Test
+<a class="jxr_linenumber" name="196" href="#196">196</a>     <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testSubtract() {
+<a class="jxr_linenumber" name="197" href="#197">197</a>         Vector3D v1 = <strong class="jxr_keyword">new</strong> Vector3D(1, 2, 3);
+<a class="jxr_linenumber" name="198" href="#198">198</a>         Vector3D v2 = <strong class="jxr_keyword">new</strong> Vector3D(-3, -2, -1);
+<a class="jxr_linenumber" name="199" href="#199">199</a>         v1 = v1.subtract(v2);
+<a class="jxr_linenumber" name="200" href="#200">200</a>         checkVector(v1, 4, 4, 4);
+<a class="jxr_linenumber" name="201" href="#201">201</a> 
+<a class="jxr_linenumber" name="202" href="#202">202</a>         checkVector(v2.subtract(v1), -7, -6, -5);
+<a class="jxr_linenumber" name="203" href="#203">203</a>         checkVector(v2.subtract(3, v1), -15, -14, -13);
+<a class="jxr_linenumber" name="204" href="#204">204</a>     }
+<a class="jxr_linenumber" name="205" href="#205">205</a> 
+<a class="jxr_linenumber" name="206" href="#206">206</a>     @Test
+<a class="jxr_linenumber" name="207" href="#207">207</a>     <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testAdd() {
+<a class="jxr_linenumber" name="208" href="#208">208</a>         Vector3D v1 = <strong class="jxr_keyword">new</strong> Vector3D(1, 2, 3);
+<a class="jxr_linenumber" name="209" href="#209">209</a>         Vector3D v2 = <strong class="jxr_keyword">new</strong> Vector3D(-3, -2, -1);
+<a class="jxr_linenumber" name="210" href="#210">210</a>         v1 = v1.add(v2);
+<a class="jxr_linenumber" name="211" href="#211">211</a>         checkVector(v1, -2, 0, 2);
+<a class="jxr_linenumber" name="212" href="#212">212</a> 
+<a class="jxr_linenumber" name="213" href="#213">213</a>         checkVector(v2.add(v1), -5, -2, 1);
+<a class="jxr_linenumber" name="214" href="#214">214</a>         checkVector(v2.add(3, v1), -9, -2, 5);
+<a class="jxr_linenumber" name="215" href="#215">215</a>     }
+<a class="jxr_linenumber" name="216" href="#216">216</a> 
+<a class="jxr_linenumber" name="217" href="#217">217</a>     @Test
+<a class="jxr_linenumber" name="218" href="#218">218</a>     <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testScalarProduct() {
+<a class="jxr_linenumber" name="219" href="#219">219</a>         Vector3D v = <strong class="jxr_keyword">new</strong> Vector3D(1, 2, 3);
+<a class="jxr_linenumber" name="220" href="#220">220</a>         v = v.scalarMultiply(3);
+<a class="jxr_linenumber" name="221" href="#221">221</a>         checkVector(v, 3, 6, 9);
+<a class="jxr_linenumber" name="222" href="#222">222</a> 
+<a class="jxr_linenumber" name="223" href="#223">223</a>         checkVector(v.scalarMultiply(0.5), 1.5, 3, 4.5);
+<a class="jxr_linenumber" name="224" href="#224">224</a>     }
+<a class="jxr_linenumber" name="225" href="#225">225</a> 
+<a class="jxr_linenumber" name="226" href="#226">226</a>     @Test
+<a class="jxr_linenumber" name="227" href="#227">227</a>     <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testVectorialProducts() {
+<a class="jxr_linenumber" name="228" href="#228">228</a>         Vector3D v1 = <strong class="jxr_keyword">new</strong> Vector3D(2, 1, -4);
+<a class="jxr_linenumber" name="229" href="#229">229</a>         Vector3D v2 = <strong class="jxr_keyword">new</strong> Vector3D(3, 1, -1);
 <a class="jxr_linenumber" name="230" href="#230">230</a> 
-<a class="jxr_linenumber" name="231" href="#231">231</a>     @Test
-<a class="jxr_linenumber" name="232" href="#232">232</a>     <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testAngle() <strong class="jxr_keyword">throws</strong> MathArithmeticException {
-<a class="jxr_linenumber" name="233" href="#233">233</a>         Assert.assertEquals(0.22572612855273393616,
-<a class="jxr_linenumber" name="234" href="#234">234</a>                             Vector3D.angle(<strong class="jxr_keyword">new</strong> Vector3D(1, 2, 3), <strong class="jxr_keyword">new</strong> Vector3D(4, 5, 6)),
-<a class="jxr_linenumber" name="235" href="#235">235</a>                             1.0e-12);
-<a class="jxr_linenumber" name="236" href="#236">236</a>         Assert.assertEquals(7.98595620686106654517199e-8,
-<a class="jxr_linenumber" name="237" href="#237">237</a>                             Vector3D.angle(<strong class="jxr_keyword">new</strong> Vector3D(1, 2, 3), <strong class="jxr_keyword">new</strong> Vector3D(2, 4, 6.000001)),
-<a class="jxr_linenumber" name="238" href="#238">238</a>                             1.0e-12);
-<a class="jxr_linenumber" name="239" href="#239">239</a>         Assert.assertEquals(3.14159257373023116985197793156,
-<a class="jxr_linenumber" name="240" href="#240">240</a>                             Vector3D.angle(<strong class="jxr_keyword">new</strong> Vector3D(1, 2, 3), <strong class="jxr_keyword">new</strong> Vector3D(-2, -4, -6.000001)),
-<a class="jxr_linenumber" name="241" href="#241">241</a>                             1.0e-12);
-<a class="jxr_linenumber" name="242" href="#242">242</a>         <strong class="jxr_keyword">try</strong> {
-<a class="jxr_linenumber" name="243" href="#243">243</a>             Vector3D.angle(Vector3D.ZERO, Vector3D.PLUS_I);
-<a class="jxr_linenumber" name="244" href="#244">244</a>             Assert.fail(<span class="jxr_string">"an exception should have been thrown"</span>);
-<a class="jxr_linenumber" name="245" href="#245">245</a>         } <strong class="jxr_keyword">catch</strong> (MathArithmeticException ae) {
-<a class="jxr_linenumber" name="246" href="#246">246</a>             <em class="jxr_comment">// expected behavior</em>
-<a class="jxr_linenumber" name="247" href="#247">247</a>         }
-<a class="jxr_linenumber" name="248" href="#248">248</a>     }
-<a class="jxr_linenumber" name="249" href="#249">249</a> 
-<a class="jxr_linenumber" name="250" href="#250">250</a>     @Test
-<a class="jxr_linenumber" name="251" href="#251">251</a>     <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testAccurateDotProduct() {
-<a class="jxr_linenumber" name="252" href="#252">252</a>         <em class="jxr_comment">// the following two vectors are nearly but not exactly orthogonal</em>
-<a class="jxr_linenumber" name="253" href="#253">253</a>         <em class="jxr_comment">// naive dot product (i.e. computing u1.x * u2.x + u1.y * u2.y + u1.z * u2.z</em>
-<a class="jxr_linenumber" name="254" href="#254">254</a>         <em class="jxr_comment">// leads to a result of 0.0, instead of the correct -1.855129...</em>
-<a class="jxr_linenumber" name="255" href="#255">255</a>         Vector3D u1 = <strong class="jxr_keyword">new</strong> Vector3D(-1321008684645961.0 /  268435456.0,
-<a class="jxr_linenumber" name="256" href="#256">256</a>                                    -5774608829631843.0 /  268435456.0,
-<a class="jxr_linenumber" name="257" href="#257">257</a>                                    -7645843051051357.0 / 8589934592.0);
-<a class="jxr_linenumber" name="258" href="#258">258</a>         Vector3D u2 = <strong class="jxr_keyword">new</strong> Vector3D(-5712344449280879.0 /    2097152.0,
-<a class="jxr_linenumber" name="259" href="#259">259</a>                                    -4550117129121957.0 /    2097152.0,
-<a class="jxr_linenumber" name="260" href="#260">260</a>                                     8846951984510141.0 /     131072.0);
-<a class="jxr_linenumber" name="261" href="#261">261</a>         <strong class="jxr_keyword">double</strong> sNaive = u1.getX() * u2.getX() + u1.getY() * u2.getY() + u1.getZ() * u2.getZ();
-<a class="jxr_linenumber" name="262" href="#262">262</a>         <strong class="jxr_keyword">double</strong> sAccurate = u1.dotProduct(u2);
-<a class="jxr_linenumber" name="263" href="#263">263</a>         Assert.assertEquals(0.0, sNaive, 1.0e-30);
-<a class="jxr_linenumber" name="264" href="#264">264</a>         Assert.assertEquals(-2088690039198397.0 / 1125899906842624.0, sAccurate, 1.0e-16);
+<a class="jxr_linenumber" name="231" href="#231">231</a>         Assert.assertTrue(FastMath.abs(Vector3D.dotProduct(v1, v2) - 11) &lt; 1.0e-12);
+<a class="jxr_linenumber" name="232" href="#232">232</a> 
+<a class="jxr_linenumber" name="233" href="#233">233</a>         Vector3D v3 = Vector3D.crossProduct(v1, v2);
+<a class="jxr_linenumber" name="234" href="#234">234</a>         checkVector(v3, 3, -10, -1);
+<a class="jxr_linenumber" name="235" href="#235">235</a> 
+<a class="jxr_linenumber" name="236" href="#236">236</a>         Assert.assertTrue(FastMath.abs(Vector3D.dotProduct(v1, v3)) &lt; 1.0e-12);
+<a class="jxr_linenumber" name="237" href="#237">237</a>         Assert.assertTrue(FastMath.abs(Vector3D.dotProduct(v2, v3)) &lt; 1.0e-12);
+<a class="jxr_linenumber" name="238" href="#238">238</a>     }
+<a class="jxr_linenumber" name="239" href="#239">239</a> 
+<a class="jxr_linenumber" name="240" href="#240">240</a>     @Test
+<a class="jxr_linenumber" name="241" href="#241">241</a>     <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testCrossProductCancellation() {
+<a class="jxr_linenumber" name="242" href="#242">242</a>         Vector3D v1 = <strong class="jxr_keyword">new</strong> Vector3D(9070467121.0, 4535233560.0, 1);
+<a class="jxr_linenumber" name="243" href="#243">243</a>         Vector3D v2 = <strong class="jxr_keyword">new</strong> Vector3D(9070467123.0, 4535233561.0, 1);
+<a class="jxr_linenumber" name="244" href="#244">244</a>         checkVector(Vector3D.crossProduct(v1, v2), -1, 2, 1);
+<a class="jxr_linenumber" name="245" href="#245">245</a> 
+<a class="jxr_linenumber" name="246" href="#246">246</a>         <strong class="jxr_keyword">double</strong> scale    = FastMath.scalb(1.0, 100);
+<a class="jxr_linenumber" name="247" href="#247">247</a>         Vector3D big1   = <strong class="jxr_keyword">new</strong> Vector3D(scale, v1);
+<a class="jxr_linenumber" name="248" href="#248">248</a>         Vector3D small2 = <strong class="jxr_keyword">new</strong> Vector3D(1 / scale, v2);
+<a class="jxr_linenumber" name="249" href="#249">249</a>         checkVector(Vector3D.crossProduct(big1, small2), -1, 2, 1);
+<a class="jxr_linenumber" name="250" href="#250">250</a> 
+<a class="jxr_linenumber" name="251" href="#251">251</a>     }
+<a class="jxr_linenumber" name="252" href="#252">252</a> 
+<a class="jxr_linenumber" name="253" href="#253">253</a>     @Test
+<a class="jxr_linenumber" name="254" href="#254">254</a>     <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testAngular() {
+<a class="jxr_linenumber" name="255" href="#255">255</a>         Assert.assertEquals(0,           Vector3D.PLUS_I.getAlpha(), 1.0e-10);
+<a class="jxr_linenumber" name="256" href="#256">256</a>         Assert.assertEquals(0,           Vector3D.PLUS_I.getDelta(), 1.0e-10);
+<a class="jxr_linenumber" name="257" href="#257">257</a>         Assert.assertEquals(FastMath.PI / 2, Vector3D.PLUS_J.getAlpha(), 1.0e-10);
+<a class="jxr_linenumber" name="258" href="#258">258</a>         Assert.assertEquals(0,           Vector3D.PLUS_J.getDelta(), 1.0e-10);
+<a class="jxr_linenumber" name="259" href="#259">259</a>         Assert.assertEquals(0,           Vector3D.PLUS_K.getAlpha(), 1.0e-10);
+<a class="jxr_linenumber" name="260" href="#260">260</a>         Assert.assertEquals(FastMath.PI / 2, Vector3D.PLUS_K.getDelta(), 1.0e-10);
+<a class="jxr_linenumber" name="261" href="#261">261</a>       
+<a class="jxr_linenumber" name="262" href="#262">262</a>         Vector3D u = <strong class="jxr_keyword">new</strong> Vector3D(-1, 1, -1);
+<a class="jxr_linenumber" name="263" href="#263">263</a>         Assert.assertEquals(3 * FastMath.PI /4, u.getAlpha(), 1.0e-10);
+<a class="jxr_linenumber" name="264" href="#264">264</a>         Assert.assertEquals(-1.0 / FastMath.sqrt(3), FastMath.sin(u.getDelta()), 1.0e-10);
 <a class="jxr_linenumber" name="265" href="#265">265</a>     }
 <a class="jxr_linenumber" name="266" href="#266">266</a> 
 <a class="jxr_linenumber" name="267" href="#267">267</a>     @Test
-<a class="jxr_linenumber" name="268" href="#268">268</a>     <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testDotProduct() {
-<a class="jxr_linenumber" name="269" href="#269">269</a>         <em class="jxr_comment">// we compare accurate versus naive dot product implementations</em>
-<a class="jxr_linenumber" name="270" href="#270">270</a>         <em class="jxr_comment">// on regular vectors (i.e. not extreme cases like in the previous test)</em>
-<a class="jxr_linenumber" name="271" href="#271">271</a>         Well1024a random = <strong class="jxr_keyword">new</strong> Well1024a(553267312521321234l);
-<a class="jxr_linenumber" name="272" href="#272">272</a>         <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">int</strong> i = 0; i &lt; 10000; ++i) {
-<a class="jxr_linenumber" name="273" href="#273">273</a>             <strong class="jxr_keyword">double</strong> ux = 10000 * random.nextDouble();
-<a class="jxr_linenumber" name="274" href="#274">274</a>             <strong class="jxr_keyword">double</strong> uy = 10000 * random.nextDouble();
-<a class="jxr_linenumber" name="275" href="#275">275</a>             <strong class="jxr_keyword">double</strong> uz = 10000 * random.nextDouble();
-<a class="jxr_linenumber" name="276" href="#276">276</a>             <strong class="jxr_keyword">double</strong> vx = 10000 * random.nextDouble();
-<a class="jxr_linenumber" name="277" href="#277">277</a>             <strong class="jxr_keyword">double</strong> vy = 10000 * random.nextDouble();
-<a class="jxr_linenumber" name="278" href="#278">278</a>             <strong class="jxr_keyword">double</strong> vz = 10000 * random.nextDouble();
-<a class="jxr_linenumber" name="279" href="#279">279</a>             <strong class="jxr_keyword">double</strong> sNaive = ux * vx + uy * vy + uz * vz;
-<a class="jxr_linenumber" name="280" href="#280">280</a>             <strong class="jxr_keyword">double</strong> sAccurate = <strong class="jxr_keyword">new</strong> Vector3D(ux, uy, uz).dotProduct(<strong class="jxr_keyword">new</strong> Vector3D(vx, vy, vz));
-<a class="jxr_linenumber" name="281" href="#281">281</a>             Assert.assertEquals(sNaive, sAccurate, 2.5e-16 * sAccurate);
-<a class="jxr_linenumber" name="282" href="#282">282</a>         }
-<a class="jxr_linenumber" name="283" href="#283">283</a>     }
-<a class="jxr_linenumber" name="284" href="#284">284</a> 
-<a class="jxr_linenumber" name="285" href="#285">285</a>     @Test
-<a class="jxr_linenumber" name="286" href="#286">286</a>     <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testAccurateCrossProduct() {
-<a class="jxr_linenumber" name="287" href="#287">287</a>         <em class="jxr_comment">// the vectors u1 and u2 are nearly but not exactly anti-parallel</em>
-<a class="jxr_linenumber" name="288" href="#288">288</a>         <em class="jxr_comment">// (7.31e-16 degrees from 180 degrees) naive cross product (i.e.</em>
-<a class="jxr_linenumber" name="289" href="#289">289</a>         <em class="jxr_comment">// computing u1.x * u2.x + u1.y * u2.y + u1.z * u2.z</em>
-<a class="jxr_linenumber" name="290" href="#290">290</a>         <em class="jxr_comment">// leads to a result of   [0.0009765, -0.0001220, -0.0039062],</em>
-<a class="jxr_linenumber" name="291" href="#291">291</a>         <em class="jxr_comment">// instead of the correct [0.0006913, -0.0001254, -0.0007909]</em>
-<a class="jxr_linenumber" name="292" href="#292">292</a>         <strong class="jxr_keyword">final</strong> Vector3D u1 = <strong class="jxr_keyword">new</strong> Vector3D(-1321008684645961.0 /   268435456.0,
-<a class="jxr_linenumber" name="293" href="#293">293</a>                                          -5774608829631843.0 /   268435456.0,
-<a class="jxr_linenumber" name="294" href="#294">294</a>                                          -7645843051051357.0 /  8589934592.0);
-<a class="jxr_linenumber" name="295" href="#295">295</a>         <strong class="jxr_keyword">final</strong> Vector3D u2 = <strong class="jxr_keyword">new</strong> Vector3D( 1796571811118507.0 /  2147483648.0,
-<a class="jxr_linenumber" name="296" href="#296">296</a>                                           7853468008299307.0 /  2147483648.0,
-<a class="jxr_linenumber" name="297" href="#297">297</a>                                           2599586637357461.0 / 17179869184.0);
-<a class="jxr_linenumber" name="298" href="#298">298</a>         <strong class="jxr_keyword">final</strong> Vector3D u3 = <strong class="jxr_keyword">new</strong> Vector3D(12753243807587107.0 / 18446744073709551616.0, 
-<a class="jxr_linenumber" name="299" href="#299">299</a>                                          -2313766922703915.0 / 18446744073709551616.0, 
-<a class="jxr_linenumber" name="300" href="#300">300</a>                                           -227970081415313.0 /   288230376151711744.0);
-<a class="jxr_linenumber" name="301" href="#301">301</a>         Vector3D cNaive = <strong class="jxr_keyword">new</strong> Vector3D(u1.getY() * u2.getZ() - u1.getZ() * u2.getY(),
-<a class="jxr_linenumber" name="302" href="#302">302</a>                                        u1.getZ() * u2.getX() - u1.getX() * u2.getZ(),
-<a class="jxr_linenumber" name="303" href="#303">303</a>                                        u1.getX() * u2.getY() - u1.getY() * u2.getX());
-<a class="jxr_linenumber" name="304" href="#304">304</a>         Vector3D cAccurate = u1.crossProduct(u2);
-<a class="jxr_linenumber" name="305" href="#305">305</a>         Assert.assertTrue(u3.distance(cNaive) &gt; 2.9 * u3.getNorm());
-<a class="jxr_linenumber" name="306" href="#306">306</a>         Assert.assertEquals(0.0, u3.distance(cAccurate), 1.0e-30 * cAccurate.getNorm());
-<a class="jxr_linenumber" name="307" href="#307">307</a>     }
-<a class="jxr_linenumber" name="308" href="#308">308</a> 
-<a class="jxr_linenumber" name="309" href="#309">309</a>     @Test
-<a class="jxr_linenumber" name="310" href="#310">310</a>     <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testCrossProduct() {
-<a class="jxr_linenumber" name="311" href="#311">311</a>         <em class="jxr_comment">// we compare accurate versus naive cross product implementations</em>
-<a class="jxr_linenumber" name="312" href="#312">312</a>         <em class="jxr_comment">// on regular vectors (i.e. not extreme cases like in the previous test)</em>
-<a class="jxr_linenumber" name="313" href="#313">313</a>         Well1024a random = <strong class="jxr_keyword">new</strong> Well1024a(885362227452043214l);
-<a class="jxr_linenumber" name="314" href="#314">314</a>         <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">int</strong> i = 0; i &lt; 10000; ++i) {
-<a class="jxr_linenumber" name="315" href="#315">315</a>             <strong class="jxr_keyword">double</strong> ux = 10000 * random.nextDouble();
-<a class="jxr_linenumber" name="316" href="#316">316</a>             <strong class="jxr_keyword">double</strong> uy = 10000 * random.nextDouble();
-<a class="jxr_linenumber" name="317" href="#317">317</a>             <strong class="jxr_keyword">double</strong> uz = 10000 * random.nextDouble();
-<a class="jxr_linenumber" name="318" href="#318">318</a>             <strong class="jxr_keyword">double</strong> vx = 10000 * random.nextDouble();
-<a class="jxr_linenumber" name="319" href="#319">319</a>             <strong class="jxr_keyword">double</strong> vy = 10000 * random.nextDouble();
-<a class="jxr_linenumber" name="320" href="#320">320</a>             <strong class="jxr_keyword">double</strong> vz = 10000 * random.nextDouble();
-<a class="jxr_linenumber" name="321" href="#321">321</a>             Vector3D cNaive = <strong class="jxr_keyword">new</strong> Vector3D(uy * vz - uz * vy, uz * vx - ux * vz, ux * vy - uy * vx);
-<a class="jxr_linenumber" name="322" href="#322">322</a>             Vector3D cAccurate = <strong class="jxr_keyword">new</strong> Vector3D(ux, uy, uz).crossProduct(<strong class="jxr_keyword">new</strong> Vector3D(vx, vy, vz));
-<a class="jxr_linenumber" name="323" href="#323">323</a>             Assert.assertEquals(0.0, cAccurate.distance(cNaive), 6.0e-15 * cAccurate.getNorm());
-<a class="jxr_linenumber" name="324" href="#324">324</a>         }
-<a class="jxr_linenumber" name="325" href="#325">325</a>     }
-<a class="jxr_linenumber" name="326" href="#326">326</a> 
-<a class="jxr_linenumber" name="327" href="#327">327</a>     <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">void</strong> checkVector(Vector3D v, <strong class="jxr_keyword">double</strong> x, <strong class="jxr_keyword">double</strong> y, <strong class="jxr_keyword">double</strong> z) {
-<a class="jxr_linenumber" name="328" href="#328">328</a>         Assert.assertEquals(x, v.getX(), 1.0e-12);
-<a class="jxr_linenumber" name="329" href="#329">329</a>         Assert.assertEquals(y, v.getY(), 1.0e-12);
-<a class="jxr_linenumber" name="330" href="#330">330</a>         Assert.assertEquals(z, v.getZ(), 1.0e-12);
-<a class="jxr_linenumber" name="331" href="#331">331</a>     }
-<a class="jxr_linenumber" name="332" href="#332">332</a> }
+<a class="jxr_linenumber" name="268" href="#268">268</a>     <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testAngularSeparation() <strong class="jxr_keyword">throws</strong> MathArithmeticException {
+<a class="jxr_linenumber" name="269" href="#269">269</a>         Vector3D v1 = <strong class="jxr_keyword">new</strong> Vector3D(2, -1, 4);
+<a class="jxr_linenumber" name="270" href="#270">270</a> 
+<a class="jxr_linenumber" name="271" href="#271">271</a>         Vector3D  k = v1.normalize();
+<a class="jxr_linenumber" name="272" href="#272">272</a>         Vector3D  i = k.orthogonal();
+<a class="jxr_linenumber" name="273" href="#273">273</a>         Vector3D v2 = k.scalarMultiply(FastMath.cos(1.2)).add(i.scalarMultiply(FastMath.sin(1.2)));
+<a class="jxr_linenumber" name="274" href="#274">274</a> 
+<a class="jxr_linenumber" name="275" href="#275">275</a>         Assert.assertTrue(FastMath.abs(Vector3D.angle(v1, v2) - 1.2) &lt; 1.0e-12);
+<a class="jxr_linenumber" name="276" href="#276">276</a>   }
+<a class="jxr_linenumber" name="277" href="#277">277</a> 
+<a class="jxr_linenumber" name="278" href="#278">278</a>     @Test
+<a class="jxr_linenumber" name="279" href="#279">279</a>     <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testNormalize() <strong class="jxr_keyword">throws</strong> MathArithmeticException {
+<a class="jxr_linenumber" name="280" href="#280">280</a>         Assert.assertEquals(1.0, <strong class="jxr_keyword">new</strong> Vector3D(5, -4, 2).normalize().getNorm(), 1.0e-12);
+<a class="jxr_linenumber" name="281" href="#281">281</a>         <strong class="jxr_keyword">try</strong> {
+<a class="jxr_linenumber" name="282" href="#282">282</a>             Vector3D.ZERO.normalize();
+<a class="jxr_linenumber" name="283" href="#283">283</a>             Assert.fail(<span class="jxr_string">"an exception should have been thrown"</span>);
+<a class="jxr_linenumber" name="284" href="#284">284</a>         } <strong class="jxr_keyword">catch</strong> (MathArithmeticException ae) {
+<a class="jxr_linenumber" name="285" href="#285">285</a>             <em class="jxr_comment">// expected behavior</em>
+<a class="jxr_linenumber" name="286" href="#286">286</a>         }
+<a class="jxr_linenumber" name="287" href="#287">287</a>     }
+<a class="jxr_linenumber" name="288" href="#288">288</a> 
+<a class="jxr_linenumber" name="289" href="#289">289</a>     @Test
+<a class="jxr_linenumber" name="290" href="#290">290</a>     <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testNegate() {
+<a class="jxr_linenumber" name="291" href="#291">291</a>         checkVector(<strong class="jxr_keyword">new</strong> Vector3D(0.1, 2.5, 1.3).negate(), -0.1, -2.5, -1.3);
+<a class="jxr_linenumber" name="292" href="#292">292</a>     }
+<a class="jxr_linenumber" name="293" href="#293">293</a> 
+<a class="jxr_linenumber" name="294" href="#294">294</a>     @Test
+<a class="jxr_linenumber" name="295" href="#295">295</a>     <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testOrthogonal() <strong class="jxr_keyword">throws</strong> MathArithmeticException {
+<a class="jxr_linenumber" name="296" href="#296">296</a>         Vector3D v1 = <strong class="jxr_keyword">new</strong> Vector3D(0.1, 2.5, 1.3);
+<a class="jxr_linenumber" name="297" href="#297">297</a>         Assert.assertEquals(0.0, Vector3D.dotProduct(v1, v1.orthogonal()), 1.0e-12);
+<a class="jxr_linenumber" name="298" href="#298">298</a>         Vector3D v2 = <strong class="jxr_keyword">new</strong> Vector3D(2.3, -0.003, 7.6);
+<a class="jxr_linenumber" name="299" href="#299">299</a>         Assert.assertEquals(0.0, Vector3D.dotProduct(v2, v2.orthogonal()), 1.0e-12);
+<a class="jxr_linenumber" name="300" href="#300">300</a>         Vector3D v3 = <strong class="jxr_keyword">new</strong> Vector3D(-1.7, 1.4, 0.2);
+<a class="jxr_linenumber" name="301" href="#301">301</a>         Assert.assertEquals(0.0, Vector3D.dotProduct(v3, v3.orthogonal()), 1.0e-12);
+<a class="jxr_linenumber" name="302" href="#302">302</a>         Vector3D v4 = <strong class="jxr_keyword">new</strong> Vector3D(4.2, 0.1, -1.8);
+<a class="jxr_linenumber" name="303" href="#303">303</a>         Assert.assertEquals(0.0, Vector3D.dotProduct(v4, v4.orthogonal()), 1.0e-12);
+<a class="jxr_linenumber" name="304" href="#304">304</a>         <strong class="jxr_keyword">try</strong> {
+<a class="jxr_linenumber" name="305" href="#305">305</a>             <strong class="jxr_keyword">new</strong> Vector3D(0, 0, 0).orthogonal();
+<a class="jxr_linenumber" name="306" href="#306">306</a>             Assert.fail(<span class="jxr_string">"an exception should have been thrown"</span>);
+<a class="jxr_linenumber" name="307" href="#307">307</a>         } <strong class="jxr_keyword">catch</strong> (MathArithmeticException ae) {
+<a class="jxr_linenumber" name="308" href="#308">308</a>             <em class="jxr_comment">// expected behavior</em>
+<a class="jxr_linenumber" name="309" href="#309">309</a>         }
+<a class="jxr_linenumber" name="310" href="#310">310</a>     }
+<a class="jxr_linenumber" name="311" href="#311">311</a>     @Test
+<a class="jxr_linenumber" name="312" href="#312">312</a>     <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testAngle() <strong class="jxr_keyword">throws</strong> MathArithmeticException {
+<a class="jxr_linenumber" name="313" href="#313">313</a>         Assert.assertEquals(0.22572612855273393616,
+<a class="jxr_linenumber" name="314" href="#314">314</a>                             Vector3D.angle(<strong class="jxr_keyword">new</strong> Vector3D(1, 2, 3), <strong class="jxr_keyword">new</strong> Vector3D(4, 5, 6)),
+<a class="jxr_linenumber" name="315" href="#315">315</a>                             1.0e-12);
+<a class="jxr_linenumber" name="316" href="#316">316</a>         Assert.assertEquals(7.98595620686106654517199e-8,
+<a class="jxr_linenumber" name="317" href="#317">317</a>                             Vector3D.angle(<strong class="jxr_keyword">new</strong> Vector3D(1, 2, 3), <strong class="jxr_keyword">new</strong> Vector3D(2, 4, 6.000001)),
+<a class="jxr_linenumber" name="318" href="#318">318</a>                             1.0e-12);
+<a class="jxr_linenumber" name="319" href="#319">319</a>         Assert.assertEquals(3.14159257373023116985197793156,
+<a class="jxr_linenumber" name="320" href="#320">320</a>                             Vector3D.angle(<strong class="jxr_keyword">new</strong> Vector3D(1, 2, 3), <strong class="jxr_keyword">new</strong> Vector3D(-2, -4, -6.000001)),
+<a class="jxr_linenumber" name="321" href="#321">321</a>                             1.0e-12);
+<a class="jxr_linenumber" name="322" href="#322">322</a>         <strong class="jxr_keyword">try</strong> {
+<a class="jxr_linenumber" name="323" href="#323">323</a>             Vector3D.angle(Vector3D.ZERO, Vector3D.PLUS_I);
+<a class="jxr_linenumber" name="324" href="#324">324</a>             Assert.fail(<span class="jxr_string">"an exception should have been thrown"</span>);
+<a class="jxr_linenumber" name="325" href="#325">325</a>         } <strong class="jxr_keyword">catch</strong> (MathArithmeticException ae) {
+<a class="jxr_linenumber" name="326" href="#326">326</a>             <em class="jxr_comment">// expected behavior</em>
+<a class="jxr_linenumber" name="327" href="#327">327</a>         }
+<a class="jxr_linenumber" name="328" href="#328">328</a>     }
+<a class="jxr_linenumber" name="329" href="#329">329</a> 
+<a class="jxr_linenumber" name="330" href="#330">330</a>     @Test
+<a class="jxr_linenumber" name="331" href="#331">331</a>     <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testAccurateDotProduct() {
+<a class="jxr_linenumber" name="332" href="#332">332</a>         <em class="jxr_comment">// the following two vectors are nearly but not exactly orthogonal</em>
+<a class="jxr_linenumber" name="333" href="#333">333</a>         <em class="jxr_comment">// naive dot product (i.e. computing u1.x * u2.x + u1.y * u2.y + u1.z * u2.z</em>
+<a class="jxr_linenumber" name="334" href="#334">334</a>         <em class="jxr_comment">// leads to a result of 0.0, instead of the correct -1.855129...</em>
+<a class="jxr_linenumber" name="335" href="#335">335</a>         Vector3D u1 = <strong class="jxr_keyword">new</strong> Vector3D(-1321008684645961.0 /  268435456.0,
+<a class="jxr_linenumber" name="336" href="#336">336</a>                                    -5774608829631843.0 /  268435456.0,
+<a class="jxr_linenumber" name="337" href="#337">337</a>                                    -7645843051051357.0 / 8589934592.0);
+<a class="jxr_linenumber" name="338" href="#338">338</a>         Vector3D u2 = <strong class="jxr_keyword">new</strong> Vector3D(-5712344449280879.0 /    2097152.0,
+<a class="jxr_linenumber" name="339" href="#339">339</a>                                    -4550117129121957.0 /    2097152.0,
+<a class="jxr_linenumber" name="340" href="#340">340</a>                                     8846951984510141.0 /     131072.0);
+<a class="jxr_linenumber" name="341" href="#341">341</a>         <strong class="jxr_keyword">double</strong> sNaive = u1.getX() * u2.getX() + u1.getY() * u2.getY() + u1.getZ() * u2.getZ();
+<a class="jxr_linenumber" name="342" href="#342">342</a>         <strong class="jxr_keyword">double</strong> sAccurate = u1.dotProduct(u2);
+<a class="jxr_linenumber" name="343" href="#343">343</a>         Assert.assertEquals(0.0, sNaive, 1.0e-30);
+<a class="jxr_linenumber" name="344" href="#344">344</a>         Assert.assertEquals(-2088690039198397.0 / 1125899906842624.0, sAccurate, 1.0e-16);
+<a class="jxr_linenumber" name="345" href="#345">345</a>     }
+<a class="jxr_linenumber" name="346" href="#346">346</a> 
+<a class="jxr_linenumber" name="347" href="#347">347</a>     @Test
+<a class="jxr_linenumber" name="348" href="#348">348</a>     <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testDotProduct() {
+<a class="jxr_linenumber" name="349" href="#349">349</a>         <em class="jxr_comment">// we compare accurate versus naive dot product implementations</em>
+<a class="jxr_linenumber" name="350" href="#350">350</a>         <em class="jxr_comment">// on regular vectors (i.e. not extreme cases like in the previous test)</em>
+<a class="jxr_linenumber" name="351" href="#351">351</a>         Well1024a random = <strong class="jxr_keyword">new</strong> Well1024a(553267312521321234l);
+<a class="jxr_linenumber" name="352" href="#352">352</a>         <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">int</strong> i = 0; i &lt; 10000; ++i) {
+<a class="jxr_linenumber" name="353" href="#353">353</a>             <strong class="jxr_keyword">double</strong> ux = 10000 * random.nextDouble();
+<a class="jxr_linenumber" name="354" href="#354">354</a>             <strong class="jxr_keyword">double</strong> uy = 10000 * random.nextDouble();
+<a class="jxr_linenumber" name="355" href="#355">355</a>             <strong class="jxr_keyword">double</strong> uz = 10000 * random.nextDouble();
+<a class="jxr_linenumber" name="356" href="#356">356</a>             <strong class="jxr_keyword">double</strong> vx = 10000 * random.nextDouble();
+<a class="jxr_linenumber" name="357" href="#357">357</a>             <strong class="jxr_keyword">double</strong> vy = 10000 * random.nextDouble();
+<a class="jxr_linenumber" name="358" href="#358">358</a>             <strong class="jxr_keyword">double</strong> vz = 10000 * random.nextDouble();
+<a class="jxr_linenumber" name="359" href="#359">359</a>             <strong class="jxr_keyword">double</strong> sNaive = ux * vx + uy * vy + uz * vz;
+<a class="jxr_linenumber" name="360" href="#360">360</a>             <strong class="jxr_keyword">double</strong> sAccurate = <strong class="jxr_keyword">new</strong> Vector3D(ux, uy, uz).dotProduct(<strong class="jxr_keyword">new</strong> Vector3D(vx, vy, vz));
+<a class="jxr_linenumber" name="361" href="#361">361</a>             Assert.assertEquals(sNaive, sAccurate, 2.5e-16 * sAccurate);
+<a class="jxr_linenumber" name="362" href="#362">362</a>         }
+<a class="jxr_linenumber" name="363" href="#363">363</a>     }
+<a class="jxr_linenumber" name="364" href="#364">364</a> 
+<a class="jxr_linenumber" name="365" href="#365">365</a>     @Test
+<a class="jxr_linenumber" name="366" href="#366">366</a>     <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testAccurateCrossProduct() {
+<a class="jxr_linenumber" name="367" href="#367">367</a>         <em class="jxr_comment">// the vectors u1 and u2 are nearly but not exactly anti-parallel</em>
+<a class="jxr_linenumber" name="368" href="#368">368</a>         <em class="jxr_comment">// (7.31e-16 degrees from 180 degrees) naive cross product (i.e.</em>
+<a class="jxr_linenumber" name="369" href="#369">369</a>         <em class="jxr_comment">// computing u1.x * u2.x + u1.y * u2.y + u1.z * u2.z</em>
+<a class="jxr_linenumber" name="370" href="#370">370</a>         <em class="jxr_comment">// leads to a result of   [0.0009765, -0.0001220, -0.0039062],</em>
+<a class="jxr_linenumber" name="371" href="#371">371</a>         <em class="jxr_comment">// instead of the correct [0.0006913, -0.0001254, -0.0007909]</em>
+<a class="jxr_linenumber" name="372" href="#372">372</a>         <strong class="jxr_keyword">final</strong> Vector3D u1 = <strong class="jxr_keyword">new</strong> Vector3D(-1321008684645961.0 /   268435456.0,
+<a class="jxr_linenumber" name="373" href="#373">373</a>                                          -5774608829631843.0 /   268435456.0,
+<a class="jxr_linenumber" name="374" href="#374">374</a>                                          -7645843051051357.0 /  8589934592.0);
+<a class="jxr_linenumber" name="375" href="#375">375</a>         <strong class="jxr_keyword">final</strong> Vector3D u2 = <strong class="jxr_keyword">new</strong> Vector3D( 1796571811118507.0 /  2147483648.0,
+<a class="jxr_linenumber" name="376" href="#376">376</a>                                           7853468008299307.0 /  2147483648.0,
+<a class="jxr_linenumber" name="377" href="#377">377</a>                                           2599586637357461.0 / 17179869184.0);
+<a class="jxr_linenumber" name="378" href="#378">378</a>         <strong class="jxr_keyword">final</strong> Vector3D u3 = <strong class="jxr_keyword">new</strong> Vector3D(12753243807587107.0 / 18446744073709551616.0, 
+<a class="jxr_linenumber" name="379" href="#379">379</a>                                          -2313766922703915.0 / 18446744073709551616.0, 
+<a class="jxr_linenumber" name="380" href="#380">380</a>                                           -227970081415313.0 /   288230376151711744.0);
+<a class="jxr_linenumber" name="381" href="#381">381</a>         Vector3D cNaive = <strong class="jxr_keyword">new</strong> Vector3D(u1.getY() * u2.getZ() - u1.getZ() * u2.getY(),
+<a class="jxr_linenumber" name="382" href="#382">382</a>                                        u1.getZ() * u2.getX() - u1.getX() * u2.getZ(),
+<a class="jxr_linenumber" name="383" href="#383">383</a>                                        u1.getX() * u2.getY() - u1.getY() * u2.getX());
+<a class="jxr_linenumber" name="384" href="#384">384</a>         Vector3D cAccurate = u1.crossProduct(u2);
+<a class="jxr_linenumber" name="385" href="#385">385</a>         Assert.assertTrue(u3.distance(cNaive) &gt; 2.9 * u3.getNorm());
+<a class="jxr_linenumber" name="386" href="#386">386</a>         Assert.assertEquals(0.0, u3.distance(cAccurate), 1.0e-30 * cAccurate.getNorm());
+<a class="jxr_linenumber" name="387" href="#387">387</a>     }
+<a class="jxr_linenumber" name="388" href="#388">388</a> 
+<a class="jxr_linenumber" name="389" href="#389">389</a>     @Test
+<a class="jxr_linenumber" name="390" href="#390">390</a>     <strong class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong> testCrossProduct() {
+<a class="jxr_linenumber" name="391" href="#391">391</a>         <em class="jxr_comment">// we compare accurate versus naive cross product implementations</em>
+<a class="jxr_linenumber" name="392" href="#392">392</a>         <em class="jxr_comment">// on regular vectors (i.e. not extreme cases like in the previous test)</em>
+<a class="jxr_linenumber" name="393" href="#393">393</a>         Well1024a random = <strong class="jxr_keyword">new</strong> Well1024a(885362227452043214l);
+<a class="jxr_linenumber" name="394" href="#394">394</a>         <strong class="jxr_keyword">for</strong> (<strong class="jxr_keyword">int</strong> i = 0; i &lt; 10000; ++i) {
+<a class="jxr_linenumber" name="395" href="#395">395</a>             <strong class="jxr_keyword">double</strong> ux = 10000 * random.nextDouble();
+<a class="jxr_linenumber" name="396" href="#396">396</a>             <strong class="jxr_keyword">double</strong> uy = 10000 * random.nextDouble();
+<a class="jxr_linenumber" name="397" href="#397">397</a>             <strong class="jxr_keyword">double</strong> uz = 10000 * random.nextDouble();
+<a class="jxr_linenumber" name="398" href="#398">398</a>             <strong class="jxr_keyword">double</strong> vx = 10000 * random.nextDouble();
+<a class="jxr_linenumber" name="399" href="#399">399</a>             <strong class="jxr_keyword">double</strong> vy = 10000 * random.nextDouble();
+<a class="jxr_linenumber" name="400" href="#400">400</a>             <strong class="jxr_keyword">double</strong> vz = 10000 * random.nextDouble();
+<a class="jxr_linenumber" name="401" href="#401">401</a>             Vector3D cNaive = <strong class="jxr_keyword">new</strong> Vector3D(uy * vz - uz * vy, uz * vx - ux * vz, ux * vy - uy * vx);
+<a class="jxr_linenumber" name="402" href="#402">402</a>             Vector3D cAccurate = <strong class="jxr_keyword">new</strong> Vector3D(ux, uy, uz).crossProduct(<strong class="jxr_keyword">new</strong> Vector3D(vx, vy, vz));
+<a class="jxr_linenumber" name="403" href="#403">403</a>             Assert.assertEquals(0.0, cAccurate.distance(cNaive), 6.0e-15 * cAccurate.getNorm());
+<a class="jxr_linenumber" name="404" href="#404">404</a>         }
+<a class="jxr_linenumber" name="405" href="#405">405</a>     }
+<a class="jxr_linenumber" name="406" href="#406">406</a> 
+<a class="jxr_linenumber" name="407" href="#407">407</a>     <strong class="jxr_keyword">private</strong> <strong class="jxr_keyword">void</strong> checkVector(Vector3D v, <strong class="jxr_keyword">double</strong> x, <strong class="jxr_keyword">double</strong> y, <strong class="jxr_keyword">double</strong> z) {
+<a class="jxr_linenumber" name="408" href="#408">408</a>         Assert.assertEquals(x, v.getX(), 1.0e-12);
+<a class="jxr_linenumber" name="409" href="#409">409</a>         Assert.assertEquals(y, v.getY(), 1.0e-12);
+<a class="jxr_linenumber" name="410" href="#410">410</a>         Assert.assertEquals(z, v.getZ(), 1.0e-12);
+<a class="jxr_linenumber" name="411" href="#411">411</a>     }
+<a class="jxr_linenumber" name="412" href="#412">412</a> }
 </pre>
 <hr/><div id="footer">This page was automatically generated by <a href="http://maven.apache.org/">Maven</a></div></body>
 </html>

Modified: websites/production/commons/content/proper/commons-math/xref-test/org/apache/commons/math3/geometry/euclidean/threed/package-frame.html
==============================================================================
--- websites/production/commons/content/proper/commons-math/xref-test/org/apache/commons/math3/geometry/euclidean/threed/package-frame.html (original)
+++ websites/production/commons/content/proper/commons-math/xref-test/org/apache/commons/math3/geometry/euclidean/threed/package-frame.html Sun Apr  7 07:52:05 2013
@@ -3,7 +3,7 @@
 <html xml:lang="en" lang="en">
 	<head>
 		<meta http-equiv="content-type" content="text/html; charset=UTF-8" />
-		<title>Commons Math 3.1.1 Reference Package org.apache.commons.math3.geometry.euclidean.threed</title>
+		<title>Commons Math 3.2 Reference Package org.apache.commons.math3.geometry.euclidean.threed</title>
 		<link rel="stylesheet" type="text/css" href="../../../../../../../stylesheet.css" title="style" />
 	</head>
 	<body>
@@ -16,6 +16,15 @@
 
       	<ul>
       		          	<li>
+            	<a href="FieldRotationDSTest.html" target="classFrame">FieldRotationDSTest</a>
+          	</li>
+          	          	<li>
+            	<a href="FieldRotationDfpTest.html" target="classFrame">FieldRotationDfpTest</a>
+          	</li>
+          	          	<li>
+            	<a href="FieldVector3DTest.html" target="classFrame">FieldVector3DTest</a>
+          	</li>
+          	          	<li>
             	<a href="FrenchVector3DFormatTest.html" target="classFrame">FrenchVector3DFormatTest</a>
           	</li>
           	          	<li>
@@ -34,6 +43,9 @@
             	<a href="RotationTest.html" target="classFrame">RotationTest</a>
           	</li>
           	          	<li>
+            	<a href="SphericalCoordinatesTest.html" target="classFrame">SphericalCoordinatesTest</a>
+          	</li>
+          	          	<li>
             	<a href="SubLineTest.html" target="classFrame">SubLineTest</a>
           	</li>
           	          	<li>