// Member of the utility package.
// Modified July 20, 2009 to handle very large arrays
// in some contexts.
package nom.tam.util;

/*
 * This code is part of the Java FITS library developed 1996-2012 by T.A. McGlynn (NASA/GSFC)
 * The code is available in the public domain and may be copied, modified and used
 * by anyone in any fashion for any purpose without restriction. 
 * 
 * No warranty regarding correctness or performance of this code is given or implied.
 * Users may contact the author if they have questions or concerns.
 * 
 * The author would like to thank many who have contributed suggestions, 
 * enhancements and bug fixes including:
 * David Glowacki, R.J. Mathar, Laurent Michel, Guillaume Belanger,
 * Laurent Bourges, Rose Early, Fred Romelfanger, Jorgo Baker, A. Kovacs, V. Forchi, J.C. Segovia,
 * Booth Hartley and Jason Weiss.  
 * I apologize to any contributors whose names may have been inadvertently omitted.
 * 
 *      Tom McGlynn
 */
import java.lang.reflect.*;
import java.util.Arrays;

/** This is a package of static functions which perform
 * computations on arrays.  Generally these routines attempt
 * to complete without throwing errors by ignoring data
 * they cannot understand.
 */
public class ArrayFuncs implements PrimitiveInfo {

    /** Compute the size of an object.  Note that this only handles
     * arrays or scalars of the primitive objects and Strings.  It
     * returns 0 for any object array element it does not understand.
     *
     * @param o The object whose size is desired.
     * @deprecated May silently underestimate the size
     *   if the size > 2 GB.
     */
    public static int computeSize(Object o) {
        return (int) computeLSize(o);
    }

    public static long computeLSize(Object o) {

        if (o == null) {
            return 0;
        }

        long size = 0;
        String classname = o.getClass().getName();
        if (classname.substring(0, 2).equals("[[")) {

            for (int i = 0; i < ((Object[]) o).length; i += 1) {
                size += computeLSize(((Object[]) o)[i]);
            }
            return size;
        }

        if (classname.charAt(0) == '[' && classname.charAt(1) != 'L') {
            char c = classname.charAt(1);

            for (int i = 0; i < PrimitiveInfo.suffixes.length; i += 1) {
                if (c == PrimitiveInfo.suffixes[i]) {
                    return (long) (Array.getLength(o)) * PrimitiveInfo.sizes[i];
                }
            }
            return 0;
        }

        // Do we have a non-primitive array?
        if (classname.charAt(0) == '[') {
            int len = 0;
            for (int i = 0; i < Array.getLength(o); i += 1) {
                len += computeLSize(Array.get(o, i));
            }
            return len;
        }

        // Now a few special scalar objects.
        if (classname.substring(0, 10).equals("java.lang.")) {
            classname = classname.substring(10, classname.length());
            if (classname.equals("Integer") || classname.equals("Float")) {
                return 4;
            } else if (classname.equals("Double") || classname.equals("Long")) {
                return 8;
            } else if (classname.equals("Short") || classname.equals("Char")) {
                return 2;
            } else if (classname.equals("Byte") || classname.equals("Boolean")) {
                return 1;
            } else if (classname.equals("String")) {
                return ((String) o).length();
            } else {
                return 0;
            }
        }
        return 0;
    }

    /** Count the number of elements in an array.
     *  @deprecated  May silently underestimate
     *               size if number is > 2 G.
     */
    public static int nElements(Object o) {
        return (int) nLElements(o);
    }

    /** Count the number of elements in an array.
     *  @deprecated  May silently underestimate
     *               size if number is > 2 G.
     */
    public static long nLElements(Object o) {

        if (o == null) {
            return 0;
        }

        String classname = o.getClass().getName();
        if (classname.charAt(1) == '[') {
            int count = 0;
            for (int i = 0; i < ((Object[]) o).length; i += 1) {
                count += nLElements(((Object[]) o)[i]);
            }
            return count;

        } else if (classname.charAt(0) == '[') {
            return Array.getLength(o);

        } else {
            return 1;
        }
    }

