diff flys-artifacts/src/main/java/org/dive4elements/river/artifacts/math/fitting/Exp.java @ 5831:bd047b71ab37

Repaired internal references
author Sascha L. Teichmann <teichmann@intevation.de>
date Thu, 25 Apr 2013 12:06:39 +0200
parents flys-artifacts/src/main/java/de/intevation/flys/artifacts/math/fitting/Exp.java@1c00ed73ccf4
children
line wrap: on
line diff
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/flys-artifacts/src/main/java/org/dive4elements/river/artifacts/math/fitting/Exp.java	Thu Apr 25 12:06:39 2013 +0200
@@ -0,0 +1,62 @@
+package org.dive4elements.river.artifacts.math.fitting;
+
+public class Exp
+extends      Function
+{
+    public static final Derivative DERIVATIVE =
+        new Derivative("W'(Q) = a^Q*log(a)*m") {
+
+        @Override
+        public org.dive4elements.river.artifacts.math.Function
+            instantiate(double [] parameters)
+        {
+                  double m     = parameters[0];
+            final double a     = parameters[1];
+            final double logam = Math.log(a)*m;
+
+            return new org.dive4elements.river.artifacts.math.Function() {
+                @Override
+                public double value(double Q) {
+                    return Math.pow(a, Q)*logam;
+                }
+            };
+        }
+    };
+
+    public static final Function INSTANCE = new Exp();
+
+    public Exp() {
+        super(
+            "exp",
+            "W(Q) = m * a^Q + b",
+            new String [] { "m", "a", "b" });
+    }
+
+    @Override
+    public double value(double x, double [] parameters) {
+        return parameters[0]*Math.pow(parameters[1], x) + parameters[2];
+    }
+
+    @Override
+    public double [] gradient(double Q, double [] parameters) {
+        double m = parameters[0];
+        double a = parameters[1];
+        double b = parameters[2];
+        return new double [] {
+            Math.pow(a, Q),
+            Math.pow(a, Q-1d)*m*Q,
+            1d
+        };
+    }
+
+    @Override
+    public Derivative getDerivative() {
+        return DERIVATIVE;
+    }
+
+    @Override
+    public Function getInverse() {
+        return InvExp.INSTANCE;
+    }
+}
+// vim:set ts=4 sw=4 si et sta sts=4 fenc=utf8 :

http://dive4elements.wald.intevation.org