nbgrid.c

nnToolKit 神经网络工具包是基于 MATLAB 神经网络工具箱自行开发的一组神经网络算法函数库

/*
 * MATLAB Compiler: 3.0
 * Date: Sun May 13 16:47:40 2007
 * 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" "-M" "-silentsetup" "-d"
 * "d:/MATLAB6p5/work/nnToolKit/src" "-B" "csglcom:nnToolKit,nnToolKit,2.0"
 * "-B" "sgl" "-m" "-W" "main" "-L" "C" "-t" "-T" "link:exe" "-h"
 * "libmmfile.mlib" "-W" "mainhg" "libmwsglm.mlib" "-t" "-W"
 * "comhg:nnToolKit,nnToolKit,2.0" "-T" "link:lib" "-h" "libmmfile.mlib" "-i"
 * "-i" "D:/MATLAB6p5/work/nnToolKit/lmnet/LmSimu.m"
 * "D:/MATLAB6p5/work/nnToolKit/lmnet/LmTrain.m"
 * "D:/MATLAB6p5/work/nnToolKit/sofm/SofmSimu.m"
 * "D:/MATLAB6p5/work/nnToolKit/sofm/SofmTrain.m" 
 */
#include "nbgrid.h"
#include "libmatlbm.h"
#include "libmmfile.h"
#include "nnmaxr.h"
#include "nntobsf.h"

static mxChar _array1_[6] = { 'n', 'b', 'g', 'r', 'i', 'd' };
static mxArray * _mxarray0_;

static mxChar _array3_[43] = { 'U', 's', 'e', ' ', 'N', 'E', 'W', 'S', 'O',
                               'M', ' ', 't', 'o', ' ', 'c', 'r', 'e', 'a',
                               't', 'e', ' ', 'a', ' ', 's', 'e', 'l', 'f',
                               '-', 'o', 'r', 'g', 'a', 'n', 'i', 'z', 'i',
                               'n', 'g', ' ', 'm', 'a', 'p', '.' };
static mxArray * _mxarray2_;

static mxChar _array5_[21] = { 'N', 'o', 't', ' ', 'e', 'n', 'o',
                               'u', 'g', 'h', ' ', 'a', 'r', 'g',
                               'u', 'm', 'e', 'n', 't', 's', '.' };
static mxArray * _mxarray4_;
static mxArray * _mxarray6_;
static mxArray * _mxarray7_;

static mxChar _array9_[3] = { 'd', '%', 'g' };
static mxArray * _mxarray8_;

void InitializeModule_nbgrid(void) {
    _mxarray0_ = mclInitializeString(6, _array1_);
    _mxarray2_ = mclInitializeString(43, _array3_);
    _mxarray4_ = mclInitializeString(21, _array5_);
    _mxarray6_ = mclInitializeDouble(1.0);
    _mxarray7_ = mclInitializeDoubleVector(0, 0, (double *)NULL);
    _mxarray8_ = mclInitializeString(3, _array9_);
}

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

static mxArray * Mnbgrid(int nargout_,
                         mxArray * d1,
                         mxArray * d2,
                         mxArray * d3,
                         mxArray * d4,
                         mxArray * d5);

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

/*
 * The function "mlfNbgrid" contains the normal interface for the "nbgrid"
 * M-function from file "d:\matlab6p5\toolbox\nnet\nnobsolete\nbgrid.m" (lines
 * 1-48). This function processes any input arguments and passes them to the
 * implementation version of the function, appearing above.
 */
mxArray * mlfNbgrid(mxArray * d1,
                    mxArray * d2,
                    mxArray * d3,
                    mxArray * d4,
                    mxArray * d5) {
    int nargout = 1;
    mxArray * m = NULL;
    mlfEnterNewContext(0, 5, d1, d2, d3, d4, d5);
    m = Mnbgrid(nargout, d1, d2, d3, d4, d5);
    mlfRestorePreviousContext(0, 5, d1, d2, d3, d4, d5);
    return mlfReturnValue(m);
}

/*
 * The function "mlxNbgrid" contains the feval interface for the "nbgrid"
 * M-function from file "d:\matlab6p5\toolbox\nnet\nnobsolete\nbgrid.m" (lines
 * 1-48). The feval function calls the implementation version of nbgrid through
 * this function. This function processes any input arguments and passes them
 * to the implementation version of the function, appearing above.
 */
