svn commit: r934426 [5/10] - in /websites/production/commons/content/proper/commons-math/apidocs: org/apache/commons/math3/analysis/interpolation/ org/apache/commons/math3/analysis/interpolation/class-use/ org/apache/commons/math3/fitting/ org/apache/c...

Tue, 30 Dec 2014 09:16:01 -0800 

Added: 
websites/production/commons/content/proper/commons-math/apidocs/org/apache/commons/math3/special/class-use/BesselJ.BesselJResult.html
==============================================================================
--- 
websites/production/commons/content/proper/commons-math/apidocs/org/apache/commons/math3/special/class-use/BesselJ.BesselJResult.html
 (added)
+++ 
websites/production/commons/content/proper/commons-math/apidocs/org/apache/commons/math3/special/class-use/BesselJ.BesselJResult.html
 Tue Dec 30 17:15:36 2014
@@ -0,0 +1,164 @@
+<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN" 
"http://www.w3.org/TR/html4/loose.dtd";>
+<!-- NewPage -->
+<html lang="en">
+<head>
+<meta http-equiv="Content-Type" content="text/html" charset="UTF-8">
+<title>Uses of Class org.apache.commons.math3.special.BesselJ.BesselJResult 
(Apache Commons Math 3.4-SNAPSHOT API)</title>
+<link rel="stylesheet" type="text/css" href="../../../../../../stylesheet.css" 
title="Style">
+</head>
+<body>
+<script type="text/javascript"><!--
+    if (location.href.indexOf('is-external=true') == -1) {
+        parent.document.title="Uses of Class 
org.apache.commons.math3.special.BesselJ.BesselJResult (Apache Commons Math 
3.4-SNAPSHOT API)";
+    }
+//-->
+</script>
+<noscript>
+<div>JavaScript is disabled on your browser.</div>
+</noscript>
+<!-- ========= START OF TOP NAVBAR ======= -->
+<div class="topNav"><a name="navbar_top">
+<!--   -->
+</a><a href="#skip-navbar_top" title="Skip navigation links"></a><a 
name="navbar_top_firstrow">
+<!--   -->
+</a>
+<ul class="navList" title="Navigation">
+<li><a href="../../../../../../overview-summary.html">Overview</a></li>
+<li><a href="../package-summary.html">Package</a></li>
+<li><a 
href="../../../../../../org/apache/commons/math3/special/BesselJ.BesselJResult.html"
 title="class in org.apache.commons.math3.special">Class</a></li>
+<li class="navBarCell1Rev">Use</li>
+<li><a href="../package-tree.html">Tree</a></li>
+<li><a href="../../../../../../deprecated-list.html">Deprecated</a></li>
+<li><a href="../../../../../../index-all.html">Index</a></li>
+<li><a href="../../../../../../help-doc.html">Help</a></li>
+</ul>
+<div class="aboutLanguage"><em><script type="text/javascript" 
src="http://cdn.mathjax.org/mathjax/latest/MathJax.js?config=TeX-AMS-MML_HTMLorMML";></script></em></div>
+</div>
+<div class="subNav">
+<ul class="navList">
+<li>Prev</li>
+<li>Next</li>
+</ul>
+<ul class="navList">
+<li><a 
href="../../../../../../index.html?org/apache/commons/math3/special/class-use/BesselJ.BesselJResult.html"
 target="_top">Frames</a></li>
+<li><a href="BesselJ.BesselJResult.html" target="_top">No Frames</a></li>
+</ul>
+<ul class="navList" id="allclasses_navbar_top">
+<li><a href="../../../../../../allclasses-noframe.html">All Classes</a></li>
+</ul>
+<div>
+<script type="text/javascript"><!--
+  allClassesLink = document.getElementById("allclasses_navbar_top");
+  if(window==top) {
+    allClassesLink.style.display = "block";
+  }
+  else {
+    allClassesLink.style.display = "none";
+  }
+  //-->
+</script>
+</div>
+<a name="skip-navbar_top">
+<!--   -->
+</a></div>
+<!-- ========= END OF TOP NAVBAR ========= -->
+<div class="header">
+<h2 title="Uses of Class 
org.apache.commons.math3.special.BesselJ.BesselJResult" class="title">Uses of 
Class<br>org.apache.commons.math3.special.BesselJ.BesselJResult</h2>
+</div>
+<div class="classUseContainer">
+<ul class="blockList">
+<li class="blockList">
+<table border="0" cellpadding="3" cellspacing="0" summary="Use table, listing 
packages, and an explanation">
+<caption><span>Packages that use <a 
href="../../../../../../org/apache/commons/math3/special/BesselJ.BesselJResult.html"
 title="class in 
org.apache.commons.math3.special">BesselJ.BesselJResult</a></span><span 
class="tabEnd">&nbsp;</span></caption>
+<tr>
+<th class="colFirst" scope="col">Package</th>
+<th class="colLast" scope="col">Description</th>
+</tr>
+<tbody>
+<tr class="altColor">
+<td class="colFirst"><a 
href="#org.apache.commons.math3.special">org.apache.commons.math3.special</a></td>
+<td class="colLast">
+<div class="block">Implementations of special functions such as Beta and 
Gamma.</div>
+</td>
+</tr>
+</tbody>
+</table>
+</li>
+<li class="blockList">
+<ul class="blockList">
+<li class="blockList"><a name="org.apache.commons.math3.special">
+<!--   -->
+</a>
+<h3>Uses of <a 
href="../../../../../../org/apache/commons/math3/special/BesselJ.BesselJResult.html"
 title="class in org.apache.commons.math3.special">BesselJ.BesselJResult</a> in 
<a 
href="../../../../../../org/apache/commons/math3/special/package-summary.html">org.apache.commons.math3.special</a></h3>
+<table border="0" cellpadding="3" cellspacing="0" summary="Use table, listing 
methods, and an explanation">
+<caption><span>Methods in <a 
href="../../../../../../org/apache/commons/math3/special/package-summary.html">org.apache.commons.math3.special</a>
 that return <a 
href="../../../../../../org/apache/commons/math3/special/BesselJ.BesselJResult.html"
 title="class in 
org.apache.commons.math3.special">BesselJ.BesselJResult</a></span><span 
class="tabEnd">&nbsp;</span></caption>
+<tr>
+<th class="colFirst" scope="col">Modifier and Type</th>
+<th class="colLast" scope="col">Method and Description</th>
+</tr>
+<tbody>
+<tr class="altColor">
+<td class="colFirst"><code>static <a 
href="../../../../../../org/apache/commons/math3/special/BesselJ.BesselJResult.html"
 title="class in 
org.apache.commons.math3.special">BesselJ.BesselJResult</a></code></td>
+<td class="colLast"><span class="strong">BesselJ.</span><code><strong><a 
href="../../../../../../org/apache/commons/math3/special/BesselJ.html#rjBesl(double,%20double,%20int)">rjBesl</a></strong>(double&nbsp;x,
+            double&nbsp;alpha,
+            int&nbsp;nb)</code>
+<div class="block">Calculates Bessel functions \(J_{n+alpha}(x)\) for
+ non-negative argument x, and non-negative order n + alpha.</div>
+</td>
+</tr>
+</tbody>
+</table>
+</li>
+</ul>
+</li>
+</ul>
+</div>
+<!-- ======= START OF BOTTOM NAVBAR ====== -->
+<div class="bottomNav"><a name="navbar_bottom">
+<!--   -->
+</a><a href="#skip-navbar_bottom" title="Skip navigation links"></a><a 
name="navbar_bottom_firstrow">
+<!--   -->
+</a>
+<ul class="navList" title="Navigation">
+<li><a href="../../../../../../overview-summary.html">Overview</a></li>
+<li><a href="../package-summary.html">Package</a></li>
+<li><a 
href="../../../../../../org/apache/commons/math3/special/BesselJ.BesselJResult.html"
 title="class in org.apache.commons.math3.special">Class</a></li>
+<li class="navBarCell1Rev">Use</li>
+<li><a href="../package-tree.html">Tree</a></li>
+<li><a href="../../../../../../deprecated-list.html">Deprecated</a></li>
+<li><a href="../../../../../../index-all.html">Index</a></li>
+<li><a href="../../../../../../help-doc.html">Help</a></li>
+</ul>
+<div class="aboutLanguage"><em><script type="text/javascript" 
src="http://cdn.mathjax.org/mathjax/latest/MathJax.js?config=TeX-AMS-MML_HTMLorMML";></script></em></div>
+</div>
+<div class="subNav">
+<ul class="navList">
+<li>Prev</li>
+<li>Next</li>
+</ul>
+<ul class="navList">
+<li><a 
href="../../../../../../index.html?org/apache/commons/math3/special/class-use/BesselJ.BesselJResult.html"
 target="_top">Frames</a></li>
+<li><a href="BesselJ.BesselJResult.html" target="_top">No Frames</a></li>
+</ul>
+<ul class="navList" id="allclasses_navbar_bottom">
+<li><a href="../../../../../../allclasses-noframe.html">All Classes</a></li>
+</ul>
+<div>
+<script type="text/javascript"><!--
+  allClassesLink = document.getElementById("allclasses_navbar_bottom");
+  if(window==top) {
+    allClassesLink.style.display = "block";
+  }
+  else {
+    allClassesLink.style.display = "none";
+  }
+  //-->
+</script>
+</div>
+<a name="skip-navbar_bottom">
+<!--   -->
+</a></div>
+<!-- ======== END OF BOTTOM NAVBAR ======= -->
+<p class="legalCopy"><small>Copyright &#169; 2003&#x2013;2014 <a 
href="http://www.apache.org/";>The Apache Software Foundation</a>. All rights 
reserved.</small></p>
+</body>
+</html>

Propchange: 
websites/production/commons/content/proper/commons-math/apidocs/org/apache/commons/math3/special/class-use/BesselJ.BesselJResult.html
------------------------------------------------------------------------------
    svn:eol-style = native

Added: 
websites/production/commons/content/proper/commons-math/apidocs/org/apache/commons/math3/special/class-use/BesselJ.html
==============================================================================
--- 
websites/production/commons/content/proper/commons-math/apidocs/org/apache/commons/math3/special/class-use/BesselJ.html
 (added)
+++ 
websites/production/commons/content/proper/commons-math/apidocs/org/apache/commons/math3/special/class-use/BesselJ.html
 Tue Dec 30 17:15:36 2014
@@ -0,0 +1,117 @@
+<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN" 
"http://www.w3.org/TR/html4/loose.dtd";>
+<!-- NewPage -->
+<html lang="en">
+<head>
+<meta http-equiv="Content-Type" content="text/html" charset="UTF-8">
+<title>Uses of Class org.apache.commons.math3.special.BesselJ (Apache Commons 
Math 3.4-SNAPSHOT API)</title>
+<link rel="stylesheet" type="text/css" href="../../../../../../stylesheet.css" 
title="Style">
+</head>
+<body>
+<script type="text/javascript"><!--
+    if (location.href.indexOf('is-external=true') == -1) {
+        parent.document.title="Uses of Class 
org.apache.commons.math3.special.BesselJ (Apache Commons Math 3.4-SNAPSHOT 
API)";
+    }
+//-->
+</script>
+<noscript>
+<div>JavaScript is disabled on your browser.</div>
+</noscript>
+<!-- ========= START OF TOP NAVBAR ======= -->
+<div class="topNav"><a name="navbar_top">
+<!--   -->
+</a><a href="#skip-navbar_top" title="Skip navigation links"></a><a 
name="navbar_top_firstrow">
+<!--   -->
+</a>
+<ul class="navList" title="Navigation">
+<li><a href="../../../../../../overview-summary.html">Overview</a></li>
+<li><a href="../package-summary.html">Package</a></li>
+<li><a href="../../../../../../org/apache/commons/math3/special/BesselJ.html" 
title="class in org.apache.commons.math3.special">Class</a></li>
+<li class="navBarCell1Rev">Use</li>
+<li><a href="../package-tree.html">Tree</a></li>
+<li><a href="../../../../../../deprecated-list.html">Deprecated</a></li>
+<li><a href="../../../../../../index-all.html">Index</a></li>
+<li><a href="../../../../../../help-doc.html">Help</a></li>
+</ul>
+<div class="aboutLanguage"><em><script type="text/javascript" 
src="http://cdn.mathjax.org/mathjax/latest/MathJax.js?config=TeX-AMS-MML_HTMLorMML";></script></em></div>
+</div>
+<div class="subNav">
+<ul class="navList">
+<li>Prev</li>
+<li>Next</li>
+</ul>
+<ul class="navList">
+<li><a 
href="../../../../../../index.html?org/apache/commons/math3/special/class-use/BesselJ.html"
 target="_top">Frames</a></li>
