/*
 * MATLAB Compiler: 2.2
 * Date: Tue Jun 25 15:00:23 2002
 * Arguments: "-B" "macro_default" "-O" "all" "-O" "fold_scalar_mxarrays:on"
 * "-O" "fold_non_scalar_mxarrays:on" "-O" "optimize_integer_for_loops:on" "-O"
 * "array_indexing:on" "-O" "optimize_conditionals:on" "-t" "-W" "lib:Templib"
 * "-T" "link:exe" "-h" "fileinc.c" "mrank" 
 */
#include "magic.h"
#include "libmatlbm.h"
#include "meshgrid.h"

static mxChar _array1_[128] = { 'R', 'u', 'n', '-', 't', 'i', 'm', 'e', ' ',
                                'E', 'r', 'r', 'o', 'r', ':', ' ', 'F', 'i',
                                'l', 'e', ':', ' ', 'm', 'a', 'g', 'i', 'c',
                                ' ', 'L', 'i', 'n', 'e', ':', ' ', '1', ' ',
                                'C', 'o', 'l', 'u', 'm', 'n', ':', ' ', '1',
                                ' ', 'T', 'h', 'e', ' ', 'f', 'u', 'n', 'c',
                                't', 'i', 'o', 'n', ' ', '"', 'm', 'a', 'g',
                                'i', 'c', '"', ' ', 'w', 'a', 's', ' ', 'c',
                                'a', 'l', 'l', 'e', 'd', ' ', 'w', 'i', 't',
                                'h', ' ', 'm', 'o', 'r', 'e', ' ', 't', 'h',
                                'a', 'n', ' ', 't', 'h', 'e', ' ', 'd', 'e',
                                'c', 'l', 'a', 'r', 'e', 'd', ' ', 'n', 'u',
                                'm', 'b', 'e', 'r', ' ', 'o', 'f', ' ', 'o',
                                'u', 't', 'p', 'u', 't', 's', ' ', '(', '1',
                                ')', '.' };
static mxArray * _mxarray0_;

static mxChar _array3_[127] = { 'R', 'u', 'n', '-', 't', 'i', 'm', 'e', ' ',
                                'E', 'r', 'r', 'o', 'r', ':', ' ', 'F', 'i',
                                'l', 'e', ':', ' ', 'm', 'a', 'g', 'i', 'c',
                                ' ', 'L', 'i', 'n', 'e', ':', ' ', '1', ' ',
                                'C', 'o', 'l', 'u', 'm', 'n', ':', ' ', '1',
                                ' ', 'T', 'h', 'e', ' ', 'f', 'u', 'n', 'c',
                                't', 'i', 'o', 'n', ' ', '"', 'm', 'a', 'g',
                                'i', 'c', '"', ' ', 'w', 'a', 's', ' ', 'c',
                                'a', 'l', 'l', 'e', 'd', ' ', 'w', 'i', 't',
                                'h', ' ', 'm', 'o', 'r', 'e', ' ', 't', 'h',
                                'a', 'n', ' ', 't', 'h', 'e', ' ', 'd', 'e',
                                'c', 'l', 'a', 'r', 'e', 'd', ' ', 'n', 'u',
                                'm', 'b', 'e', 'r', ' ', 'o', 'f', ' ', 'i',
                                'n', 'p', 'u', 't', 's', ' ', '(', '1', ')',
                                '.' };
static mxArray * _mxarray2_;
static mxArray * _mxarray4_;
static mxArray * _mxarray5_;
static mxArray * _mxarray6_;
static mxArray * _mxarray7_;
static mxArray * _mxarray8_;

void InitializeModule_magic(void) {
    _mxarray0_ = mclInitializeString(128, _array1_);
    _mxarray2_ = mclInitializeString(127, _array3_);
    _mxarray4_ = mclInitializeDouble(2.0);
    _mxarray5_ = mclInitializeDouble(1.0);
    _mxarray6_ = mclInitializeDouble(3.0);
    _mxarray7_ = mclInitializeDouble(4.0);
    _mxarray8_ = mclInitializeDouble(0.0);
}