void mlxNbgrid(int nlhs, mxArray * plhs[], int nrhs, mxArray * prhs[]) {
    mxArray * mprhs[5];
    mxArray * mplhs[1];
    int i;
    if (nlhs > 1) {
        mlfError(
          mxCreateString(
            "Run-time Error: File: nbgrid Line: 1 Column: "
            "1 The function \"nbgrid\" was called with mor"
            "e than the declared number of outputs (1)."),
          NULL);
    }
    if (nrhs > 5) {
        mlfError(
          mxCreateString(
            "Run-time Error: File: nbgrid Line: 1 Column: "
            "1 The function \"nbgrid\" was called with mor"
            "e than the declared number of inputs (5)."),
          NULL);
    }
    for (i = 0; i < 1; ++i) {
        mplhs[i] = NULL;
    }
    for (i = 0; i < 5 && i < nrhs; ++i) {
        mprhs[i] = prhs[i];
    }
    for (; i < 5; ++i) {
        mprhs[i] = NULL;
    }
    mlfEnterNewContext(0, 5, mprhs[0], mprhs[1], mprhs[2], mprhs[3], mprhs[4]);
    mplhs[0] = Mnbgrid(nlhs, mprhs[0], mprhs[1], mprhs[2], mprhs[3], mprhs[4]);
    mlfRestorePreviousContext(
      0, 5, mprhs[0], mprhs[1], mprhs[2], mprhs[3], mprhs[4]);
    plhs[0] = mplhs[0];
}

/*
 * The function "Mnbgrid" is the implementation version of the "nbgrid"
 * M-function from file "d:\matlab6p5\toolbox\nnet\nnobsolete\nbgrid.m" (lines
 * 1-48). 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=nbgrid(d1,d2,d3,d4,d5)
 */