+<li><a href="BesselJ.html" target="_top">No Frames</a></li>
+</ul>
+<ul class="navList" id="allclasses_navbar_top">
+<li><a href="../../../../../../allclasses-noframe.html">All Classes</a></li>
+</ul>
+<div>
+<script type="text/javascript"><!--
+  allClassesLink = document.getElementById("allclasses_navbar_top");
+  if(window==top) {
+    allClassesLink.style.display = "block";
+  }
+  else {
+    allClassesLink.style.display = "none";
+  }
+  //-->
+</script>
+</div>
+<a name="skip-navbar_top">
+<!--   -->
+</a></div>
+<!-- ========= END OF TOP NAVBAR ========= -->
+<div class="header">
+<h2 title="Uses of Class org.apache.commons.math3.special.BesselJ" 
class="title">Uses of Class<br>org.apache.commons.math3.special.BesselJ</h2>
+</div>
+<div class="classUseContainer">No usage of 
org.apache.commons.math3.special.BesselJ</div>
+<!-- ======= START OF BOTTOM NAVBAR ====== -->
+<div class="bottomNav"><a name="navbar_bottom">
+<!--   -->
+</a><a href="#skip-navbar_bottom" title="Skip navigation links"></a><a 
name="navbar_bottom_firstrow">
+<!--   -->
+</a>
+<ul class="navList" title="Navigation">
+<li><a href="../../../../../../overview-summary.html">Overview</a></li>
+<li><a href="../package-summary.html">Package</a></li>
+<li><a href="../../../../../../org/apache/commons/math3/special/BesselJ.html" 
title="class in org.apache.commons.math3.special">Class</a></li>
+<li class="navBarCell1Rev">Use</li>
+<li><a href="../package-tree.html">Tree</a></li>
+<li><a href="../../../../../../deprecated-list.html">Deprecated</a></li>
+<li><a href="../../../../../../index-all.html">Index</a></li>
+<li><a href="../../../../../../help-doc.html">Help</a></li>
+</ul>
+<div class="aboutLanguage"><em><script type="text/javascript" 
src="http://cdn.mathjax.org/mathjax/latest/MathJax.js?config=TeX-AMS-MML_HTMLorMML";></script></em></div>
+</div>
+<div class="subNav">
+<ul class="navList">
+<li>Prev</li>
+<li>Next</li>
+</ul>
+<ul class="navList">
+<li><a 
href="../../../../../../index.html?org/apache/commons/math3/special/class-use/BesselJ.html"
 target="_top">Frames</a></li>
+<li><a href="BesselJ.html" target="_top">No Frames</a></li>
+</ul>
+<ul class="navList" id="allclasses_navbar_bottom">
+<li><a href="../../../../../../allclasses-noframe.html">All Classes</a></li>
+</ul>
+<div>
+<script type="text/javascript"><!--
+  allClassesLink = document.getElementById("allclasses_navbar_bottom");
+  if(window==top) {
+    allClassesLink.style.display = "block";
+  }
+  else {
+    allClassesLink.style.display = "none";
+  }
+  //-->
+</script>
+</div>
+<a name="skip-navbar_bottom">
+<!--   -->
+</a></div>
+<!-- ======== END OF BOTTOM NAVBAR ======= -->
+<p class="legalCopy"><small>Copyright &#169; 2003&#x2013;2014 <a 
href="http://www.apache.org/";>The Apache Software Foundation</a>. All rights 
reserved.</small></p>
+</body>
+</html>

Propchange: 
websites/production/commons/content/proper/commons-math/apidocs/org/apache/commons/math3/special/class-use/BesselJ.html
------------------------------------------------------------------------------
    svn:eol-style = native

Added: 
websites/production/commons/content/proper/commons-math/apidocs/src-html/org/apache/commons/math3/analysis/interpolation/BicubicInterpolatingFunction.html
==============================================================================
--- 
websites/production/commons/content/proper/commons-math/apidocs/src-html/org/apache/commons/math3/analysis/interpolation/BicubicInterpolatingFunction.html
 (added)
+++ 
websites/production/commons/content/proper/commons-math/apidocs/src-html/org/apache/commons/math3/analysis/interpolation/BicubicInterpolatingFunction.html
 Tue Dec 30 17:15:36 2014