void TerminateModule_magic(void) {
    mxDestroyArray(_mxarray8_);
    mxDestroyArray(_mxarray7_);
    mxDestroyArray(_mxarray6_);
    mxDestroyArray(_mxarray5_);
    mxDestroyArray(_mxarray4_);
    mxDestroyArray(_mxarray2_);
    mxDestroyArray(_mxarray0_);
}

static mxArray * Mmagic(int nargout_, mxArray * n);

_mexLocalFunctionTable _local_function_table_magic
  = { 0, (mexFunctionTableEntry *)NULL };

/*
 * The function "mlfMagic" contains the normal interface for the "magic"
 * M-function from file "d:\MATLAB6p1\toolbox\matlab\elmat\magic.m" (lines
 * 1-43). This function processes any input arguments and passes them to the
 * implementation version of the function, appearing above.
 */
mxArray * mlfMagic(mxArray * n) {
    int nargout = 1;
    mxArray * M = mclGetUninitializedArray();
    mlfEnterNewContext(0, 1, n);
    M = Mmagic(nargout, n);
    mlfRestorePreviousContext(0, 1, n);
    return mlfReturnValue(M);
}

/*
 * The function "mlxMagic" contains the feval interface for the "magic"
 * M-function from file "d:\MATLAB6p1\toolbox\matlab\elmat\magic.m" (lines
 * 1-43). The feval function calls the implementation version of magic through
 * this function. This function processes any input arguments and passes them
 * to the implementation version of the function, appearing above.
 */
void mlxMagic(int nlhs, mxArray * plhs[], int nrhs, mxArray * prhs[]) {
    mxArray * mprhs[1];
    mxArray * mplhs[1];
    int i;
    if (nlhs > 1) {
        mlfError(_mxarray0_);
    }
    if (nrhs > 1) {
        mlfError(_mxarray2_);
    }
    for (i = 0; i < 1; ++i) {
        mplhs[i] = mclGetUninitializedArray();
    }
    for (i = 0; i < 1 && i < nrhs; ++i) {
        mprhs[i] = prhs[i];
    }
    for (; i < 1; ++i) {
        mprhs[i] = NULL;
    }
    mlfEnterNewContext(0, 1, mprhs[0]);
    mplhs[0] = Mmagic(nlhs, mprhs[0]);
    mlfRestorePreviousContext(0, 1, mprhs[0]);
    plhs[0] = mplhs[0];
}

/*
 * The function "Mmagic" is the implementation version of the "magic"
 * M-function from file "d:\MATLAB6p1\toolbox\matlab\elmat\magic.m" (lines
 * 1-43). It contains the actual compiled code for that M-function. It is a
 * static function and must only be called from one of the interface functions,
 * appearing below.
 */
/*
 * function M = magic(n)
 */
