Mercurial > hg > ishara

diff general/numerical/scalar/safe_exp.m @ 17:cb249ba61e9e

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Moved from parent directory.
author samer
date Thu, 17 Jan 2013 13:21:46 +0000
parents
children
line wrap: on
line diff
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/general/numerical/scalar/safe_exp.m	Thu Jan 17 13:21:46 2013 +0000
@@ -0,0 +1,11 @@
+function [y,k]=safe_exp(x,l)
+% safe_exp - high dynamic range exponential
+%
+% safe_exp :: X:[[N,M]] -> [[N,M]], [[1,M]].
+%
+% returns y and k such that exp(X) = y * exp(k) and
+% maximum value in y is 1.
+
+k=max(x,[],1);
+y=exp(x-repmat(k,size(x,1),1));
+