From f618466c25d43f3bae9e40920273bf77de1e1149 Mon Sep 17 00:00:00 2001
From: leochanj105 <leochanj@live.unc.edu>
Date: Mon, 19 Oct 2020 23:09:30 -0400
Subject: initial sd-vbs

initial sd-vbs

add sd-vbs

sd-vbs
---
 .../src/matlabPyrTools/MEX/pointOp.c               | 126 +++++++++++++++++++++
 1 file changed, 126 insertions(+)
 create mode 100755 SD-VBS/benchmarks/texture_synthesis/src/matlabPyrTools/MEX/pointOp.c

(limited to 'SD-VBS/benchmarks/texture_synthesis/src/matlabPyrTools/MEX/pointOp.c')

diff --git a/SD-VBS/benchmarks/texture_synthesis/src/matlabPyrTools/MEX/pointOp.c b/SD-VBS/benchmarks/texture_synthesis/src/matlabPyrTools/MEX/pointOp.c
new file mode 100755
index 0000000..6ffcb45
--- /dev/null
+++ b/SD-VBS/benchmarks/texture_synthesis/src/matlabPyrTools/MEX/pointOp.c
@@ -0,0 +1,126 @@
+/* 
+RES = pointOp(IM, LUT, ORIGIN, INCREMENT, WARNINGS)
+  >>> See pointOp.m for documentation <<<
+  EPS, ported from OBVIUS, 7/96.
+*/
+
+#define V4_COMPAT
+#include <matrix.h>  /* Matlab matrices */
+#include <mex.h>
+
+#include <stddef.h>  /* NULL */
+
+#define notDblMtx(it) (!mxIsNumeric(it) || !mxIsDouble(it) || mxIsSparse(it) || mxIsComplex(it))
+
+void internal_pointop();
+
+void mexFunction(int nlhs,	     /* Num return vals on lhs */
+		 mxArray *plhs[],    /* Matrices on lhs      */
+		 int nrhs,	     /* Num args on rhs    */
+		 const mxArray *prhs[]     /* Matrices on rhs */
+		 )
+  {
+  double *image, *lut, *res;
+  double origin, increment;
+  int x_dim, y_dim, lx_dim, ly_dim;
+  int warnings = 1;
+  mxArray *arg;
+  double *mxMat;
+
+  if (nrhs < 4 ) mexErrMsgTxt("requres  at least 4 args.");
+
+  /* ARG 1: IMAGE  */
+  arg = prhs[0];
+  if notDblMtx(arg) mexErrMsgTxt("IMAGE arg must be a real non-sparse matrix.");
+  image = mxGetPr(arg);
+  x_dim = (int) mxGetM(arg); /* X is inner index! */
+  y_dim = (int) mxGetN(arg);
+
+  /* ARG 2: Lookup table */
+  arg = prhs[1];
+  if notDblMtx(arg) mexErrMsgTxt("LUT arg must be a real non-sparse matrix.");
+  lut = mxGetPr(arg);
+  lx_dim = (int) mxGetM(arg); /* X is inner index! */
+  ly_dim = (int) mxGetN(arg);
+  if ( (lx_dim != 1) && (ly_dim != 1) )
+    mexErrMsgTxt("Lookup table must be a row or column vector.");
+
+  /* ARG 3: ORIGIN */
+  arg = prhs[2];
+  if notDblMtx(arg) mexErrMsgTxt("ORIGIN arg must be a real scalar.");
+  if (mxGetM(arg) * mxGetN(arg) != 1)
+     mexErrMsgTxt("ORIGIN arg must be a real scalar.");
+  mxMat = mxGetPr(arg);
+  origin = *mxMat;
+
+  /* ARG 4: INCREMENT */
+  arg = prhs[3];
+  if notDblMtx(arg) mexErrMsgTxt("INCREMENT arg must be a real scalar.");
+  if (mxGetM(arg) * mxGetN(arg) != 1)
+     mexErrMsgTxt("INCREMENT arg must be a real scalar.");
+  mxMat = mxGetPr(arg);
+  increment = *mxMat;
+
+  /* ARG 5: WARNINGS */
+  if (nrhs>4)
+    {
+    arg = prhs[4];
+    if notDblMtx(arg) mexErrMsgTxt("WARINGS arg must be a real scalar.");
+    if (mxGetM(arg) * mxGetN(arg) != 1)
+      mexErrMsgTxt("WARNINGS arg must be a real scalar.");
+    mxMat = mxGetPr(arg);
+    warnings = (int) *mxMat;
+    }
+
+  plhs[0] = (mxArray *)  mxCreateDoubleMatrix(x_dim,y_dim,mxREAL);
+  if (plhs[0] == NULL) mexErrMsgTxt("Cannot allocate result matrix");
+  res = mxGetPr(plhs[0]);
+      
+  internal_pointop(image, res, x_dim*y_dim, lut, lx_dim*ly_dim, 
+		   origin, increment, warnings);
+  return;
+  }      
+
+
+/* Use linear interpolation on a lookup table.
+   Taken from OBVIUS.  EPS, Spring, 1987.
+ */
+void internal_pointop (im, res, size, lut, lutsize, origin, increment, warnings)
+  register double *im, *res, *lut;
+  register double origin, increment; 
+  register int size, lutsize, warnings;
+  {
+  register int i, index;
+  register double pos;
+  register int l_unwarned = warnings;
+  register int r_unwarned = warnings;
+
+  lutsize = lutsize - 2;	/* Maximum index value */
+  if (increment > 0)
+    for (i=0; i<size; i++)
+	{
+	pos = (im[i] - origin) / increment;
+	index = (int) pos;   /* Floor */
+	if (index < 0)
+	    {
+	    index = 0;
+	    if (l_unwarned)
+		{
+		mexPrintf("Warning: Extrapolating to left of lookup table...\n");
+		l_unwarned = 0;
+		}
+	    }
+	else if (index > lutsize)
+	    {
+	    index = lutsize;
+	    if (r_unwarned)
+		{
+		mexPrintf("Warning: Extrapolating to right of lookup table...\n");
+		r_unwarned = 0;
+		}
+	    }
+	res[i] = lut[index] + (lut[index+1] - lut[index]) * (pos - index);
+	}
+  else
+    for (i=0; i<size; i++) res[i] = *lut;
+  }
-- 
cgit v1.2.2