/*      This example of the use of both BiCubicInterpolation and BiCubicSpline demonstrates
        interpolation within a data set of the refractive indices of
        water tabulated as a function of both wavelength and temperature.

        Michael Thomas Flanagan
        Created:    May 2003
        Updated:    1 May 2005, 10 January 2011

*/

import flanagan.interpolation.*;

public class BiCubicExampleTwo{

        public static void main(String arg[]){

        // Array of wavelengths (m)
        double[] wavelength = {404.6e-9, 589.32e-9, 706.52e-9};
        // Array of temperatures (Celsius)
        double[] temperature = {0,10,20,30,40,50,60,70,80,90,100};
        // Two dimensional array of corresponding refractive indices
        double[][] refrindex = {{1.34359, 1.34351, 1.34287, 1.3418, 1.34039, 1.33867, 1.33669, 1.33447, 1.33204, 1.32942, 1.32663}, {1.33346, 1.33341, 1.33283, 1.33184, 1.33052, 1.32892, 1.32707, 1.325, 1.32274, 1.32029, 1.31766}, {1.33086, 1.33073, 1.33007, 1.32903, 1.32766, 1.32603, 1.32417, 1.32209, 1.31983, 1.31739, 1.31481}};
        // Interpolation variables
        double x1, x2, y1, y2;

        // Create a BiCubicSpline instance and initialise it to the data stored in the arrays wavelength, temperature and refrindex
        // also calculates the required derivatives for susequent calls to this method
        BiCubicSpline bcs = new BiCubicSpline(wavelength, temperature, refrindex);

        // Create a BiCubicInterpolation instance and initialise it to the data stored in the arrays wavelength, temperature and refrindex
        // Derivatives calculated internally using supplied data points
        BiCubicInterpolation bci = new BiCubicInterpolation(wavelength, temperature, refrindex, 0);


        // First interpolation at a wavelength of 404.6 nm and a temperature of 20 degrees Celsius
        x1 = 4.046e-7;
        x2 = 20.0;
        y1 = bcs.interpolate(x1, x2);
        System.out.println("BicubicSpline:        The refractive index of water at " + x1*1.0e9 + " nm and " + x2 + " C is "+ y1);
        y2 = bci.interpolate(x1, x2);
        System.out.println("BicubicInterpolation: The refractive index of water at " + x1*1.0e9 + " nm and " + x2 + " C is "+ y2);

        // Second interpolation at a wavelength of 590 nm and a temperature of 33 degrees Celsius
        x1 = 5.9e-7;
        x2 = 33.0;
        y1 = bcs.interpolate(x1, x2);
        System.out.println("BicubicSpline:        The refractive index of water at " + x1*1.0e9 + " nm and " + x2 + " C is "+ y1);
        y2 = bci.interpolate(x1, x2);
        System.out.println("BicubicInterpolation: The refractive index of water at " + x1*1.0e9 + " nm and " + x2 + " C is "+ y2);

    }
}