@@ -0,0 +1,397 @@
+<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN" 
"http://www.w3.org/TR/html4/loose.dtd";>
+<html lang="en">
+<head>
+<title>Source code</title>
+<link rel="stylesheet" type="text/css" 
href="../../../../../../../stylesheet.css" title="Style">
+</head>
+<body>
+<div class="sourceContainer">
+<pre><span class="sourceLineNo">001</span>/*<a name="line.1"></a>
+<span class="sourceLineNo">002</span> * Licensed to the Apache Software 
Foundation (ASF) under one or more<a name="line.2"></a>
+<span class="sourceLineNo">003</span> * contributor license agreements.  See 
the NOTICE file distributed with<a name="line.3"></a>
+<span class="sourceLineNo">004</span> * this work for additional information 
regarding copyright ownership.<a name="line.4"></a>
+<span class="sourceLineNo">005</span> * The ASF licenses this file to You 
under the Apache License, Version 2.0<a name="line.5"></a>
+<span class="sourceLineNo">006</span> * (the "License"); you may not use this 
file except in compliance with<a name="line.6"></a>
+<span class="sourceLineNo">007</span> * the License.  You may obtain a copy of 
the License at<a name="line.7"></a>
+<span class="sourceLineNo">008</span> *<a name="line.8"></a>
+<span class="sourceLineNo">009</span> *      
http://www.apache.org/licenses/LICENSE-2.0<a name="line.9"></a>
+<span class="sourceLineNo">010</span> *<a name="line.10"></a>
+<span class="sourceLineNo">011</span> * Unless required by applicable law or 
agreed to in writing, software<a name="line.11"></a>
+<span class="sourceLineNo">012</span> * distributed under the License is 
distributed on an "AS IS" BASIS,<a name="line.12"></a>
+<span class="sourceLineNo">013</span> * WITHOUT WARRANTIES OR CONDITIONS OF 
ANY KIND, either express or implied.<a name="line.13"></a>
+<span class="sourceLineNo">014</span> * See the License for the specific 
language governing permissions and<a name="line.14"></a>
+<span class="sourceLineNo">015</span> * limitations under the License.<a 
name="line.15"></a>
+<span class="sourceLineNo">016</span> */<a name="line.16"></a>
+<span class="sourceLineNo">017</span>package 
org.apache.commons.math3.analysis.interpolation;<a name="line.17"></a>
+<span class="sourceLineNo">018</span><a name="line.18"></a>
+<span class="sourceLineNo">019</span>import java.util.Arrays;<a 
name="line.19"></a>
+<span class="sourceLineNo">020</span>import 
org.apache.commons.math3.analysis.BivariateFunction;<a name="line.20"></a>
+<span class="sourceLineNo">021</span>import 
org.apache.commons.math3.exception.DimensionMismatchException;<a 
name="line.21"></a>
+<span class="sourceLineNo">022</span>import 
org.apache.commons.math3.exception.NoDataException;<a name="line.22"></a>
+<span class="sourceLineNo">023</span>import 
org.apache.commons.math3.exception.OutOfRangeException;<a name="line.23"></a>
+<span class="sourceLineNo">024</span>import 
org.apache.commons.math3.exception.NonMonotonicSequenceException;<a 
name="line.24"></a>
+<span class="sourceLineNo">025</span>import 
org.apache.commons.math3.util.MathArrays;<a name="line.25"></a>
+<span class="sourceLineNo">026</span><a name="line.26"></a>
+<span class="sourceLineNo">027</span>/**<a name="line.27"></a>
+<span class="sourceLineNo">028</span> * Function that implements the<a 
name="line.28"></a>
+<span class="sourceLineNo">029</span> * &lt;a 
href="http://en.wikipedia.org/wiki/Bicubic_interpolation"&gt;<a 
name="line.29"></a>
+<span class="sourceLineNo">030</span> * bicubic spline 
interpolation&lt;/a&gt;.<a name="line.30"></a>
+<span class="sourceLineNo">031</span> *<a name="line.31"></a>
+<span class="sourceLineNo">032</span> * @since 3.4<a name="line.32"></a>
+<span class="sourceLineNo">033</span> */<a name="line.33"></a>
+<span class="sourceLineNo">034</span>public class 
BicubicInterpolatingFunction<a name="line.34"></a>
+<span class="sourceLineNo">035</span>    implements BivariateFunction {<a 
name="line.35"></a>
+<span class="sourceLineNo">036</span>    /** Number of coefficients. */<a 
name="line.36"></a>
+<span class="sourceLineNo">037</span>    private static final int NUM_COEFF = 
16;<a name="line.37"></a>
+<span class="sourceLineNo">038</span>    /**<a name="line.38"></a>
+<span class="sourceLineNo">039</span>     * Matrix to compute the spline 
coefficients from the function values<a name="line.39"></a>
+<span class="sourceLineNo">040</span>     * and function derivatives values<a 
name="line.40"></a>
+<span class="sourceLineNo">041</span>     */<a name="line.41"></a>
+<span class="sourceLineNo">042</span>    private static final double[][] AINV 
= {<a name="line.42"></a>
+<span class="sourceLineNo">043</span>        { 1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0 
},<a name="line.43"></a>
+<span class="sourceLineNo">044</span>        { 0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0 
},<a name="line.44"></a>
+<span class="sourceLineNo">045</span>        { 
-3,3,0,0,-2,-1,0,0,0,0,0,0,0,0,0,0 },<a name="line.45"></a>
+<span class="sourceLineNo">046</span>        { 
2,-2,0,0,1,1,0,0,0,0,0,0,0,0,0,0 },<a name="line.46"></a>
+<span class="sourceLineNo">047</span>        { 0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0 
},<a name="line.47"></a>
+<span class="sourceLineNo">048</span>        { 0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0 
},<a name="line.48"></a>
+<span class="sourceLineNo">049</span>        { 
0,0,0,0,0,0,0,0,-3,3,0,0,-2,-1,0,0 },<a name="line.49"></a>
+<span class="sourceLineNo">050</span>        { 
0,0,0,0,0,0,0,0,2,-2,0,0,1,1,0,0 },<a name="line.50"></a>
+<span class="sourceLineNo">051</span>        { 
-3,0,3,0,0,0,0,0,-2,0,-1,0,0,0,0,0 },<a name="line.51"></a>
+<span class="sourceLineNo">052</span>        { 
0,0,0,0,-3,0,3,0,0,0,0,0,-2,0,-1,0 },<a name="line.52"></a>
+<span class="sourceLineNo">053</span>        { 
9,-9,-9,9,6,3,-6,-3,6,-6,3,-3,4,2,2,1 },<a name="line.53"></a>
+<span class="sourceLineNo">054</span>        { 
-6,6,6,-6,-3,-3,3,3,-4,4,-2,2,-2,-2,-1,-1 },<a name="line.54"></a>
+<span class="sourceLineNo">055</span>        { 
2,0,-2,0,0,0,0,0,1,0,1,0,0,0,0,0 },<a name="line.55"></a>
+<span class="sourceLineNo">056</span>        { 
0,0,0,0,2,0,-2,0,0,0,0,0,1,0,1,0 },<a name="line.56"></a>
+<span class="sourceLineNo">057</span>        { 
-6,6,6,-6,-4,-2,4,2,-3,3,-3,3,-2,-1,-2,-1 },<a name="line.57"></a>
+<span class="sourceLineNo">058</span>        { 
4,-4,-4,4,2,2,-2,-2,2,-2,2,-2,1,1,1,1 }<a name="line.58"></a>
+<span class="sourceLineNo">059</span>    };<a name="line.59"></a>
+<span class="sourceLineNo">060</span><a name="line.60"></a>
+<span class="sourceLineNo">061</span>    /** Samples x-coordinates */<a 
name="line.61"></a>
+<span class="sourceLineNo">062</span>    private final double[] xval;<a 
name="line.62"></a>
+<span class="sourceLineNo">063</span>    /** Samples y-coordinates */<a 
name="line.63"></a>
+<span class="sourceLineNo">064</span>    private final double[] yval;<a 
name="line.64"></a>
+<span class="sourceLineNo">065</span>    /** Set of cubic splines patching the 
whole data grid */<a name="line.65"></a>
+<span class="sourceLineNo">066</span>    private final BicubicFunction[][] 
splines;<a name="line.66"></a>
+<span class="sourceLineNo">067</span><a name="line.67"></a>
+<span class="sourceLineNo">068</span>    /**<a name="line.68"></a>
+<span class="sourceLineNo">069</span>     * @param x Sample values of the 
x-coordinate, in increasing order.<a name="line.69"></a>
+<span class="sourceLineNo">070</span>     * @param y Sample values of the 
y-coordinate, in increasing order.<a name="line.70"></a>
+<span class="sourceLineNo">071</span>     * @param f Values of the function on 
every grid point.<a name="line.71"></a>
+<span class="sourceLineNo">072</span>     * @param dFdX Values of the partial 
derivative of function with respect<a name="line.72"></a>
+<span class="sourceLineNo">073</span>     * to x on every grid point.<a 
name="line.73"></a>
+<span class="sourceLineNo">074</span>     * @param dFdY Values of the partial 
derivative of function with respect<a name="line.74"></a>
+<span class="sourceLineNo">075</span>     * to y on every grid point.<a 
name="line.75"></a>
+<span class="sourceLineNo">076</span>     * @param d2FdXdY Values of the cross 
partial derivative of function on<a name="line.76"></a>
+<span class="sourceLineNo">077</span>     * every grid point.<a 
name="line.77"></a>
+<span class="sourceLineNo">078</span>     * @throws DimensionMismatchException 
if the various arrays do not contain<a name="line.78"></a>
+<span class="sourceLineNo">079</span>     * the expected number of elements.<a 
name="line.79"></a>
+<span class="sourceLineNo">080</span>     * @throws 
NonMonotonicSequenceException if {@code x} or {@code y} are<a 
name="line.80"></a>
+<span class="sourceLineNo">081</span>     * not strictly increasing.<a 
name="line.81"></a>
+<span class="sourceLineNo">082</span>     * @throws NoDataException if any of 
the arrays has zero length.<a name="line.82"></a>
+<span class="sourceLineNo">083</span>     */<a name="line.83"></a>
+<span class="sourceLineNo">084</span>    public 
BicubicInterpolatingFunction(double[] x,<a name="line.84"></a>
+<span class="sourceLineNo">085</span>                                        
double[] y,<a name="line.85"></a>
+<span class="sourceLineNo">086</span>                                        
double[][] f,<a name="line.86"></a>
+<span class="sourceLineNo">087</span>                                        
double[][] dFdX,<a name="line.87"></a>
+<span class="sourceLineNo">088</span>                                        
double[][] dFdY,<a name="line.88"></a>
+<span class="sourceLineNo">089</span>                                        
double[][] d2FdXdY)<a name="line.89"></a>
+<span class="sourceLineNo">090</span>        throws 
DimensionMismatchException,<a name="line.90"></a>
+<span class="sourceLineNo">091</span>               NoDataException,<a 
name="line.91"></a>
+<span class="sourceLineNo">092</span>               
NonMonotonicSequenceException {<a name="line.92"></a>
+<span class="sourceLineNo">093</span>        final int xLen = x.length;<a 
name="line.93"></a>
+<span class="sourceLineNo">094</span>        final int yLen = y.length;<a 
name="line.94"></a>
+<span class="sourceLineNo">095</span><a name="line.95"></a>
+<span class="sourceLineNo">096</span>        if (xLen == 0 || yLen == 0 || 
f.length == 0 || f[0].length == 0) {<a name="line.96"></a>
+<span class="sourceLineNo">097</span>            throw new 
NoDataException();<a name="line.97"></a>
+<span class="sourceLineNo">098</span>        }<a name="line.98"></a>
+<span class="sourceLineNo">099</span>        if (xLen != f.length) {<a 
name="line.99"></a>
+<span class="sourceLineNo">100</span>            throw new 
DimensionMismatchException(xLen, f.length);<a name="line.100"></a>
+<span class="sourceLineNo">101</span>        }<a name="line.101"></a>
+<span class="sourceLineNo">102</span>        if (xLen != dFdX.length) {<a 
name="line.102"></a>
+<span class="sourceLineNo">103</span>            throw new 
DimensionMismatchException(xLen, dFdX.length);<a name="line.103"></a>
+<span class="sourceLineNo">104</span>        }<a name="line.104"></a>
+<span class="sourceLineNo">105</span>        if (xLen != dFdY.length) {<a 
name="line.105"></a>
+<span class="sourceLineNo">106</span>            throw new 
DimensionMismatchException(xLen, dFdY.length);<a name="line.106"></a>
+<span class="sourceLineNo">107</span>        }<a name="line.107"></a>
+<span class="sourceLineNo">108</span>        if (xLen != d2FdXdY.length) {<a 
name="line.108"></a>
+<span class="sourceLineNo">109</span>            throw new 
DimensionMismatchException(xLen, d2FdXdY.length);<a name="line.109"></a>
+<span class="sourceLineNo">110</span>        }<a name="line.110"></a>
+<span class="sourceLineNo">111</span><a name="line.111"></a>
+<span class="sourceLineNo">112</span>        MathArrays.checkOrder(x);<a 
name="line.112"></a>
+<span class="sourceLineNo">113</span>        MathArrays.checkOrder(y);<a 
name="line.113"></a>
+<span class="sourceLineNo">114</span><a name="line.114"></a>
+<span class="sourceLineNo">115</span>        xval = x.clone();<a 
name="line.115"></a>
+<span class="sourceLineNo">116</span>        yval = y.clone();<a 
name="line.116"></a>
+<span class="sourceLineNo">117</span><a name="line.117"></a>
+<span class="sourceLineNo">118</span>        final int lastI = xLen - 1;<a 
name="line.118"></a>
+<span class="sourceLineNo">119</span>        final int lastJ = yLen - 1;<a 
name="line.119"></a>
+<span class="sourceLineNo">120</span>        splines = new 
BicubicFunction[lastI][lastJ];<a name="line.120"></a>
+<span class="sourceLineNo">121</span><a name="line.121"></a>
+<span class="sourceLineNo">122</span>        for (int i = 0; i &lt; lastI; 
i++) {<a name="line.122"></a>
+<span class="sourceLineNo">123</span>            if (f[i].length != yLen) {<a 
name="line.123"></a>
+<span class="sourceLineNo">124</span>                throw new 
DimensionMismatchException(f[i].length, yLen);<a name="line.124"></a>
+<span class="sourceLineNo">125</span>            }<a name="line.125"></a>
+<span class="sourceLineNo">126</span>            if (dFdX[i].length != yLen) 
{<a name="line.126"></a>
+<span class="sourceLineNo">127</span>                throw new 
DimensionMismatchException(dFdX[i].length, yLen);<a name="line.127"></a>
+<span class="sourceLineNo">128</span>            }<a name="line.128"></a>
+<span class="sourceLineNo">129</span>            if (dFdY[i].length != yLen) 
{<a name="line.129"></a>
+<span class="sourceLineNo">130</span>                throw new 
DimensionMismatchException(dFdY[i].length, yLen);<a name="line.130"></a>
+<span class="sourceLineNo">131</span>            }<a name="line.131"></a>
+<span class="sourceLineNo">132</span>            if (d2FdXdY[i].length != 
yLen) {<a name="line.132"></a>
+<span class="sourceLineNo">133</span>                throw new 
DimensionMismatchException(d2FdXdY[i].length, yLen);<a name="line.133"></a>
+<span class="sourceLineNo">134</span>            }<a name="line.134"></a>
+<span class="sourceLineNo">135</span>            final int ip1 = i + 1;<a 
name="line.135"></a>
+<span class="sourceLineNo">136</span>            final double xR = xval[ip1] - 
xval[i];<a name="line.136"></a>
+<span class="sourceLineNo">137</span>            for (int j = 0; j &lt; lastJ; 
j++) {<a name="line.137"></a>
+<span class="sourceLineNo">138</span>                final int jp1 = j + 1;<a 
name="line.138"></a>
+<span class="sourceLineNo">139</span>                final double yR = 
yval[jp1] - yval[j];<a name="line.139"></a>
+<span class="sourceLineNo">140</span>                final double xRyR = xR * 
yR;<a name="line.140"></a>
+<span class="sourceLineNo">141</span>                final double[] beta = new 
double[] {<a name="line.141"></a>
+<span class="sourceLineNo">142</span>                    f[i][j], f[ip1][j], 
f[i][jp1], f[ip1][jp1],<a name="line.142"></a>
+<span class="sourceLineNo">143</span>                    dFdX[i][j] * xR, 
dFdX[ip1][j] * xR, dFdX[i][jp1] * xR, dFdX[ip1][jp1] * xR,<a 
name="line.143"></a>
+<span class="sourceLineNo">144</span>                    dFdY[i][j] * yR, 
dFdY[ip1][j] * yR, dFdY[i][jp1] * yR, dFdY[ip1][jp1] * yR,<a 
name="line.144"></a>
+<span class="sourceLineNo">145</span>                    d2FdXdY[i][j] * xRyR, 
d2FdXdY[ip1][j] * xRyR, d2FdXdY[i][jp1] * xRyR, d2FdXdY[ip1][jp1] * xRyR<a 
name="line.145"></a>
+<span class="sourceLineNo">146</span>                };<a name="line.146"></a>
+<span class="sourceLineNo">147</span><a name="line.147"></a>
+<span class="sourceLineNo">148</span>                splines[i][j] = new 
BicubicFunction(computeSplineCoefficients(beta));<a name="line.148"></a>
+<span class="sourceLineNo">149</span>            }<a name="line.149"></a>
+<span class="sourceLineNo">150</span>        }<a name="line.150"></a>
+<span class="sourceLineNo">151</span>    }<a name="line.151"></a>
+<span class="sourceLineNo">152</span><a name="line.152"></a>
+<span class="sourceLineNo">153</span>    /**<a name="line.153"></a>
+<span class="sourceLineNo">154</span>     * {@inheritDoc}<a 
name="line.154"></a>
+<span class="sourceLineNo">155</span>     */<a name="line.155"></a>
+<span class="sourceLineNo">156</span>    public double value(double x, double 
y)<a name="line.156"></a>
+<span class="sourceLineNo">157</span>        throws OutOfRangeException {<a 
name="line.157"></a>
+<span class="sourceLineNo">158</span>        final int i = searchIndex(x, 
xval);<a name="line.158"></a>
+<span class="sourceLineNo">159</span>        final int j = searchIndex(y, 
yval);<a name="line.159"></a>
+<span class="sourceLineNo">160</span><a name="line.160"></a>
+<span class="sourceLineNo">161</span>        final double xN = (x - xval[i]) / 
(xval[i + 1] - xval[i]);<a name="line.161"></a>
+<span class="sourceLineNo">162</span>        final double yN = (y - yval[j]) / 
(yval[j + 1] - yval[j]);<a name="line.162"></a>
+<span class="sourceLineNo">163</span><a name="line.163"></a>
+<span class="sourceLineNo">164</span>        return splines[i][j].value(xN, 
yN);<a name="line.164"></a>
+<span class="sourceLineNo">165</span>    }<a name="line.165"></a>
+<span class="sourceLineNo">166</span><a name="line.166"></a>
+<span class="sourceLineNo">167</span>    /**<a name="line.167"></a>
+<span class="sourceLineNo">168</span>     * Indicates whether a point is 
within the interpolation range.<a name="line.168"></a>
+<span class="sourceLineNo">169</span>     *<a name="line.169"></a>
+<span class="sourceLineNo">170</span>     * @param x First coordinate.<a 
name="line.170"></a>
+<span class="sourceLineNo">171</span>     * @param y Second coordinate.<a 
name="line.171"></a>
+<span class="sourceLineNo">172</span>     * @return {@code true} if (x, y) is 
a valid point.<a name="line.172"></a>
+<span class="sourceLineNo">173</span>     */<a name="line.173"></a>
+<span class="sourceLineNo">174</span>    public boolean isValidPoint(double x, 
double y) {<a name="line.174"></a>
+<span class="sourceLineNo">175</span>        if (x &lt; xval[0] ||<a 
name="line.175"></a>
+<span class="sourceLineNo">176</span>            x &gt; xval[xval.length - 1] 
||<a name="line.176"></a>
+<span class="sourceLineNo">177</span>            y &lt; yval[0] ||<a 
name="line.177"></a>
+<span class="sourceLineNo">178</span>            y &gt; yval[yval.length - 1]) 
{<a name="line.178"></a>
+<span class="sourceLineNo">179</span>            return false;<a 
name="line.179"></a>
+<span class="sourceLineNo">180</span>        } else {<a name="line.180"></a>
+<span class="sourceLineNo">181</span>            return true;<a 
name="line.181"></a>
+<span class="sourceLineNo">182</span>        }<a name="line.182"></a>
+<span class="sourceLineNo">183</span>    }<a name="line.183"></a>
+<span class="sourceLineNo">184</span><a name="line.184"></a>
+<span class="sourceLineNo">185</span>    /**<a name="line.185"></a>
+<span class="sourceLineNo">186</span>     * @param c Coordinate.<a 
name="line.186"></a>
+<span class="sourceLineNo">187</span>     * @param val Coordinate samples.<a 
name="line.187"></a>
+<span class="sourceLineNo">188</span>     * @return the index in {@code val} 
corresponding to the interval<a name="line.188"></a>
+<span class="sourceLineNo">189</span>     * containing {@code c}.<a 
name="line.189"></a>
+<span class="sourceLineNo">190</span>     * @throws OutOfRangeException if 
{@code c} is out of the<a name="line.190"></a>
+<span class="sourceLineNo">191</span>     * range defined by the boundary 
values of {@code val}.<a name="line.191"></a>
+<span class="sourceLineNo">192</span>     */<a name="line.192"></a>
+<span class="sourceLineNo">193</span>    private int searchIndex(double c, 
double[] val) {<a name="line.193"></a>
+<span class="sourceLineNo">194</span>        final int r = 
Arrays.binarySearch(val, c);<a name="line.194"></a>
+<span class="sourceLineNo">195</span><a name="line.195"></a>
+<span class="sourceLineNo">196</span>        if (r == -1 ||<a 
name="line.196"></a>
+<span class="sourceLineNo">197</span>            r == -val.length - 1) {<a 
name="line.197"></a>
+<span class="sourceLineNo">198</span>            throw new 
OutOfRangeException(c, val[0], val[val.length - 1]);<a name="line.198"></a>
+<span class="sourceLineNo">199</span>        }<a name="line.199"></a>
+<span class="sourceLineNo">200</span><a name="line.200"></a>
+<span class="sourceLineNo">201</span>        if (r &lt; 0) {<a 
name="line.201"></a>
+<span class="sourceLineNo">202</span>            // "c" in within an 
interpolation sub-interval: Return the<a name="line.202"></a>
+<span class="sourceLineNo">203</span>            // index of the sample at the 
lower end of the sub-interval.<a name="line.203"></a>
+<span class="sourceLineNo">204</span>            return -r - 2;<a 
name="line.204"></a>
+<span class="sourceLineNo">205</span>        }<a name="line.205"></a>
+<span class="sourceLineNo">206</span>        final int last = val.length - 
1;<a name="line.206"></a>
+<span class="sourceLineNo">207</span>        if (r == last) {<a 
name="line.207"></a>
+<span class="sourceLineNo">208</span>            // "c" is the last sample of 
the range: Return the index<a name="line.208"></a>
+<span class="sourceLineNo">209</span>            // of the sample at the lower 
end of the last sub-interval.<a name="line.209"></a>
+<span class="sourceLineNo">210</span>            return last - 1;<a 
name="line.210"></a>
+<span class="sourceLineNo">211</span>        }<a name="line.211"></a>
+<span class="sourceLineNo">212</span><a name="line.212"></a>
+<span class="sourceLineNo">213</span>        // "c" is another sample point.<a 
name="line.213"></a>
+<span class="sourceLineNo">214</span>        return r;<a name="line.214"></a>
+<span class="sourceLineNo">215</span>    }<a name="line.215"></a>
+<span class="sourceLineNo">216</span><a name="line.216"></a>
+<span class="sourceLineNo">217</span>    /**<a name="line.217"></a>
+<span class="sourceLineNo">218</span>     * Compute the spline coefficients 
from the list of function values and<a name="line.218"></a>
+<span class="sourceLineNo">219</span>     * function partial derivatives 
values at the four corners of a grid<a name="line.219"></a>
+<span class="sourceLineNo">220</span>     * element. They must be specified in 
the following order:<a name="line.220"></a>
+<span class="sourceLineNo">221</span>     * &lt;ul&gt;<a name="line.221"></a>
+<span class="sourceLineNo">222</span>     *  &lt;li&gt;f(0,0)&lt;/li&gt;<a 
name="line.222"></a>
+<span class="sourceLineNo">223</span>     *  &lt;li&gt;f(1,0)&lt;/li&gt;<a 
name="line.223"></a>
+<span class="sourceLineNo">224</span>     *  &lt;li&gt;f(0,1)&lt;/li&gt;<a 
name="line.224"></a>
+<span class="sourceLineNo">225</span>     *  &lt;li&gt;f(1,1)&lt;/li&gt;<a 
name="line.225"></a>
+<span class="sourceLineNo">226</span>     *  
&lt;li&gt;f&lt;sub&gt;x&lt;/sub&gt;(0,0)&lt;/li&gt;<a name="line.226"></a>
+<span class="sourceLineNo">227</span>     *  
&lt;li&gt;f&lt;sub&gt;x&lt;/sub&gt;(1,0)&lt;/li&gt;<a name="line.227"></a>
+<span class="sourceLineNo">228</span>     *  
&lt;li&gt;f&lt;sub&gt;x&lt;/sub&gt;(0,1)&lt;/li&gt;<a name="line.228"></a>
+<span class="sourceLineNo">229</span>     *  
&lt;li&gt;f&lt;sub&gt;x&lt;/sub&gt;(1,1)&lt;/li&gt;<a name="line.229"></a>
+<span class="sourceLineNo">230</span>     *  
&lt;li&gt;f&lt;sub&gt;y&lt;/sub&gt;(0,0)&lt;/li&gt;<a name="line.230"></a>
+<span class="sourceLineNo">231</span>     *  
&lt;li&gt;f&lt;sub&gt;y&lt;/sub&gt;(1,0)&lt;/li&gt;<a name="line.231"></a>
+<span class="sourceLineNo">232</span>     *  
&lt;li&gt;f&lt;sub&gt;y&lt;/sub&gt;(0,1)&lt;/li&gt;<a name="line.232"></a>
+<span class="sourceLineNo">233</span>     *  
&lt;li&gt;f&lt;sub&gt;y&lt;/sub&gt;(1,1)&lt;/li&gt;<a name="line.233"></a>
+<span class="sourceLineNo">234</span>     *  
&lt;li&gt;f&lt;sub&gt;xy&lt;/sub&gt;(0,0)&lt;/li&gt;<a name="line.234"></a>
+<span class="sourceLineNo">235</span>     *  
&lt;li&gt;f&lt;sub&gt;xy&lt;/sub&gt;(1,0)&lt;/li&gt;<a name="line.235"></a>
+<span class="sourceLineNo">236</span>     *  
&lt;li&gt;f&lt;sub&gt;xy&lt;/sub&gt;(0,1)&lt;/li&gt;<a name="line.236"></a>
+<span class="sourceLineNo">237</span>     *  
&lt;li&gt;f&lt;sub&gt;xy&lt;/sub&gt;(1,1)&lt;/li&gt;<a name="line.237"></a>
+<span class="sourceLineNo">238</span>     * &lt;/ul&gt;<a name="line.238"></a>
+<span class="sourceLineNo">239</span>     * where the subscripts indicate the 
partial derivative with respect to<a name="line.239"></a>
+<span class="sourceLineNo">240</span>     * the corresponding variable(s).<a 
name="line.240"></a>
+<span class="sourceLineNo">241</span>     *<a name="line.241"></a>
+<span class="sourceLineNo">242</span>     * @param beta List of function 
values and function partial derivatives<a name="line.242"></a>
+<span class="sourceLineNo">243</span>     * values.<a name="line.243"></a>
+<span class="sourceLineNo">244</span>     * @return the spline coefficients.<a 
name="line.244"></a>
+<span class="sourceLineNo">245</span>     */<a name="line.245"></a>
+<span class="sourceLineNo">246</span>    private double[] 
computeSplineCoefficients(double[] beta) {<a name="line.246"></a>
+<span class="sourceLineNo">247</span>        final double[] a = new 
double[NUM_COEFF];<a name="line.247"></a>
+<span class="sourceLineNo">248</span><a name="line.248"></a>
+<span class="sourceLineNo">249</span>        for (int i = 0; i &lt; NUM_COEFF; 
i++) {<a name="line.249"></a>
+<span class="sourceLineNo">250</span>            double result = 0;<a 
name="line.250"></a>
+<span class="sourceLineNo">251</span>            final double[] row = 
AINV[i];<a name="line.251"></a>
+<span class="sourceLineNo">252</span>            for (int j = 0; j &lt; 
NUM_COEFF; j++) {<a name="line.252"></a>
+<span class="sourceLineNo">253</span>                result += row[j] * 
beta[j];<a name="line.253"></a>
+<span class="sourceLineNo">254</span>            }<a name="line.254"></a>
+<span class="sourceLineNo">255</span>            a[i] = result;<a 
name="line.255"></a>
+<span class="sourceLineNo">256</span>        }<a name="line.256"></a>
+<span class="sourceLineNo">257</span><a name="line.257"></a>
+<span class="sourceLineNo">258</span>        return a;<a name="line.258"></a>
+<span class="sourceLineNo">259</span>    }<a name="line.259"></a>
+<span class="sourceLineNo">260</span>}<a name="line.260"></a>
+<span class="sourceLineNo">261</span><a name="line.261"></a>
+<span class="sourceLineNo">262</span>/**<a name="line.262"></a>
+<span class="sourceLineNo">263</span> * Bicubic function.<a 
name="line.263"></a>
+<span class="sourceLineNo">264</span> */<a name="line.264"></a>
+<span class="sourceLineNo">265</span>class BicubicFunction implements 
BivariateFunction {<a name="line.265"></a>
+<span class="sourceLineNo">266</span>    /** Number of points. */<a 
name="line.266"></a>
+<span class="sourceLineNo">267</span>    private static final short N = 4;<a 
name="line.267"></a>
+<span class="sourceLineNo">268</span>    /** Coefficients */<a 
name="line.268"></a>
+<span class="sourceLineNo">269</span>    private final double[][] a;<a 
name="line.269"></a>
+<span class="sourceLineNo">270</span><a name="line.270"></a>
+<span class="sourceLineNo">271</span>    /**<a name="line.271"></a>
+<span class="sourceLineNo">272</span>     * Simple constructor.<a 
name="line.272"></a>
+<span class="sourceLineNo">273</span>     *<a name="line.273"></a>
+<span class="sourceLineNo">274</span>     * @param coeff Spline 
coefficients.<a name="line.274"></a>
+<span class="sourceLineNo">275</span>     */<a name="line.275"></a>
+<span class="sourceLineNo">276</span>    public BicubicFunction(double[] 
coeff) {<a name="line.276"></a>
+<span class="sourceLineNo">277</span>        a = new double[N][N];<a 
name="line.277"></a>
+<span class="sourceLineNo">278</span>        for (int j = 0; j &lt; N; j++) 
{<a name="line.278"></a>
+<span class="sourceLineNo">279</span>            final double[] aJ = a[j];<a 
name="line.279"></a>
+<span class="sourceLineNo">280</span>            for (int i = 0; i &lt; N; 
i++) {<a name="line.280"></a>
+<span class="sourceLineNo">281</span>                aJ[i] = coeff[i * N + 
j];<a name="line.281"></a>
+<span class="sourceLineNo">282</span>            }<a name="line.282"></a>
+<span class="sourceLineNo">283</span>        }<a name="line.283"></a>
+<span class="sourceLineNo">284</span>    }<a name="line.284"></a>
+<span class="sourceLineNo">285</span><a name="line.285"></a>
+<span class="sourceLineNo">286</span>    /**<a name="line.286"></a>
+<span class="sourceLineNo">287</span>     * {@inheritDoc}<a 
name="line.287"></a>
+<span class="sourceLineNo">288</span>     */<a name="line.288"></a>
+<span class="sourceLineNo">289</span>    public double value(double x, double 
y) {<a name="line.289"></a>
+<span class="sourceLineNo">290</span>        if (x &lt; 0 || x &gt; 1) {<a 
name="line.290"></a>
+<span class="sourceLineNo">291</span>            throw new 
OutOfRangeException(x, 0, 1);<a name="line.291"></a>
+<span class="sourceLineNo">292</span>        }<a name="line.292"></a>
+<span class="sourceLineNo">293</span>        if (y &lt; 0 || y &gt; 1) {<a 
name="line.293"></a>
+<span class="sourceLineNo">294</span>            throw new 
OutOfRangeException(y, 0, 1);<a name="line.294"></a>
+<span class="sourceLineNo">295</span>        }<a name="line.295"></a>
+<span class="sourceLineNo">296</span><a name="line.296"></a>
+<span class="sourceLineNo">297</span>        final double x2 = x * x;<a 
name="line.297"></a>
+<span class="sourceLineNo">298</span>        final double x3 = x2 * x;<a 
name="line.298"></a>
+<span class="sourceLineNo">299</span>        final double[] pX = {1, x, x2, 
x3};<a name="line.299"></a>
+<span class="sourceLineNo">300</span><a name="line.300"></a>
+<span class="sourceLineNo">301</span>        final double y2 = y * y;<a 
name="line.301"></a>
+<span class="sourceLineNo">302</span>        final double y3 = y2 * y;<a 
name="line.302"></a>
+<span class="sourceLineNo">303</span>        final double[] pY = {1, y, y2, 
y3};<a name="line.303"></a>
+<span class="sourceLineNo">304</span><a name="line.304"></a>
+<span class="sourceLineNo">305</span>        return apply(pX, pY, a);<a 
name="line.305"></a>
+<span class="sourceLineNo">306</span>    }<a name="line.306"></a>
+<span class="sourceLineNo">307</span><a name="line.307"></a>
+<span class="sourceLineNo">308</span>    /**<a name="line.308"></a>
+<span class="sourceLineNo">309</span>     * Compute the value of the bicubic 
polynomial.<a name="line.309"></a>
+<span class="sourceLineNo">310</span>     *<a name="line.310"></a>
+<span class="sourceLineNo">311</span>     * @param pX Powers of the 
x-coordinate.<a name="line.311"></a>
+<span class="sourceLineNo">312</span>     * @param pY Powers of the 
y-coordinate.<a name="line.312"></a>
+<span class="sourceLineNo">313</span>     * @param coeff Spline 
coefficients.<a name="line.313"></a>
+<span class="sourceLineNo">314</span>     * @return the interpolated value.<a 
name="line.314"></a>
+<span class="sourceLineNo">315</span>     */<a name="line.315"></a>
+<span class="sourceLineNo">316</span>    private double apply(double[] pX, 
double[] pY, double[][] coeff) {<a name="line.316"></a>
+<span class="sourceLineNo">317</span>        double result = 0;<a 
name="line.317"></a>
+<span class="sourceLineNo">318</span>        for (int i = 0; i &lt; N; i++) 
{<a name="line.318"></a>
+<span class="sourceLineNo">319</span>            final double r = 
MathArrays.linearCombination(coeff[i], pY);<a name="line.319"></a>
+<span class="sourceLineNo">320</span>            result += r * pX[i];<a 
name="line.320"></a>
+<span class="sourceLineNo">321</span>        }<a name="line.321"></a>
+<span class="sourceLineNo">322</span><a name="line.322"></a>
+<span class="sourceLineNo">323</span>        return result;<a 
name="line.323"></a>
+<span class="sourceLineNo">324</span>    }<a name="line.324"></a>
+<span class="sourceLineNo">325</span>}<a name="line.325"></a>
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+</pre>
+</div>
+</body>
+</html>

