Added:
websites/production/commons/content/proper/commons-math/apidocs/src-html/org/apache/commons/math3/analysis/interpolation/TricubicInterpolatingFunction.html
==============================================================================
---
websites/production/commons/content/proper/commons-math/apidocs/src-html/org/apache/commons/math3/analysis/interpolation/TricubicInterpolatingFunction.html
(added)
+++
websites/production/commons/content/proper/commons-math/apidocs/src-html/org/apache/commons/math3/analysis/interpolation/TricubicInterpolatingFunction.html
Tue Dec 30 17:15:36 2014
@@ -0,0 +1,582 @@
+<!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.analysis.TrivariateFunction;<a name="line.19"></a>
+<span class="sourceLineNo">020</span>import
org.apache.commons.math3.exception.DimensionMismatchException;<a
name="line.20"></a>
+<span class="sourceLineNo">021</span>import
org.apache.commons.math3.exception.NoDataException;<a name="line.21"></a>
+<span class="sourceLineNo">022</span>import
org.apache.commons.math3.exception.OutOfRangeException;<a name="line.22"></a>
+<span class="sourceLineNo">023</span>import
org.apache.commons.math3.exception.NonMonotonicSequenceException;<a
name="line.23"></a>
+<span class="sourceLineNo">024</span>import
org.apache.commons.math3.util.MathArrays;<a name="line.24"></a>
+<span class="sourceLineNo">025</span><a name="line.25"></a>
+<span class="sourceLineNo">026</span>/**<a name="line.26"></a>
+<span class="sourceLineNo">027</span> * Function that implements the<a
name="line.27"></a>
+<span class="sourceLineNo">028</span> * <a
href="http://en.wikipedia.org/wiki/Tricubic_interpolation"><a
name="line.28"></a>
+<span class="sourceLineNo">029</span> * tricubic spline
interpolation</a>, as proposed in<a name="line.29"></a>
+<span class="sourceLineNo">030</span> * <quote><a name="line.30"></a>
+<span class="sourceLineNo">031</span> * Tricubic interpolation in three
dimensions<br/><a name="line.31"></a>
+<span class="sourceLineNo">032</span> * F. Lekien and J. Marsden<br/><a
name="line.32"></a>
+<span class="sourceLineNo">033</span> * <em>Int. J. Numer. Meth.
Eng</em> 2005; <b>63</b>:455-471<a name="line.33"></a>
+<span class="sourceLineNo">034</span> * </quote><a name="line.34"></a>
+<span class="sourceLineNo">035</span> *<a name="line.35"></a>
+<span class="sourceLineNo">036</span> * @since 2.2<a name="line.36"></a>
+<span class="sourceLineNo">037</span> * @deprecated in 3.4.<a
name="line.37"></a>
+<span class="sourceLineNo">038</span> */<a name="line.38"></a>
+<span class="sourceLineNo">039</span>@Deprecated<a name="line.39"></a>
+<span class="sourceLineNo">040</span>public class
TricubicInterpolatingFunction<a name="line.40"></a>
+<span class="sourceLineNo">041</span> implements TrivariateFunction {<a
name="line.41"></a>
+<span class="sourceLineNo">042</span> /**<a name="line.42"></a>
+<span class="sourceLineNo">043</span> * Matrix to compute the spline
coefficients from the function values<a name="line.43"></a>
+<span class="sourceLineNo">044</span> * and function derivatives values<a
name="line.44"></a>
+<span class="sourceLineNo">045</span> */<a name="line.45"></a>
+<span class="sourceLineNo">046</span> private static final double[][] AINV
= {<a name="line.46"></a>
+<span class="sourceLineNo">047</span> {
1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,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,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0
},<a name="line.48"></a>
+<span class="sourceLineNo">049</span> {
-3,3,0,0,0,0,0,0,-2,-1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0
},<a name="line.49"></a>
+<span class="sourceLineNo">050</span> {
2,-2,0,0,0,0,0,0,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0
},<a name="line.50"></a>
+<span class="sourceLineNo">051</span> {
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0
},<a name="line.51"></a>
+<span class="sourceLineNo">052</span> {
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0
},<a name="line.52"></a>
+<span class="sourceLineNo">053</span> {
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-3,3,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-2,-1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0
},<a name="line.53"></a>
+<span class="sourceLineNo">054</span> {
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,2,-2,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0
},<a name="line.54"></a>
+<span class="sourceLineNo">055</span> {
-3,0,3,0,0,0,0,0,0,0,0,0,0,0,0,0,-2,0,-1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0
},<a name="line.55"></a>
+<span class="sourceLineNo">056</span> {
0,0,0,0,0,0,0,0,-3,0,3,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-2,0,-1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0
},<a name="line.56"></a>
+<span class="sourceLineNo">057</span> {
9,-9,-9,9,0,0,0,0,6,3,-6,-3,0,0,0,0,6,-6,3,-3,0,0,0,0,0,0,0,0,0,0,0,0,4,2,2,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0
},<a name="line.57"></a>
+<span class="sourceLineNo">058</span> {
-6,6,6,-6,0,0,0,0,-3,-3,3,3,0,0,0,0,-4,4,-2,2,0,0,0,0,0,0,0,0,0,0,0,0,-2,-2,-1,-1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0
},<a name="line.58"></a>
+<span class="sourceLineNo">059</span> {
2,0,-2,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0
},<a name="line.59"></a>
+<span class="sourceLineNo">060</span> {
0,0,0,0,0,0,0,0,2,0,-2,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0
},<a name="line.60"></a>
+<span class="sourceLineNo">061</span> {
-6,6,6,-6,0,0,0,0,-4,-2,4,2,0,0,0,0,-3,3,-3,3,0,0,0,0,0,0,0,0,0,0,0,0,-2,-1,-2,-1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0
},<a name="line.61"></a>
+<span class="sourceLineNo">062</span> {
4,-4,-4,4,0,0,0,0,2,2,-2,-2,0,0,0,0,2,-2,2,-2,0,0,0,0,0,0,0,0,0,0,0,0,1,1,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0
},<a name="line.62"></a>
+<span class="sourceLineNo">063</span> {
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0
},<a name="line.63"></a>
+<span class="sourceLineNo">064</span> {
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0
},<a name="line.64"></a>
+<span class="sourceLineNo">065</span> {
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-3,3,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-2,-1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0
},<a name="line.65"></a>
+<span class="sourceLineNo">066</span> {
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,2,-2,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0
},<a name="line.66"></a>
+<span class="sourceLineNo">067</span> {
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0
},<a name="line.67"></a>
+<span class="sourceLineNo">068</span> {
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0
},<a name="line.68"></a>
+<span class="sourceLineNo">069</span> {
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-3,3,0,0,0,0,0,0,-2,-1,0,0,0,0,0,0
},<a name="line.69"></a>
+<span class="sourceLineNo">070</span> {
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,2,-2,0,0,0,0,0,0,1,1,0,0,0,0,0,0
},<a name="line.70"></a>
+<span class="sourceLineNo">071</span> {
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-3,0,3,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-2,0,-1,0,0,0,0,0,0,0,0,0,0,0,0,0
},<a name="line.71"></a>
+<span class="sourceLineNo">072</span> {
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-3,0,3,0,0,0,0,0,0,0,0,0,0,0,0,0,-2,0,-1,0,0,0,0,0
},<a name="line.72"></a>
+<span class="sourceLineNo">073</span> {
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,9,-9,-9,9,0,0,0,0,0,0,0,0,0,0,0,0,6,3,-6,-3,0,0,0,0,6,-6,3,-3,0,0,0,0,4,2,2,1,0,0,0,0
},<a name="line.73"></a>
+<span class="sourceLineNo">074</span> {
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-6,6,6,-6,0,0,0,0,0,0,0,0,0,0,0,0,-3,-3,3,3,0,0,0,0,-4,4,-2,2,0,0,0,0,-2,-2,-1,-1,0,0,0,0
},<a name="line.74"></a>
+<span class="sourceLineNo">075</span> {
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,2,0,-2,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0
},<a name="line.75"></a>
+<span class="sourceLineNo">076</span> {
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,2,0,-2,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,1,0,0,0,0,0
},<a name="line.76"></a>
+<span class="sourceLineNo">077</span> {
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-6,6,6,-6,0,0,0,0,0,0,0,0,0,0,0,0,-4,-2,4,2,0,0,0,0,-3,3,-3,3,0,0,0,0,-2,-1,-2,-1,0,0,0,0
},<a name="line.77"></a>
+<span class="sourceLineNo">078</span> {
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,4,-4,-4,4,0,0,0,0,0,0,0,0,0,0,0,0,2,2,-2,-2,0,0,0,0,2,-2,2,-2,0,0,0,0,1,1,1,1,0,0,0,0
},<a name="line.78"></a>
+<span class="sourceLineNo">079</span>
{-3,0,0,0,3,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-2,0,0,0,-1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0
},<a name="line.79"></a>
+<span class="sourceLineNo">080</span> {
0,0,0,0,0,0,0,0,-3,0,0,0,3,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-2,0,0,0,-1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0
},<a name="line.80"></a>
+<span class="sourceLineNo">081</span> {
9,-9,0,0,-9,9,0,0,6,3,0,0,-6,-3,0,0,0,0,0,0,0,0,0,0,6,-6,0,0,3,-3,0,0,0,0,0,0,0,0,0,0,4,2,0,0,2,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0
},<a name="line.81"></a>
+<span class="sourceLineNo">082</span> {
-6,6,0,0,6,-6,0,0,-3,-3,0,0,3,3,0,0,0,0,0,0,0,0,0,0,-4,4,0,0,-2,2,0,0,0,0,0,0,0,0,0,0,-2,-2,0,0,-1,-1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0
},<a name="line.82"></a>
+<span class="sourceLineNo">083</span> {
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-3,0,0,0,3,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-2,0,0,0,-1,0,0,0,0,0,0,0,0,0,0,0
},<a name="line.83"></a>
+<span class="sourceLineNo">084</span> {
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-3,0,0,0,3,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-2,0,0,0,-1,0,0,0
},<a name="line.84"></a>
+<span class="sourceLineNo">085</span> {
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,9,-9,0,0,-9,9,0,0,0,0,0,0,0,0,0,0,6,3,0,0,-6,-3,0,0,0,0,0,0,0,0,0,0,6,-6,0,0,3,-3,0,0,4,2,0,0,2,1,0,0
},<a name="line.85"></a>
+<span class="sourceLineNo">086</span> {
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-6,6,0,0,6,-6,0,0,0,0,0,0,0,0,0,0,-3,-3,0,0,3,3,0,0,0,0,0,0,0,0,0,0,-4,4,0,0,-2,2,0,0,-2,-2,0,0,-1,-1,0,0
},<a name="line.86"></a>
+<span class="sourceLineNo">087</span> {
9,0,-9,0,-9,0,9,0,0,0,0,0,0,0,0,0,6,0,3,0,-6,0,-3,0,6,0,-6,0,3,0,-3,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,4,0,2,0,2,0,1,0,0,0,0,0,0,0,0,0
},<a name="line.87"></a>
+<span class="sourceLineNo">088</span> {
0,0,0,0,0,0,0,0,9,0,-9,0,-9,0,9,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,6,0,3,0,-6,0,-3,0,6,0,-6,0,3,0,-3,0,0,0,0,0,0,0,0,0,4,0,2,0,2,0,1,0
},<a name="line.88"></a>
+<span class="sourceLineNo">089</span> {
-27,27,27,-27,27,-27,-27,27,-18,-9,18,9,18,9,-18,-9,-18,18,-9,9,18,-18,9,-9,-18,18,18,-18,-9,9,9,-9,-12,-6,-6,-3,12,6,6,3,-12,-6,12,6,-6,-3,6,3,-12,12,-6,6,-6,6,-3,3,-8,-4,-4,-2,-4,-2,-2,-1
},<a name="line.89"></a>
+<span class="sourceLineNo">090</span> {
18,-18,-18,18,-18,18,18,-18,9,9,-9,-9,-9,-9,9,9,12,-12,6,-6,-12,12,-6,6,12,-12,-12,12,6,-6,-6,6,6,6,3,3,-6,-6,-3,-3,6,6,-6,-6,3,3,-3,-3,8,-8,4,-4,4,-4,2,-2,4,4,2,2,2,2,1,1
},<a name="line.90"></a>
+<span class="sourceLineNo">091</span> {
-6,0,6,0,6,0,-6,0,0,0,0,0,0,0,0,0,-3,0,-3,0,3,0,3,0,-4,0,4,0,-2,0,2,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-2,0,-2,0,-1,0,-1,0,0,0,0,0,0,0,0,0
},<a name="line.91"></a>
+<span class="sourceLineNo">092</span> {
0,0,0,0,0,0,0,0,-6,0,6,0,6,0,-6,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-3,0,-3,0,3,0,3,0,-4,0,4,0,-2,0,2,0,0,0,0,0,0,0,0,0,-2,0,-2,0,-1,0,-1,0
},<a name="line.92"></a>
+<span class="sourceLineNo">093</span> {
18,-18,-18,18,-18,18,18,-18,12,6,-12,-6,-12,-6,12,6,9,-9,9,-9,-9,9,-9,9,12,-12,-12,12,6,-6,-6,6,6,3,6,3,-6,-3,-6,-3,8,4,-8,-4,4,2,-4,-2,6,-6,6,-6,3,-3,3,-3,4,2,4,2,2,1,2,1
},<a name="line.93"></a>
+<span class="sourceLineNo">094</span> {
-12,12,12,-12,12,-12,-12,12,-6,-6,6,6,6,6,-6,-6,-6,6,-6,6,6,-6,6,-6,-8,8,8,-8,-4,4,4,-4,-3,-3,-3,-3,3,3,3,3,-4,-4,4,4,-2,-2,2,2,-4,4,-4,4,-2,2,-2,2,-2,-2,-2,-2,-1,-1,-1,-1
},<a name="line.94"></a>
+<span class="sourceLineNo">095</span> {
2,0,0,0,-2,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0
},<a name="line.95"></a>
+<span class="sourceLineNo">096</span> {
0,0,0,0,0,0,0,0,2,0,0,0,-2,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0
},<a name="line.96"></a>
+<span class="sourceLineNo">097</span> {
-6,6,0,0,6,-6,0,0,-4,-2,0,0,4,2,0,0,0,0,0,0,0,0,0,0,-3,3,0,0,-3,3,0,0,0,0,0,0,0,0,0,0,-2,-1,0,0,-2,-1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0
},<a name="line.97"></a>
+<span class="sourceLineNo">098</span> {
4,-4,0,0,-4,4,0,0,2,2,0,0,-2,-2,0,0,0,0,0,0,0,0,0,0,2,-2,0,0,2,-2,0,0,0,0,0,0,0,0,0,0,1,1,0,0,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0
},<a name="line.98"></a>
+<span class="sourceLineNo">099</span> {
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,2,0,0,0,-2,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0
},<a name="line.99"></a>
+<span class="sourceLineNo">100</span> {
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,2,0,0,0,-2,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,1,0,0,0
},<a name="line.100"></a>
+<span class="sourceLineNo">101</span> {
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-6,6,0,0,6,-6,0,0,0,0,0,0,0,0,0,0,-4,-2,0,0,4,2,0,0,0,0,0,0,0,0,0,0,-3,3,0,0,-3,3,0,0,-2,-1,0,0,-2,-1,0,0
},<a name="line.101"></a>
+<span class="sourceLineNo">102</span> {
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,4,-4,0,0,-4,4,0,0,0,0,0,0,0,0,0,0,2,2,0,0,-2,-2,0,0,0,0,0,0,0,0,0,0,2,-2,0,0,2,-2,0,0,1,1,0,0,1,1,0,0
},<a name="line.102"></a>
+<span class="sourceLineNo">103</span> {
-6,0,6,0,6,0,-6,0,0,0,0,0,0,0,0,0,-4,0,-2,0,4,0,2,0,-3,0,3,0,-3,0,3,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-2,0,-1,0,-2,0,-1,0,0,0,0,0,0,0,0,0
},<a name="line.103"></a>
+<span class="sourceLineNo">104</span> {
0,0,0,0,0,0,0,0,-6,0,6,0,6,0,-6,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-4,0,-2,0,4,0,2,0,-3,0,3,0,-3,0,3,0,0,0,0,0,0,0,0,0,-2,0,-1,0,-2,0,-1,0
},<a name="line.104"></a>
+<span class="sourceLineNo">105</span> {
18,-18,-18,18,-18,18,18,-18,12,6,-12,-6,-12,-6,12,6,12,-12,6,-6,-12,12,-6,6,9,-9,-9,9,9,-9,-9,9,8,4,4,2,-8,-4,-4,-2,6,3,-6,-3,6,3,-6,-3,6,-6,3,-3,6,-6,3,-3,4,2,2,1,4,2,2,1
},<a name="line.105"></a>
+<span class="sourceLineNo">106</span> {
-12,12,12,-12,12,-12,-12,12,-6,-6,6,6,6,6,-6,-6,-8,8,-4,4,8,-8,4,-4,-6,6,6,-6,-6,6,6,-6,-4,-4,-2,-2,4,4,2,2,-3,-3,3,3,-3,-3,3,3,-4,4,-2,2,-4,4,-2,2,-2,-2,-1,-1,-2,-2,-1,-1
},<a name="line.106"></a>
+<span class="sourceLineNo">107</span> {
4,0,-4,0,-4,0,4,0,0,0,0,0,0,0,0,0,2,0,2,0,-2,0,-2,0,2,0,-2,0,2,0,-2,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,1,0,1,0,1,0,0,0,0,0,0,0,0,0
},<a name="line.107"></a>
+<span class="sourceLineNo">108</span> {
0,0,0,0,0,0,0,0,4,0,-4,0,-4,0,4,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,2,0,2,0,-2,0,-2,0,2,0,-2,0,2,0,-2,0,0,0,0,0,0,0,0,0,1,0,1,0,1,0,1,0
},<a name="line.108"></a>
+<span class="sourceLineNo">109</span> {
-12,12,12,-12,12,-12,-12,12,-8,-4,8,4,8,4,-8,-4,-6,6,-6,6,6,-6,6,-6,-6,6,6,-6,-6,6,6,-6,-4,-2,-4,-2,4,2,4,2,-4,-2,4,2,-4,-2,4,2,-3,3,-3,3,-3,3,-3,3,-2,-1,-2,-1,-2,-1,-2,-1
},<a name="line.109"></a>
+<span class="sourceLineNo">110</span> {
8,-8,-8,8,-8,8,8,-8,4,4,-4,-4,-4,-4,4,4,4,-4,4,-4,-4,4,-4,4,4,-4,-4,4,4,-4,-4,4,2,2,2,2,-2,-2,-2,-2,2,2,-2,-2,2,2,-2,-2,2,-2,2,-2,2,-2,2,-2,1,1,1,1,1,1,1,1
}<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>
+<span class="sourceLineNo">113</span> /** Samples x-coordinates */<a
name="line.113"></a>
+<span class="sourceLineNo">114</span> private final double[] xval;<a
name="line.114"></a>
+<span class="sourceLineNo">115</span> /** Samples y-coordinates */<a
name="line.115"></a>
+<span class="sourceLineNo">116</span> private final double[] yval;<a
name="line.116"></a>
+<span class="sourceLineNo">117</span> /** Samples z-coordinates */<a
name="line.117"></a>
+<span class="sourceLineNo">118</span> private final double[] zval;<a
name="line.118"></a>
+<span class="sourceLineNo">119</span> /** Set of cubic splines pacthing the
whole data grid */<a name="line.119"></a>
+<span class="sourceLineNo">120</span> private final TricubicFunction[][][]
splines;<a name="line.120"></a>
+<span class="sourceLineNo">121</span><a name="line.121"></a>
+<span class="sourceLineNo">122</span> /**<a name="line.122"></a>
+<span class="sourceLineNo">123</span> * @param x Sample values of the
x-coordinate, in increasing order.<a name="line.123"></a>
+<span class="sourceLineNo">124</span> * @param y Sample values of the
y-coordinate, in increasing order.<a name="line.124"></a>
+<span class="sourceLineNo">125</span> * @param z Sample values of the
y-coordinate, in increasing order.<a name="line.125"></a>
+<span class="sourceLineNo">126</span> * @param f Values of the function on
every grid point.<a name="line.126"></a>
+<span class="sourceLineNo">127</span> * @param dFdX Values of the partial
derivative of function with respect to x on every grid point.<a
name="line.127"></a>
+<span class="sourceLineNo">128</span> * @param dFdY Values of the partial
derivative of function with respect to y on every grid point.<a
name="line.128"></a>
+<span class="sourceLineNo">129</span> * @param dFdZ Values of the partial
derivative of function with respect to z on every grid point.<a
name="line.129"></a>
+<span class="sourceLineNo">130</span> * @param d2FdXdY Values of the cross
partial derivative of function on every grid point.<a name="line.130"></a>
+<span class="sourceLineNo">131</span> * @param d2FdXdZ Values of the cross
partial derivative of function on every grid point.<a name="line.131"></a>
+<span class="sourceLineNo">132</span> * @param d2FdYdZ Values of the cross
partial derivative of function on every grid point.<a name="line.132"></a>
+<span class="sourceLineNo">133</span> * @param d3FdXdYdZ Values of the
cross partial derivative of function on every grid point.<a name="line.133"></a>
+<span class="sourceLineNo">134</span> * @throws NoDataException if any of
the arrays has zero length.<a name="line.134"></a>
+<span class="sourceLineNo">135</span> * @throws DimensionMismatchException
if the various arrays do not contain the expected number of elements.<a
name="line.135"></a>
+<span class="sourceLineNo">136</span> * @throws
NonMonotonicSequenceException if {@code x}, {@code y} or {@code z} are not
strictly increasing.<a name="line.136"></a>
+<span class="sourceLineNo">137</span> */<a name="line.137"></a>
+<span class="sourceLineNo">138</span> public
TricubicInterpolatingFunction(double[] x,<a name="line.138"></a>
+<span class="sourceLineNo">139</span>
double[] y,<a name="line.139"></a>
+<span class="sourceLineNo">140</span>
double[] z,<a name="line.140"></a>
+<span class="sourceLineNo">141</span>
double[][][] f,<a name="line.141"></a>
+<span class="sourceLineNo">142</span>
double[][][] dFdX,<a name="line.142"></a>
+<span class="sourceLineNo">143</span>
double[][][] dFdY,<a name="line.143"></a>
+<span class="sourceLineNo">144</span>
double[][][] dFdZ,<a name="line.144"></a>
+<span class="sourceLineNo">145</span>
double[][][] d2FdXdY,<a name="line.145"></a>
+<span class="sourceLineNo">146</span>
double[][][] d2FdXdZ,<a name="line.146"></a>
+<span class="sourceLineNo">147</span>
double[][][] d2FdYdZ,<a name="line.147"></a>
+<span class="sourceLineNo">148</span>
double[][][] d3FdXdYdZ)<a name="line.148"></a>
+<span class="sourceLineNo">149</span> throws NoDataException,<a
name="line.149"></a>
+<span class="sourceLineNo">150</span>
DimensionMismatchException,<a name="line.150"></a>
+<span class="sourceLineNo">151</span>
NonMonotonicSequenceException {<a name="line.151"></a>
+<span class="sourceLineNo">152</span> final int xLen = x.length;<a
name="line.152"></a>
+<span class="sourceLineNo">153</span> final int yLen = y.length;<a
name="line.153"></a>
+<span class="sourceLineNo">154</span> final int zLen = z.length;<a
name="line.154"></a>
+<span class="sourceLineNo">155</span><a name="line.155"></a>
+<span class="sourceLineNo">156</span> if (xLen == 0 || yLen == 0 ||
z.length == 0 || f.length == 0 || f[0].length == 0) {<a name="line.156"></a>
+<span class="sourceLineNo">157</span> throw new
NoDataException();<a name="line.157"></a>
+<span class="sourceLineNo">158</span> }<a name="line.158"></a>
+<span class="sourceLineNo">159</span> if (xLen != f.length) {<a
name="line.159"></a>
+<span class="sourceLineNo">160</span> throw new
DimensionMismatchException(xLen, f.length);<a name="line.160"></a>
+<span class="sourceLineNo">161</span> }<a name="line.161"></a>
+<span class="sourceLineNo">162</span> if (xLen != dFdX.length) {<a
name="line.162"></a>
+<span class="sourceLineNo">163</span> throw new
DimensionMismatchException(xLen, dFdX.length);<a name="line.163"></a>
+<span class="sourceLineNo">164</span> }<a name="line.164"></a>
+<span class="sourceLineNo">165</span> if (xLen != dFdY.length) {<a
name="line.165"></a>
+<span class="sourceLineNo">166</span> throw new
DimensionMismatchException(xLen, dFdY.length);<a name="line.166"></a>
+<span class="sourceLineNo">167</span> }<a name="line.167"></a>
+<span class="sourceLineNo">168</span> if (xLen != dFdZ.length) {<a
name="line.168"></a>
+<span class="sourceLineNo">169</span> throw new
DimensionMismatchException(xLen, dFdZ.length);<a name="line.169"></a>
+<span class="sourceLineNo">170</span> }<a name="line.170"></a>
+<span class="sourceLineNo">171</span> if (xLen != d2FdXdY.length) {<a
name="line.171"></a>
+<span class="sourceLineNo">172</span> throw new
DimensionMismatchException(xLen, d2FdXdY.length);<a name="line.172"></a>
+<span class="sourceLineNo">173</span> }<a name="line.173"></a>
+<span class="sourceLineNo">174</span> if (xLen != d2FdXdZ.length) {<a
name="line.174"></a>
+<span class="sourceLineNo">175</span> throw new
DimensionMismatchException(xLen, d2FdXdZ.length);<a name="line.175"></a>
+<span class="sourceLineNo">176</span> }<a name="line.176"></a>
+<span class="sourceLineNo">177</span> if (xLen != d2FdYdZ.length) {<a
name="line.177"></a>
+<span class="sourceLineNo">178</span> throw new
DimensionMismatchException(xLen, d2FdYdZ.length);<a name="line.178"></a>
+<span class="sourceLineNo">179</span> }<a name="line.179"></a>
+<span class="sourceLineNo">180</span> if (xLen != d3FdXdYdZ.length) {<a
name="line.180"></a>
+<span class="sourceLineNo">181</span> throw new
DimensionMismatchException(xLen, d3FdXdYdZ.length);<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> MathArrays.checkOrder(x);<a
name="line.184"></a>
+<span class="sourceLineNo">185</span> MathArrays.checkOrder(y);<a
name="line.185"></a>
+<span class="sourceLineNo">186</span> MathArrays.checkOrder(z);<a
name="line.186"></a>
+<span class="sourceLineNo">187</span><a name="line.187"></a>
+<span class="sourceLineNo">188</span> xval = x.clone();<a
name="line.188"></a>
+<span class="sourceLineNo">189</span> yval = y.clone();<a
name="line.189"></a>
+<span class="sourceLineNo">190</span> zval = z.clone();<a
name="line.190"></a>
+<span class="sourceLineNo">191</span><a name="line.191"></a>
+<span class="sourceLineNo">192</span> final int lastI = xLen - 1;<a
name="line.192"></a>
+<span class="sourceLineNo">193</span> final int lastJ = yLen - 1;<a
name="line.193"></a>
+<span class="sourceLineNo">194</span> final int lastK = zLen - 1;<a
name="line.194"></a>
+<span class="sourceLineNo">195</span> splines = new
TricubicFunction[lastI][lastJ][lastK];<a name="line.195"></a>
+<span class="sourceLineNo">196</span><a name="line.196"></a>
+<span class="sourceLineNo">197</span> for (int i = 0; i < lastI;
i++) {<a name="line.197"></a>
+<span class="sourceLineNo">198</span> if (f[i].length != yLen) {<a
name="line.198"></a>
+<span class="sourceLineNo">199</span> throw new
DimensionMismatchException(f[i].length, yLen);<a name="line.199"></a>
+<span class="sourceLineNo">200</span> }<a name="line.200"></a>
+<span class="sourceLineNo">201</span> if (dFdX[i].length != yLen)
{<a name="line.201"></a>
+<span class="sourceLineNo">202</span> throw new
DimensionMismatchException(dFdX[i].length, yLen);<a name="line.202"></a>
+<span class="sourceLineNo">203</span> }<a name="line.203"></a>
+<span class="sourceLineNo">204</span> if (dFdY[i].length != yLen)
{<a name="line.204"></a>
+<span class="sourceLineNo">205</span> throw new
DimensionMismatchException(dFdY[i].length, yLen);<a name="line.205"></a>
+<span class="sourceLineNo">206</span> }<a name="line.206"></a>
+<span class="sourceLineNo">207</span> if (dFdZ[i].length != yLen)
{<a name="line.207"></a>
+<span class="sourceLineNo">208</span> throw new
DimensionMismatchException(dFdZ[i].length, yLen);<a name="line.208"></a>
+<span class="sourceLineNo">209</span> }<a name="line.209"></a>
+<span class="sourceLineNo">210</span> if (d2FdXdY[i].length !=
yLen) {<a name="line.210"></a>
+<span class="sourceLineNo">211</span> throw new
DimensionMismatchException(d2FdXdY[i].length, yLen);<a name="line.211"></a>
+<span class="sourceLineNo">212</span> }<a name="line.212"></a>
+<span class="sourceLineNo">213</span> if (d2FdXdZ[i].length !=
yLen) {<a name="line.213"></a>
+<span class="sourceLineNo">214</span> throw new
DimensionMismatchException(d2FdXdZ[i].length, yLen);<a name="line.214"></a>
+<span class="sourceLineNo">215</span> }<a name="line.215"></a>
+<span class="sourceLineNo">216</span> if (d2FdYdZ[i].length !=
yLen) {<a name="line.216"></a>
+<span class="sourceLineNo">217</span> throw new
DimensionMismatchException(d2FdYdZ[i].length, yLen);<a name="line.217"></a>
+<span class="sourceLineNo">218</span> }<a name="line.218"></a>
+<span class="sourceLineNo">219</span> if (d3FdXdYdZ[i].length !=
yLen) {<a name="line.219"></a>
+<span class="sourceLineNo">220</span> throw new
DimensionMismatchException(d3FdXdYdZ[i].length, yLen);<a name="line.220"></a>
+<span class="sourceLineNo">221</span> }<a name="line.221"></a>
+<span class="sourceLineNo">222</span><a name="line.222"></a>
+<span class="sourceLineNo">223</span> final int ip1 = i + 1;<a
name="line.223"></a>
+<span class="sourceLineNo">224</span> final double xR = xval[ip1] -
xval[i];<a name="line.224"></a>
+<span class="sourceLineNo">225</span> for (int j = 0; j < lastJ;
j++) {<a name="line.225"></a>
+<span class="sourceLineNo">226</span> if (f[i][j].length !=
zLen) {<a name="line.226"></a>
+<span class="sourceLineNo">227</span> throw new
DimensionMismatchException(f[i][j].length, zLen);<a name="line.227"></a>
+<span class="sourceLineNo">228</span> }<a name="line.228"></a>
+<span class="sourceLineNo">229</span> if (dFdX[i][j].length !=
zLen) {<a name="line.229"></a>
+<span class="sourceLineNo">230</span> throw new
DimensionMismatchException(dFdX[i][j].length, zLen);<a name="line.230"></a>
+<span class="sourceLineNo">231</span> }<a name="line.231"></a>
+<span class="sourceLineNo">232</span> if (dFdY[i][j].length !=
zLen) {<a name="line.232"></a>
+<span class="sourceLineNo">233</span> throw new
DimensionMismatchException(dFdY[i][j].length, zLen);<a name="line.233"></a>
+<span class="sourceLineNo">234</span> }<a name="line.234"></a>
+<span class="sourceLineNo">235</span> if (dFdZ[i][j].length !=
zLen) {<a name="line.235"></a>
+<span class="sourceLineNo">236</span> throw new
DimensionMismatchException(dFdZ[i][j].length, zLen);<a name="line.236"></a>
+<span class="sourceLineNo">237</span> }<a name="line.237"></a>
+<span class="sourceLineNo">238</span> if (d2FdXdY[i][j].length
!= zLen) {<a name="line.238"></a>
+<span class="sourceLineNo">239</span> throw new
DimensionMismatchException(d2FdXdY[i][j].length, zLen);<a name="line.239"></a>
+<span class="sourceLineNo">240</span> }<a name="line.240"></a>
+<span class="sourceLineNo">241</span> if (d2FdXdZ[i][j].length
!= zLen) {<a name="line.241"></a>
+<span class="sourceLineNo">242</span> throw new
DimensionMismatchException(d2FdXdZ[i][j].length, zLen);<a name="line.242"></a>
+<span class="sourceLineNo">243</span> }<a name="line.243"></a>
+<span class="sourceLineNo">244</span> if (d2FdYdZ[i][j].length
!= zLen) {<a name="line.244"></a>
+<span class="sourceLineNo">245</span> throw new
DimensionMismatchException(d2FdYdZ[i][j].length, zLen);<a name="line.245"></a>
+<span class="sourceLineNo">246</span> }<a name="line.246"></a>
+<span class="sourceLineNo">247</span> if
(d3FdXdYdZ[i][j].length != zLen) {<a name="line.247"></a>
+<span class="sourceLineNo">248</span> throw new
DimensionMismatchException(d3FdXdYdZ[i][j].length, zLen);<a name="line.248"></a>
+<span class="sourceLineNo">249</span> }<a name="line.249"></a>
+<span class="sourceLineNo">250</span><a name="line.250"></a>
+<span class="sourceLineNo">251</span> final int jp1 = j + 1;<a
name="line.251"></a>
+<span class="sourceLineNo">252</span> final double yR =
yval[jp1] - yval[j];<a name="line.252"></a>
+<span class="sourceLineNo">253</span> final double xRyR = xR *
yR;<a name="line.253"></a>
+<span class="sourceLineNo">254</span> for (int k = 0; k <
lastK; k++) {<a name="line.254"></a>
+<span class="sourceLineNo">255</span> final int kp1 = k +
1;<a name="line.255"></a>
+<span class="sourceLineNo">256</span> final double zR =
zval[kp1] - zval[k];<a name="line.256"></a>
+<span class="sourceLineNo">257</span> final double xRzR =
xR * zR;<a name="line.257"></a>
+<span class="sourceLineNo">258</span> final double yRzR =
yR * zR;<a name="line.258"></a>
+<span class="sourceLineNo">259</span> final double xRyRzR =
xR * yRzR;<a name="line.259"></a>
+<span class="sourceLineNo">260</span><a name="line.260"></a>
+<span class="sourceLineNo">261</span> final double[] beta =
new double[] {<a name="line.261"></a>
+<span class="sourceLineNo">262</span> f[i][j][k],
f[ip1][j][k],<a name="line.262"></a>
+<span class="sourceLineNo">263</span> f[i][jp1][k],
f[ip1][jp1][k],<a name="line.263"></a>
+<span class="sourceLineNo">264</span> f[i][j][kp1],
f[ip1][j][kp1],<a name="line.264"></a>
+<span class="sourceLineNo">265</span> f[i][jp1][kp1],
f[ip1][jp1][kp1],<a name="line.265"></a>
+<span class="sourceLineNo">266</span><a name="line.266"></a>
+<span class="sourceLineNo">267</span> dFdX[i][j][k] *
xR, dFdX[ip1][j][k] * xR,<a name="line.267"></a>
+<span class="sourceLineNo">268</span> dFdX[i][jp1][k] *
xR, dFdX[ip1][jp1][k] * xR,<a name="line.268"></a>
+<span class="sourceLineNo">269</span> dFdX[i][j][kp1] *
xR, dFdX[ip1][j][kp1] * xR,<a name="line.269"></a>
+<span class="sourceLineNo">270</span> dFdX[i][jp1][kp1]
* xR, dFdX[ip1][jp1][kp1] * xR,<a name="line.270"></a>
+<span class="sourceLineNo">271</span><a name="line.271"></a>
+<span class="sourceLineNo">272</span> dFdY[i][j][k] *
yR, dFdY[ip1][j][k] * yR,<a name="line.272"></a>
+<span class="sourceLineNo">273</span> dFdY[i][jp1][k] *
yR, dFdY[ip1][jp1][k] * yR,<a name="line.273"></a>
+<span class="sourceLineNo">274</span> dFdY[i][j][kp1] *
yR, dFdY[ip1][j][kp1] * yR,<a name="line.274"></a>
+<span class="sourceLineNo">275</span> dFdY[i][jp1][kp1]
* yR, dFdY[ip1][jp1][kp1] * yR,<a name="line.275"></a>
+<span class="sourceLineNo">276</span><a name="line.276"></a>
+<span class="sourceLineNo">277</span> dFdZ[i][j][k] *
zR, dFdZ[ip1][j][k] * zR,<a name="line.277"></a>
+<span class="sourceLineNo">278</span> dFdZ[i][jp1][k] *
zR, dFdZ[ip1][jp1][k] * zR,<a name="line.278"></a>
+<span class="sourceLineNo">279</span> dFdZ[i][j][kp1] *
zR, dFdZ[ip1][j][kp1] * zR,<a name="line.279"></a>
+<span class="sourceLineNo">280</span> dFdZ[i][jp1][kp1]
* zR, dFdZ[ip1][jp1][kp1] * zR,<a name="line.280"></a>
+<span class="sourceLineNo">281</span><a name="line.281"></a>
+<span class="sourceLineNo">282</span> d2FdXdY[i][j][k]
* xRyR, d2FdXdY[ip1][j][k] * xRyR,<a name="line.282"></a>
+<span class="sourceLineNo">283</span>
d2FdXdY[i][jp1][k] * xRyR, d2FdXdY[ip1][jp1][k] * xRyR,<a name="line.283"></a>
+<span class="sourceLineNo">284</span>
d2FdXdY[i][j][kp1] * xRyR, d2FdXdY[ip1][j][kp1] * xRyR,<a name="line.284"></a>
+<span class="sourceLineNo">285</span>
d2FdXdY[i][jp1][kp1] * xRyR, d2FdXdY[ip1][jp1][kp1] * xRyR,<a
name="line.285"></a>
+<span class="sourceLineNo">286</span><a name="line.286"></a>
+<span class="sourceLineNo">287</span> d2FdXdZ[i][j][k]
* xRzR, d2FdXdZ[ip1][j][k] * xRzR,<a name="line.287"></a>
+<span class="sourceLineNo">288</span>
d2FdXdZ[i][jp1][k] * xRzR, d2FdXdZ[ip1][jp1][k] * xRzR,<a name="line.288"></a>
+<span class="sourceLineNo">289</span>
d2FdXdZ[i][j][kp1] * xRzR, d2FdXdZ[ip1][j][kp1] * xRzR,<a name="line.289"></a>
+<span class="sourceLineNo">290</span>
d2FdXdZ[i][jp1][kp1] * xRzR, d2FdXdZ[ip1][jp1][kp1] * xRzR,<a
name="line.290"></a>
+<span class="sourceLineNo">291</span><a name="line.291"></a>
+<span class="sourceLineNo">292</span> d2FdYdZ[i][j][k]
* yRzR, d2FdYdZ[ip1][j][k] * yRzR,<a name="line.292"></a>
+<span class="sourceLineNo">293</span>
d2FdYdZ[i][jp1][k] * yRzR, d2FdYdZ[ip1][jp1][k] * yRzR,<a name="line.293"></a>
+<span class="sourceLineNo">294</span>
d2FdYdZ[i][j][kp1] * yRzR, d2FdYdZ[ip1][j][kp1] * yRzR,<a name="line.294"></a>
+<span class="sourceLineNo">295</span>
d2FdYdZ[i][jp1][kp1] * yRzR, d2FdYdZ[ip1][jp1][kp1] * yRzR,<a
name="line.295"></a>
+<span class="sourceLineNo">296</span><a name="line.296"></a>
+<span class="sourceLineNo">297</span>
d3FdXdYdZ[i][j][k] * xRyRzR, d3FdXdYdZ[ip1][j][k] * xRyRzR,<a
name="line.297"></a>
+<span class="sourceLineNo">298</span>
d3FdXdYdZ[i][jp1][k] * xRyRzR, d3FdXdYdZ[ip1][jp1][k] * xRyRzR,<a
name="line.298"></a>
+<span class="sourceLineNo">299</span>
d3FdXdYdZ[i][j][kp1] * xRyRzR, d3FdXdYdZ[ip1][j][kp1] * xRyRzR,<a
name="line.299"></a>
+<span class="sourceLineNo">300</span>
d3FdXdYdZ[i][jp1][kp1] * xRyRzR, d3FdXdYdZ[ip1][jp1][kp1] * xRyRzR,<a
name="line.300"></a>
+<span class="sourceLineNo">301</span> };<a
name="line.301"></a>
+<span class="sourceLineNo">302</span><a name="line.302"></a>
+<span class="sourceLineNo">303</span> splines[i][j][k] =
new TricubicFunction(computeCoefficients(beta));<a name="line.303"></a>
+<span class="sourceLineNo">304</span> }<a name="line.304"></a>
+<span class="sourceLineNo">305</span> }<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> /**<a name="line.309"></a>
+<span class="sourceLineNo">310</span> * {@inheritDoc}<a
name="line.310"></a>
+<span class="sourceLineNo">311</span> *<a name="line.311"></a>
+<span class="sourceLineNo">312</span> * @throws OutOfRangeException if any
of the variables is outside its interpolation range.<a name="line.312"></a>
+<span class="sourceLineNo">313</span> */<a name="line.313"></a>
+<span class="sourceLineNo">314</span> public double value(double x, double
y, double z)<a name="line.314"></a>
+<span class="sourceLineNo">315</span> throws OutOfRangeException {<a
name="line.315"></a>
+<span class="sourceLineNo">316</span> final int i = searchIndex(x,
xval);<a name="line.316"></a>
+<span class="sourceLineNo">317</span> if (i == -1) {<a
name="line.317"></a>
+<span class="sourceLineNo">318</span> throw new
OutOfRangeException(x, xval[0], xval[xval.length - 1]);<a name="line.318"></a>
+<span class="sourceLineNo">319</span> }<a name="line.319"></a>
+<span class="sourceLineNo">320</span> final int j = searchIndex(y,
yval);<a name="line.320"></a>
+<span class="sourceLineNo">321</span> if (j == -1) {<a
name="line.321"></a>
+<span class="sourceLineNo">322</span> throw new
OutOfRangeException(y, yval[0], yval[yval.length - 1]);<a name="line.322"></a>
+<span class="sourceLineNo">323</span> }<a name="line.323"></a>
+<span class="sourceLineNo">324</span> final int k = searchIndex(z,
zval);<a name="line.324"></a>
+<span class="sourceLineNo">325</span> if (k == -1) {<a
name="line.325"></a>
+<span class="sourceLineNo">326</span> throw new
OutOfRangeException(z, zval[0], zval[zval.length - 1]);<a name="line.326"></a>
+<span class="sourceLineNo">327</span> }<a name="line.327"></a>
+<span class="sourceLineNo">328</span><a name="line.328"></a>
+<span class="sourceLineNo">329</span> final double xN = (x - xval[i]) /
(xval[i + 1] - xval[i]);<a name="line.329"></a>
+<span class="sourceLineNo">330</span> final double yN = (y - yval[j]) /
(yval[j + 1] - yval[j]);<a name="line.330"></a>
+<span class="sourceLineNo">331</span> final double zN = (z - zval[k]) /
(zval[k + 1] - zval[k]);<a name="line.331"></a>
+<span class="sourceLineNo">332</span><a name="line.332"></a>
+<span class="sourceLineNo">333</span> return splines[i][j][k].value(xN,
yN, zN);<a name="line.333"></a>
+<span class="sourceLineNo">334</span> }<a name="line.334"></a>
+<span class="sourceLineNo">335</span><a name="line.335"></a>
+<span class="sourceLineNo">336</span> /**<a name="line.336"></a>
+<span class="sourceLineNo">337</span> * Indicates whether a point is
within the interpolation range.<a name="line.337"></a>
+<span class="sourceLineNo">338</span> *<a name="line.338"></a>
+<span class="sourceLineNo">339</span> * @param x First coordinate.<a
name="line.339"></a>
+<span class="sourceLineNo">340</span> * @param y Second coordinate.<a
name="line.340"></a>
+<span class="sourceLineNo">341</span> * @param z Third coordinate.<a
name="line.341"></a>
+<span class="sourceLineNo">342</span> * @return {@code true} if (x, y, z)
is a valid point.<a name="line.342"></a>
+<span class="sourceLineNo">343</span> */<a name="line.343"></a>
+<span class="sourceLineNo">344</span> public boolean isValidPoint(double x,
double y, double z) {<a name="line.344"></a>
+<span class="sourceLineNo">345</span> if (x < xval[0] ||<a
name="line.345"></a>
+<span class="sourceLineNo">346</span> x > xval[xval.length - 1]
||<a name="line.346"></a>
+<span class="sourceLineNo">347</span> y < yval[0] ||<a
name="line.347"></a>
+<span class="sourceLineNo">348</span> y > yval[yval.length - 1]
||<a name="line.348"></a>
+<span class="sourceLineNo">349</span> z < zval[0] ||<a
name="line.349"></a>
+<span class="sourceLineNo">350</span> z > zval[zval.length - 1])
{<a name="line.350"></a>
+<span class="sourceLineNo">351</span> return false;<a
name="line.351"></a>
+<span class="sourceLineNo">352</span> } else {<a name="line.352"></a>
+<span class="sourceLineNo">353</span> return true;<a
name="line.353"></a>
+<span class="sourceLineNo">354</span> }<a name="line.354"></a>
+<span class="sourceLineNo">355</span> }<a name="line.355"></a>
+<span class="sourceLineNo">356</span><a name="line.356"></a>
+<span class="sourceLineNo">357</span> /**<a name="line.357"></a>
+<span class="sourceLineNo">358</span> * @param c Coordinate.<a
name="line.358"></a>
+<span class="sourceLineNo">359</span> * @param val Coordinate samples.<a
name="line.359"></a>
+<span class="sourceLineNo">360</span> * @return the index in {@code val}
corresponding to the interval containing {@code c}, or {@code -1}<a
name="line.360"></a>
+<span class="sourceLineNo">361</span> * if {@code c} is out of the range
defined by the end values of {@code val}.<a name="line.361"></a>
+<span class="sourceLineNo">362</span> */<a name="line.362"></a>
+<span class="sourceLineNo">363</span> private int searchIndex(double c,
double[] val) {<a name="line.363"></a>
+<span class="sourceLineNo">364</span> if (c < val[0]) {<a
name="line.364"></a>
+<span class="sourceLineNo">365</span> return -1;<a
name="line.365"></a>
+<span class="sourceLineNo">366</span> }<a name="line.366"></a>
+<span class="sourceLineNo">367</span><a name="line.367"></a>
+<span class="sourceLineNo">368</span> final int max = val.length;<a
name="line.368"></a>
+<span class="sourceLineNo">369</span> for (int i = 1; i < max; i++)
{<a name="line.369"></a>
+<span class="sourceLineNo">370</span> if (c <= val[i]) {<a
name="line.370"></a>
+<span class="sourceLineNo">371</span> return i - 1;<a
name="line.371"></a>
+<span class="sourceLineNo">372</span> }<a name="line.372"></a>
+<span class="sourceLineNo">373</span> }<a name="line.373"></a>
+<span class="sourceLineNo">374</span><a name="line.374"></a>
+<span class="sourceLineNo">375</span> return -1;<a name="line.375"></a>
+<span class="sourceLineNo">376</span> }<a name="line.376"></a>
+<span class="sourceLineNo">377</span><a name="line.377"></a>
+<span class="sourceLineNo">378</span> /**<a name="line.378"></a>
+<span class="sourceLineNo">379</span> * Compute the spline coefficients
from the list of function values and<a name="line.379"></a>
+<span class="sourceLineNo">380</span> * function partial derivatives
values at the four corners of a grid<a name="line.380"></a>
+<span class="sourceLineNo">381</span> * element. They must be specified in
the following order:<a name="line.381"></a>
+<span class="sourceLineNo">382</span> * <ul><a name="line.382"></a>
+<span class="sourceLineNo">383</span> * <li>f(0,0,0)</li><a
name="line.383"></a>
+<span class="sourceLineNo">384</span> * <li>f(1,0,0)</li><a
name="line.384"></a>
+<span class="sourceLineNo">385</span> * <li>f(0,1,0)</li><a
name="line.385"></a>
+<span class="sourceLineNo">386</span> * <li>f(1,1,0)</li><a
name="line.386"></a>
+<span class="sourceLineNo">387</span> * <li>f(0,0,1)</li><a
name="line.387"></a>
+<span class="sourceLineNo">388</span> * <li>f(1,0,1)</li><a
name="line.388"></a>
+<span class="sourceLineNo">389</span> * <li>f(0,1,1)</li><a
name="line.389"></a>
+<span class="sourceLineNo">390</span> * <li>f(1,1,1)</li><a
name="line.390"></a>
+<span class="sourceLineNo">391</span> *<a name="line.391"></a>
+<span class="sourceLineNo">392</span> *
<li>f<sub>x</sub>(0,0,0)</li><a name="line.392"></a>
+<span class="sourceLineNo">393</span> * <li>... <em>(same
order as above)</em></li><a name="line.393"></a>
+<span class="sourceLineNo">394</span> *
<li>f<sub>x</sub>(1,1,1)</li><a name="line.394"></a>
+<span class="sourceLineNo">395</span> *<a name="line.395"></a>
+<span class="sourceLineNo">396</span> *
<li>f<sub>y</sub>(0,0,0)</li><a name="line.396"></a>
+<span class="sourceLineNo">397</span> * <li>... <em>(same
order as above)</em></li><a name="line.397"></a>
+<span class="sourceLineNo">398</span> *
<li>f<sub>y</sub>(1,1,1)</li><a name="line.398"></a>
+<span class="sourceLineNo">399</span> *<a name="line.399"></a>
+<span class="sourceLineNo">400</span> *
<li>f<sub>z</sub>(0,0,0)</li><a name="line.400"></a>
+<span class="sourceLineNo">401</span> * <li>... <em>(same
order as above)</em></li><a name="line.401"></a>
+<span class="sourceLineNo">402</span> *
<li>f<sub>z</sub>(1,1,1)</li><a name="line.402"></a>
+<span class="sourceLineNo">403</span> *<a name="line.403"></a>
+<span class="sourceLineNo">404</span> *
<li>f<sub>xy</sub>(0,0,0)</li><a name="line.404"></a>
+<span class="sourceLineNo">405</span> * <li>... <em>(same
order as above)</em></li><a name="line.405"></a>
+<span class="sourceLineNo">406</span> *
<li>f<sub>xy</sub>(1,1,1)</li><a name="line.406"></a>
+<span class="sourceLineNo">407</span> *<a name="line.407"></a>
+<span class="sourceLineNo">408</span> *
<li>f<sub>xz</sub>(0,0,0)</li><a name="line.408"></a>
+<span class="sourceLineNo">409</span> * <li>... <em>(same
order as above)</em></li><a name="line.409"></a>
+<span class="sourceLineNo">410</span> *
<li>f<sub>xz</sub>(1,1,1)</li><a name="line.410"></a>
+<span class="sourceLineNo">411</span> *<a name="line.411"></a>
+<span class="sourceLineNo">412</span> *
<li>f<sub>yz</sub>(0,0,0)</li><a name="line.412"></a>
+<span class="sourceLineNo">413</span> * <li>... <em>(same
order as above)</em></li><a name="line.413"></a>
+<span class="sourceLineNo">414</span> *
<li>f<sub>yz</sub>(1,1,1)</li><a name="line.414"></a>
+<span class="sourceLineNo">415</span> *<a name="line.415"></a>
+<span class="sourceLineNo">416</span> *
<li>f<sub>xyz</sub>(0,0,0)</li><a name="line.416"></a>
+<span class="sourceLineNo">417</span> * <li>... <em>(same
order as above)</em></li><a name="line.417"></a>
+<span class="sourceLineNo">418</span> *
<li>f<sub>xyz</sub>(1,1,1)</li><a name="line.418"></a>
+<span class="sourceLineNo">419</span> * </ul><a name="line.419"></a>
+<span class="sourceLineNo">420</span> * where the subscripts indicate the
partial derivative with respect to<a name="line.420"></a>
+<span class="sourceLineNo">421</span> * the corresponding variable(s).<a
name="line.421"></a>
+<span class="sourceLineNo">422</span> *<a name="line.422"></a>
+<span class="sourceLineNo">423</span> * @param beta List of function
values and function partial derivatives values.<a name="line.423"></a>
+<span class="sourceLineNo">424</span> * @return the spline coefficients.<a
name="line.424"></a>
+<span class="sourceLineNo">425</span> */<a name="line.425"></a>
+<span class="sourceLineNo">426</span> private double[]
computeCoefficients(double[] beta) {<a name="line.426"></a>
+<span class="sourceLineNo">427</span> final int sz = 64;<a
name="line.427"></a>
+<span class="sourceLineNo">428</span> final double[] a = new
double[sz];<a name="line.428"></a>
+<span class="sourceLineNo">429</span><a name="line.429"></a>
+<span class="sourceLineNo">430</span> for (int i = 0; i < sz; i++)
{<a name="line.430"></a>
+<span class="sourceLineNo">431</span> double result = 0;<a
name="line.431"></a>
+<span class="sourceLineNo">432</span> final double[] row =
AINV[i];<a name="line.432"></a>
+<span class="sourceLineNo">433</span> for (int j = 0; j < sz;
j++) {<a name="line.433"></a>
+<span class="sourceLineNo">434</span> result += row[j] *
beta[j];<a name="line.434"></a>
+<span class="sourceLineNo">435</span> }<a name="line.435"></a>
+<span class="sourceLineNo">436</span> a[i] = result;<a
name="line.436"></a>
+<span class="sourceLineNo">437</span> }<a name="line.437"></a>
+<span class="sourceLineNo">438</span><a name="line.438"></a>
+<span class="sourceLineNo">439</span> return a;<a name="line.439"></a>
+<span class="sourceLineNo">440</span> }<a name="line.440"></a>
+<span class="sourceLineNo">441</span>}<a name="line.441"></a>
+<span class="sourceLineNo">442</span><a name="line.442"></a>
+<span class="sourceLineNo">443</span>/**<a name="line.443"></a>
+<span class="sourceLineNo">444</span> * 3D-spline function.<a
name="line.444"></a>
+<span class="sourceLineNo">445</span> *<a name="line.445"></a>
+<span class="sourceLineNo">446</span> */<a name="line.446"></a>
+<span class="sourceLineNo">447</span>class TricubicFunction<a
name="line.447"></a>
+<span class="sourceLineNo">448</span> implements TrivariateFunction {<a
name="line.448"></a>
+<span class="sourceLineNo">449</span> /** Number of points. */<a
name="line.449"></a>
+<span class="sourceLineNo">450</span> private static final short N = 4;<a
name="line.450"></a>
+<span class="sourceLineNo">451</span> /** Coefficients */<a
name="line.451"></a>
+<span class="sourceLineNo">452</span> private final double[][][] a = new
double[N][N][N];<a name="line.452"></a>
+<span class="sourceLineNo">453</span><a name="line.453"></a>
+<span class="sourceLineNo">454</span> /**<a name="line.454"></a>
+<span class="sourceLineNo">455</span> * @param aV List of spline
coefficients.<a name="line.455"></a>
+<span class="sourceLineNo">456</span> */<a name="line.456"></a>
+<span class="sourceLineNo">457</span> public TricubicFunction(double[] aV)
{<a name="line.457"></a>
+<span class="sourceLineNo">458</span> for (int i = 0; i < N; i++)
{<a name="line.458"></a>
+<span class="sourceLineNo">459</span> for (int j = 0; j < N;
j++) {<a name="line.459"></a>
+<span class="sourceLineNo">460</span> for (int k = 0; k < N;
k++) {<a name="line.460"></a>
+<span class="sourceLineNo">461</span> a[i][j][k] = aV[i + N
* (j + N * k)];<a name="line.461"></a>
+<span class="sourceLineNo">462</span> }<a name="line.462"></a>
+<span class="sourceLineNo">463</span> }<a name="line.463"></a>
+<span class="sourceLineNo">464</span> }<a name="line.464"></a>
+<span class="sourceLineNo">465</span> }<a name="line.465"></a>
+<span class="sourceLineNo">466</span><a name="line.466"></a>
+<span class="sourceLineNo">467</span> /**<a name="line.467"></a>
+<span class="sourceLineNo">468</span> * @param x x-coordinate of the
interpolation point.<a name="line.468"></a>
+<span class="sourceLineNo">469</span> * @param y y-coordinate of the
interpolation point.<a name="line.469"></a>
+<span class="sourceLineNo">470</span> * @param z z-coordinate of the
interpolation point.<a name="line.470"></a>
+<span class="sourceLineNo">471</span> * @return the interpolated value.<a
name="line.471"></a>
+<span class="sourceLineNo">472</span> * @throws OutOfRangeException if
{@code x}, {@code y} or<a name="line.472"></a>
+<span class="sourceLineNo">473</span> * {@code z} are not in the interval
{@code [0, 1]}.<a name="line.473"></a>
+<span class="sourceLineNo">474</span> */<a name="line.474"></a>
+<span class="sourceLineNo">475</span> public double value(double x, double
y, double z)<a name="line.475"></a>
+<span class="sourceLineNo">476</span> throws OutOfRangeException {<a
name="line.476"></a>
+<span class="sourceLineNo">477</span> if (x < 0 || x > 1) {<a
name="line.477"></a>
+<span class="sourceLineNo">478</span> throw new
OutOfRangeException(x, 0, 1);<a name="line.478"></a>
+<span class="sourceLineNo">479</span> }<a name="line.479"></a>
+<span class="sourceLineNo">480</span> if (y < 0 || y > 1) {<a
name="line.480"></a>
+<span class="sourceLineNo">481</span> throw new
OutOfRangeException(y, 0, 1);<a name="line.481"></a>
+<span class="sourceLineNo">482</span> }<a name="line.482"></a>
+<span class="sourceLineNo">483</span> if (z < 0 || z > 1) {<a
name="line.483"></a>
+<span class="sourceLineNo">484</span> throw new
OutOfRangeException(z, 0, 1);<a name="line.484"></a>
+<span class="sourceLineNo">485</span> }<a name="line.485"></a>
+<span class="sourceLineNo">486</span><a name="line.486"></a>
+<span class="sourceLineNo">487</span> final double x2 = x * x;<a
name="line.487"></a>
+<span class="sourceLineNo">488</span> final double x3 = x2 * x;<a
name="line.488"></a>
+<span class="sourceLineNo">489</span> final double[] pX = { 1, x, x2,
x3 };<a name="line.489"></a>
+<span class="sourceLineNo">490</span><a name="line.490"></a>
+<span class="sourceLineNo">491</span> final double y2 = y * y;<a
name="line.491"></a>
+<span class="sourceLineNo">492</span> final double y3 = y2 * y;<a
name="line.492"></a>
+<span class="sourceLineNo">493</span> final double[] pY = { 1, y, y2,
y3 };<a name="line.493"></a>
+<span class="sourceLineNo">494</span><a name="line.494"></a>
+<span class="sourceLineNo">495</span> final double z2 = z * z;<a
name="line.495"></a>
+<span class="sourceLineNo">496</span> final double z3 = z2 * z;<a
name="line.496"></a>
+<span class="sourceLineNo">497</span> final double[] pZ = { 1, z, z2,
z3 };<a name="line.497"></a>
+<span class="sourceLineNo">498</span><a name="line.498"></a>
+<span class="sourceLineNo">499</span> double result = 0;<a
name="line.499"></a>
+<span class="sourceLineNo">500</span> for (int i = 0; i < N; i++)
{<a name="line.500"></a>
+<span class="sourceLineNo">501</span> for (int j = 0; j < N;
j++) {<a name="line.501"></a>
+<span class="sourceLineNo">502</span> for (int k = 0; k < N;
k++) {<a name="line.502"></a>
+<span class="sourceLineNo">503</span> result += a[i][j][k]
* pX[i] * pY[j] * pZ[k];<a name="line.503"></a>
+<span class="sourceLineNo">504</span> }<a name="line.504"></a>
+<span class="sourceLineNo">505</span> }<a name="line.505"></a>
+<span class="sourceLineNo">506</span> }<a name="line.506"></a>
+<span class="sourceLineNo">507</span><a name="line.507"></a>
+<span class="sourceLineNo">508</span> return result;<a
name="line.508"></a>
+<span class="sourceLineNo">509</span> }<a name="line.509"></a>
+<span class="sourceLineNo">510</span>}<a name="line.510"></a>
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+</pre>
+</div>
+</body>
+</html>
Propchange:
websites/production/commons/content/proper/commons-math/apidocs/src-html/org/apache/commons/math3/analysis/interpolation/TricubicInterpolatingFunction.html
------------------------------------------------------------------------------
svn:eol-style = native
Added:
websites/production/commons/content/proper/commons-math/apidocs/src-html/org/apache/commons/math3/analysis/interpolation/TricubicInterpolator.html
==============================================================================
---
websites/production/commons/content/proper/commons-math/apidocs/src-html/org/apache/commons/math3/analysis/interpolation/TricubicInterpolator.html
(added)
+++
websites/production/commons/content/proper/commons-math/apidocs/src-html/org/apache/commons/math3/analysis/interpolation/TricubicInterpolator.html
Tue Dec 30 17:15:36 2014
@@ -0,0 +1,214 @@
+<!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 tricubic interpolating
function.<a name="line.26"></a>
+<span class="sourceLineNo">027</span> *<a name="line.27"></a>
+<span class="sourceLineNo">028</span> * @since 3.4<a name="line.28"></a>
+<span class="sourceLineNo">029</span> */<a name="line.29"></a>
+<span class="sourceLineNo">030</span>public class TricubicInterpolator<a
name="line.30"></a>
+<span class="sourceLineNo">031</span> implements TrivariateGridInterpolator
{<a name="line.31"></a>
+<span class="sourceLineNo">032</span> /**<a name="line.32"></a>
+<span class="sourceLineNo">033</span> * {@inheritDoc}<a name="line.33"></a>
+<span class="sourceLineNo">034</span> */<a name="line.34"></a>
+<span class="sourceLineNo">035</span> public TricubicInterpolatingFunction
interpolate(final double[] xval,<a name="line.35"></a>
+<span class="sourceLineNo">036</span>
final double[] yval,<a name="line.36"></a>
+<span class="sourceLineNo">037</span>
final double[] zval,<a name="line.37"></a>
+<span class="sourceLineNo">038</span>
final double[][][] fval)<a name="line.38"></a>
+<span class="sourceLineNo">039</span> throws NoDataException,
NumberIsTooSmallException,<a name="line.39"></a>
+<span class="sourceLineNo">040</span>
DimensionMismatchException, NonMonotonicSequenceException {<a
name="line.40"></a>
+<span class="sourceLineNo">041</span> if (xval.length == 0 ||
yval.length == 0 || zval.length == 0 || fval.length == 0) {<a
name="line.41"></a>
+<span class="sourceLineNo">042</span> throw new
NoDataException();<a name="line.42"></a>
+<span class="sourceLineNo">043</span> }<a name="line.43"></a>
+<span class="sourceLineNo">044</span> if (xval.length != fval.length)
{<a name="line.44"></a>
+<span class="sourceLineNo">045</span> throw new
DimensionMismatchException(xval.length, fval.length);<a name="line.45"></a>
+<span class="sourceLineNo">046</span> }<a name="line.46"></a>
+<span class="sourceLineNo">047</span><a name="line.47"></a>
+<span class="sourceLineNo">048</span> MathArrays.checkOrder(xval);<a
name="line.48"></a>
+<span class="sourceLineNo">049</span> MathArrays.checkOrder(yval);<a
name="line.49"></a>
+<span class="sourceLineNo">050</span> MathArrays.checkOrder(zval);<a
name="line.50"></a>
+<span class="sourceLineNo">051</span><a name="line.51"></a>
+<span class="sourceLineNo">052</span> final int xLen = xval.length;<a
name="line.52"></a>
+<span class="sourceLineNo">053</span> final int yLen = yval.length;<a
name="line.53"></a>
+<span class="sourceLineNo">054</span> final int zLen = zval.length;<a
name="line.54"></a>
+<span class="sourceLineNo">055</span><a name="line.55"></a>
+<span class="sourceLineNo">056</span> // Approximation to the partial
derivatives using finite differences.<a name="line.56"></a>
+<span class="sourceLineNo">057</span> final double[][][] dFdX = new
double[xLen][yLen][zLen];<a name="line.57"></a>
+<span class="sourceLineNo">058</span> final double[][][] dFdY = new
double[xLen][yLen][zLen];<a name="line.58"></a>
+<span class="sourceLineNo">059</span> final double[][][] dFdZ = new
double[xLen][yLen][zLen];<a name="line.59"></a>
+<span class="sourceLineNo">060</span> final double[][][] d2FdXdY = new
double[xLen][yLen][zLen];<a name="line.60"></a>
+<span class="sourceLineNo">061</span> final double[][][] d2FdXdZ = new
double[xLen][yLen][zLen];<a name="line.61"></a>
+<span class="sourceLineNo">062</span> final double[][][] d2FdYdZ = new
double[xLen][yLen][zLen];<a name="line.62"></a>
+<span class="sourceLineNo">063</span> final double[][][] d3FdXdYdZ =
new double[xLen][yLen][zLen];<a name="line.63"></a>
+<span class="sourceLineNo">064</span><a name="line.64"></a>
+<span class="sourceLineNo">065</span> for (int i = 1; i < xLen - 1;
i++) {<a name="line.65"></a>
+<span class="sourceLineNo">066</span> if (yval.length !=
fval[i].length) {<a name="line.66"></a>
+<span class="sourceLineNo">067</span> throw new
DimensionMismatchException(yval.length, fval[i].length);<a name="line.67"></a>
+<span class="sourceLineNo">068</span> }<a name="line.68"></a>
+<span class="sourceLineNo">069</span><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> if (zval.length !=
fval[i][j].length) {<a name="line.79"></a>
+<span class="sourceLineNo">080</span> throw new
DimensionMismatchException(zval.length, fval[i][j].length);<a
name="line.80"></a>
+<span class="sourceLineNo">081</span> }<a name="line.81"></a>
+<span class="sourceLineNo">082</span><a name="line.82"></a>
+<span class="sourceLineNo">083</span> final int nJ = j + 1;<a
name="line.83"></a>
+<span class="sourceLineNo">084</span> final int pJ = j - 1;<a
name="line.84"></a>
+<span class="sourceLineNo">085</span><a name="line.85"></a>
+<span class="sourceLineNo">086</span> final double nY =
yval[nJ];<a name="line.86"></a>
+<span class="sourceLineNo">087</span> final double pY =
yval[pJ];<a name="line.87"></a>
+<span class="sourceLineNo">088</span><a name="line.88"></a>
+<span class="sourceLineNo">089</span> final double deltaY = nY
- pY;<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> for (int k = 1; k <
zLen - 1; k++) {<a name="line.92"></a>
+<span class="sourceLineNo">093</span> final int nK = k +
1;<a name="line.93"></a>
+<span class="sourceLineNo">094</span> final int pK = k -
1;<a name="line.94"></a>
+<span class="sourceLineNo">095</span><a name="line.95"></a>
+<span class="sourceLineNo">096</span> final double nZ =
zval[nK];<a name="line.96"></a>
+<span class="sourceLineNo">097</span> final double pZ =
zval[pK];<a name="line.97"></a>
+<span class="sourceLineNo">098</span><a name="line.98"></a>
+<span class="sourceLineNo">099</span> final double deltaZ =
nZ - pZ;<a name="line.99"></a>
+<span class="sourceLineNo">100</span><a name="line.100"></a>
+<span class="sourceLineNo">101</span> dFdX[i][j][k] =
(fval[nI][j][k] - fval[pI][j][k]) / deltaX;<a name="line.101"></a>
+<span class="sourceLineNo">102</span> dFdY[i][j][k] =
(fval[i][nJ][k] - fval[i][pJ][k]) / deltaY;<a name="line.102"></a>
+<span class="sourceLineNo">103</span> dFdZ[i][j][k] =
(fval[i][j][nK] - fval[i][j][pK]) / deltaZ;<a name="line.103"></a>
+<span class="sourceLineNo">104</span><a name="line.104"></a>
+<span class="sourceLineNo">105</span> final double deltaXZ
= deltaX * deltaZ;<a name="line.105"></a>
+<span class="sourceLineNo">106</span> final double deltaYZ
= deltaY * deltaZ;<a name="line.106"></a>
+<span class="sourceLineNo">107</span><a name="line.107"></a>
+<span class="sourceLineNo">108</span> d2FdXdY[i][j][k] =
(fval[nI][nJ][k] - fval[nI][pJ][k] - fval[pI][nJ][k] + fval[pI][pJ][k]) /
deltaXY;<a name="line.108"></a>
+<span class="sourceLineNo">109</span> d2FdXdZ[i][j][k] =
(fval[nI][j][nK] - fval[nI][j][pK] - fval[pI][j][nK] + fval[pI][j][pK]) /
deltaXZ;<a name="line.109"></a>
+<span class="sourceLineNo">110</span> d2FdYdZ[i][j][k] =
(fval[i][nJ][nK] - fval[i][nJ][pK] - fval[i][pJ][nK] + fval[i][pJ][pK]) /
deltaYZ;<a name="line.110"></a>
+<span class="sourceLineNo">111</span><a name="line.111"></a>
+<span class="sourceLineNo">112</span> final double deltaXYZ
= deltaXY * deltaZ;<a name="line.112"></a>
+<span class="sourceLineNo">113</span><a name="line.113"></a>
+<span class="sourceLineNo">114</span> d3FdXdYdZ[i][j][k] =
(fval[nI][nJ][nK] - fval[nI][pJ][nK] -<a name="line.114"></a>
+<span class="sourceLineNo">115</span>
fval[pI][nJ][nK] + fval[pI][pJ][nK] -<a name="line.115"></a>
+<span class="sourceLineNo">116</span>
fval[nI][nJ][pK] + fval[nI][pJ][pK] +<a name="line.116"></a>
+<span class="sourceLineNo">117</span>
fval[pI][nJ][pK] - fval[pI][pJ][pK]) / deltaXYZ;<a name="line.117"></a>
+<span class="sourceLineNo">118</span> }<a name="line.118"></a>
+<span class="sourceLineNo">119</span> }<a name="line.119"></a>
+<span class="sourceLineNo">120</span> }<a name="line.120"></a>
+<span class="sourceLineNo">121</span><a name="line.121"></a>
+<span class="sourceLineNo">122</span> // Create the interpolating
function.<a name="line.122"></a>
+<span class="sourceLineNo">123</span> return new
TricubicInterpolatingFunction(xval, yval, zval, fval,<a name="line.123"></a>
+<span class="sourceLineNo">124</span>
dFdX, dFdY, dFdZ,<a name="line.124"></a>
+<span class="sourceLineNo">125</span>
d2FdXdY, d2FdXdZ, d2FdYdZ,<a name="line.125"></a>
+<span class="sourceLineNo">126</span>
d3FdXdYdZ) {<a name="line.126"></a>
+<span class="sourceLineNo">127</span> @Override<a
name="line.127"></a>
+<span class="sourceLineNo">128</span> public boolean
isValidPoint(double x, double y, double z) {<a name="line.128"></a>
+<span class="sourceLineNo">129</span> if (x < xval[1] ||<a
name="line.129"></a>
+<span class="sourceLineNo">130</span> x >
xval[xval.length - 2] ||<a name="line.130"></a>
+<span class="sourceLineNo">131</span> y < yval[1] ||<a
name="line.131"></a>
+<span class="sourceLineNo">132</span> y >
yval[yval.length - 2] ||<a name="line.132"></a>
+<span class="sourceLineNo">133</span> z < zval[1] ||<a
name="line.133"></a>
+<span class="sourceLineNo">134</span> z >
zval[zval.length - 2]) {<a name="line.134"></a>
+<span class="sourceLineNo">135</span> return false;<a
name="line.135"></a>
+<span class="sourceLineNo">136</span> } else {<a
name="line.136"></a>
+<span class="sourceLineNo">137</span> return true;<a
name="line.137"></a>
+<span class="sourceLineNo">138</span> }<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> }<a name="line.141"></a>
+<span class="sourceLineNo">142</span>}<a name="line.142"></a>
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+</pre>
+</div>
+</body>
+</html>
Propchange:
websites/production/commons/content/proper/commons-math/apidocs/src-html/org/apache/commons/math3/analysis/interpolation/TricubicInterpolator.html
------------------------------------------------------------------------------
svn:eol-style = native
Added:
websites/production/commons/content/proper/commons-math/apidocs/src-html/org/apache/commons/math3/fitting/BasicCurveFitter.html
==============================================================================
---
websites/production/commons/content/proper/commons-math/apidocs/src-html/org/apache/commons/math3/fitting/BasicCurveFitter.html
(added)
+++
websites/production/commons/content/proper/commons-math/apidocs/src-html/org/apache/commons/math3/fitting/BasicCurveFitter.html
Tue Dec 30 17:15:36 2014
@@ -0,0 +1,198 @@
+<!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.fitting;<a name="line.17"></a>
+<span class="sourceLineNo">018</span><a name="line.18"></a>
+<span class="sourceLineNo">019</span>import java.util.Collection;<a
name="line.19"></a>
+<span class="sourceLineNo">020</span><a name="line.20"></a>
+<span class="sourceLineNo">021</span>import
org.apache.commons.math3.analysis.ParametricUnivariateFunction;<a
name="line.21"></a>
+<span class="sourceLineNo">022</span>import
org.apache.commons.math3.fitting.leastsquares.LeastSquaresBuilder;<a
name="line.22"></a>
+<span class="sourceLineNo">023</span>import
org.apache.commons.math3.fitting.leastsquares.LeastSquaresProblem;<a
name="line.23"></a>
+<span class="sourceLineNo">024</span>import
org.apache.commons.math3.linear.DiagonalMatrix;<a name="line.24"></a>
+<span class="sourceLineNo">025</span><a name="line.25"></a>
+<span class="sourceLineNo">026</span>/**<a name="line.26"></a>
+<span class="sourceLineNo">027</span> * Fits points to a user-defined {@link
ParametricUnivariateFunction function}.<a name="line.27"></a>
+<span class="sourceLineNo">028</span> *<a name="line.28"></a>
+<span class="sourceLineNo">029</span> * @since 3.4<a name="line.29"></a>
+<span class="sourceLineNo">030</span> */<a name="line.30"></a>
+<span class="sourceLineNo">031</span>public class BasicCurveFitter extends
AbstractCurveFitter {<a name="line.31"></a>
+<span class="sourceLineNo">032</span> /** Function to fit. */<a
name="line.32"></a>
+<span class="sourceLineNo">033</span> private final
ParametricUnivariateFunction function;<a name="line.33"></a>
+<span class="sourceLineNo">034</span> /** Initial guess for the parameters.
*/<a name="line.34"></a>
+<span class="sourceLineNo">035</span> private final double[]
initialGuess;<a name="line.35"></a>
+<span class="sourceLineNo">036</span> /** Maximum number of iterations of
the optimization algorithm. */<a name="line.36"></a>
+<span class="sourceLineNo">037</span> private final int maxIter;<a
name="line.37"></a>
+<span class="sourceLineNo">038</span><a name="line.38"></a>
+<span class="sourceLineNo">039</span> /**<a name="line.39"></a>
+<span class="sourceLineNo">040</span> * Contructor used by the factory
methods.<a name="line.40"></a>
+<span class="sourceLineNo">041</span> *<a name="line.41"></a>
+<span class="sourceLineNo">042</span> * @param function Function to fit.<a
name="line.42"></a>
+<span class="sourceLineNo">043</span> * @param initialGuess Initial guess.
Cannot be {@code null}. Its length must<a name="line.43"></a>
+<span class="sourceLineNo">044</span> * be consistent with the number of
parameters of the {@code function} to fit.<a name="line.44"></a>
+<span class="sourceLineNo">045</span> * @param maxIter Maximum number of
iterations of the optimization algorithm.<a name="line.45"></a>
+<span class="sourceLineNo">046</span> */<a name="line.46"></a>
+<span class="sourceLineNo">047</span> private
BasicCurveFitter(ParametricUnivariateFunction function,<a name="line.47"></a>
+<span class="sourceLineNo">048</span> double[]
initialGuess,<a name="line.48"></a>
+<span class="sourceLineNo">049</span> int maxIter)
{<a name="line.49"></a>
+<span class="sourceLineNo">050</span> this.function = function;<a
name="line.50"></a>
+<span class="sourceLineNo">051</span> this.initialGuess =
initialGuess;<a name="line.51"></a>
+<span class="sourceLineNo">052</span> this.maxIter = maxIter;<a
name="line.52"></a>
+<span class="sourceLineNo">053</span> }<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> * Creates a curve fitter.<a
name="line.56"></a>
+<span class="sourceLineNo">057</span> * The initial guess for the
parameters will be {@link ParameterGuesser}<a name="line.57"></a>
+<span class="sourceLineNo">058</span> * computed automatically, and the
maximum number of iterations of the<a name="line.58"></a>
+<span class="sourceLineNo">059</span> * optimization algorithm is set to
{@link Integer#MAX_VALUE}.<a name="line.59"></a>
+<span class="sourceLineNo">060</span> *<a name="line.60"></a>
+<span class="sourceLineNo">061</span> * @param f Function to fit.<a
name="line.61"></a>
+<span class="sourceLineNo">062</span> * @param start Initial guess for the
parameters. Cannot be {@code null}.<a name="line.62"></a>
+<span class="sourceLineNo">063</span> * Its length must be consistent with
the number of parameters of the<a name="line.63"></a>
+<span class="sourceLineNo">064</span> * function to fit.<a
name="line.64"></a>
+<span class="sourceLineNo">065</span> * @return a curve fitter.<a
name="line.65"></a>
+<span class="sourceLineNo">066</span> *<a name="line.66"></a>
+<span class="sourceLineNo">067</span> * @see #withStartPoint(double[])<a
name="line.67"></a>
+<span class="sourceLineNo">068</span> * @see #withMaxIterations(int)<a
name="line.68"></a>
+<span class="sourceLineNo">069</span> */<a name="line.69"></a>
+<span class="sourceLineNo">070</span> public static BasicCurveFitter
create(ParametricUnivariateFunction f,<a name="line.70"></a>
+<span class="sourceLineNo">071</span>
double[] start) {<a name="line.71"></a>
+<span class="sourceLineNo">072</span> return new BasicCurveFitter(f,
start, Integer.MAX_VALUE);<a name="line.72"></a>
+<span class="sourceLineNo">073</span> }<a name="line.73"></a>
+<span class="sourceLineNo">074</span><a name="line.74"></a>
+<span class="sourceLineNo">075</span> /**<a name="line.75"></a>
+<span class="sourceLineNo">076</span> * Configure the start point (initial
guess).<a name="line.76"></a>
+<span class="sourceLineNo">077</span> * @param newStart new start point
(initial guess)<a name="line.77"></a>
+<span class="sourceLineNo">078</span> * @return a new instance.<a
name="line.78"></a>
+<span class="sourceLineNo">079</span> */<a name="line.79"></a>
+<span class="sourceLineNo">080</span> public BasicCurveFitter
withStartPoint(double[] newStart) {<a name="line.80"></a>
+<span class="sourceLineNo">081</span> return new
BasicCurveFitter(function,<a name="line.81"></a>
+<span class="sourceLineNo">082</span>
newStart.clone(),<a name="line.82"></a>
+<span class="sourceLineNo">083</span>
maxIter);<a name="line.83"></a>
+<span class="sourceLineNo">084</span> }<a name="line.84"></a>
+<span class="sourceLineNo">085</span><a name="line.85"></a>
+<span class="sourceLineNo">086</span> /**<a name="line.86"></a>
+<span class="sourceLineNo">087</span> * Configure the maximum number of
iterations.<a name="line.87"></a>
+<span class="sourceLineNo">088</span> * @param newMaxIter maximum number
of iterations<a name="line.88"></a>
+<span class="sourceLineNo">089</span> * @return a new instance.<a
name="line.89"></a>
+<span class="sourceLineNo">090</span> */<a name="line.90"></a>
+<span class="sourceLineNo">091</span> public BasicCurveFitter
withMaxIterations(int newMaxIter) {<a name="line.91"></a>
+<span class="sourceLineNo">092</span> return new
BasicCurveFitter(function,<a name="line.92"></a>
+<span class="sourceLineNo">093</span>
initialGuess,<a name="line.93"></a>
+<span class="sourceLineNo">094</span>
newMaxIter);<a name="line.94"></a>
+<span class="sourceLineNo">095</span> }<a name="line.95"></a>
+<span class="sourceLineNo">096</span><a name="line.96"></a>
+<span class="sourceLineNo">097</span> /** {@inheritDoc} */<a
name="line.97"></a>
+<span class="sourceLineNo">098</span> @Override<a name="line.98"></a>
+<span class="sourceLineNo">099</span> protected LeastSquaresProblem
getProblem(Collection<WeightedObservedPoint> observations) {<a
name="line.99"></a>
+<span class="sourceLineNo">100</span> // Prepare least-squares
problem.<a name="line.100"></a>
+<span class="sourceLineNo">101</span> final int len =
observations.size();<a name="line.101"></a>
+<span class="sourceLineNo">102</span> final double[] target = new
double[len];<a name="line.102"></a>
+<span class="sourceLineNo">103</span> final double[] weights = new
double[len];<a name="line.103"></a>
+<span class="sourceLineNo">104</span><a name="line.104"></a>
+<span class="sourceLineNo">105</span> int count = 0;<a
name="line.105"></a>
+<span class="sourceLineNo">106</span> for (WeightedObservedPoint obs :
observations) {<a name="line.106"></a>
+<span class="sourceLineNo">107</span> target[count] =
obs.getY();<a name="line.107"></a>
+<span class="sourceLineNo">108</span> weights[count] =
obs.getWeight();<a name="line.108"></a>
+<span class="sourceLineNo">109</span> ++count;<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> final
AbstractCurveFitter.TheoreticalValuesFunction model<a name="line.112"></a>
+<span class="sourceLineNo">113</span> = new
AbstractCurveFitter.TheoreticalValuesFunction(function,<a name="line.113"></a>
+<span class="sourceLineNo">114</span>
observations);<a name="line.114"></a>
+<span class="sourceLineNo">115</span><a name="line.115"></a>
+<span class="sourceLineNo">116</span> // Create an optimizer for
fitting the curve to the observed points.<a name="line.116"></a>
+<span class="sourceLineNo">117</span> return new
LeastSquaresBuilder().<a name="line.117"></a>
+<span class="sourceLineNo">118</span>
maxEvaluations(Integer.MAX_VALUE).<a name="line.118"></a>
+<span class="sourceLineNo">119</span> maxIterations(maxIter).<a
name="line.119"></a>
+<span class="sourceLineNo">120</span> start(initialGuess).<a
name="line.120"></a>
+<span class="sourceLineNo">121</span> target(target).<a
name="line.121"></a>
+<span class="sourceLineNo">122</span> weight(new
DiagonalMatrix(weights)).<a name="line.122"></a>
+<span class="sourceLineNo">123</span>
model(model.getModelFunction(), model.getModelFunctionJacobian()).<a
name="line.123"></a>
+<span class="sourceLineNo">124</span> build();<a
name="line.124"></a>
+<span class="sourceLineNo">125</span> }<a name="line.125"></a>
+<span class="sourceLineNo">126</span>}<a name="line.126"></a>
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+</pre>
+</div>
+</body>
+</html>
Propchange:
websites/production/commons/content/proper/commons-math/apidocs/src-html/org/apache/commons/math3/fitting/BasicCurveFitter.html
------------------------------------------------------------------------------
svn:eol-style = native