static mxArray * Mmagic(int nargout_, mxArray * n) {
    mexLocalFunctionTable save_local_function_table_
      = mclSetCurrentLocalFunctionTable(&_local_function_table_magic);
    mxArray * M = mclGetUninitializedArray();
    mxArray * j = mclGetUninitializedArray();
    mxArray * k = mclGetUninitializedArray();
    mxArray * i = mclGetUninitializedArray();
    mxArray * p = mclGetUninitializedArray();
    mxArray * K = mclGetUninitializedArray();
    mxArray * B = mclGetUninitializedArray();
    mxArray * A = mclGetUninitializedArray();
    mxArray * I = mclGetUninitializedArray();
    mxArray * J = mclGetUninitializedArray();
    mclCopyArray(&n);
    /*
     * %MAGIC  Magic square.
     * %   MAGIC(N) is an N-by-N matrix constructed from the integers
     * %   1 through N^2 with equal row, column, and diagonal sums.
     * %   Produces valid magic squares for all N > 0 except N = 2.
     * 
     * %   Copyright 1984-2001 The MathWorks, Inc. 
     * %   $Revision: 5.13 $  $Date: 2001/04/15 12:02:43 $
     * 
     * % Historically, MATLAB's magic was a built-in function.
     * % This M-file uses a new algorithm to generate the same matrices.
     * 
     * n = floor(real(double(n(1))));
     */
    mlfAssign(
      &n,
      mlfFloor(
        mclVe(
          mlfReal(
            mclVe(mlfDouble(mclVe(mclIntArrayRef1(mclVsa(n, "n"), 1))))))));
    /*
     * 
     * % Odd order.
     * if mod(n,2) == 1
     */
    if (mclEqBool(mclVe(mlfMod(mclVa(n, "n"), _mxarray4_)), _mxarray5_)) {
        /*
         * [J,I] = meshgrid(1:n);
         */
        mlfAssign(
          &J,
          mlfNMeshgrid(
            2,
            &I,
            NULL,
            mlfColon(_mxarray5_, mclVa(n, "n"), NULL),
            NULL,
            NULL));
        /*
         * A = mod(I+J-(n+3)/2,n);
         */
        mlfAssign(
          &A,
          mlfMod(
            mclMinus(
              mclPlus(mclVv(I, "I"), mclVv(J, "J")),
              mclMrdivide(mclPlus(mclVa(n, "n"), _mxarray6_), _mxarray4_)),
            mclVa(n, "n")));
        /*
         * B = mod(I+2*J-2,n);
         */
        mlfAssign(
          &B,
          mlfMod(
            mclMinus(
              mclPlus(mclVv(I, "I"), mclMtimes(_mxarray4_, mclVv(J, "J"))),
              _mxarray4_),
            mclVa(n, "n")));
        /*
         * M = n*A + B + 1;
         */
        mlfAssign(
          &M,
          mclPlus(
            mclPlus(mclMtimes(mclVa(n, "n"), mclVv(A, "A")), mclVv(B, "B")),
            _mxarray5_));
    /*
     * 
     * % Doubly even order.
     * elseif mod(n,4) == 0
     */
    } else if (mclEqBool(
                 mclVe(mlfMod(mclVa(n, "n"), _mxarray7_)), _mxarray8_)) {
        /*
         * [J,I] = meshgrid(1:n);
         */
        mlfAssign(
          &J,
          mlfNMeshgrid(
            2,
            &I,
            NULL,
            mlfColon(_mxarray5_, mclVa(n, "n"), NULL),
            NULL,
            NULL));
        /*
         * K = fix(mod(I,4)/2) == fix(mod(J,4)/2);
         */
        mlfAssign(
          &K,
          mclEq(
            mclVe(
              mlfFix(
                mclMrdivide(
                  mclVe(mlfMod(mclVv(I, "I"), _mxarray7_)), _mxarray4_))),
            mclVe(
              mlfFix(
                mclMrdivide(
                  mclVe(mlfMod(mclVv(J, "J"), _mxarray7_)), _mxarray4_)))));
        /*
         * M = reshape(1:n*n,n,n)';
         */
        mlfAssign(
          &M,
          mlfCtranspose(
            mclVe(
              mlfReshape(
                mlfColon(
                  _mxarray5_, mclMtimes(mclVa(n, "n"), mclVa(n, "n")), NULL),
                mclVa(n, "n"), mclVa(n, "n"), NULL))));
        /*
         * M(K) = n*n+1 - M(K);
         */
        mclArrayAssign1(
          &M,
          mclMinus(
            mclPlus(mclMtimes(mclVa(n, "n"), mclVa(n, "n")), _mxarray5_),
            mclVe(mclArrayRef1(mclVsv(M, "M"), mclVsv(K, "K")))),
          mclVsv(K, "K"));
    /*
     * 
     * % Singly even order.
     * else
     */
    } else {
        /*
         * p = n/2;
         */
        mlfAssign(&p, mclMrdivide(mclVa(n, "n"), _mxarray4_));
        /*
         * M = magic(p);
         */
        mlfAssign(&M, mlfMagic(mclVv(p, "p")));
        /*
         * M = [M M+2*p^2; M+3*p^2 M+p^2];
         */
        mlfAssign(
          &M,
          mlfVertcat(
            mlfHorzcat(
              mclVv(M, "M"),
              mclPlus(
                mclVv(M, "M"),
                mclMtimes(_mxarray4_, mclMpower(mclVv(p, "p"), _mxarray4_))),
              NULL),
            mlfHorzcat(
              mclPlus(
                mclVv(M, "M"),
                mclMtimes(_mxarray6_, mclMpower(mclVv(p, "p"), _mxarray4_))),
              mclPlus(mclVv(M, "M"), mclMpower(mclVv(p, "p"), _mxarray4_)),
              NULL),
            NULL));
        /*
         * if n == 2, break, end
         */
        if (mclEqBool(mclVa(n, "n"), _mxarray4_)) {
            goto return_;
        }
        /*
         * i = (1:p)';
         */
        mlfAssign(&i, mlfCtranspose(mlfColon(_mxarray5_, mclVv(p, "p"), NULL)));
        /*
         * k = (n-2)/4;
         */
        mlfAssign(
          &k, mclMrdivide(mclMinus(mclVa(n, "n"), _mxarray4_), _mxarray7_));
        /*
         * j = [1:k (n-k+2):n];
         */
        mlfAssign(
          &j,
          mlfHorzcat(
            mlfColon(_mxarray5_, mclVv(k, "k"), NULL),
            mlfColon(
              mclPlus(mclMinus(mclVa(n, "n"), mclVv(k, "k")), _mxarray4_),
              mclVa(n, "n"),
              NULL),
            NULL));
        /*
         * M([i; i+p],j) = M([i+p; i],j);
         */
        mclArrayAssign2(
          &M,
          mclArrayRef2(
            mclVsv(M, "M"),
            mlfVertcat(
              mclPlus(mclVv(i, "i"), mclVv(p, "p")), mclVv(i, "i"), NULL),
            mclVsv(j, "j")),
          mlfVertcat(
            mclVv(i, "i"), mclPlus(mclVv(i, "i"), mclVv(p, "p")), NULL),
          mclVsv(j, "j"));
        /*
         * i = k+1;
         */
        mlfAssign(&i, mclPlus(mclVv(k, "k"), _mxarray5_));
        /*
         * j = [1 i];
         */
        mlfAssign(&j, mlfHorzcat(_mxarray5_, mclVv(i, "i"), NULL));
        /*
         * M([i; i+p],j) = M([i+p; i],j);
         */
        mclArrayAssign2(
          &M,
          mclArrayRef2(
            mclVsv(M, "M"),
            mlfVertcat(
              mclPlus(mclVv(i, "i"), mclVv(p, "p")), mclVv(i, "i"), NULL),
            mclVsv(j, "j")),
          mlfVertcat(
            mclVv(i, "i"), mclPlus(mclVv(i, "i"), mclVv(p, "p")), NULL),
          mclVsv(j, "j"));
    /*
     * end
     */
    }
    return_:
    mclValidateOutput(M, 1, nargout_, "M", "magic");
    mxDestroyArray(J);
    mxDestroyArray(I);
    mxDestroyArray(A);
    mxDestroyArray(B);
    mxDestroyArray(K);
    mxDestroyArray(p);
    mxDestroyArray(i);
    mxDestroyArray(k);
    mxDestroyArray(j);
    mxDestroyArray(n);
    mclSetCurrentLocalFunctionTable(save_local_function_table_);
    return M;
}