//******************************************************************************
//** SCATMECH: Polarized Light Scattering C++ Class Library
//** 
//** File: allrough.cpp
//**
//** Thomas A. Germer
//** Optical Technology Division, National Institute of Standards and Technology
//** 100 Bureau Dr. Stop 8443; Gaithersburg, MD 20899-8443
//** Phone: (301) 975-2876; FAX: (301) 975-6991
//** Email: thomas.germer@nist.gov
//**
//** Version: 6.00 (February 2008)
//**
//******************************************************************************
#include "allrough.h"

using namespace std;

namespace SCATMECH {

    JonesMatrix 
    Correlated_Roughness_Stack_BRDF_Model::
    jones()
    {
        SETUP();
        
        JonesMatrix J=JonesZero(); // Added initialization for version 4.01
        for (int layer=0;layer<=stack.get_n();++layer) {
            set_this_layer(layer);
            J = J + Roughness_Stack_BRDF_Model::jones();
        }
        return J;
    }

    MuellerMatrix 
    Uncorrelated_Roughness_Stack_BRDF_Model::
    mueller()
    {
        SETUP();
        
        MuellerMatrix M=MuellerZero();  // Added initialization for version 4.01
        for (int layer=0;layer<=stack.get_n();++layer) {
            set_this_layer(layer);
            M = M + Roughness_Stack_BRDF_Model::jones();
        }
        return M;
    }

    DEFINE_MODEL(Correlated_Roughness_Stack_BRDF_Model,Roughness_Stack_BRDF_Model,
                "Correlated_Roughness_Stack_BRDF_Model",
                "All films in a stack equally rough and correlated");

    DEFINE_MODEL(Uncorrelated_Roughness_Stack_BRDF_Model,Roughness_Stack_BRDF_Model,
                "Uncorrelated_Roughness_Stack_BRDF_Model",
                "All films in a stack equally rough but uncorrelated");


} // namespace SCATMECH;