Image Component Library (ICL)
|
00001 /******************************************************************** 00002 ** Image Component Library (ICL) ** 00003 ** ** 00004 ** Copyright (C) 2006-2013 CITEC, University of Bielefeld ** 00005 ** Neuroinformatics Group ** 00006 ** Website: www.iclcv.org and ** 00007 ** http://opensource.cit-ec.de/projects/icl ** 00008 ** ** 00009 ** File : ICLCore/src/ICLCore/Parable.h ** 00010 ** Module : ICLCore ** 00011 ** Authors: Christof Elbrechter ** 00012 ** ** 00013 ** ** 00014 ** GNU LESSER GENERAL PUBLIC LICENSE ** 00015 ** This file may be used under the terms of the GNU Lesser General ** 00016 ** Public License version 3.0 as published by the ** 00017 ** ** 00018 ** Free Software Foundation and appearing in the file LICENSE.LGPL ** 00019 ** included in the packaging of this file. Please review the ** 00020 ** following information to ensure the license requirements will ** 00021 ** be met: http://www.gnu.org/licenses/lgpl-3.0.txt ** 00022 ** ** 00023 ** The development of this software was supported by the ** 00024 ** Excellence Cluster EXC 277 Cognitive Interaction Technology. ** 00025 ** The Excellence Cluster EXC 277 is a grant of the Deutsche ** 00026 ** Forschungsgemeinschaft (DFG) in the context of the German ** 00027 ** Excellence Initiative. ** 00028 ** ** 00029 ********************************************************************/ 00030 00031 #pragma once 00032 00033 #include <ICLUtils/CompatMacros.h> 00034 #include <ICLUtils/Point32f.h> 00035 #include <stdio.h> 00036 00037 namespace icl{ 00038 namespace core{ 00039 00041 00048 struct Parable{ 00049 00051 Parable():a(0),b(0),c(0){} 00052 00054 Parable(float a, float b,float c):a(a),b(b),c(c){} 00055 00057 00064 Parable(utils::Point32f p1,utils::Point32f p2,utils::Point32f p3){ 00065 p1.x+=0.00000001; // avoid numerical problems 00066 p2.x+=0.00000002; 00067 p3.x+=0.00000003; 00068 00069 float xx1 = p1.x * p1.x; 00070 float xx2 = p2.x * p2.x; 00071 float xx3 = p3.x * p3.x; 00072 b = ((p3.y-p1.y)/(xx3-xx1)-(p2.y-p1.y)/(xx2-xx1)) / 00073 ((p1.x-p2.x)/(xx2-xx1)-(p1.x-p3.x)/(xx3-xx1)); 00074 00075 a = (b*(p1.x-p2.x)+(p2.y-p1.y)) / (xx2-xx1); 00076 00077 c = p1.y-a*xx1-b*p1.x; 00078 } 00079 00080 float a; 00081 float b; 00082 float c; 00084 00085 float operator()(float x) const{ return a*x*x+b*x+c; } 00086 00088 void show()const{ 00089 printf("Parable:\n"); 00090 printf("f(x)=%f*x²+%fx+%f\n",a,b,c); 00091 } 00092 }; 00093 00094 } // namespace core 00095 } 00096