    /** Try to create a deep clone of an Array or a standard clone of a scalar.
     * The object may comprise arrays of
     * any primitive type or any Object type which implements Cloneable.
     * However, if the Object is some kind of collection, e.g., a Vector
     * then only a shallow copy of that object is made.  I.e., deep refers
     * only to arrays.
     *
     * @param o The object to be copied.
     */
    public static Object deepClone(Object o) {

        if (o == null) {
            return null;
        }

        String classname = o.getClass().getName();

        // Is this an array?
        if (classname.charAt(0) != '[') {
            return genericClone(o);
        }

        // Check if this is a 1D primitive array.
        if (classname.charAt(1) != '[' && classname.charAt(1) != 'L') {
            try {
                // Some compilers (SuperCede, e.g.) still
                // think you have to catch this...
                if (false) {
                    throw new CloneNotSupportedException();
                }
                switch (classname.charAt(1)) {
                    case 'B':
                        return ((byte[]) o).clone();
                    case 'Z':
                        return ((boolean[]) o).clone();
                    case 'C':
                        return ((char[]) o).clone();
                    case 'S':
                        return ((short[]) o).clone();
                    case 'I':
                        return ((int[]) o).clone();
                    case 'J':
                        return ((long[]) o).clone();
                    case 'F':
                        return ((float[]) o).clone();
                    case 'D':
                        return ((double[]) o).clone();
                    default:
                        System.err.println("Unknown primtive array class:" + classname);
                        return null;

                }
            } catch (CloneNotSupportedException e) {
            }
        }

        // Get the base type.
        int ndim = 1;
        while (classname.charAt(ndim) == '[') {
            ndim += 1;
        }
        Class baseClass;
        if (classname.charAt(ndim) != 'L') {
            baseClass = getBaseClass(o);
        } else {
            try {
                baseClass = Class.forName(classname.substring(ndim + 1, classname.length() - 1));
            } catch (ClassNotFoundException e) {
                System.err.println("Internal error: class definition inconsistency: " + classname);
                return null;
            }
        }

        // Allocate the array but make all but the first dimension 0.
        int[] dims = new int[ndim];
        dims[0] = Array.getLength(o);
        for (int i = 1; i < ndim; i += 1) {
            dims[i] = 0;
        }

        Object copy = ArrayFuncs.newInstance(baseClass, dims);

        // Now fill in the next level down by recursion.
        for (int i = 0; i < dims[0]; i += 1) {
            Array.set(copy, i, deepClone(Array.get(o, i)));
        }

        return copy;
    }

    /** Clone an Object if possible.
     *
     * This method returns an Object which is a clone of the
     * input object.  It checks if the method implements the
     * Cloneable interface and then uses reflection to invoke
     * the clone method.  This can't be done directly since
     * as far as the compiler is concerned the clone method for
     * Object is protected and someone could implement Cloneable but
     * leave the clone method protected.  The cloning can fail in a
     * variety of ways which are trapped so that it returns null instead.
     * This method will generally create a shallow clone.  If you
     * wish a deep copy of an array the method deepClone should be used.
     *
     * @param o The object to be cloned.
     */
    public static Object genericClone(Object o) {

        if (!(o instanceof Cloneable)) {
            return null;
        }

        Class[] argTypes = new Class[0];
        Object[] args = new Object[0];
        Class type = o.getClass();

        try {
            return type.getMethod("clone", argTypes).invoke(o, args);
        } catch (Exception e) {
            if (type.isArray()) {
                return deepClone(o);
            }
            // Implements cloneable, but does not
            // apparently make clone public.
            return null;
        }
    }

    /** Copy one array into another.
     * This function copies the contents of one array
     * into a previously allocated array.
     * The arrays must agree in type and size.
     * @param original The array to be copied.
     * @param copy     The array to be copied into.  This
     *                 array must already be fully allocated.
     */
    public static void copyArray(Object original, Object copy) {
        String oname = original.getClass().getName();
        String cname = copy.getClass().getName();

        if (!oname.equals(cname)) {
            return;
        }

        if (oname.charAt(0) != '[') {
            return;
        }

        if (oname.charAt(1) == '[') {
            Object[] x = (Object[]) original;
            Object[] y = (Object[]) copy;
            if (x.length != y.length) {
                return;
            }
            for (int i = 0; i < x.length; i += 1) {
                copyArray(x, y);
            }
        }
        int len = Array.getLength(original);

        System.arraycopy(original, 0, copy, 0, len);
    }

    /** Find the dimensions of an object.
     *
     * This method returns an integer array with the dimensions
     * of the object o which should usually be an array.
     *
     * It returns an array of dimension 0 for scalar objects
     * and it returns -1 for dimension which have not been allocated,
     *  e.g., int[][][] x = new int[100][][]; should return [100,-1,-1].
     *
     * @param o The object to get the dimensions of.
     */
    public static int[] getDimensions(Object o) {

        if (o == null) {
            return null;
        }

        String classname = o.getClass().getName();

        int ndim = 0;

        while (classname.charAt(ndim) == '[') {
            ndim += 1;
        }

        int[] dimens = new int[ndim];

        for (int i = 0; i < ndim; i += 1) {
            dimens[i] = -1;  // So that we can distinguish a null from a 0 length.
        }

        for (int i = 0; i < ndim; i += 1) {
            dimens[i] = java.lang.reflect.Array.getLength(o);
            if (dimens[i] == 0) {
                return dimens;
            }
            if (i != ndim - 1) {
                o = ((Object[]) o)[0];
                if (o == null) {
                    return dimens;
                }
            }
        }
        return dimens;
    }

