#include <iostream>
#include <iomanip>
#include <cstdio>
#include "libnormal.h"    //This is produced later by mcc (MATLAB compiler) that translates your .m files to a library with its header.

void matlab_init()
{
    //1. Initialize MCR and library function of Matlab 7 R14.
    //
    /* Call the mclInitializeApplication routine. Make sure that the application
     * was initialized properly by checking the return status. This initialization
     * has to be done before calling any MATLAB API's or MATLAB Compiler generated
     * shared library functions.
     */
    if( !mclInitializeApplication(NULL,0) )
    {
        fprintf(stderr, "Could not initialize the application.\n");
        exit(1);
    }
    /* Call the library intialization routine and make sure that the
     * library was initialized properly.
     */

    if (!libnormalInitialize()){
        fprintf(stderr,"Could not initialize the library.\n");
        exit(1);
    }
}
int compute_normcdf(double* data_in, double* data_out, double mu, double sigma, int num)
{
    try
    {
        // Create input data
		double* mus = new double[num];
		double* sigmas = new double[num];
		for(int i=0; i<num; i++)
		{
			mus[i] = mu;
			sigmas[i] = sigma;
		}
        mwArray array_in(num, 1, mxDOUBLE_CLASS, mxREAL);
		mwArray array_mus(num, 1, mxDOUBLE_CLASS, mxREAL);
		mwArray array_sigmas(num, 1, mxDOUBLE_CLASS, mxREAL);
        array_in.SetData(data_in, num);
		array_mus.SetData(mus, num);
		array_sigmas.SetData(sigmas, num);
		delete [] mus;
		delete [] sigmas;

        ////Show inputs.
        //std::cout << "Input matrix 1 is:" << std::endl;
        //std::cout << array_in << std::endl;

        //Create output array
        mwArray array_out;
        ////Call the library function
        matlab_normcdf(1, array_out, array_in, array_mus, array_sigmas);

		array_out.GetData(data_out, num);
        //Display the return value of the library function
        //std::cout << "The value of added matrix is:" << std::endl;
        //std::cout << array_out << std::endl;

    }
    catch (const mwException& e)
    {
        std::cerr << e.what() << std::endl;
        return -1;
    }
    catch (...) 
    {
        std::cerr << "Unexpected error thrown" << std::endl;
        return -1;
    }

    return 0;
}
void matlab_terminate()
{
    //3. Terminate MCR and library function and free resources.
    //
    libnormalTerminate();
    mclTerminateApplication();

}