Solving issue 1583: "Some tools always collapsed after restart"

rtengine/LUT.h

@@ -0,0 +1,213 @@
+/*
+ * LUT.h
+ *  This file is part of RawTherapee.
+ *
+ *  Copyright (c) 2011 Jan Rinze Peterzon (janrinze@gmail.com)
+ *
+ *  RawTherapee is free software: you can redistribute it and/or modify
+ *  it under the terms of the GNU General Public License as published by
+ *  the Free Software Foundation, either version 3 of the License, or
+ *  (at your option) any later version.
+ *
+ *  RawTherapee 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 General Public License for more details.
+ *
+ *  You should have received a copy of the GNU General Public License
+ *  along with RawTherapee.  If not, see <http://www.gnu.org/licenses/>.
+ */
+
+/*
+ *  Declaration of flexible Lookup Tables
+ *
+ *  Usage:
+ *
+ *  	LUT<type> name (size);
+ *		LUT<type> name (size, flags);
+ *
+ *		creates an array which is valid within the normal C/C++ scope "{ ... }"
+ *
+ *      access to elements is a simple as:
+ *
+ *      	LUT<float> my_lut (10);
+ *      	float value = my_lut[3];
+ *          float value = my_lut[2.5]; // this will interpolate
+ *
+ *      when using a float type index it will interpolate the lookup values
+ *
+ *      extra setting in flags: (clipping is set by default)
+ *      LUT_CLIP_ABOVE
+ *      LUT_CLIP_BELOW
+ *
+ *      example:
+ *      	LUT<float> my_lut (10,LUT_CLIP_BELOW);
+ *          float value = my_lut[22.5];  // this will extrapolate
+ *          float value = my_lut[-22.5]; // this will not extrapolate
+ *
+ *          LUT<float> my_lut (10,0); // this will extrapolate on either side
+ *
+ *      shotcuts:
+ *
+ *      	LUTf stands for LUT<float>
+ *          LUTi stands for LUT<int>
+ *          LUTu stands for LUT<unsigned int>
+ */
+
+#ifndef LUT_H_
+#define LUT_H_
+
+// bit representations of flags
+#define LUT_CLIP_BELOW 1
+#define LUT_CLIP_ABOVE 2
+
+#define LUTf LUT<float>
+#define LUTi LUT<int>
+#define LUTu LUT<unsigned int>
+#define LUTd LUT<double>
+
+#include <cstring>
+#ifndef NDEBUG
+#include <glibmm.h>
+#include <fstream>
+#endif
+
+template<typename T>
+class LUT {
+private:
+	// list of variables ordered to improve cache speed
+	unsigned int maxs; 
+	T * data;
+	unsigned int clip, size, owner;
+public:
+	LUT(int s, int flags = 0xfffffff) {
+		clip = flags;
+		data = new T[s];
+		owner = 1;
+		size = s;
+		maxs=size-2;
+	}
+	void operator ()(int s, int flags = 0xfffffff) {
+		if (owner&&data)
+			delete[] data;
+		clip = flags;
+		data = new T[s];
+		owner = 1;
+		size = s;
+		maxs=size-2;
+	}
+
+	LUT(int s, T * source, int flags = 0xfffffff) {
+		clip = flags;
+		data = new T[s];
+		owner = 1;
+		size = s;
+		maxs=size-2;
+		for (int i = 0; i < s; i++) {
+			data[i] = source[i];
+		}
+	}
+
+	LUT(void) {
+		data = NULL;
+		reset();
+	}
+
+	~LUT() {
+		if (owner)
+			delete[] data;
+	}
+
+	void setClip(int flags) {
+		clip = flags;
+	}
+
+	LUT<T> & operator=(const LUT<T> &rhs) {
+	    if (this != &rhs) {
+	      if (rhs.size>this->size)
+	      {
+	    	delete [] this->data;
+	    	this->data=NULL;
+	      }
+	      if (this->data==NULL) this->data=new T[rhs.size];
+	      this->clip=rhs.clip;
+	      this->owner=1;
+	      memcpy(this->data,rhs.data,rhs.size*sizeof(T));
+	      this->size=rhs.size;
+	      this->maxs=this->size-2;
+	    }
+
+	    return *this;
+	  }
+	// use with integer indices
+	T& operator[](int index) const {
+		if (((unsigned int)index)<size) return data[index];
+		else
+		{
+			if (index < 0)
+				return data[0];
+			else
+				return data[size - 1];
+		}
+		
+	}
+	// use with float indices
+	T operator[](float index) const {
+		int idx = (int)index;  // don't use floor! The difference in negative space is no problems here
+		if (((unsigned int)idx) > maxs) {
+			if (idx<0)
+			{
+				if (clip & LUT_CLIP_BELOW)
+					return data[0];
+				idx=0;
+			}
+			else
+			{
+				if (clip & LUT_CLIP_ABOVE)
+					return data[size - 1];
+				idx =maxs;
+			}
+		}
+		float diff = index - (float) idx;
+		T p1 = data[idx];
+		T p2 = data[idx + 1]-p1;
+		return (p1 + p2*diff);
+	}
+	
+	
+#ifndef NDEBUG
+	// Debug facility ; dump the content of the LUT in a file. No control of the filename is done
+	void dump(Glib::ustring fname) {
+		if (size) {
+		    Glib::ustring fname_ = fname + ".xyz"; // TopSolid'Design "plot" file format
+			std::ofstream f (fname_.c_str());
+			f << "$" << std::endl;;
+			for (unsigned int iter=0; iter<size; iter++) {
+				f << iter << ", " << data[iter] << ", 0." << std::endl;;
+			}
+			f << "$" << std::endl;;
+			f.close ();
+		}
+	}
+#endif
+
+
+	operator bool (void) const
+		{
+			return size>0;
+		}
+
+	void clear(void) {
+		memset(data, 0, size * sizeof(T));
+	}
+
+	void reset(void) {
+		delete[] data;
+		data = NULL;
+		owner = 1;
+		size = 0;
+		maxs=0;
+    }
+};
+
+#endif /* LUT_H_ */