/*=================================================================
 *
 * TRIANGLE.CPP
 *
 * Calculate the points in Sierpinski's triangle, an elementary
 * fractal. 
 *
 * This function calls a 
 * using MATLAB Compiler.
 *
 * Copyright 1984-2004 The MathWorks, Inc.
 *
 *=================================================================*/

// Include the MCR header file and the library specific header file 
// as generated by MATLAB Compiler
#include "libtriangle.h"

void usage(const char *name)
{
    std::cout << "Usage: %s [number of points]" << std::endl; 
    exit (-1);
}

int main(int ac, char *av[])
{
    // Default number of iterations
    int num_points = 7500;

    // Validate the number of inputs
    if (ac < 1 || ac > 2)
    {
	fprintf(stderr, "Expecting 0 or 1 input(s). Found %d\n", ac);
	usage(av[0]);
    }

    // If we have the right number of inputs (1), try to convert the input
    // string to an integer.
    //
    if (ac == 2)
        num_points = atoi(av[1]);

    // Type check on input argument -- atoi() will fail if the input is
    // not an integer.
    //
    if (num_points == 0)
    {
	std::cerr << "First argument must be an integer." << std::endl;
	usage(av[0]);
    }

    // Call the library intialization routine and make sure that the
    // library was initialized properly
    //
    if (!mclInitializeApplication(NULL,0) || !libtriangleInitialize())
    {
        std::cerr << "could not initialize the triangle library properly" << std::endl;
        return -1;
    }

    // Create the input data
    // Input parameters:
    //
    // iterations: Number of points to draw in the triangle
    // draw: If true, draw the triangle in a figure window before returning.
    //
    mwArray iterations((mxDouble)num_points);
    mwArray draw(1.0);
    // Create the output variables
    // The Sierpinski function returns the X and Y coordinates of the points
    // forming the pattern in the triangle.
    mwArray x;
    mwArray y;

    // Call the library function
    sierpinski(2, x, y, iterations, draw);

    // Display the return value of the library function
    std::cout << "Calculated "  << x.NumberOfElements() << " points" << std::endl;

    // Block until user dismisses the figure
    mclWaitForFiguresToDie(NULL);

    // Call the library termination routine
    libtriangleTerminate();
    // Shut down all MCR instances
    mclTerminateApplication();

    /* Success */
    return 0;
}