static mxArray * Mnbgrid(int nargout_,
                         mxArray * d1,
                         mxArray * d2,
                         mxArray * d3,
                         mxArray * d4,
                         mxArray * d5) {
    mexLocalFunctionTable save_local_function_table_
      = mclSetCurrentLocalFunctionTable(&_local_function_table_nbgrid);
    int nargin_ = mclNargin(5, d1, d2, d3, d4, d5, NULL);
    mxArray * m = NULL;
    mxArray * distij = NULL;
    mxArray * j = NULL;
    mxArray * y = NULL;
    mxArray * x = NULL;
    mxArray * di = NULL;
    mxArray * i = NULL;
    mxArray * S = NULL;
    mxArray * coord = NULL;
    mxArray * n = NULL;
    mxArray * ans = NULL;
    mclCopyArray(&d1);
    mclCopyArray(&d2);
    mclCopyArray(&d3);
    mclCopyArray(&d4);
    mclCopyArray(&d5);
    /*
     * %NBGRID Neighborhood matrix using grid distance.
     * %  
     * %  This function is obselete.
     * %  Use NEWSOM to create a self-organizing map.
     * 
     * nntobsf('nbgrid','Use NEWSOM to create a self-organizing map.')
     */
    mlfNntobsf(_mxarray0_, _mxarray2_, NULL);
    /*
     * 
     * %  NBGRID(D1,D2,...,DN)
     * %    Di - Size of layers ith dimension.
     * %  Returns neighborhood matrix M for (D1*D2*...*DN)
     * %    neurons arranged in an N dimensional grid, where
     * %    each element (i,j) is equal to the grid distance
     * %    between neurons i and j.
     * %  
     * %  The N-neighborhood matrix can be found by: M <= N.
     * %  
     * %  EXAMPLE: m = nbgrid(2,3)
     * %  
     * %  See also NBMAN, NBDIST.
     * 
     * % Mark Beale, 1-31-92
     * % Revised 12-15-93, MB
     * % Copyright 1992-2002 The MathWorks, Inc.
     * % $Revision: 1.11 $  $Date: 2002/03/25 16:53:57 $
     * 
     * if nargin < 1, error('Not enough arguments.'); end
     */
    if (nargin_ < 1) {
        mlfError(_mxarray4_, NULL);
    }
    /*
     * if nargin == 1, n = 1; end
     */
    if (nargin_ == 1) {
        mlfAssign(&n, _mxarray6_);
    }
    /*
     * 
     * % COORDINATES
     * coord = (1:d1)';
     */
    mlfAssign(
      &coord, mlfCtranspose(mlfColon(_mxarray6_, mclVa(d1, "d1"), NULL)));
    /*
     * S = d1;
     */
    mlfAssign(&S, mclVa(d1, "d1"));
    {
        /*
         * for i=2:nargin
         */
        int v_ = mclForIntStart(2);
        int e_ = nargin_;
        if (v_ > e_) {
            mlfAssign(&i, _mxarray7_);
        } else {
            /*
             * di = eval(sprintf('d%g',i));
             * [x,y] = meshgrid(1:S,1:di);
             * S = S*di;
             * coord = [reshape(coord(x),S,1) reshape(y,S,1)];
             * end
             */
            for (; ; ) {
                mlfAssign(
                  &di,
                  mlfEval(
                    mclValueVarargout(),
                    mlfSprintf(NULL, _mxarray8_, mlfScalar(v_), NULL),
                    NULL));
                mlfAssign(
                  &x,
                  mlfNMeshgrid(
                    2,
                    &y,
                    NULL,
                    mlfColon(_mxarray6_, mclVv(S, "S"), NULL),
                    mlfColon(_mxarray6_, mclVv(di, "di"), NULL),
                    NULL));
                mlfAssign(&S, mclMtimes(mclVv(S, "S"), mclVv(di, "di")));
                mlfAssign(
                  &coord,
                  mlfHorzcat(
                    mlfReshape(
                      mclArrayRef1(mclVv(coord, "coord"), mclVv(x, "x")),
                      mclVv(S, "S"), _mxarray6_, NULL),
                    mlfReshape(mclVv(y, "y"), mclVv(S, "S"), _mxarray6_, NULL),
                    NULL));
                if (v_ == e_) {
                    break;
                }
                ++v_;
            }
            mlfAssign(&i, mlfScalar(v_));
        }
    }
    /*
     * 
     * % DISTANCES
     * m = zeros(S,S);
     */
    mlfAssign(&m, mlfZeros(mclVv(S, "S"), mclVv(S, "S"), NULL));
    /*
     * for i=1:(S-1)
     */
    {
        int v_ = mclForIntStart(1);
        int e_ = mclForIntEnd(mclMinus(mclVv(S, "S"), _mxarray6_));
        if (v_ > e_) {
            mlfAssign(&i, _mxarray7_);
        } else {
            /*
             * j=(i+1):S;
             * distij = nnmaxr((abs(ones(length(j),1)*coord(i,:)-coord(j,:))));
             * m(i,j) = distij';
             * m(j,i) = distij;
             * end
             */
            for (; ; ) {
                mlfAssign(&j, mlfColon(mlfScalar(v_ + 1), mclVv(S, "S"), NULL));
                mlfAssign(
                  &distij,
                  mlfNnmaxr(
                    mlfAbs(
                      mclMinus(
                        mclMtimes(
                          mlfOnes(
                            mlfScalar(mclLengthInt(mclVv(j, "j"))),
                            _mxarray6_,
                            NULL),
                          mclArrayRef2(
                            mclVv(coord, "coord"),
                            mlfScalar(v_),
                            mlfCreateColonIndex())),
                        mclArrayRef2(
                          mclVv(coord, "coord"),
                          mclVv(j, "j"),
                          mlfCreateColonIndex())))));
                mclArrayAssign2(
                  &m,
                  mlfCtranspose(mclVv(distij, "distij")),
                  mlfScalar(v_),
                  mclVv(j, "j"));
                mclArrayAssign2(
                  &m, mclVv(distij, "distij"), mclVv(j, "j"), mlfScalar(v_));
                if (v_ == e_) {
                    break;
                }
                ++v_;
            }
            mlfAssign(&i, mlfScalar(v_));
        }
    }
    mclValidateOutput(m, 1, nargout_, "m", "nbgrid");
    mxDestroyArray(ans);
    mxDestroyArray(n);
    mxDestroyArray(coord);
    mxDestroyArray(S);
    mxDestroyArray(i);
    mxDestroyArray(di);
    mxDestroyArray(x);
    mxDestroyArray(y);
    mxDestroyArray(j);
    mxDestroyArray(distij);
    mxDestroyArray(d5);
    mxDestroyArray(d4);
    mxDestroyArray(d3);
    mxDestroyArray(d2);
    mxDestroyArray(d1);
    mclSetCurrentLocalFunctionTable(save_local_function_table_);
    return m;
}