import com.imsl.math.*;

public abstract class ReflectivityModel extends Object {

    public static final int DELTAQ = 0;
    public static final int DELTAQOVERQ = 1;

    // Profile array is alternating rho real, rho imag, deltaZ
    public void setProfile( double[] inProfile) {
        numlayers = inProfile.length/3;
        deltaZ = new double[numlayers];
        rho = new Complex[numlayers];
        for(int i=0;i<numlayers;i++) {
            rho[i] = new Complex(inProfile[3*i], inProfile[3*i + 1]);
            deltaZ[i] = inProfile[3*i + 2];
        }
    }

    public void setDeltaZ( double[] inDeltaZ) {
         numlayers = inDeltaZ.length;
         deltaZ = new double[numlayers];
         for(int i=0;i<numlayers;i++) deltaZ[i]=inDeltaZ[i];
    }
 
    public void setRho( Complex[] inRho) {
         numlayers = inRho.length;
         rho = new Complex[numlayers];
         for(int i=0;i<numlayers;i++) rho[i] = new Complex(inRho[i]);
    }
 
    public double[] getDeltaZ() {
         return deltaZ;
    }
 
    public Complex[] getRho() {
         return rho;
    }
 
    public void setResolution( double inResolution) {
         resolution = inResolution;
    }

    public void setFronting( Complex inRho) {
         fronting = new Complex(inRho);
    }

    public void setBacking( Complex inRho) {
         backing = new Complex(inRho);
    }

    public double[] getProfile() {
        double[] outProfile = new double[numlayers*3];
        int j = 0;
        for(int i=0;i<numlayers;i++) {
            outProfile[j++] = rho[i].real();
            outProfile[j++] = rho[i].imag();
            outProfile[j++] = deltaZ[i];
        }
        return outProfile;
    }

    public double getResolution() {
        return resolution;
    }

    public Complex getFronting() {
        return new Complex(fronting);
    }

    public Complex getBacking() {
        return new Complex(backing);
    }

    public int getNumLayers() {
        return numlayers;
    }

    public abstract Complex getReflectivity( double inQ);
    public abstract double getSmearedReflectivity(double inQ, int resolution);

    private double resolution = 0.000001;
    private double[] deltaZ;
    private Complex fronting, backing;
    private Complex[] rho;
    private int numlayers;

}