# HG changeset patch
# User samer
# Date 1328543708 0
# Node ID 0f16c069423bec2eb36bb0cbfc8b6db18539d757
# Parent  80e9103b057f891d5b39723e9fcf9dd833f1476f
Added randomise.m to provide random init; also added minimum radius to mt_get_transmat_near

diff -r 80e9103b057f -r 0f16c069423b mt_get_transmat_near.m
--- a/mt_get_transmat_near.m	Sun Feb 05 18:20:17 2012 +0000
+++ b/mt_get_transmat_near.m	Mon Feb 06 15:55:08 2012 +0000
@@ -14,17 +14,18 @@
 % The selected transition matrix will be shown using mt_show_transmat,
 
 
-function [T,P,I]=mt_get_transmat_at(Sys,Id,K,X,Y)
+function [T,P,I]=mt_get_transmat_near(Sys,Id,K,X,Y,MinRad)
 	logK=log(K);
 	normpos = Sys.map([X;Y])';
 	target = normpos*logK;
 	tmats = Sys.transmats{K};
 	info = Sys.info{K};
 	L = size(info,1);
+	if nargin<6, MinRad=0.01; end
 
 	% distance from target
 	d2 = sum((info(:,1:2) - repmat(target,L,1)).^2,2);
-	Mask = d2<=max(min(d2)*2,0.01);
+	Mask = d2<=max(min(d2)*2,MinRad);
 
 	M = find(Mask);
 %	J = M(argmax(info(M,3)));
diff -r 80e9103b057f -r 0f16c069423b private/randomise.m
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/private/randomise.m	Mon Feb 06 15:55:08 2012 +0000
@@ -0,0 +1,11 @@
+function s=randomise
+% randomise - get entropy from time and reset random number generators
+
+
+m=[5000,700000,1000000,1000000,1000000,10000000];
+rand('state',sum(m.*clock));
+randn('state',sum(m.*clock));
+if nargout>0,
+	s=getrndstate;
+end
+