wolffd@0: function A = softmax_activation( softmax_group)
wolffd@0: A = exp(softmax_group);
wolffd@0: A = A./repmat(sum(A,2),1,size(A,2));
wolffd@0: for i=1:size(A,1)
wolffd@0:     A(i,:) = cumsum(A(i,:))>rand();
wolffd@0:     m = max(A(i,:));
wolffd@0:     inx = find(A(i,:)==m,1,'first');
wolffd@0:     A(i,:) = 0;
wolffd@0:     A(i,inx) = 1;
wolffd@0: end
wolffd@0: end
wolffd@0: