From 995fc1b91f005c8ca430dd8f2bcca0d676907f06 Mon Sep 17 00:00:00 2001
From: HOSOKAWA Kenchi <hskwk@users.sourceforge.jp>
Date: Wed, 4 Apr 2012 00:22:53 +0900
Subject: [PATCH] 12

---
 dev5/psychlops/core/math/matrix.cs |  15 ++++
 test5/PsychlopsMain.cs             | 143 ++++++++++++++++++++++++++++++++++++-
 2 files changed, 157 insertions(+), 1 deletion(-)

diff --git a/dev5/psychlops/core/math/matrix.cs b/dev5/psychlops/core/math/matrix.cs
index de4a836..4d011e5 100644
--- a/dev5/psychlops/core/math/matrix.cs
+++ b/dev5/psychlops/core/math/matrix.cs
@@ -6,6 +6,11 @@ namespace Psychlops
 
 	public abstract class Matrix
 	{
+		public static Matrix gen(int row, int col)
+		{
+			return new Internal.MatrixImplementation(row, col);
+		}
+
 		public abstract double this[int row, int col]
 		{
 			get;
@@ -26,6 +31,16 @@ namespace Psychlops
 			return new Internal.MatrixExpression();
 		}
 		*/
+		public void each(Func<double, double> f)
+		{
+			for (int i = 1; i <= getRows(); i++)
+			{
+				for (int j = 1; j <= getCols(); j++)
+				{
+					this[i, j] = f(this[i, j]);
+				}
+			}
+		}
 
 		public abstract int rows { get; }
 		public abstract int cols { get; }
diff --git a/test5/PsychlopsMain.cs b/test5/PsychlopsMain.cs
index f1d7bb8..4242429 100644
--- a/test5/PsychlopsMain.cs
+++ b/test5/PsychlopsMain.cs
@@ -1,4 +1,143 @@
-//Psychlops Code Template
+// Motion with glass patterns.
+// Ross, J., Badcock, D. R., and Hayes, A. (2000)
+// Coherent global motion in the absence of coherent velocity signals.
+// Current Biology, 10, 679-682.
+///+ Prefix
+//// Include Psychlops Package
+using Psychlops;
+
+namespace PsychlopsSilverlightApp
+{
+
+	public class PsychlopsMain
+	{
+		Psychlops.Widgets.Slider distance, duration, stimulus_type;
+		///- Prefix
+
+		public void psychlops_main()
+		{
+			GlassPattern();
+		}
+
+		readonly int DOTNUM = 512;
+
+		///+ Stimulus drawing function
+		//// A function for stimulus drawing (main body)
+		void GlassPattern() {
+
+			///+ Preperation
+			//// Declare and initialize local variables
+			Canvas cnvs = new Canvas(400,400); //Prepare drawing window
+			double dotsize = 2;
+			double dot_lum  = 0.5;
+			double bg_lum    = 0.2;
+			double fieldsize = 300;
+			double[] orientation = new double[DOTNUM];
+			Rectangle dots = new Rectangle(dotsize, dotsize);
+			///+ prepare dots position
+			//// prepare dots position and paired orientation
+			Matrix positionmat = Matrix.gen(DOTNUM,2);
+			//Math.random(positionmat, -0.5*fieldsize, 0.5*fieldsize);
+			positionmat.each( (v) => Math.random(-0.5*fieldsize, 0.5*fieldsize) );
+
+			
+			for(int i=0; i<DOTNUM; i++){orientation[i] = Math.atan2(positionmat[i+1,2], positionmat[i+1,1]);}
+			///- prepare dots position
+
+
+			///+ Draw Gaussian Envelope
+			////Draw offscreen Gaussian Envelope
+			Image envelope = new Image(fieldsize+100, fieldsize+100);
+			envelope.clear(new Color(bg_lum)); //clear offscreen image
+			double _x, _y;
+			for(int i=0; i<fieldsize+100; i++){
+				_x=i-0.5*fieldsize-50;
+				for(int j=0; j<fieldsize+100; j++){
+					_y=j-0.5*fieldsize-50;
+					envelope.alpha(i,j,1.0-Math.exp(-((_x*_x+_y*_y)/ (2.0*Math.pow(fieldsize/6.0, 2.0) ))));
+				}
+			}
+			envelope.cache();
+			///- Draw Gaussian Envelope
+			///+ set Independent
+			//// set variables and value ranges for interaction
+			Interval rng = new Interval();
+			stimulus_type = new Psychlops.Widgets.Slider("Stimulus Type", 0.0 <= rng <= 2.0, 1, 1);
+			stimulus_type.setValue(0);
+			distance = new Psychlops.Widgets.Slider("distance between pairs", 1 <= rng <= 100, 1, 1);
+			distance.setValue(3.0);
+			duration = new Psychlops.Widgets.Slider("Refresh", 0.0 <= rng <= 10.0, 1, 1);
+			duration.setValue(10);
+			///- set Independent
+
+			int frame = 0;
+			while(!Keyboard.esc.pushed()) {
+				///+ position change
+				//// Re-randomize position with a desingated interval
+				if(frame > duration)
+				{
+					//Math.random(positionmat, -0.5 * fieldsize, 0.5 * fieldsize);
+					positionmat.each((v) => Math.random(-0.5 * fieldsize, 0.5 * fieldsize));
+					for (int i = 0; i < DOTNUM; i++) { orientation[i] = Math.atan2(positionmat[i + 1, 2], positionmat[i + 1, 1]); }
+					frame = 0;
+				}
+				///- position change
+				Display.clear(new Color(bg_lum)); //Clear window
+
+				//// Draw DOTNUM pairs of dots
+				for(int i=0; i<DOTNUM; i++){
+					///+ draw first dots
+
+					///+ switch stimulus type
+					//// draw the second dots to generate selected stimulu type.
+					switch((int)stimulus_type){
+						///+ case rotation
+						//// Rotation
+						case 0:
+						dots.centering().shift(positionmat[i+1,1],positionmat[i+1,2]); // move to center position
+						dots.shift(-0.5*distance*Math.cos(orientation[i]+0.5*Math.PI),-0.5*distance*Math.sin(orientation[i]+0.5*Math.PI));// move to the first position
+						dots.draw(new Color(dot_lum)); // draw first dot
+						dots.shift(distance*Math.cos(orientation[i]+0.5*Math.PI),distance*Math.sin(orientation[i]+0.5*Math.PI));// move to the second position
+						//cnvs.msg("Rotation", cnvs.getCenter().x, cnvs.getHeight() - 20);
+						break;
+						///- case rotation
+						///+ case diversion
+						////Diversion
+						case 1:
+						dots.centering().shift(positionmat[i+1,1],positionmat[i+1,2]);// move to center position
+						dots.shift(-0.5*distance*Math.cos(orientation[i]+0.5*Math.PI),-0.5*distance*Math.sin(orientation[i]+0.5*Math.PI));// move to the first position
+						dots.draw(new Color(dot_lum));// draw first dot
+						dots.shift(distance*Math.cos(orientation[i]),distance*Math.sin(orientation[i]));// move to the second position
+						break;
+						///- case diversion
+						///+ case spiral
+						//// Spiral
+						case 2:
+						dots.centering().shift(positionmat[i+1,1],positionmat[i+1,2]);// move to center position
+						dots.shift(-0.5*distance*Math.cos(orientation[i]+0.5*Math.PI),-0.5*distance*Math.sin(orientation[i]+0.5*Math.PI));// move to the first position
+						dots.draw(new Color(dot_lum));// draw first dot
+						dots.shift(distance*Math.cos(orientation[i]+0.25*Math.PI),distance*Math.sin(orientation[i]+0.25*Math.PI));// move to the second position
+						break;
+						///- case spiral
+					}
+					///- switch stimulus type
+
+					dots.draw(new Color(dot_lum)); // Draw dots at a desinated position.
+				}
+				///- Draw dots
+				envelope.centering().draw(); //Draw offscreen Gaussian onto screen
+				Display.flip(); // Flip frame buffers
+				frame++;
+			}
+		}
+	}
+}
+
+
+
+
+/*
+ * //Psychlops Code Template
 //    Please visit following web site to get sample codes.
 //    http://psychlops.sourceforge.jp/ja/?StartCode
 //    CodeDresser at following address is also available to view the code.
@@ -117,6 +256,8 @@ namespace PsychlopsSilverlightApp
 
 }
 
+*/
+
 /*
 		public void psychlops_main() {
 
-- 
2.11.0