1 files changed, 45 insertions, 0 deletions

diff --git a/SD-VBS/benchmarks/texture_synthesis/src/matlab/rcosFn.m b/SD-VBS/benchmarks/texture_synthesis/src/matlab/rcosFn.m
new file mode 100755
index 0000000..5dac344
--- /dev/null
+++ b/SD-VBS/benchmarks/texture_synthesis/src/matlab/rcosFn.m
@@ -0,0 +1,45 @@
+% [X, Y] = rcosFn(WIDTH, POSITION, VALUES)
+%
+% Return a lookup table (suitable for use by INTERP1) 
+% containing a "raised cosine" soft threshold function:
+% 
+%    Y =  VALUES(1) + (VALUES(2)-VALUES(1)) *
+%              cos^2( PI/2 * (X - POSITION + WIDTH)/WIDTH )
+%
+% WIDTH is the width of the region over which the transition occurs
+% (default = 1). POSITION is the location of the center of the
+% threshold (default = 0).  VALUES (default = [0,1]) specifies the
+% values to the left and right of the transition.
+
+% Eero Simoncelli, 7/96.
+
+function [X, Y] = rcosFn(width,position,values)
+
+%------------------------------------------------------------
+% OPTIONAL ARGS:
+
+if (exist('width') ~= 1)
+  width = 1;
+end
+
+if (exist('position') ~= 1)
+  position = 0;
+end
+
+if (exist('values') ~= 1)
+  values = [0,1];
+end
+
+%------------------------------------------------------------
+
+sz = 256;  %% arbitrary!
+
+X    = pi * [-sz-1:1] / (2*sz);
+
+Y = values(1) + (values(2)-values(1)) * cos(X).^2;
+
+%    Make sure end values are repeated, for extrapolation...
+Y(1) = Y(2);
+Y(sz+3) = Y(sz+2);
+
+X = position + (2*width/pi) * (X + pi/4);