1 files changed, 77 insertions, 0 deletions
diff --git a/fuzzylite/src/defuzzifier/SmallestOfMaximum.cpp b/fuzzylite/src/defuzzifier/SmallestOfMaximum.cpp
new file mode 100644
index 0000000..1e67395
--- /dev/null
+++ b/fuzzylite/src/defuzzifier/SmallestOfMaximum.cpp
@@ -0,0 +1,77 @@
+/*
+ Author: Juan Rada-Vilela, Ph.D.
+ Copyright (C) 2010-2014 FuzzyLite Limited
+ All rights reserved
+
+ This file is part of fuzzylite.
+
+ fuzzylite is free software: you can redistribute it and/or modify it under
+ the terms of the GNU Lesser General Public License as published by the Free
+ Software Foundation, either version 3 of the License, or (at your option)
+ any later version.
+
+ fuzzylite is distributed in the hope that it will be useful, but WITHOUT
+ ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+ FITNESS FOR A PARTICULAR PURPOSE.  See the GNU Lesser General Public License
+ for more details.
+
+ You should have received a copy of the GNU Lesser General Public License
+ along with fuzzylite.  If not, see <http://www.gnu.org/licenses/>.
+
+ fuzzylite™ is a trademark of FuzzyLite Limited.
+
+ */
+
+#include "fl/defuzzifier/SmallestOfMaximum.h"
+
+#include "fl/Exception.h"
+#include "fl/term/Term.h"
+
+
+namespace fl {
+
+    SmallestOfMaximum::SmallestOfMaximum(int resolution)
+    : IntegralDefuzzifier(resolution) {
+    }
+
+    SmallestOfMaximum::~SmallestOfMaximum() {
+    }
+
+    std::string SmallestOfMaximum::className() const {
+        return "SmallestOfMaximum";
+    }
+
+    scalar SmallestOfMaximum::defuzzify(const Term* term, scalar minimum, scalar maximum) const {
+        if (not fl::Op::isFinite(minimum + maximum)) {
+            return fl::nan;
+        }
+        if (maximum - minimum > _resolution) {
+            FL_DBG("[accuracy warning] the resolution <" << _resolution << "> "
+                    "is smaller than the range <" << minimum << ", " << maximum << ">. In order to "
+                    "improve the accuracy, the resolution should be at least equal to the range.");
+        }
+        scalar dx = (maximum - minimum) / _resolution;
+        scalar x, y;
+        scalar ymax = -1.0, xsmallest = minimum;
+        for (int i = 0; i < _resolution; ++i) {
+            x = minimum + (i + 0.5) * dx;
+            y = term->membership(x);
+
+            if (Op::isGt(y, ymax)) {
+                xsmallest = x;
+                ymax = y;
+            }
+        }
+        return xsmallest;
+    }
+
+    SmallestOfMaximum* SmallestOfMaximum::clone() const {
+        return new SmallestOfMaximum(*this);
+    }
+
+    Defuzzifier* SmallestOfMaximum::constructor() {
+        return new SmallestOfMaximum;
+    }
+
+}
+