    /** This routine returns the base array of a multi-dimensional
     *  array.  I.e., a one-d array of whatever the array is composed
     *  of.  Note that arrays are not guaranteed to be rectangular,
     *  so this returns o[0][0]....
     */
    public static Object getBaseArray(Object o) {
        String cname = o.getClass().getName();
        if (cname.charAt(1) == '[') {
            return getBaseArray(((Object[]) o)[0]);
        } else {
            return o;
        }
    }

    /** This routine returns the base class of an object.  This is just
     * the class of the object for non-arrays.
     */
    public static Class getBaseClass(Object o) {

        if (o == null) {
            return Void.TYPE;
        }

        String className = o.getClass().getName();

        int dims = 0;
        while (className.charAt(dims) == '[') {
            dims += 1;
        }

        if (dims == 0) {
            return o.getClass();
        }

        char c = className.charAt(dims);
        for (int i = 0; i < PrimitiveInfo.suffixes.length; i += 1) {
            if (c == PrimitiveInfo.suffixes[i]) {
                return PrimitiveInfo.classes[i];
            }
        }

        if (c == 'L') {
            try {
                return Class.forName(className.substring(dims + 1, className.length() - 1));
            } catch (ClassNotFoundException e) {
                return null;
            }
        }
        return null;
    }

    /** This routine returns the size of the base element of an array.
     * @param o The array object whose base length is desired.
     * @return the size of the object in bytes, 0 if null, or
     * -1 if not a primitive array.
     */
    public static int getBaseLength(Object o) {

        if (o == null) {
            return 0;
        }

        String className = o.getClass().getName();

        int dims = 0;

        while (className.charAt(dims) == '[') {
            dims += 1;
        }

        if (dims == 0) {
            return -1;
        }

        char c = className.charAt(dims);
        for (int i = 0; i < PrimitiveInfo.suffixes.length; i += 1) {
            if (c == PrimitiveInfo.suffixes[i]) {
                return PrimitiveInfo.sizes[i];
            }
        }
        return -1;
    }

    /** Create an array and populate it with a test pattern.
     *
     * @param baseType  The base type of the array.  This is expected to
     *                  be a numeric type, but this is not checked.
     * @param dims      The desired dimensions.
     * @return An array object populated with a simple test pattern.
     */
    public static Object generateArray(Class baseType, int[] dims) {

        // Generate an array and populate it with a test pattern of
        // data.

        Object x = ArrayFuncs.newInstance(baseType, dims);
        testPattern(x, (byte) 0);
        return x;
    }

    /** Just create a simple pattern cycling through valid byte values.
     * We use bytes because they can be cast to any other numeric type.
     * @param o      The array in which the test pattern is to be set.
     * @param start  The value for the first element.
     */
    public static byte testPattern(Object o, byte start) {

        int[] dims = getDimensions(o);
        if (dims.length > 1) {
            for (int i = 0; i < ((Object[]) o).length; i += 1) {
                start = testPattern(((Object[]) o)[i], start);
            }

        } else if (dims.length == 1) {
            for (int i = 0; i < dims[0]; i += 1) {
                java.lang.reflect.Array.setByte(o, i, start);
                start += 1;
            }
        }
        return start;
    }

    /** Generate a description of an array (presumed rectangular).
     * @param o The array to be described.
     */
    public static String arrayDescription(Object o) {

        Class base = getBaseClass(o);
        if (base == Void.TYPE) {
            return "NULL";
        }

        int[] dims = getDimensions(o);

        StringBuffer desc = new StringBuffer();

        // Note that all instances Class describing a given class are
        // the same so we can use == here.
        boolean found = false;

        for (int i = 0; i < PrimitiveInfo.classes.length; i += 1) {
            if (base == PrimitiveInfo.classes[i]) {
                found = true;
                desc.append(PrimitiveInfo.types[i]);
                break;
            }
        }

        if (!found) {
            desc.append(base.getName());
        }

        if (dims != null) {
            desc.append("[");
            for (int i = 0; i < dims.length; i += 1) {
                desc.append("" + dims[i]);
                if (i < dims.length - 1) {
                    desc.append("][");
                }
            }
            desc.append("]");
        }
        return new String(desc);
    }