Propchange: 
websites/production/commons/content/proper/commons-math/apidocs/src-html/org/apache/commons/math3/analysis/interpolation/BicubicInterpolatingFunction.html
------------------------------------------------------------------------------
    svn:eol-style = native

Added: 
websites/production/commons/content/proper/commons-math/apidocs/src-html/org/apache/commons/math3/analysis/interpolation/BicubicInterpolator.html
==============================================================================
--- 
websites/production/commons/content/proper/commons-math/apidocs/src-html/org/apache/commons/math3/analysis/interpolation/BicubicInterpolator.html
 (added)
+++ 
websites/production/commons/content/proper/commons-math/apidocs/src-html/org/apache/commons/math3/analysis/interpolation/BicubicInterpolator.html
 Tue Dec 30 17:15:36 2014
@@ -0,0 +1,184 @@
+<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN" 
"http://www.w3.org/TR/html4/loose.dtd";>
+<html lang="en">
+<head>
+<title>Source code</title>
+<link rel="stylesheet" type="text/css" 
href="../../../../../../../stylesheet.css" title="Style">
+</head>
+<body>
+<div class="sourceContainer">
+<pre><span class="sourceLineNo">001</span>/*<a name="line.1"></a>
+<span class="sourceLineNo">002</span> * Licensed to the Apache Software 
Foundation (ASF) under one or more<a name="line.2"></a>
+<span class="sourceLineNo">003</span> * contributor license agreements.  See 
the NOTICE file distributed with<a name="line.3"></a>
+<span class="sourceLineNo">004</span> * this work for additional information 
regarding copyright ownership.<a name="line.4"></a>
+<span class="sourceLineNo">005</span> * The ASF licenses this file to You 
under the Apache License, Version 2.0<a name="line.5"></a>
+<span class="sourceLineNo">006</span> * (the "License"); you may not use this 
file except in compliance with<a name="line.6"></a>
+<span class="sourceLineNo">007</span> * the License.  You may obtain a copy of 
the License at<a name="line.7"></a>
+<span class="sourceLineNo">008</span> *<a name="line.8"></a>
+<span class="sourceLineNo">009</span> *      
http://www.apache.org/licenses/LICENSE-2.0<a name="line.9"></a>
+<span class="sourceLineNo">010</span> *<a name="line.10"></a>
+<span class="sourceLineNo">011</span> * Unless required by applicable law or 
agreed to in writing, software<a name="line.11"></a>
+<span class="sourceLineNo">012</span> * distributed under the License is 
distributed on an "AS IS" BASIS,<a name="line.12"></a>
+<span class="sourceLineNo">013</span> * WITHOUT WARRANTIES OR CONDITIONS OF 
ANY KIND, either express or implied.<a name="line.13"></a>
+<span class="sourceLineNo">014</span> * See the License for the specific 
language governing permissions and<a name="line.14"></a>
+<span class="sourceLineNo">015</span> * limitations under the License.<a 
name="line.15"></a>
+<span class="sourceLineNo">016</span> */<a name="line.16"></a>
+<span class="sourceLineNo">017</span>package 
org.apache.commons.math3.analysis.interpolation;<a name="line.17"></a>
+<span class="sourceLineNo">018</span><a name="line.18"></a>
+<span class="sourceLineNo">019</span>import 
org.apache.commons.math3.exception.DimensionMismatchException;<a 
name="line.19"></a>
+<span class="sourceLineNo">020</span>import 
org.apache.commons.math3.exception.NoDataException;<a name="line.20"></a>
+<span class="sourceLineNo">021</span>import 
org.apache.commons.math3.exception.NonMonotonicSequenceException;<a 
name="line.21"></a>
+<span class="sourceLineNo">022</span>import 
org.apache.commons.math3.exception.NumberIsTooSmallException;<a 
name="line.22"></a>
+<span class="sourceLineNo">023</span>import 
org.apache.commons.math3.util.MathArrays;<a name="line.23"></a>
+<span class="sourceLineNo">024</span><a name="line.24"></a>
+<span class="sourceLineNo">025</span>/**<a name="line.25"></a>
+<span class="sourceLineNo">026</span> * Generates a {@link 
BicubicInterpolatingFunction bicubic interpolating<a name="line.26"></a>
+<span class="sourceLineNo">027</span> * function}.<a name="line.27"></a>
+<span class="sourceLineNo">028</span> * &lt;p&gt;<a name="line.28"></a>
+<span class="sourceLineNo">029</span> *  Caveat: Because the interpolation 
scheme requires that derivatives be<a name="line.29"></a>
+<span class="sourceLineNo">030</span> *  specified at the sample points, those 
are approximated with finite<a name="line.30"></a>
+<span class="sourceLineNo">031</span> *  differences (using the 2-points 
symmetric formulae).<a name="line.31"></a>
+<span class="sourceLineNo">032</span> *  Since their values are undefined at 
the borders of the provided<a name="line.32"></a>
+<span class="sourceLineNo">033</span> *  interpolation ranges, the 
interpolated values will be wrong at the<a name="line.33"></a>
+<span class="sourceLineNo">034</span> *  edges of the patch.<a 
name="line.34"></a>
+<span class="sourceLineNo">035</span> *  The {@code interpolate} method will 
return a function that overrides<a name="line.35"></a>
+<span class="sourceLineNo">036</span> *  {@link 
BicubicInterpolatingFunction#isValidPoint(double,double)} to<a 
name="line.36"></a>
+<span class="sourceLineNo">037</span> *  indicate points where the 
interpolation will be inaccurate.<a name="line.37"></a>
+<span class="sourceLineNo">038</span> * &lt;/p&gt;<a name="line.38"></a>
+<span class="sourceLineNo">039</span> *<a name="line.39"></a>
+<span class="sourceLineNo">040</span> * @since 3.4<a name="line.40"></a>
+<span class="sourceLineNo">041</span> */<a name="line.41"></a>
+<span class="sourceLineNo">042</span>public class BicubicInterpolator<a 
name="line.42"></a>
+<span class="sourceLineNo">043</span>    implements BivariateGridInterpolator 
{<a name="line.43"></a>
+<span class="sourceLineNo">044</span>    /**<a name="line.44"></a>
+<span class="sourceLineNo">045</span>     * {@inheritDoc}<a name="line.45"></a>
+<span class="sourceLineNo">046</span>     */<a name="line.46"></a>
+<span class="sourceLineNo">047</span>    public BicubicInterpolatingFunction 
interpolate(final double[] xval,<a name="line.47"></a>
+<span class="sourceLineNo">048</span>                                          
          final double[] yval,<a name="line.48"></a>
+<span class="sourceLineNo">049</span>                                          
          final double[][] fval)<a name="line.49"></a>
+<span class="sourceLineNo">050</span>        throws NoDataException, 
DimensionMismatchException,<a name="line.50"></a>
+<span class="sourceLineNo">051</span>               
NonMonotonicSequenceException, NumberIsTooSmallException {<a name="line.51"></a>
+<span class="sourceLineNo">052</span>        if (xval.length == 0 || 
yval.length == 0 || fval.length == 0) {<a name="line.52"></a>
+<span class="sourceLineNo">053</span>            throw new 
NoDataException();<a name="line.53"></a>
+<span class="sourceLineNo">054</span>        }<a name="line.54"></a>
+<span class="sourceLineNo">055</span>        if (xval.length != fval.length) 
{<a name="line.55"></a>
+<span class="sourceLineNo">056</span>            throw new 
DimensionMismatchException(xval.length, fval.length);<a name="line.56"></a>
+<span class="sourceLineNo">057</span>        }<a name="line.57"></a>
+<span class="sourceLineNo">058</span><a name="line.58"></a>
+<span class="sourceLineNo">059</span>        MathArrays.checkOrder(xval);<a 
name="line.59"></a>
+<span class="sourceLineNo">060</span>        MathArrays.checkOrder(yval);<a 
name="line.60"></a>
+<span class="sourceLineNo">061</span><a name="line.61"></a>
+<span class="sourceLineNo">062</span>        final int xLen = xval.length;<a 
name="line.62"></a>
+<span class="sourceLineNo">063</span>        final int yLen = yval.length;<a 
name="line.63"></a>
+<span class="sourceLineNo">064</span><a name="line.64"></a>
+<span class="sourceLineNo">065</span>        // Approximation to the partial 
derivatives using finite differences.<a name="line.65"></a>
+<span class="sourceLineNo">066</span>        final double[][] dFdX = new 
double[xLen][yLen];<a name="line.66"></a>
+<span class="sourceLineNo">067</span>        final double[][] dFdY = new 
double[xLen][yLen];<a name="line.67"></a>
+<span class="sourceLineNo">068</span>        final double[][] d2FdXdY = new 
double[xLen][yLen];<a name="line.68"></a>
+<span class="sourceLineNo">069</span>        for (int i = 1; i &lt; xLen - 1; 
i++) {<a name="line.69"></a>
+<span class="sourceLineNo">070</span>            final int nI = i + 1;<a 
name="line.70"></a>
+<span class="sourceLineNo">071</span>            final int pI = i - 1;<a 
name="line.71"></a>
+<span class="sourceLineNo">072</span><a name="line.72"></a>
+<span class="sourceLineNo">073</span>            final double nX = xval[nI];<a 
name="line.73"></a>
+<span class="sourceLineNo">074</span>            final double pX = xval[pI];<a 
name="line.74"></a>
+<span class="sourceLineNo">075</span><a name="line.75"></a>
+<span class="sourceLineNo">076</span>            final double deltaX = nX - 
pX;<a name="line.76"></a>
+<span class="sourceLineNo">077</span><a name="line.77"></a>
+<span class="sourceLineNo">078</span>            for (int j = 1; j &lt; yLen - 
1; j++) {<a name="line.78"></a>
+<span class="sourceLineNo">079</span>                final int nJ = j + 1;<a 
name="line.79"></a>
+<span class="sourceLineNo">080</span>                final int pJ = j - 1;<a 
name="line.80"></a>
+<span class="sourceLineNo">081</span><a name="line.81"></a>
+<span class="sourceLineNo">082</span>                final double nY = 
yval[nJ];<a name="line.82"></a>
+<span class="sourceLineNo">083</span>                final double pY = 
yval[pJ];<a name="line.83"></a>
+<span class="sourceLineNo">084</span><a name="line.84"></a>
+<span class="sourceLineNo">085</span>                final double deltaY = nY 
- pY;<a name="line.85"></a>
+<span class="sourceLineNo">086</span><a name="line.86"></a>
+<span class="sourceLineNo">087</span>                dFdX[i][j] = (fval[nI][j] 
- fval[pI][j]) / deltaX;<a name="line.87"></a>
+<span class="sourceLineNo">088</span>                dFdY[i][j] = (fval[i][nJ] 
- fval[i][pJ]) / deltaY;<a name="line.88"></a>
+<span class="sourceLineNo">089</span><a name="line.89"></a>
+<span class="sourceLineNo">090</span>                final double deltaXY = 
deltaX * deltaY;<a name="line.90"></a>
+<span class="sourceLineNo">091</span><a name="line.91"></a>
+<span class="sourceLineNo">092</span>                d2FdXdY[i][j] = 
(fval[nI][nJ] - fval[nI][pJ] - fval[pI][nJ] + fval[pI][pJ]) / deltaXY;<a 
name="line.92"></a>
+<span class="sourceLineNo">093</span>            }<a name="line.93"></a>
+<span class="sourceLineNo">094</span>        }<a name="line.94"></a>
+<span class="sourceLineNo">095</span><a name="line.95"></a>
+<span class="sourceLineNo">096</span>        // Create the interpolating 
function.<a name="line.96"></a>
+<span class="sourceLineNo">097</span>        return new 
BicubicInterpolatingFunction(xval, yval, fval,<a name="line.97"></a>
+<span class="sourceLineNo">098</span>                                          
      dFdX, dFdY, d2FdXdY) {<a name="line.98"></a>
+<span class="sourceLineNo">099</span>            @Override<a 
name="line.99"></a>
+<span class="sourceLineNo">100</span>            public boolean 
isValidPoint(double x, double y) {<a name="line.100"></a>
+<span class="sourceLineNo">101</span>                if (x &lt; xval[1] ||<a 
name="line.101"></a>
+<span class="sourceLineNo">102</span>                    x &gt; 
xval[xval.length - 2] ||<a name="line.102"></a>
+<span class="sourceLineNo">103</span>                    y &lt; yval[1] ||<a 
name="line.103"></a>
+<span class="sourceLineNo">104</span>                    y &gt; 
yval[yval.length - 2]) {<a name="line.104"></a>
+<span class="sourceLineNo">105</span>                    return false;<a 
name="line.105"></a>
+<span class="sourceLineNo">106</span>                } else {<a 
name="line.106"></a>
+<span class="sourceLineNo">107</span>                    return true;<a 
name="line.107"></a>
+<span class="sourceLineNo">108</span>                }<a name="line.108"></a>
+<span class="sourceLineNo">109</span>            }<a name="line.109"></a>
+<span class="sourceLineNo">110</span>        };<a name="line.110"></a>
+<span class="sourceLineNo">111</span>    }<a name="line.111"></a>
+<span class="sourceLineNo">112</span>}<a name="line.112"></a>
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+</pre>
+</div>
+</body>
+</html>

