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"> </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"> </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 x,
+ double alpha,
+ int 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 © 2003–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 © 2003–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> * <a
href="http://en.wikipedia.org/wiki/Bicubic_interpolation"><a
name="line.29"></a>
+<span class="sourceLineNo">030</span> * bicubic spline
interpolation</a>.<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 < 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 < 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 < xval[0] ||<a
name="line.175"></a>
+<span class="sourceLineNo">176</span> x > xval[xval.length - 1]
||<a name="line.176"></a>
+<span class="sourceLineNo">177</span> y < yval[0] ||<a
name="line.177"></a>
+<span class="sourceLineNo">178</span> y > 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 < 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> * <ul><a name="line.221"></a>
+<span class="sourceLineNo">222</span> * <li>f(0,0)</li><a
name="line.222"></a>
+<span class="sourceLineNo">223</span> * <li>f(1,0)</li><a
name="line.223"></a>
+<span class="sourceLineNo">224</span> * <li>f(0,1)</li><a
name="line.224"></a>
+<span class="sourceLineNo">225</span> * <li>f(1,1)</li><a
name="line.225"></a>
+<span class="sourceLineNo">226</span> *
<li>f<sub>x</sub>(0,0)</li><a name="line.226"></a>
+<span class="sourceLineNo">227</span> *
<li>f<sub>x</sub>(1,0)</li><a name="line.227"></a>
+<span class="sourceLineNo">228</span> *
<li>f<sub>x</sub>(0,1)</li><a name="line.228"></a>
+<span class="sourceLineNo">229</span> *
<li>f<sub>x</sub>(1,1)</li><a name="line.229"></a>
+<span class="sourceLineNo">230</span> *
<li>f<sub>y</sub>(0,0)</li><a name="line.230"></a>
+<span class="sourceLineNo">231</span> *
<li>f<sub>y</sub>(1,0)</li><a name="line.231"></a>
+<span class="sourceLineNo">232</span> *
<li>f<sub>y</sub>(0,1)</li><a name="line.232"></a>
+<span class="sourceLineNo">233</span> *
<li>f<sub>y</sub>(1,1)</li><a name="line.233"></a>
+<span class="sourceLineNo">234</span> *
<li>f<sub>xy</sub>(0,0)</li><a name="line.234"></a>
+<span class="sourceLineNo">235</span> *
<li>f<sub>xy</sub>(1,0)</li><a name="line.235"></a>
+<span class="sourceLineNo">236</span> *
<li>f<sub>xy</sub>(0,1)</li><a name="line.236"></a>
+<span class="sourceLineNo">237</span> *
<li>f<sub>xy</sub>(1,1)</li><a name="line.237"></a>
+<span class="sourceLineNo">238</span> * </ul><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 < 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 <
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 < 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 < 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 < 0 || x > 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 < 0 || y > 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 < 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> * <p><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> * </p><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 < 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 < 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 < xval[1] ||<a
name="line.101"></a>
+<span class="sourceLineNo">102</span> x >
xval[xval.length - 2] ||<a name="line.102"></a>
+<span class="sourceLineNo">103</span> y < yval[1] ||<a
name="line.103"></a>
+<span class="sourceLineNo">104</span> y >
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> * <a
href="http://www.paulinternet.nl/?page=bicubic">bicubic spline</a><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 < MIN_NUM_POINTS ||<a
name="line.88"></a>
+<span class="sourceLineNo">089</span> yLen < MIN_NUM_POINTS ||<a
name="line.89"></a>
+<span class="sourceLineNo">090</span> f.length < MIN_NUM_POINTS
||<a name="line.90"></a>
+<span class="sourceLineNo">091</span> f[0].length <
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 <
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 <
count; zIndex++) {<a name="line.133"></a>
+<span class="sourceLineNo">134</span> for (int index = 0; index
< 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 < xval[0] ||<a
name="line.158"></a>
+<span class="sourceLineNo">159</span> x > xval[xval.length - 1]
||<a name="line.159"></a>
+<span class="sourceLineNo">160</span> y < yval[0] ||<a
name="line.160"></a>
+<span class="sourceLineNo">161</span> y > 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 < 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 < 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) >= 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