    /** Examine the structure of an array in detail.
     * @param o The array to be examined.
     */
    public static void examinePrimitiveArray(Object o) {
        String className = o.getClass().getName();

        // If we have a two-d array, or if the array is a one-d array
        // of Objects, then recurse over the next dimension.  We handle
        // Object specially because each element could itself be an array.
        if (className.substring(0, 2).equals("[[")
                || className.equals("[Ljava.lang.Object;")) {
            System.out.println("[");
            for (int i = 0; i < ((Object[]) o).length; i += 1) {
                examinePrimitiveArray(((Object[]) o)[i]);
            }
            System.out.print("]");
        } else if (className.charAt(0) != '[') {
            System.out.println(className);
        } else {
            System.out.println("[" + java.lang.reflect.Array.getLength(o) + "]"
                    + className.substring(1));
        }
    }

    /** Given an array of arbitrary dimensionality return
     * the array flattened into a single dimension.
     * @param input The input array.
     */
    public static Object flatten(Object input) {

        int[] dimens = getDimensions(input);
        if (dimens.length <= 1) {
            return input;
        }
        int size = 1;
        for (int i = 0; i < dimens.length; i += 1) {
            size *= dimens[i];
        }

        Object flat = ArrayFuncs.newInstance(getBaseClass(input), size);

        if (size == 0) {
            return flat;
        }

        int offset = 0;
        doFlatten(input, flat, offset);
        return flat;
    }

    /** This routine does the actually flattening of multi-dimensional
     * arrays.
     * @param input  The input array to be flattened.
     * @param output The flattened array.
     * @param offset The current offset within the output array.
     * @return       The number of elements within the array.
     */
    protected static int doFlatten(Object input, Object output, int offset) {

        String classname = input.getClass().getName();
        if (classname.charAt(0) != '[') {
            throw new RuntimeException("Attempt to flatten non-array");
        }
        int size = Array.getLength(input);

        if (classname.charAt(1) != '[') {
            System.arraycopy(input, 0, output, offset, size);
            return size;
        }
        int total = 0;
        Object[] xx = (Object[]) input;
        for (int i = 0; i < size; i += 1) {
            int len = doFlatten(xx[i], output, offset + total);
            total += len;
        }
        return total;
    }

    /** Curl an input array up into a multi-dimensional array.
     *
     * @param input The one dimensional array to be curled.
     * @param dimens The desired dimensions
     * @return The curled array.
     */
    public static Object curl(Object input, int[] dimens) {

        if (input == null) {
            return null;
        }
        String classname = input.getClass().getName();
        if (classname.charAt(0) != '[' || classname.charAt(1) == '[') {
            throw new RuntimeException("Attempt to curl non-1D array");
        }

        int size = Array.getLength(input);

        int test = 1;
        for (int i = 0; i < dimens.length; i += 1) {
            test *= dimens[i];
        }

        if (test != size) {
            throw new RuntimeException("Curled array does not fit desired dimensions");
        }

        Class base = getBaseClass(input);

        Object newArray = ArrayFuncs.newInstance(base, dimens);

        int offset = 0;

        doCurl(input, newArray, dimens, offset);
        return newArray;

    }

    /** Do the curling of the 1-d to multi-d array.
     * @param input  The 1-d array to be curled.
     * @param output The multi-dimensional array to be filled.
     * @param dimens The desired output dimensions.
     * @param offset The current offset in the input array.
     * @return       The number of elements curled.
     */
    protected static int doCurl(Object input, Object output,
            int[] dimens, int offset) {

        if (dimens.length == 1) {
            System.arraycopy(input, offset, output, 0, dimens[0]);
            return dimens[0];
        }

        int total = 0;
        int[] xdimens = new int[dimens.length - 1];
        for (int i = 1; i < dimens.length; i += 1) {
            xdimens[i - 1] = dimens[i];
        }

        for (int i = 0; i < dimens[0]; i += 1) {
            total += doCurl(input, ((Object[]) output)[i], xdimens, offset + total);
        }
        return total;
    }

