Propchange:
websites/production/commons/content/proper/commons-math/xref-test/org/apache/commons/math3/geometry/euclidean/threed/FieldVector3DTest.html
------------------------------------------------------------------------------
svn:eol-style = native
Propchange:
websites/production/commons/content/proper/commons-math/xref-test/org/apache/commons/math3/geometry/euclidean/threed/FieldVector3DTest.html
------------------------------------------------------------------------------
svn:keywords = Author Date Id Revision
Modified:
websites/production/commons/content/proper/commons-math/xref-test/org/apache/commons/math3/geometry/euclidean/threed/LineTest.html
==============================================================================
---
websites/production/commons/content/proper/commons-math/xref-test/org/apache/commons/math3/geometry/euclidean/threed/LineTest.html
(original)
+++
websites/production/commons/content/proper/commons-math/xref-test/org/apache/commons/math3/geometry/euclidean/threed/LineTest.html
Sun Apr 7 07:52:05 2013
@@ -139,7 +139,23 @@
<a class="jxr_linenumber" name="129" href="#129">129</a>
Assert.assertNull(l.intersection(<strong class="jxr_keyword">new</strong>
Line(<strong class="jxr_keyword">new</strong> Vector3D(0, -4, 0), <strong
class="jxr_keyword">new</strong> Vector3D(1, -4, 0))));
<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> }
+<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>
testRevert() {
+<a class="jxr_linenumber" name="134" href="#134">134</a>
+<a class="jxr_linenumber" name="135" href="#135">135</a> <em
class="jxr_comment">// setup</em>
+<a class="jxr_linenumber" name="136" href="#136">136</a> Line line =
<strong class="jxr_keyword">new</strong> Line(<strong
class="jxr_keyword">new</strong> Vector3D(1653345.6696423641,
6170370.041579291, 90000),
+<a class="jxr_linenumber" name="137" href="#137">137</a>
<strong class="jxr_keyword">new</strong> Vector3D(1650757.5050732433,
6160710.879908984, 0.9));
+<a class="jxr_linenumber" name="138" href="#138">138</a> Vector3D
expected = line.getDirection().negate();
+<a class="jxr_linenumber" name="139" href="#139">139</a>
+<a class="jxr_linenumber" name="140" href="#140">140</a> <em
class="jxr_comment">// action</em>
+<a class="jxr_linenumber" name="141" href="#141">141</a> Line reverted
= line.revert();
+<a class="jxr_linenumber" name="142" href="#142">142</a>
+<a class="jxr_linenumber" name="143" href="#143">143</a> <em
class="jxr_comment">// verify</em>
+<a class="jxr_linenumber" name="144" href="#144">144</a>
Assert.assertArrayEquals(expected.toArray(), reverted.getDirection().toArray(),
0);
+<a class="jxr_linenumber" name="145" href="#145">145</a>
+<a class="jxr_linenumber" name="146" href="#146">146</a> }
+<a class="jxr_linenumber" name="147" href="#147">147</a>
+<a class="jxr_linenumber" name="148" href="#148">148</a> }
</pre>
<hr/><div id="footer">This page was automatically generated by <a
href="http://maven.apache.org/";>Maven</a></div></body>
</html>
Added:
websites/production/commons/content/proper/commons-math/xref-test/org/apache/commons/math3/geometry/euclidean/threed/SphericalCoordinatesTest.html
==============================================================================
---
websites/production/commons/content/proper/commons-math/xref-test/org/apache/commons/math3/geometry/euclidean/threed/SphericalCoordinatesTest.html
(added)
+++
websites/production/commons/content/proper/commons-math/xref-test/org/apache/commons/math3/geometry/euclidean/threed/SphericalCoordinatesTest.html
Sun Apr 7 07:52:05 2013
@@ -0,0 +1,200 @@
+<!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>SphericalCoordinatesTest xref</title>
+<link type="text/css" rel="stylesheet"
href="../../../../../../../stylesheet.css" />
+</head>
+<body>
+<div id="overview"><a
href="../../../../../../../../testapidocs/org/apache/commons/math3/geometry/euclidean/threed/SphericalCoordinatesTest.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">
* Licensed to the Apache Software Foundation (ASF) under one or more</em>
+<a class="jxr_linenumber" name="3" href="#3">3</a> <em class="jxr_comment">
* contributor license agreements. See the NOTICE file distributed with</em>
+<a class="jxr_linenumber" name="4" href="#4">4</a> <em class="jxr_comment">
* this work for additional information regarding copyright ownership.</em>
+<a class="jxr_linenumber" name="5" href="#5">5</a> <em class="jxr_comment">
* The ASF licenses this file to You under the Apache License, Version 2.0</em>
+<a class="jxr_linenumber" name="6" href="#6">6</a> <em class="jxr_comment">
* (the "License"); you may not use this file except in compliance with</em>
+<a class="jxr_linenumber" name="7" href="#7">7</a> <em class="jxr_comment">
* the License. You may obtain a copy of the License at</em>
+<a class="jxr_linenumber" name="8" href="#8">8</a> <em class="jxr_comment">
*</em>
+<a class="jxr_linenumber" name="9" href="#9">9</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="10" href="#10">10</a> <em
class="jxr_comment"> *</em>
+<a class="jxr_linenumber" name="11" href="#11">11</a> <em
class="jxr_comment"> * Unless required by applicable law or agreed to in
writing, software</em>
+<a class="jxr_linenumber" name="12" href="#12">12</a> <em
class="jxr_comment"> * distributed under the License is distributed on an "AS
IS" BASIS,</em>
+<a class="jxr_linenumber" name="13" href="#13">13</a> <em
class="jxr_comment"> * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either
express or implied.</em>
+<a class="jxr_linenumber" name="14" href="#14">14</a> <em
class="jxr_comment"> * See the License for the specific language governing
permissions and</em>
+<a class="jxr_linenumber" name="15" href="#15">15</a> <em
class="jxr_comment"> * limitations under the License.</em>
+<a class="jxr_linenumber" name="16" href="#16">16</a> <em
class="jxr_comment"> */</em>
+<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.TestUtils;
+<a class="jxr_linenumber" name="21" href="#21">21</a> <strong
class="jxr_keyword">import</strong>
org.apache.commons.math3.analysis.differentiation.DerivativeStructure;
+<a class="jxr_linenumber" name="22" href="#22">22</a> <strong
class="jxr_keyword">import</strong>
org.apache.commons.math3.exception.DimensionMismatchException;
+<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/SphericalCoordinatesTest.html">SphericalCoordinatesTest</a>
{
+<a class="jxr_linenumber" name="28" href="#28">28</a>
+<a class="jxr_linenumber" name="29" href="#29">29</a> @Test
+<a class="jxr_linenumber" name="30" href="#30">30</a> <strong
class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong>
testCoordinatesStoC() <strong class="jxr_keyword">throws</strong>
DimensionMismatchException {
+<a class="jxr_linenumber" name="31" href="#31">31</a> <strong
class="jxr_keyword">double</strong> piO2 = 0.5 * FastMath.PI;
+<a class="jxr_linenumber" name="32" href="#32">32</a>
SphericalCoordinates sc1 = <strong class="jxr_keyword">new</strong>
SphericalCoordinates(2.0, 0, piO2);
+<a class="jxr_linenumber" name="33" href="#33">33</a>
Assert.assertEquals(0, sc1.getCartesian().distance(<strong
class="jxr_keyword">new</strong> Vector3D(2, 0, 0)), 1.0e-10);
+<a class="jxr_linenumber" name="34" href="#34">34</a>
SphericalCoordinates sc2 = <strong class="jxr_keyword">new</strong>
SphericalCoordinates(2.0, piO2, piO2);
+<a class="jxr_linenumber" name="35" href="#35">35</a>
Assert.assertEquals(0, sc2.getCartesian().distance(<strong
class="jxr_keyword">new</strong> Vector3D(0, 2, 0)), 1.0e-10);
+<a class="jxr_linenumber" name="36" href="#36">36</a>
SphericalCoordinates sc3 = <strong class="jxr_keyword">new</strong>
SphericalCoordinates(2.0, FastMath.PI, piO2);
+<a class="jxr_linenumber" name="37" href="#37">37</a>
Assert.assertEquals(0, sc3.getCartesian().distance(<strong
class="jxr_keyword">new</strong> Vector3D(-2, 0, 0)), 1.0e-10);
+<a class="jxr_linenumber" name="38" href="#38">38</a>
SphericalCoordinates sc4 = <strong class="jxr_keyword">new</strong>
SphericalCoordinates(2.0, -piO2, piO2);
+<a class="jxr_linenumber" name="39" href="#39">39</a>
Assert.assertEquals(0, sc4.getCartesian().distance(<strong
class="jxr_keyword">new</strong> Vector3D(0, -2, 0)), 1.0e-10);
+<a class="jxr_linenumber" name="40" href="#40">40</a>
SphericalCoordinates sc5 = <strong class="jxr_keyword">new</strong>
SphericalCoordinates(2.0, 1.23456, 0);
+<a class="jxr_linenumber" name="41" href="#41">41</a>
Assert.assertEquals(0, sc5.getCartesian().distance(<strong
class="jxr_keyword">new</strong> Vector3D(0, 0, 2)), 1.0e-10);
+<a class="jxr_linenumber" name="42" href="#42">42</a>
SphericalCoordinates sc6 = <strong class="jxr_keyword">new</strong>
SphericalCoordinates(2.0, 6.54321, FastMath.PI);
+<a class="jxr_linenumber" name="43" href="#43">43</a>
Assert.assertEquals(0, sc6.getCartesian().distance(<strong
class="jxr_keyword">new</strong> Vector3D(0, 0, -2)), 1.0e-10);
+<a class="jxr_linenumber" name="44" href="#44">44</a> }
+<a class="jxr_linenumber" name="45" href="#45">45</a>
+<a class="jxr_linenumber" name="46" href="#46">46</a> @Test
+<a class="jxr_linenumber" name="47" href="#47">47</a> <strong
class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong>
testCoordinatesCtoS() <strong class="jxr_keyword">throws</strong>
DimensionMismatchException {
+<a class="jxr_linenumber" name="48" href="#48">48</a> <strong
class="jxr_keyword">double</strong> piO2 = 0.5 * FastMath.PI;
+<a class="jxr_linenumber" name="49" href="#49">49</a>
SphericalCoordinates sc1 = <strong class="jxr_keyword">new</strong>
SphericalCoordinates(<strong class="jxr_keyword">new</strong> Vector3D(2, 0,
0));
+<a class="jxr_linenumber" name="50" href="#50">50</a>
Assert.assertEquals(2, sc1.getR(), 1.0e-10);
+<a class="jxr_linenumber" name="51" href="#51">51</a>
Assert.assertEquals(0, sc1.getTheta(), 1.0e-10);
+<a class="jxr_linenumber" name="52" href="#52">52</a>
Assert.assertEquals(piO2, sc1.getPhi(), 1.0e-10);
+<a class="jxr_linenumber" name="53" href="#53">53</a>
SphericalCoordinates sc2 = <strong class="jxr_keyword">new</strong>
SphericalCoordinates(<strong class="jxr_keyword">new</strong> Vector3D(0, 2,
0));
+<a class="jxr_linenumber" name="54" href="#54">54</a>
Assert.assertEquals(2, sc2.getR(), 1.0e-10);
+<a class="jxr_linenumber" name="55" href="#55">55</a>
Assert.assertEquals(piO2, sc2.getTheta(), 1.0e-10);
+<a class="jxr_linenumber" name="56" href="#56">56</a>
Assert.assertEquals(piO2, sc2.getPhi(), 1.0e-10);
+<a class="jxr_linenumber" name="57" href="#57">57</a>
SphericalCoordinates sc3 = <strong class="jxr_keyword">new</strong>
SphericalCoordinates(<strong class="jxr_keyword">new</strong> Vector3D(-2, 0,
0));
+<a class="jxr_linenumber" name="58" href="#58">58</a>
Assert.assertEquals(2, sc3.getR(), 1.0e-10);
+<a class="jxr_linenumber" name="59" href="#59">59</a>
Assert.assertEquals(FastMath.PI, sc3.getTheta(), 1.0e-10);
+<a class="jxr_linenumber" name="60" href="#60">60</a>
Assert.assertEquals(piO2, sc3.getPhi(), 1.0e-10);
+<a class="jxr_linenumber" name="61" href="#61">61</a>
SphericalCoordinates sc4 = <strong class="jxr_keyword">new</strong>
SphericalCoordinates(<strong class="jxr_keyword">new</strong> Vector3D(0, -2,
0));
+<a class="jxr_linenumber" name="62" href="#62">62</a>
Assert.assertEquals(2, sc4.getR(), 1.0e-10);
+<a class="jxr_linenumber" name="63" href="#63">63</a>
Assert.assertEquals(-piO2, sc4.getTheta(), 1.0e-10);
+<a class="jxr_linenumber" name="64" href="#64">64</a>
Assert.assertEquals(piO2, sc4.getPhi(), 1.0e-10);
+<a class="jxr_linenumber" name="65" href="#65">65</a>
SphericalCoordinates sc5 = <strong class="jxr_keyword">new</strong>
SphericalCoordinates(<strong class="jxr_keyword">new</strong> Vector3D(0, 0,
2));
+<a class="jxr_linenumber" name="66" href="#66">66</a>
Assert.assertEquals(2, sc5.getR(), 1.0e-10);
+<a class="jxr_linenumber" name="67" href="#67">67</a> <em
class="jxr_comment">// don't check theta on poles, as it is singular</em>
+<a class="jxr_linenumber" name="68" href="#68">68</a>
Assert.assertEquals(0, sc5.getPhi(), 1.0e-10);
+<a class="jxr_linenumber" name="69" href="#69">69</a>
SphericalCoordinates sc6 = <strong class="jxr_keyword">new</strong>
SphericalCoordinates(<strong class="jxr_keyword">new</strong> Vector3D(0, 0,
-2));
+<a class="jxr_linenumber" name="70" href="#70">70</a>
Assert.assertEquals(2, sc6.getR(), 1.0e-10);
+<a class="jxr_linenumber" name="71" href="#71">71</a> <em
class="jxr_comment">// don't check theta on poles, as it is singular</em>
+<a class="jxr_linenumber" name="72" href="#72">72</a>
Assert.assertEquals(FastMath.PI, sc6.getPhi(), 1.0e-10);
+<a class="jxr_linenumber" name="73" href="#73">73</a> }
+<a class="jxr_linenumber" name="74" href="#74">74</a>
+<a class="jxr_linenumber" name="75" href="#75">75</a> @Test
+<a class="jxr_linenumber" name="76" href="#76">76</a> <strong
class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong>
testGradient() {
+<a class="jxr_linenumber" name="77" href="#77">77</a> <strong
class="jxr_keyword">for</strong> (<strong class="jxr_keyword">double</strong> r
= 0.2; r < 10; r += 0.5) {
+<a class="jxr_linenumber" name="78" href="#78">78</a> <strong
class="jxr_keyword">for</strong> (<strong class="jxr_keyword">double</strong>
theta = 0; theta < 2 * FastMath.PI; theta += 0.1) {
+<a class="jxr_linenumber" name="79" href="#79">79</a> <strong
class="jxr_keyword">for</strong> (<strong class="jxr_keyword">double</strong>
phi = 0.1; phi < FastMath.PI; phi += 0.1) {
+<a class="jxr_linenumber" name="80" href="#80">80</a>
SphericalCoordinates sc = <strong class="jxr_keyword">new</strong>
SphericalCoordinates(r, theta, phi);
+<a class="jxr_linenumber" name="81" href="#81">81</a>
+<a class="jxr_linenumber" name="82" href="#82">82</a>
DerivativeStructure svalue = valueSpherical(<strong
class="jxr_keyword">new</strong> DerivativeStructure(3, 1, 0, r),
+<a class="jxr_linenumber" name="83" href="#83">83</a>
<strong
class="jxr_keyword">new</strong> DerivativeStructure(3, 1, 1, theta),
+<a class="jxr_linenumber" name="84" href="#84">84</a>
<strong
class="jxr_keyword">new</strong> DerivativeStructure(3, 1, 2, phi));
+<a class="jxr_linenumber" name="85" href="#85">85</a>
<strong class="jxr_keyword">double</strong>[] sGradient = <strong
class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[] {
+<a class="jxr_linenumber" name="86" href="#86">86</a>
svalue.getPartialDerivative(1, 0, 0),
+<a class="jxr_linenumber" name="87" href="#87">87</a>
svalue.getPartialDerivative(0, 1, 0),
+<a class="jxr_linenumber" name="88" href="#88">88</a>
svalue.getPartialDerivative(0, 0, 1),
+<a class="jxr_linenumber" name="89" href="#89">89</a> };
+<a class="jxr_linenumber" name="90" href="#90">90</a>
+<a class="jxr_linenumber" name="91" href="#91">91</a>
DerivativeStructure cvalue = valueCartesian(<strong
class="jxr_keyword">new</strong> DerivativeStructure(3, 1, 0,
sc.getCartesian().getX()),
+<a class="jxr_linenumber" name="92" href="#92">92</a>
<strong
class="jxr_keyword">new</strong> DerivativeStructure(3, 1, 1,
sc.getCartesian().getY()),
+<a class="jxr_linenumber" name="93" href="#93">93</a>
<strong
class="jxr_keyword">new</strong> DerivativeStructure(3, 1, 2,
sc.getCartesian().getZ()));
+<a class="jxr_linenumber" name="94" href="#94">94</a>
Vector3D refCGradient = <strong class="jxr_keyword">new</strong>
Vector3D(cvalue.getPartialDerivative(1, 0, 0),
+<a class="jxr_linenumber" name="95" href="#95">95</a>
cvalue.getPartialDerivative(0, 1, 0),
+<a class="jxr_linenumber" name="96" href="#96">96</a>
cvalue.getPartialDerivative(0, 0, 1));
+<a class="jxr_linenumber" name="97" href="#97">97</a>
+<a class="jxr_linenumber" name="98" href="#98">98</a>
Vector3D testCGradient = <strong class="jxr_keyword">new</strong>
Vector3D(sc.toCartesianGradient(sGradient));
+<a class="jxr_linenumber" name="99" href="#99">99</a>
+<a class="jxr_linenumber" name="100" href="#100">100</a>
Assert.assertEquals(0, testCGradient.distance(refCGradient) /
refCGradient.getNorm(), 5.0e-14);
+<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> }
+<a class="jxr_linenumber" name="104" href="#104">104</a> }
+<a class="jxr_linenumber" name="105" href="#105">105</a> }
+<a class="jxr_linenumber" name="106" href="#106">106</a>
+<a class="jxr_linenumber" name="107" href="#107">107</a> @Test
+<a class="jxr_linenumber" name="108" href="#108">108</a> <strong
class="jxr_keyword">public</strong> <strong class="jxr_keyword">void</strong>
testHessian() {
+<a class="jxr_linenumber" name="109" href="#109">109</a> <strong
class="jxr_keyword">for</strong> (<strong class="jxr_keyword">double</strong> r
= 0.2; r < 10; r += 0.5) {
+<a class="jxr_linenumber" name="110" href="#110">110</a> <strong
class="jxr_keyword">for</strong> (<strong class="jxr_keyword">double</strong>
theta = 0; theta < 2 * FastMath.PI; theta += 0.2) {
+<a class="jxr_linenumber" name="111" href="#111">111</a>
<strong class="jxr_keyword">for</strong> (<strong
class="jxr_keyword">double</strong> phi = 0.1; phi < FastMath.PI; phi +=
0.2) {
+<a class="jxr_linenumber" name="112" href="#112">112</a>
SphericalCoordinates sc = <strong class="jxr_keyword">new</strong>
SphericalCoordinates(r, theta, phi);
+<a class="jxr_linenumber" name="113" href="#113">113</a>
+<a class="jxr_linenumber" name="114" href="#114">114</a>
DerivativeStructure svalue = valueSpherical(<strong
class="jxr_keyword">new</strong> DerivativeStructure(3, 2, 0, r),
+<a class="jxr_linenumber" name="115" href="#115">115</a>
<strong
class="jxr_keyword">new</strong> DerivativeStructure(3, 2, 1, theta),
+<a class="jxr_linenumber" name="116" href="#116">116</a>
<strong
class="jxr_keyword">new</strong> DerivativeStructure(3, 2, 2, phi));
+<a class="jxr_linenumber" name="117" href="#117">117</a>
<strong class="jxr_keyword">double</strong>[] sGradient = <strong
class="jxr_keyword">new</strong> <strong class="jxr_keyword">double</strong>[] {
+<a class="jxr_linenumber" name="118" href="#118">118</a>
svalue.getPartialDerivative(1, 0, 0),
+<a class="jxr_linenumber" name="119" href="#119">119</a>
svalue.getPartialDerivative(0, 1, 0),
+<a class="jxr_linenumber" name="120" href="#120">120</a>
svalue.getPartialDerivative(0, 0, 1),
+<a class="jxr_linenumber" name="121" href="#121">121</a> };
+<a class="jxr_linenumber" name="122" href="#122">122</a>
<strong class="jxr_keyword">double</strong>[][] sHessian = <strong
class="jxr_keyword">new</strong> <strong
class="jxr_keyword">double</strong>[3][3];
+<a class="jxr_linenumber" name="123" href="#123">123</a>
sHessian[0][0] = svalue.getPartialDerivative(2, 0, 0); <em
class="jxr_comment">// d2F/dR2</em>
+<a class="jxr_linenumber" name="124" href="#124">124</a>
sHessian[1][0] = svalue.getPartialDerivative(1, 1, 0); <em
class="jxr_comment">// d2F/dRdTheta</em>
+<a class="jxr_linenumber" name="125" href="#125">125</a>
sHessian[2][0] = svalue.getPartialDerivative(1, 0, 1); <em
class="jxr_comment">// d2F/dRdPhi</em>
+<a class="jxr_linenumber" name="126" href="#126">126</a>
sHessian[0][1] = Double.NaN; <em class="jxr_comment">// just to check
upper-right part is not used</em>
+<a class="jxr_linenumber" name="127" href="#127">127</a>
sHessian[1][1] = svalue.getPartialDerivative(0, 2, 0); <em
class="jxr_comment">// d2F/dTheta2</em>
+<a class="jxr_linenumber" name="128" href="#128">128</a>
sHessian[2][1] = svalue.getPartialDerivative(0, 1, 1); <em
class="jxr_comment">// d2F/dThetadPhi</em>
+<a class="jxr_linenumber" name="129" href="#129">129</a>
sHessian[0][2] = Double.NaN; <em class="jxr_comment">// just to check
upper-right part is not used</em>
+<a class="jxr_linenumber" name="130" href="#130">130</a>
sHessian[1][2] = Double.NaN; <em class="jxr_comment">// just to check
upper-right part is not used</em>
+<a class="jxr_linenumber" name="131" href="#131">131</a>
sHessian[2][2] = svalue.getPartialDerivative(0, 0, 2); <em
class="jxr_comment">// d2F/dPhi2</em>
+<a class="jxr_linenumber" name="132" href="#132">132</a>
+<a class="jxr_linenumber" name="133" href="#133">133</a>
DerivativeStructure cvalue = valueCartesian(<strong
class="jxr_keyword">new</strong> DerivativeStructure(3, 2, 0,
sc.getCartesian().getX()),
+<a class="jxr_linenumber" name="134" href="#134">134</a>
<strong
class="jxr_keyword">new</strong> DerivativeStructure(3, 2, 1,
sc.getCartesian().getY()),
+<a class="jxr_linenumber" name="135" href="#135">135</a>
<strong
class="jxr_keyword">new</strong> DerivativeStructure(3, 2, 2,
sc.getCartesian().getZ()));
+<a class="jxr_linenumber" name="136" href="#136">136</a>
<strong class="jxr_keyword">double</strong>[][] refCHessian = <strong
class="jxr_keyword">new</strong> <strong
class="jxr_keyword">double</strong>[3][3];
+<a class="jxr_linenumber" name="137" href="#137">137</a>
refCHessian[0][0] = cvalue.getPartialDerivative(2, 0, 0); <em
class="jxr_comment">// d2F/dX2</em>
+<a class="jxr_linenumber" name="138" href="#138">138</a>
refCHessian[1][0] = cvalue.getPartialDerivative(1, 1, 0); <em
class="jxr_comment">// d2F/dXdY</em>
+<a class="jxr_linenumber" name="139" href="#139">139</a>
refCHessian[2][0] = cvalue.getPartialDerivative(1, 0, 1); <em
class="jxr_comment">// d2F/dXdZ</em>
+<a class="jxr_linenumber" name="140" href="#140">140</a>
refCHessian[0][1] = refCHessian[1][0];
+<a class="jxr_linenumber" name="141" href="#141">141</a>
refCHessian[1][1] = cvalue.getPartialDerivative(0, 2, 0); <em
class="jxr_comment">// d2F/dY2</em>
+<a class="jxr_linenumber" name="142" href="#142">142</a>
refCHessian[2][1] = cvalue.getPartialDerivative(0, 1, 1); <em
class="jxr_comment">// d2F/dYdZ</em>
+<a class="jxr_linenumber" name="143" href="#143">143</a>
refCHessian[0][2] = refCHessian[2][0];
+<a class="jxr_linenumber" name="144" href="#144">144</a>
refCHessian[1][2] = refCHessian[2][1];
+<a class="jxr_linenumber" name="145" href="#145">145</a>
refCHessian[2][2] = cvalue.getPartialDerivative(0, 0, 2); <em
class="jxr_comment">// d2F/dZ2</em>
+<a class="jxr_linenumber" name="146" href="#146">146</a>
<strong class="jxr_keyword">double</strong> norm = 0;
+<a class="jxr_linenumber" name="147" href="#147">147</a>
<strong class="jxr_keyword">for</strong> (<strong
class="jxr_keyword">int</strong> i = 0; i < 3; ++i) {
+<a class="jxr_linenumber" name="148" href="#148">148</a>
<strong class="jxr_keyword">for</strong> (<strong
class="jxr_keyword">int</strong> j = 0; j < 3; ++j) {
+<a class="jxr_linenumber" name="149" href="#149">149</a>
norm = FastMath.max(norm, FastMath.abs(refCHessian[i][j]));
+<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>
+<a class="jxr_linenumber" name="153" href="#153">153</a>
<strong class="jxr_keyword">double</strong>[][] testCHessian =
sc.toCartesianHessian(sHessian, sGradient);
+<a class="jxr_linenumber" name="154" href="#154">154</a>
<strong class="jxr_keyword">for</strong> (<strong
class="jxr_keyword">int</strong> i = 0; i < 3; ++i) {
+<a class="jxr_linenumber" name="155" href="#155">155</a>
<strong class="jxr_keyword">for</strong> (<strong
class="jxr_keyword">int</strong> j = 0; j < 3; ++j) {
+<a class="jxr_linenumber" name="156" href="#156">156</a>
Assert.assertEquals(<span class="jxr_string">""</span> +
FastMath.abs((refCHessian[i][j] - testCHessian[i][j]) / norm),
+<a class="jxr_linenumber" name="157" href="#157">157</a>
refCHessian[i][j], testCHessian[i][j], 1.0e-14 *
norm);
+<a class="jxr_linenumber" name="158" href="#158">158</a>
}
+<a class="jxr_linenumber" name="159" href="#159">159</a> }
+<a class="jxr_linenumber" name="160" href="#160">160</a>
+<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> }
+<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> <strong
class="jxr_keyword">public</strong> DerivativeStructure
valueCartesian(DerivativeStructure x, DerivativeStructure y,
DerivativeStructure z) {
+<a class="jxr_linenumber" name="167" href="#167">167</a> <strong
class="jxr_keyword">return</strong>
x.divide(y.multiply(5).add(10)).multiply(z.pow(3));
+<a class="jxr_linenumber" name="168" href="#168">168</a> }
+<a class="jxr_linenumber" name="169" href="#169">169</a>
+<a class="jxr_linenumber" name="170" href="#170">170</a> <strong
class="jxr_keyword">public</strong> DerivativeStructure
valueSpherical(DerivativeStructure r, DerivativeStructure theta,
DerivativeStructure phi) {
+<a class="jxr_linenumber" name="171" href="#171">171</a> <strong
class="jxr_keyword">return</strong>
valueCartesian(r.multiply(theta.cos()).multiply(phi.sin()),
+<a class="jxr_linenumber" name="172" href="#172">172</a>
r.multiply(theta.sin()).multiply(phi.sin()),
+<a class="jxr_linenumber" name="173" href="#173">173</a>
r.multiply(phi.cos()));
+<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>
testSerialization() {
+<a class="jxr_linenumber" name="178" href="#178">178</a>
SphericalCoordinates a = <strong class="jxr_keyword">new</strong>
SphericalCoordinates(3, 2, 1);
+<a class="jxr_linenumber" name="179" href="#179">179</a>
SphericalCoordinates b = (SphericalCoordinates)
TestUtils.serializeAndRecover(a);
+<a class="jxr_linenumber" name="180" href="#180">180</a>
Assert.assertEquals(0, a.getCartesian().distance(b.getCartesian()), 1.0e-10);
+<a class="jxr_linenumber" name="181" href="#181">181</a>
Assert.assertEquals(a.getR(), b.getR(), 1.0e-10);
+<a class="jxr_linenumber" name="182" href="#182">182</a>
Assert.assertEquals(a.getTheta(), b.getTheta(), 1.0e-10);
+<a class="jxr_linenumber" name="183" href="#183">183</a>
Assert.assertEquals(a.getPhi(), b.getPhi(), 1.0e-10);
+<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> }
+</pre>
+<hr/><div id="footer">This page was automatically generated by <a
href="http://maven.apache.org/";>Maven</a></div></body>
+</html>
+
Propchange:
websites/production/commons/content/proper/commons-math/xref-test/org/apache/commons/math3/geometry/euclidean/threed/SphericalCoordinatesTest.html
------------------------------------------------------------------------------
svn:eol-style = native
Propchange:
websites/production/commons/content/proper/commons-math/xref-test/org/apache/commons/math3/geometry/euclidean/threed/SphericalCoordinatesTest.html
------------------------------------------------------------------------------
svn:keywords = Author Date Id Revision