Propchange: 
websites/production/commons/content/proper/commons-math/apidocs/src-html/org/apache/commons/math3/analysis/interpolation/BicubicInterpolator.html
------------------------------------------------------------------------------
    svn:eol-style = native

Added: 
websites/production/commons/content/proper/commons-math/apidocs/src-html/org/apache/commons/math3/analysis/interpolation/PiecewiseBicubicSplineInterpolatingFunction.html
==============================================================================
--- 
websites/production/commons/content/proper/commons-math/apidocs/src-html/org/apache/commons/math3/analysis/interpolation/PiecewiseBicubicSplineInterpolatingFunction.html
 (added)
+++ 
websites/production/commons/content/proper/commons-math/apidocs/src-html/org/apache/commons/math3/analysis/interpolation/PiecewiseBicubicSplineInterpolatingFunction.html
 Tue Dec 30 17:15:36 2014
@@ -0,0 +1,282 @@
+<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN" 
"http://www.w3.org/TR/html4/loose.dtd";>
+<html lang="en">
+<head>
+<title>Source code</title>
+<link rel="stylesheet" type="text/css" 
href="../../../../../../../stylesheet.css" title="Style">
+</head>
+<body>
+<div class="sourceContainer">
+<pre><span class="sourceLineNo">001</span>/*<a name="line.1"></a>
+<span class="sourceLineNo">002</span> * Licensed to the Apache Software 
Foundation (ASF) under one or more<a name="line.2"></a>
+<span class="sourceLineNo">003</span> * contributor license agreements.  See 
the NOTICE file distributed with<a name="line.3"></a>
+<span class="sourceLineNo">004</span> * this work for additional information 
regarding copyright ownership.<a name="line.4"></a>
+<span class="sourceLineNo">005</span> * The ASF licenses this file to You 
under the Apache License, Version 2.0<a name="line.5"></a>
+<span class="sourceLineNo">006</span> * (the "License"); you may not use this 
file except in compliance with<a name="line.6"></a>
+<span class="sourceLineNo">007</span> * the License.  You may obtain a copy of 
the License at<a name="line.7"></a>
+<span class="sourceLineNo">008</span> *<a name="line.8"></a>
+<span class="sourceLineNo">009</span> *      
http://www.apache.org/licenses/LICENSE-2.0<a name="line.9"></a>
+<span class="sourceLineNo">010</span> *<a name="line.10"></a>
+<span class="sourceLineNo">011</span> * Unless required by applicable law or 
agreed to in writing, software<a name="line.11"></a>
+<span class="sourceLineNo">012</span> * distributed under the License is 
distributed on an "AS IS" BASIS,<a name="line.12"></a>
+<span class="sourceLineNo">013</span> * WITHOUT WARRANTIES OR CONDITIONS OF 
ANY KIND, either express or implied.<a name="line.13"></a>
+<span class="sourceLineNo">014</span> * See the License for the specific 
language governing permissions and<a name="line.14"></a>
+<span class="sourceLineNo">015</span> * limitations under the License.<a 
name="line.15"></a>
+<span class="sourceLineNo">016</span> */<a name="line.16"></a>
+<span class="sourceLineNo">017</span>package 
org.apache.commons.math3.analysis.interpolation;<a name="line.17"></a>
+<span class="sourceLineNo">018</span><a name="line.18"></a>
+<span class="sourceLineNo">019</span>import java.util.Arrays;<a 
name="line.19"></a>
+<span class="sourceLineNo">020</span>import 
org.apache.commons.math3.analysis.BivariateFunction;<a name="line.20"></a>
+<span class="sourceLineNo">021</span>import 
org.apache.commons.math3.analysis.polynomials.PolynomialSplineFunction;<a 
name="line.21"></a>
+<span class="sourceLineNo">022</span>import 
org.apache.commons.math3.exception.DimensionMismatchException;<a 
name="line.22"></a>
+<span class="sourceLineNo">023</span>import 
org.apache.commons.math3.exception.InsufficientDataException;<a 
name="line.23"></a>
+<span class="sourceLineNo">024</span>import 
org.apache.commons.math3.exception.NoDataException;<a name="line.24"></a>
+<span class="sourceLineNo">025</span>import 
org.apache.commons.math3.exception.NullArgumentException;<a name="line.25"></a>
+<span class="sourceLineNo">026</span>import 
org.apache.commons.math3.exception.OutOfRangeException;<a name="line.26"></a>
+<span class="sourceLineNo">027</span>import 
org.apache.commons.math3.exception.NonMonotonicSequenceException;<a 
name="line.27"></a>
+<span class="sourceLineNo">028</span>import 
org.apache.commons.math3.util.MathArrays;<a name="line.28"></a>
+<span class="sourceLineNo">029</span><a name="line.29"></a>
+<span class="sourceLineNo">030</span>/**<a name="line.30"></a>
+<span class="sourceLineNo">031</span> * Function that implements the<a 
name="line.31"></a>
+<span class="sourceLineNo">032</span> * &lt;a 
href="http://www.paulinternet.nl/?page=bicubic"&gt;bicubic spline&lt;/a&gt;<a 
name="line.32"></a>
+<span class="sourceLineNo">033</span> * interpolation.<a name="line.33"></a>
+<span class="sourceLineNo">034</span> * This implementation currently uses 
{@link AkimaSplineInterpolator} as the<a name="line.34"></a>
+<span class="sourceLineNo">035</span> * underlying one-dimensional 
interpolator, which requires 5 sample points;<a name="line.35"></a>
+<span class="sourceLineNo">036</span> * insufficient data will raise an 
exception when the<a name="line.36"></a>
+<span class="sourceLineNo">037</span> * {@link #value(double,double) value} 
method is called.<a name="line.37"></a>
+<span class="sourceLineNo">038</span> *<a name="line.38"></a>
+<span class="sourceLineNo">039</span> * @since 3.4<a name="line.39"></a>
+<span class="sourceLineNo">040</span> */<a name="line.40"></a>
+<span class="sourceLineNo">041</span>public class 
PiecewiseBicubicSplineInterpolatingFunction<a name="line.41"></a>
+<span class="sourceLineNo">042</span>    implements BivariateFunction {<a 
name="line.42"></a>
+<span class="sourceLineNo">043</span>    /** The minimum number of points that 
are needed to compute the function. */<a name="line.43"></a>
+<span class="sourceLineNo">044</span>    private static final int 
MIN_NUM_POINTS = 5;<a name="line.44"></a>
+<span class="sourceLineNo">045</span>    /** Samples x-coordinates */<a 
name="line.45"></a>
+<span class="sourceLineNo">046</span>    private final double[] xval;<a 
name="line.46"></a>
+<span class="sourceLineNo">047</span>    /** Samples y-coordinates */<a 
name="line.47"></a>
+<span class="sourceLineNo">048</span>    private final double[] yval;<a 
name="line.48"></a>
+<span class="sourceLineNo">049</span>    /** Set of cubic splines patching the 
whole data grid */<a name="line.49"></a>
+<span class="sourceLineNo">050</span>    private final double[][] fval;<a 
name="line.50"></a>
+<span class="sourceLineNo">051</span><a name="line.51"></a>
+<span class="sourceLineNo">052</span>    /**<a name="line.52"></a>
+<span class="sourceLineNo">053</span>     * @param x Sample values of the 
x-coordinate, in increasing order.<a name="line.53"></a>
+<span class="sourceLineNo">054</span>     * @param y Sample values of the 
y-coordinate, in increasing order.<a name="line.54"></a>
+<span class="sourceLineNo">055</span>     * @param f Values of the function on 
every grid point. the expected number<a name="line.55"></a>
+<span class="sourceLineNo">056</span>     *        of elements.<a 
name="line.56"></a>
+<span class="sourceLineNo">057</span>     * @throws 
NonMonotonicSequenceException if {@code x} or {@code y} are not<a 
name="line.57"></a>
+<span class="sourceLineNo">058</span>     *         strictly increasing.<a 
name="line.58"></a>
+<span class="sourceLineNo">059</span>     * @throws NullArgumentException if 
any of the arguments are null<a name="line.59"></a>
+<span class="sourceLineNo">060</span>     * @throws NoDataException if any of 
the arrays has zero length.<a name="line.60"></a>
+<span class="sourceLineNo">061</span>     * @throws DimensionMismatchException 
if the length of x and y don't match the row, column<a name="line.61"></a>
+<span class="sourceLineNo">062</span>     *         height of f<a 
name="line.62"></a>
+<span class="sourceLineNo">063</span>     */<a name="line.63"></a>
+<span class="sourceLineNo">064</span>    public 
PiecewiseBicubicSplineInterpolatingFunction(double[] x,<a name="line.64"></a>
+<span class="sourceLineNo">065</span>                                          
             double[] y,<a name="line.65"></a>
+<span class="sourceLineNo">066</span>                                          
             double[][] f)<a name="line.66"></a>
+<span class="sourceLineNo">067</span>        throws 
DimensionMismatchException,<a name="line.67"></a>
+<span class="sourceLineNo">068</span>               NullArgumentException,<a 
name="line.68"></a>
+<span class="sourceLineNo">069</span>               NoDataException,<a 
name="line.69"></a>
+<span class="sourceLineNo">070</span>               
NonMonotonicSequenceException {<a name="line.70"></a>
+<span class="sourceLineNo">071</span>        if (x == null ||<a 
name="line.71"></a>
+<span class="sourceLineNo">072</span>            y == null ||<a 
name="line.72"></a>
+<span class="sourceLineNo">073</span>            f == null ||<a 
name="line.73"></a>
+<span class="sourceLineNo">074</span>            f[0] == null) {<a 
name="line.74"></a>
+<span class="sourceLineNo">075</span>            throw new 
NullArgumentException();<a name="line.75"></a>
+<span class="sourceLineNo">076</span>        }<a name="line.76"></a>
+<span class="sourceLineNo">077</span><a name="line.77"></a>
+<span class="sourceLineNo">078</span>        final int xLen = x.length;<a 
name="line.78"></a>
+<span class="sourceLineNo">079</span>        final int yLen = y.length;<a 
name="line.79"></a>
+<span class="sourceLineNo">080</span><a name="line.80"></a>
+<span class="sourceLineNo">081</span>        if (xLen == 0 ||<a 
name="line.81"></a>
+<span class="sourceLineNo">082</span>            yLen == 0 ||<a 
name="line.82"></a>
+<span class="sourceLineNo">083</span>            f.length == 0 ||<a 
name="line.83"></a>
+<span class="sourceLineNo">084</span>            f[0].length == 0) {<a 
name="line.84"></a>
+<span class="sourceLineNo">085</span>            throw new 
NoDataException();<a name="line.85"></a>
+<span class="sourceLineNo">086</span>        }<a name="line.86"></a>
+<span class="sourceLineNo">087</span><a name="line.87"></a>
+<span class="sourceLineNo">088</span>        if (xLen &lt; MIN_NUM_POINTS ||<a 
name="line.88"></a>
+<span class="sourceLineNo">089</span>            yLen &lt; MIN_NUM_POINTS ||<a 
name="line.89"></a>
+<span class="sourceLineNo">090</span>            f.length &lt; MIN_NUM_POINTS 
||<a name="line.90"></a>
+<span class="sourceLineNo">091</span>            f[0].length &lt; 
MIN_NUM_POINTS) {<a name="line.91"></a>
+<span class="sourceLineNo">092</span>            throw new 
InsufficientDataException();<a name="line.92"></a>
+<span class="sourceLineNo">093</span>        }<a name="line.93"></a>
+<span class="sourceLineNo">094</span><a name="line.94"></a>
+<span class="sourceLineNo">095</span>        if (xLen != f.length) {<a 
name="line.95"></a>
+<span class="sourceLineNo">096</span>            throw new 
DimensionMismatchException(xLen, f.length);<a name="line.96"></a>
+<span class="sourceLineNo">097</span>        }<a name="line.97"></a>
+<span class="sourceLineNo">098</span><a name="line.98"></a>
+<span class="sourceLineNo">099</span>        if (yLen != f[0].length) {<a 
name="line.99"></a>
+<span class="sourceLineNo">100</span>            throw new 
DimensionMismatchException(yLen, f[0].length);<a name="line.100"></a>
+<span class="sourceLineNo">101</span>        }<a name="line.101"></a>
+<span class="sourceLineNo">102</span><a name="line.102"></a>
+<span class="sourceLineNo">103</span>        MathArrays.checkOrder(x);<a 
name="line.103"></a>
+<span class="sourceLineNo">104</span>        MathArrays.checkOrder(y);<a 
name="line.104"></a>
+<span class="sourceLineNo">105</span><a name="line.105"></a>
+<span class="sourceLineNo">106</span>        xval = x.clone();<a 
name="line.106"></a>
+<span class="sourceLineNo">107</span>        yval = y.clone();<a 
name="line.107"></a>
+<span class="sourceLineNo">108</span>        fval = f.clone();<a 
name="line.108"></a>
+<span class="sourceLineNo">109</span>    }<a name="line.109"></a>
+<span class="sourceLineNo">110</span><a name="line.110"></a>
+<span class="sourceLineNo">111</span>    /**<a name="line.111"></a>
+<span class="sourceLineNo">112</span>     * {@inheritDoc}<a 
name="line.112"></a>
+<span class="sourceLineNo">113</span>     */<a name="line.113"></a>
+<span class="sourceLineNo">114</span>    public double value(double x,<a 
name="line.114"></a>
+<span class="sourceLineNo">115</span>                        double y)<a 
name="line.115"></a>
+<span class="sourceLineNo">116</span>        throws OutOfRangeException {<a 
name="line.116"></a>
+<span class="sourceLineNo">117</span>        final AkimaSplineInterpolator 
interpolator = new AkimaSplineInterpolator();<a name="line.117"></a>
+<span class="sourceLineNo">118</span>        final int offset = 2;<a 
name="line.118"></a>
+<span class="sourceLineNo">119</span>        final int count = offset + 3;<a 
name="line.119"></a>
+<span class="sourceLineNo">120</span>        final int i = searchIndex(x, 
xval, offset, count);<a name="line.120"></a>
+<span class="sourceLineNo">121</span>        final int j = searchIndex(y, 
yval, offset, count);<a name="line.121"></a>
+<span class="sourceLineNo">122</span><a name="line.122"></a>
+<span class="sourceLineNo">123</span>        final double xArray[] = new 
double[count];<a name="line.123"></a>
+<span class="sourceLineNo">124</span>        final double yArray[] = new 
double[count];<a name="line.124"></a>
+<span class="sourceLineNo">125</span>        final double zArray[] = new 
double[count];<a name="line.125"></a>
+<span class="sourceLineNo">126</span>        final double interpArray[] = new 
double[count];<a name="line.126"></a>
+<span class="sourceLineNo">127</span><a name="line.127"></a>
+<span class="sourceLineNo">128</span>        for (int index = 0; index &lt; 
count; index++) {<a name="line.128"></a>
+<span class="sourceLineNo">129</span>            xArray[index] = xval[i + 
index];<a name="line.129"></a>
+<span class="sourceLineNo">130</span>            yArray[index] = yval[j + 
index];<a name="line.130"></a>
+<span class="sourceLineNo">131</span>        }<a name="line.131"></a>
+<span class="sourceLineNo">132</span><a name="line.132"></a>
+<span class="sourceLineNo">133</span>        for (int zIndex = 0; zIndex &lt; 
count; zIndex++) {<a name="line.133"></a>
+<span class="sourceLineNo">134</span>            for (int index = 0; index 
&lt; count; index++) {<a name="line.134"></a>
+<span class="sourceLineNo">135</span>                zArray[index] = fval[i + 
index][j + zIndex];<a name="line.135"></a>
+<span class="sourceLineNo">136</span>            }<a name="line.136"></a>
+<span class="sourceLineNo">137</span>            final 
PolynomialSplineFunction spline = interpolator.interpolate(xArray, zArray);<a 
name="line.137"></a>
+<span class="sourceLineNo">138</span>            interpArray[zIndex] = 
spline.value(x);<a name="line.138"></a>
+<span class="sourceLineNo">139</span>        }<a name="line.139"></a>
+<span class="sourceLineNo">140</span><a name="line.140"></a>
+<span class="sourceLineNo">141</span>        final PolynomialSplineFunction 
spline = interpolator.interpolate(yArray, interpArray);<a name="line.141"></a>
+<span class="sourceLineNo">142</span><a name="line.142"></a>
+<span class="sourceLineNo">143</span>        double returnValue = 
spline.value(y);<a name="line.143"></a>
+<span class="sourceLineNo">144</span><a name="line.144"></a>
+<span class="sourceLineNo">145</span>        return returnValue;<a 
name="line.145"></a>
+<span class="sourceLineNo">146</span>    }<a name="line.146"></a>
+<span class="sourceLineNo">147</span><a name="line.147"></a>
+<span class="sourceLineNo">148</span>    /**<a name="line.148"></a>
+<span class="sourceLineNo">149</span>     * Indicates whether a point is 
within the interpolation range.<a name="line.149"></a>
+<span class="sourceLineNo">150</span>     *<a name="line.150"></a>
+<span class="sourceLineNo">151</span>     * @param x First coordinate.<a 
name="line.151"></a>
+<span class="sourceLineNo">152</span>     * @param y Second coordinate.<a 
name="line.152"></a>
+<span class="sourceLineNo">153</span>     * @return {@code true} if (x, y) is 
a valid point.<a name="line.153"></a>
+<span class="sourceLineNo">154</span>     * @since 3.3<a name="line.154"></a>
+<span class="sourceLineNo">155</span>     */<a name="line.155"></a>
+<span class="sourceLineNo">156</span>    public boolean isValidPoint(double 
x,<a name="line.156"></a>
+<span class="sourceLineNo">157</span>                                double y) 
{<a name="line.157"></a>
+<span class="sourceLineNo">158</span>        if (x &lt; xval[0] ||<a 
name="line.158"></a>
+<span class="sourceLineNo">159</span>            x &gt; xval[xval.length - 1] 
||<a name="line.159"></a>
+<span class="sourceLineNo">160</span>            y &lt; yval[0] ||<a 
name="line.160"></a>
+<span class="sourceLineNo">161</span>            y &gt; yval[yval.length - 1]) 
{<a name="line.161"></a>
+<span class="sourceLineNo">162</span>            return false;<a 
name="line.162"></a>
+<span class="sourceLineNo">163</span>        } else {<a name="line.163"></a>
+<span class="sourceLineNo">164</span>            return true;<a 
name="line.164"></a>
+<span class="sourceLineNo">165</span>        }<a name="line.165"></a>
+<span class="sourceLineNo">166</span>    }<a name="line.166"></a>
+<span class="sourceLineNo">167</span><a name="line.167"></a>
+<span class="sourceLineNo">168</span>    /**<a name="line.168"></a>
+<span class="sourceLineNo">169</span>     * @param c Coordinate.<a 
name="line.169"></a>
+<span class="sourceLineNo">170</span>     * @param val Coordinate samples.<a 
name="line.170"></a>
+<span class="sourceLineNo">171</span>     * @param offset how far back from 
found value to offset for querying<a name="line.171"></a>
+<span class="sourceLineNo">172</span>     * @param count total number of 
elements forward from beginning that will be<a name="line.172"></a>
+<span class="sourceLineNo">173</span>     *        queried<a 
name="line.173"></a>
+<span class="sourceLineNo">174</span>     * @return the index in {@code val} 
corresponding to the interval containing<a name="line.174"></a>
+<span class="sourceLineNo">175</span>     *         {@code c}.<a 
name="line.175"></a>
+<span class="sourceLineNo">176</span>     * @throws OutOfRangeException if 
{@code c} is out of the range defined by<a name="line.176"></a>
+<span class="sourceLineNo">177</span>     *         the boundary values of 
{@code val}.<a name="line.177"></a>
+<span class="sourceLineNo">178</span>     */<a name="line.178"></a>
+<span class="sourceLineNo">179</span>    private int searchIndex(double c,<a 
name="line.179"></a>
+<span class="sourceLineNo">180</span>                            double[] 
val,<a name="line.180"></a>
+<span class="sourceLineNo">181</span>                            int offset,<a 
name="line.181"></a>
+<span class="sourceLineNo">182</span>                            int count) 
{<a name="line.182"></a>
+<span class="sourceLineNo">183</span>        int r = Arrays.binarySearch(val, 
c);<a name="line.183"></a>
+<span class="sourceLineNo">184</span><a name="line.184"></a>
+<span class="sourceLineNo">185</span>        if (r == -1 || r == -val.length - 
1) {<a name="line.185"></a>
+<span class="sourceLineNo">186</span>            throw new 
OutOfRangeException(c, val[0], val[val.length - 1]);<a name="line.186"></a>
+<span class="sourceLineNo">187</span>        }<a name="line.187"></a>
+<span class="sourceLineNo">188</span><a name="line.188"></a>
+<span class="sourceLineNo">189</span>        if (r &lt; 0) {<a 
name="line.189"></a>
+<span class="sourceLineNo">190</span>            // "c" in within an 
interpolation sub-interval, which returns<a name="line.190"></a>
+<span class="sourceLineNo">191</span>            // negative<a 
name="line.191"></a>
+<span class="sourceLineNo">192</span>            // need to remove the 
negative sign for consistency<a name="line.192"></a>
+<span class="sourceLineNo">193</span>            r = -r - offset - 1;<a 
name="line.193"></a>
+<span class="sourceLineNo">194</span>        } else {<a name="line.194"></a>
+<span class="sourceLineNo">195</span>            r -= offset;<a 
name="line.195"></a>
+<span class="sourceLineNo">196</span>        }<a name="line.196"></a>
+<span class="sourceLineNo">197</span><a name="line.197"></a>
+<span class="sourceLineNo">198</span>        if (r &lt; 0) {<a 
name="line.198"></a>
+<span class="sourceLineNo">199</span>            r = 0;<a name="line.199"></a>
+<span class="sourceLineNo">200</span>        }<a name="line.200"></a>
+<span class="sourceLineNo">201</span><a name="line.201"></a>
+<span class="sourceLineNo">202</span>        if ((r + count) &gt;= val.length) 
{<a name="line.202"></a>
+<span class="sourceLineNo">203</span>            // "c" is the last sample of 
the range: Return the index<a name="line.203"></a>
+<span class="sourceLineNo">204</span>            // of the sample at the lower 
end of the last sub-interval.<a name="line.204"></a>
+<span class="sourceLineNo">205</span>            r = val.length - count;<a 
name="line.205"></a>
+<span class="sourceLineNo">206</span>        }<a name="line.206"></a>
+<span class="sourceLineNo">207</span><a name="line.207"></a>
+<span class="sourceLineNo">208</span>        return r;<a name="line.208"></a>
+<span class="sourceLineNo">209</span>    }<a name="line.209"></a>
+<span class="sourceLineNo">210</span>}<a name="line.210"></a>
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+</pre>
+</div>
+</body>
+</html>