    /** Create an array of a type given by new type with
     * the dimensionality given in array.
     * @param array   A possibly multidimensional array to be converted.
     * @param newType The desired output type.  This should be one of the
     *                class descriptors for primitive numeric data, e.g., double.type.
     */
    public static Object mimicArray(Object array, Class newType) {

        String classname = array.getClass().getName();
        if (classname.charAt(0) != '[') {
            return null;
        }

        int dims = 1;

        while (classname.charAt(dims) == '[') {
            dims += 1;
        }

        Object mimic;

        if (dims > 1) {

            Object[] xarray = (Object[]) array;
            int[] dimens = new int[dims];
            dimens[0] = xarray.length;  // Leave other dimensions at 0.


            mimic = ArrayFuncs.newInstance(newType, dimens);

            for (int i = 0; i < xarray.length; i += 1) {
                Object temp = mimicArray(xarray[i], newType);
                ((Object[]) mimic)[i] = temp;
            }

        } else {
            mimic = ArrayFuncs.newInstance(newType, Array.getLength(array));
        }

        return mimic;
    }

    /** Convert an array to a specified type.  This method supports conversions
     * only among the primitive numeric types.
     * @param array   A possibly multidimensional array to be converted.
     * @param newType The desired output type.  This should be one of the
     *                class descriptors for primitive numeric data, e.g., double.type.
     * @param reuse   If set, and the requested type is the same as the
     *                  original, then the original is returned.
     */
    public static Object convertArray(Object array, Class newType, boolean reuse) {

        if (getBaseClass(array) == newType && reuse) {
            return array;
        } else {
            return convertArray(array, newType);
        }
    }

    /** Convert an array to a specified type.  This method supports conversions
     * only among the primitive numeric types.
     * @param array   A possibly multidimensional array to be converted.
     * @param newType The desired output type.  This should be one of the
     *                class descriptors for primitive numeric data, e.g., double.type.
     */
    public static Object convertArray(Object array, Class newType) {

        /* We break this up into two steps so that users
         * can reuse an array many times and only allocate a
         * new array when needed.
         */

        /* First create the full new array. */
        Object mimic = mimicArray(array, newType);
        if (mimic == null) {
            return mimic;
        }

        /* Now copy the info into the new array */
        copyInto(array, mimic);

        return mimic;
    }

