divide_by_table_ndxd.c

来自「麻省理工学院的人工智能工具箱,很珍贵,希望对大家有用!」· C语言 代码 · 共 59 行

/* C mex version for divide_by_table_ndxD.m in potential/Tables directory  */

#include "mex.h"

void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[]){
	double  *pb, *ps, value;
	int     i, j, I, J, N, pointer, temp;
	int     *ndx;
	char    *used;

	plhs[0] = mxDuplicateArray(prhs[0]);
	N = mxGetNumberOfElements(plhs[0]);
	pb = mxGetPr(plhs[0]);
	ps = mxGetPr(prhs[1]);
	ndx = mxGetData(prhs[2]);
	J = mxGetNumberOfElements(prhs[2]);
	I = N / J;

	if(J == 1){
		for(i=0; i<N; i++){
			value = *ps++;
			if(value != 0) *pb /= value;
			pb++;
		}
		return;
	}

	if(I == 1){
		value = *ps;
		if(value == 0) value = 1;
		for(i=0; i<N; i++){
			*pb++ /= value;
		}
		return;
	}

	used = (char *)malloc(N * sizeof(char));
	for(i=0; i<N; i++){
		used[i] = 0;
	}
	pointer = 0;
	for(i=0; i<I; i++){
		while(used[pointer]){
			pointer++;
		}
		temp = pointer;
		used[pointer] = 1;
		for(j=0; j<J; j++){
			temp = pointer + ndx[j];
			value = *ps;
			if(value != 0) pb[temp] /= value;
			used[temp] = 1;
		}
		ps++;
		pointer++;
	}
	free(used);
}