Propchange: 
websites/production/commons/content/proper/commons-math/apidocs/src-html/org/apache/commons/math3/analysis/interpolation/PiecewiseBicubicSplineInterpolatingFunction.html
------------------------------------------------------------------------------
    svn:eol-style = native

Added: 
websites/production/commons/content/proper/commons-math/apidocs/src-html/org/apache/commons/math3/analysis/interpolation/PiecewiseBicubicSplineInterpolator.html
==============================================================================
--- 
websites/production/commons/content/proper/commons-math/apidocs/src-html/org/apache/commons/math3/analysis/interpolation/PiecewiseBicubicSplineInterpolator.html
 (added)
+++ 
websites/production/commons/content/proper/commons-math/apidocs/src-html/org/apache/commons/math3/analysis/interpolation/PiecewiseBicubicSplineInterpolator.html
 Tue Dec 30 17:15:36 2014
@@ -0,0 +1,133 @@
+<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN" 
"http://www.w3.org/TR/html4/loose.dtd";>
+<html lang="en">
+<head>
+<title>Source code</title>
+<link rel="stylesheet" type="text/css" 
href="../../../../../../../stylesheet.css" title="Style">
+</head>
+<body>
+<div class="sourceContainer">
+<pre><span class="sourceLineNo">001</span>/*<a name="line.1"></a>
+<span class="sourceLineNo">002</span> * Licensed to the Apache Software 
Foundation (ASF) under one or more<a name="line.2"></a>
+<span class="sourceLineNo">003</span> * contributor license agreements.  See 
the NOTICE file distributed with<a name="line.3"></a>
+<span class="sourceLineNo">004</span> * this work for additional information 
regarding copyright ownership.<a name="line.4"></a>
+<span class="sourceLineNo">005</span> * The ASF licenses this file to You 
under the Apache License, Version 2.0<a name="line.5"></a>
+<span class="sourceLineNo">006</span> * (the "License"); you may not use this 
file except in compliance with<a name="line.6"></a>
+<span class="sourceLineNo">007</span> * the License.  You may obtain a copy of 
the License at<a name="line.7"></a>
+<span class="sourceLineNo">008</span> *<a name="line.8"></a>
+<span class="sourceLineNo">009</span> *      
http://www.apache.org/licenses/LICENSE-2.0<a name="line.9"></a>
+<span class="sourceLineNo">010</span> *<a name="line.10"></a>
+<span class="sourceLineNo">011</span> * Unless required by applicable law or 
agreed to in writing, software<a name="line.11"></a>
+<span class="sourceLineNo">012</span> * distributed under the License is 
distributed on an "AS IS" BASIS,<a name="line.12"></a>
+<span class="sourceLineNo">013</span> * WITHOUT WARRANTIES OR CONDITIONS OF 
ANY KIND, either express or implied.<a name="line.13"></a>
+<span class="sourceLineNo">014</span> * See the License for the specific 
language governing permissions and<a name="line.14"></a>
+<span class="sourceLineNo">015</span> * limitations under the License.<a 
name="line.15"></a>
+<span class="sourceLineNo">016</span> */<a name="line.16"></a>
+<span class="sourceLineNo">017</span>package 
org.apache.commons.math3.analysis.interpolation;<a name="line.17"></a>
+<span class="sourceLineNo">018</span><a name="line.18"></a>
+<span class="sourceLineNo">019</span>import 
org.apache.commons.math3.exception.DimensionMismatchException;<a 
name="line.19"></a>
+<span class="sourceLineNo">020</span>import 
org.apache.commons.math3.exception.NoDataException;<a name="line.20"></a>
+<span class="sourceLineNo">021</span>import 
org.apache.commons.math3.exception.NonMonotonicSequenceException;<a 
name="line.21"></a>
+<span class="sourceLineNo">022</span>import 
org.apache.commons.math3.exception.NullArgumentException;<a name="line.22"></a>
+<span class="sourceLineNo">023</span>import 
org.apache.commons.math3.util.MathArrays;<a name="line.23"></a>
+<span class="sourceLineNo">024</span><a name="line.24"></a>
+<span class="sourceLineNo">025</span>/**<a name="line.25"></a>
+<span class="sourceLineNo">026</span> * Generates a piecewise-bicubic 
interpolating function.<a name="line.26"></a>
+<span class="sourceLineNo">027</span> *<a name="line.27"></a>
+<span class="sourceLineNo">028</span> * @since 2.2<a name="line.28"></a>
+<span class="sourceLineNo">029</span> */<a name="line.29"></a>
+<span class="sourceLineNo">030</span>public class 
PiecewiseBicubicSplineInterpolator<a name="line.30"></a>
+<span class="sourceLineNo">031</span>    implements BivariateGridInterpolator 
{<a name="line.31"></a>
+<span class="sourceLineNo">032</span><a name="line.32"></a>
+<span class="sourceLineNo">033</span>    /**<a name="line.33"></a>
+<span class="sourceLineNo">034</span>     * {@inheritDoc}<a name="line.34"></a>
+<span class="sourceLineNo">035</span>     */<a name="line.35"></a>
+<span class="sourceLineNo">036</span>    public 
PiecewiseBicubicSplineInterpolatingFunction interpolate( final double[] xval,<a 
name="line.36"></a>
+<span class="sourceLineNo">037</span>                                          
                          final double[] yval,<a name="line.37"></a>
+<span class="sourceLineNo">038</span>                                          
                          final double[][] fval)<a name="line.38"></a>
+<span class="sourceLineNo">039</span>        throws 
DimensionMismatchException,<a name="line.39"></a>
+<span class="sourceLineNo">040</span>               NullArgumentException,<a 
name="line.40"></a>
+<span class="sourceLineNo">041</span>               NoDataException,<a 
name="line.41"></a>
+<span class="sourceLineNo">042</span>               
NonMonotonicSequenceException {<a name="line.42"></a>
+<span class="sourceLineNo">043</span>        if ( xval == null ||<a 
name="line.43"></a>
+<span class="sourceLineNo">044</span>             yval == null ||<a 
name="line.44"></a>
+<span class="sourceLineNo">045</span>             fval == null ||<a 
name="line.45"></a>
+<span class="sourceLineNo">046</span>             fval[0] == null ) {<a 
name="line.46"></a>
+<span class="sourceLineNo">047</span>            throw new 
NullArgumentException();<a name="line.47"></a>
+<span class="sourceLineNo">048</span>        }<a name="line.48"></a>
+<span class="sourceLineNo">049</span><a name="line.49"></a>
+<span class="sourceLineNo">050</span>        if ( xval.length == 0 ||<a 
name="line.50"></a>
+<span class="sourceLineNo">051</span>             yval.length == 0 ||<a 
name="line.51"></a>
+<span class="sourceLineNo">052</span>             fval.length == 0 ) {<a 
name="line.52"></a>
+<span class="sourceLineNo">053</span>            throw new 
NoDataException();<a name="line.53"></a>
+<span class="sourceLineNo">054</span>        }<a name="line.54"></a>
+<span class="sourceLineNo">055</span><a name="line.55"></a>
+<span class="sourceLineNo">056</span>        MathArrays.checkOrder(xval);<a 
name="line.56"></a>
+<span class="sourceLineNo">057</span>        MathArrays.checkOrder(yval);<a 
name="line.57"></a>
+<span class="sourceLineNo">058</span><a name="line.58"></a>
+<span class="sourceLineNo">059</span>        return new 
PiecewiseBicubicSplineInterpolatingFunction( xval, yval, fval );<a 
name="line.59"></a>
+<span class="sourceLineNo">060</span>    }<a name="line.60"></a>
+<span class="sourceLineNo">061</span>}<a name="line.61"></a>
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+</pre>
+</div>
+</body>
+</html>

Propchange: 
websites/production/commons/content/proper/commons-math/apidocs/src-html/org/apache/commons/math3/analysis/interpolation/PiecewiseBicubicSplineInterpolator.html
------------------------------------------------------------------------------
    svn:eol-style = native

Reply via email to