    /** Copy an array into an array of a different type.
     *  The dimensions and dimensionalities of the two
     *  arrays should be the same.
     *  @param array  The original array.
     *  @param mimic  The array mimicking the original.
     */
    public static void copyInto(Object array, Object mimic) {

        String classname = array.getClass().getName();
        if (classname.charAt(0) != '[') {
            return;
        }

        /* Do multidimensional arrays recursively */
        if (classname.charAt(1) == '[') {

            for (int i = 0; i < ((Object[]) array).length; i += 1) {
                copyInto(((Object[]) array)[i], ((Object[]) mimic)[i]);
            }

        } else {

            byte[] xbarr;
            short[] xsarr;
            char[] xcarr;
            int[] xiarr;
            long[] xlarr;
            float[] xfarr;
            double[] xdarr;

            Class base = getBaseClass(array);
            Class newType = getBaseClass(mimic);

            if (base == byte.class) {
                byte[] barr = (byte[]) array;

                if (newType == byte.class) {
                    System.arraycopy(array, 0, mimic, 0, barr.length);

                } else if (newType == short.class) {
                    xsarr = (short[]) mimic;
                    for (int i = 0; i < barr.length; i += 1) {
                        xsarr[i] = barr[i];
                    }

                } else if (newType == char.class) {
                    xcarr = (char[]) mimic;
                    for (int i = 0; i < barr.length; i += 1) {
                        xcarr[i] = (char) barr[i];
                    }

                } else if (newType == int.class) {
                    xiarr = (int[]) mimic;
                    for (int i = 0; i < barr.length; i += 1) {
                        xiarr[i] = barr[i];
                    }

                } else if (newType == long.class) {
                    xlarr = (long[]) mimic;
                    for (int i = 0; i < barr.length; i += 1) {
                        xlarr[i] = barr[i];
                    }

                } else if (newType == float.class) {
                    xfarr = (float[]) mimic;
                    for (int i = 0; i < barr.length; i += 1) {
                        xfarr[i] = barr[i];
                    }

                } else if (newType == double.class) {
                    xdarr = (double[]) mimic;
                    for (int i = 0; i < barr.length; i += 1) {
                        xdarr[i] = barr[i];
                    }
                }

            } else if (base == short.class) {
                short[] sarr = (short[]) array;

                if (newType == byte.class) {
                    xbarr = (byte[]) mimic;
                    for (int i = 0; i < sarr.length; i += 1) {
                        xbarr[i] = (byte) sarr[i];
                    }

                } else if (newType == short.class) {
                    System.arraycopy(array, 0, mimic, 0, sarr.length);

                } else if (newType == char.class) {
                    xcarr = (char[]) mimic;
                    for (int i = 0; i < sarr.length; i += 1) {
                        xcarr[i] = (char) sarr[i];
                    }

                } else if (newType == int.class) {
                    xiarr = (int[]) mimic;
                    for (int i = 0; i < sarr.length; i += 1) {
                        xiarr[i] = sarr[i];
                    }

                } else if (newType == long.class) {
                    xlarr = (long[]) mimic;
                    for (int i = 0; i < sarr.length; i += 1) {
                        xlarr[i] = sarr[i];
                    }

                } else if (newType == float.class) {
                    xfarr = (float[]) mimic;
                    for (int i = 0; i < sarr.length; i += 1) {
                        xfarr[i] = sarr[i];
                    }

                } else if (newType == double.class) {
                    xdarr = (double[]) mimic;
                    for (int i = 0; i < sarr.length; i += 1) {
                        xdarr[i] = sarr[i];
                    }
                }

            } else if (base == char.class) {
                char[] carr = (char[]) array;

                if (newType == byte.class) {
                    xbarr = (byte[]) mimic;
                    for (int i = 0; i < carr.length; i += 1) {
                        xbarr[i] = (byte) carr[i];
                    }

                } else if (newType == short.class) {
                    xsarr = (short[]) mimic;
                    for (int i = 0; i < carr.length; i += 1) {
                        xsarr[i] = (short) carr[i];
                    }

                } else if (newType == char.class) {
                    System.arraycopy(array, 0, mimic, 0, carr.length);

                } else if (newType == int.class) {
                    xiarr = (int[]) mimic;
                    for (int i = 0; i < carr.length; i += 1) {
                        xiarr[i] = carr[i];
                    }

                } else if (newType == long.class) {
                    xlarr = (long[]) mimic;
                    for (int i = 0; i < carr.length; i += 1) {
                        xlarr[i] = carr[i];
                    }

                } else if (newType == float.class) {
                    xfarr = (float[]) mimic;
                    for (int i = 0; i < carr.length; i += 1) {
                        xfarr[i] = carr[i];
                    }

                } else if (newType == double.class) {
                    xdarr = (double[]) mimic;
                    for (int i = 0; i < carr.length; i += 1) {
                        xdarr[i] = carr[i];
                    }
                }

            } else if (base == int.class) {
                int[] iarr = (int[]) array;

                if (newType == byte.class) {
                    xbarr = (byte[]) mimic;
                    for (int i = 0; i < iarr.length; i += 1) {
                        xbarr[i] = (byte) iarr[i];
                    }

                } else if (newType == short.class) {
                    xsarr = (short[]) mimic;
                    for (int i = 0; i < iarr.length; i += 1) {
                        xsarr[i] = (short) iarr[i];
                    }

                } else if (newType == char.class) {
                    xcarr = (char[]) mimic;
                    for (int i = 0; i < iarr.length; i += 1) {
                        xcarr[i] = (char) iarr[i];
                    }

                } else if (newType == int.class) {
                    System.arraycopy(array, 0, mimic, 0, iarr.length);

                } else if (newType == long.class) {
                    xlarr = (long[]) mimic;
                    for (int i = 0; i < iarr.length; i += 1) {
                        xlarr[i] = iarr[i];
                    }

                } else if (newType == float.class) {
                    xfarr = (float[]) mimic;
                    for (int i = 0; i < iarr.length; i += 1) {
                        xfarr[i] = iarr[i];
                    }

                } else if (newType == double.class) {
                    xdarr = (double[]) mimic;
                    for (int i = 0; i < iarr.length; i += 1) {
                        xdarr[i] = iarr[i];
                    }
                }


            } else if (base == long.class) {
                long[] larr = (long[]) array;

                if (newType == byte.class) {
                    xbarr = (byte[]) mimic;
                    for (int i = 0; i < larr.length; i += 1) {
                        xbarr[i] = (byte) larr[i];
                    }

                } else if (newType == short.class) {
                    xsarr = (short[]) mimic;
                    for (int i = 0; i < larr.length; i += 1) {
                        xsarr[i] = (short) larr[i];
                    }

                } else if (newType == char.class) {
                    xcarr = (char[]) mimic;
                    for (int i = 0; i < larr.length; i += 1) {
                        xcarr[i] = (char) larr[i];
                    }

                } else if (newType == int.class) {
                    xiarr = (int[]) mimic;
                    for (int i = 0; i < larr.length; i += 1) {
                        xiarr[i] = (int) larr[i];
                    }

                } else if (newType == long.class) {
                    System.arraycopy(array, 0, mimic, 0, larr.length);

                } else if (newType == float.class) {
                    xfarr = (float[]) mimic;
                    for (int i = 0; i < larr.length; i += 1) {
                        xfarr[i] = (float) larr[i];
                    }

                } else if (newType == double.class) {
                    xdarr = (double[]) mimic;
                    for (int i = 0; i < larr.length; i += 1) {
                        xdarr[i] = (double) larr[i];
                    }
                }

            } else if (base == float.class) {
                float[] farr = (float[]) array;

                if (newType == byte.class) {
                    xbarr = (byte[]) mimic;
                    for (int i = 0; i < farr.length; i += 1) {
                        xbarr[i] = (byte) farr[i];
                    }

                } else if (newType == short.class) {
                    xsarr = (short[]) mimic;
                    for (int i = 0; i < farr.length; i += 1) {
                        xsarr[i] = (short) farr[i];
                    }

                } else if (newType == char.class) {
                    xcarr = (char[]) mimic;
                    for (int i = 0; i < farr.length; i += 1) {
                        xcarr[i] = (char) farr[i];
                    }

                } else if (newType == int.class) {
                    xiarr = (int[]) mimic;
                    for (int i = 0; i < farr.length; i += 1) {
                        xiarr[i] = (int) farr[i];
                    }

                } else if (newType == long.class) {
                    xlarr = (long[]) mimic;
                    for (int i = 0; i < farr.length; i += 1) {
                        xlarr[i] = (long) farr[i];
                    }

                } else if (newType == float.class) {
                    System.arraycopy(array, 0, mimic, 0, farr.length);

                } else if (newType == double.class) {
                    xdarr = (double[]) mimic;
                    for (int i = 0; i < farr.length; i += 1) {
                        xdarr[i] = farr[i];
                    }
                }


            } else if (base == double.class) {
                double[] darr = (double[]) array;

                if (newType == byte.class) {
                    xbarr = (byte[]) mimic;
                    for (int i = 0; i < darr.length; i += 1) {
                        xbarr[i] = (byte) darr[i];
                    }

                } else if (newType == short.class) {
                    xsarr = (short[]) mimic;
                    for (int i = 0; i < darr.length; i += 1) {
                        xsarr[i] = (short) darr[i];
                    }

                } else if (newType == char.class) {
                    xcarr = (char[]) mimic;
                    for (int i = 0; i < darr.length; i += 1) {
                        xcarr[i] = (char) darr[i];
                    }

                } else if (newType == int.class) {
                    xiarr = (int[]) mimic;
                    for (int i = 0; i < darr.length; i += 1) {
                        xiarr[i] = (int) darr[i];
                    }

                } else if (newType == long.class) {
                    xlarr = (long[]) mimic;
                    for (int i = 0; i < darr.length; i += 1) {
                        xlarr[i] = (long) darr[i];
                    }

                } else if (newType == float.class) {
                    xfarr = (float[]) mimic;
                    for (int i = 0; i < darr.length; i += 1) {
                        xfarr[i] = (float) darr[i];
                    }

                } else if (newType == double.class) {
                    System.arraycopy(array, 0, mimic, 0, darr.length);
                }
            }
        }

        return;

    }

