Add a useful package (from Source forge) for octave

[CreaPhase.git] / octave_packages / secs2d-0.0.8 / Utilities / Utemplogm.m

diff --git a/octave_packages/secs2d-0.0.8/Utilities/Utemplogm.m b/octave_packages/secs2d-0.0.8/Utilities/Utemplogm.m
new file mode 100644 (file)
index 0000000..06a1e4b
--- /dev/null
+++ b/octave_packages/secs2d-0.0.8/Utilities/Utemplogm.m
@@ -0,0 +1,8 @@
+function T = Utemplogm(t1,t2)
+
+T = zeros(size(t2));
+
+sing     = abs(t2-t1)< 100*eps ;
+T(sing)  = (t2(sing)+t1(sing))/2;
+T(~sing) = (t2(~sing)-t1(~sing))./log(t2(~sing)./t1(~sing));
+