    /** Allocate an array dynamically. The Array.newInstance method
     *  does not throw an error when there is insufficient memory
     *  and silently returns a null.
     *  @param cl	The class of the array.
     *  @param dim      The dimension of the array.
     *  @return The allocated array.
     *  @throws An OutOfMemoryError if insufficient space is available.
     */
    public static Object newInstance(Class cl, int dim) {

        Object o = Array.newInstance(cl, dim);
        if (o == null) {
            String desc = cl + "[" + dim + "]";
            throw new OutOfMemoryError("Unable to allocate array: " + desc);
        }
        return o;
    }

    /** Allocate an array dynamically. The Array.newInstance method
     *  does not throw an error and silently returns a null.
     *
     *  @param cl	The class of the array.
     *  @param dims     The dimensions of the array.
     *  @return The allocated array.
     *  @throws An OutOfMemoryError if insufficient space is available.
     */
    public static Object newInstance(Class cl, int[] dims) {

        if (dims.length == 0) {
            // Treat a scalar as a 1-d array of length 1
            dims = new int[]{1};
        }

        Object o = Array.newInstance(cl, dims);
        if (o == null) {
            String desc = cl + "[";
            String comma = "";
            for (int i = 0; i < dims.length; i += 1) {
                desc += comma + dims[i];
                comma = ",";
            }
            desc += "]";
            throw new OutOfMemoryError("Unable to allocate array: " + desc);
        }
        return o;
    }

    /** Are two objects equal?  Arrays have the standard object equals
     *  method which only returns true if the two object are the same.
     *  This method returns true if every element of the arrays match.
     *  The inputs may be of any dimensionality.  The dimensionality
     *  and dimensions of the arrays must match as well as any elements.
     *  If the elements are non-primitive. non-array objects, then the
     *  equals method is called for each element.
     *  If both elements are multi-dimensional arrays, then
     *  the method recurses.
     */
    public static boolean arrayEquals(Object x, Object y) {
        return arrayEquals(x, y, 0, 0);
    }

    /** Are two objects equal?  Arrays have the standard object equals
     *  method which only returns true if the two object are the same.
     *  This method returns true if every element of the arrays match.
     *  The inputs may be of any dimensionality.  The dimensionality
     *  and dimensions of the arrays must match as well as any elements.
     *  If the elements are non-primitive. non-array objects, then the
     *  equals method is called for each element.
     *  If both elements are multi-dimensional arrays, then
     *  the method recurses.
     */
    public static boolean arrayEquals(Object x, Object y, double tolf, double told) {

        // Handle the special cases first.
        // We treat null == null so that two object arrays
        // can match if they have matching null elements.
        if (x == null && y == null) {
            return true;
        }

        if (x == null || y == null) {
            return false;
        }

        Class xClass = x.getClass();
        Class yClass = y.getClass();

        if (xClass != yClass) {
            return false;
        }

        if (!xClass.isArray()) {
            return x.equals(y);

        } else {
            if (xClass.equals(int[].class)) {
                return Arrays.equals((int[]) x, (int[]) y);

            } else if (xClass.equals(double[].class)) {
                if (told == 0) {
                    return Arrays.equals((double[]) x, (double[]) y);
                } else {
                    return doubleArrayEquals((double[]) x, (double[]) y, told);
                }

            } else if (xClass.equals(long[].class)) {
                return Arrays.equals((long[]) x, (long[]) y);

            } else if (xClass.equals(float[].class)) {
                if (tolf == 0) {
                    return Arrays.equals((float[]) x, (float[]) y);
                } else {
                    return floatArrayEquals((float[]) x, (float[]) y, (float) tolf);
                }

            } else if (xClass.equals(byte[].class)) {
                return Arrays.equals((byte[]) x, (byte[]) y);

            } else if (xClass.equals(short[].class)) {
                return Arrays.equals((short[]) x, (short[]) y);

            } else if (xClass.equals(char[].class)) {
                return Arrays.equals((char[]) x, (char[]) y);

            } else if (xClass.equals(boolean[].class)) {
                return Arrays.equals((boolean[]) x, (boolean[]) y);

            } else {
                // Non-primitive and multidimensional arrays can be
                // cast to Object[]
                Object[] xo = (Object[]) x;
                Object[] yo = (Object[]) y;
                if (xo.length != yo.length) {
                    return false;
                }
                for (int i = 0; i < xo.length; i += 1) {
                    if (!arrayEquals(xo[i], yo[i], tolf, told)) {
                        return false;
                    }
                }

                return true;

            }
        }
    }

    /** Compare two double arrays using a given tolerance */
    public static boolean doubleArrayEquals(double[] x, double[] y, double tol) {

        for (int i = 0; i < x.length; i += 1) {
            if (x[i] == 0) {
                return y[i] == 0;
            }
            if (Math.abs((y[i] - x[i]) / x[i]) > tol) {
                return false;
            }
        }
        return true;
    }

    /** Compare two float arrays using a given tolerance */
    public static boolean floatArrayEquals(float[] x, float[] y, float tol) {

        for (int i = 0; i < x.length; i += 1) {
            if (x[i] == 0) {
                return y[i] == 0;
            }
            if (Math.abs((y[i] - x[i]) / x[i]) > tol) {
                return false;
            }
        }
        return true;
    }

    /** Dump an array on the given print steam */
    public static void dumpArray(java.io.PrintStream p, Object arr) {
        // Get the dimensionality and then dump.
        if (arr == null) {
            p.print("null ");
        } else {
            Class nm = arr.getClass();
            if (nm.isArray()) {
                p.print("[");
                for (int i = 0; i < java.lang.reflect.Array.getLength(arr); i += 1) {
                    dumpArray(p, java.lang.reflect.Array.get(arr, i));
                }
                p.print("]\n");
            } else {
                p.print(" " + arr.toString() + " ");
            }
        }
    }
}