/*! ========================================================================

** Extended Template and Library Test Suite

** Handle Template Class Test

** $Id$

**

** Copyright (c) 2002 Robert B. Quattlebaum Jr.

**

** This package 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 2 of

** the License, or (at your option) any later version.

**

** This package 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.

**

** === N O T E S ===========================================================

**

** ========================================================================= */

/* === H E A D E R S ======================================================= */

#include <list></list>

#include <cstdio></cstdio>

#include <cstdlib></cstdlib>

#include <string></string>

#include <utility></utility>

#include <memory></memory>

#include <map></map>

#include <etl pen=""></etl>

#include <etl boxblur=""></etl>

//#include <etl gaussian=""></etl>

#include <cmath></cmath>

/* === M A C R O S ========================================================= */

using namespace std;

using namespace etl;

/* === C L A S S E S ======================================================= */

int generic_pen_test(int w, int h)

{

	printf("generic_pen(w:%d,h:%d): ",w,h);

	auto_ptr<float> data(new float[w*h]);</float>

	//unique_ptr<float> data(new float[w*h]);</float>

	if(!data.get())

	{

		printf("Um..... malloc failure on line %d of " __FILE__ "...\n", __LINE__);

		abort();

	}

	generic_pen<float> pen(data.get(),w,h);</float>

	generic_pen<float> pen2;</float>

	if(!pen)

	{

		printf("FAILURE! " __FILE__ "@%d: On pen bool test\n", __LINE__);

		return 1;

	}

	if(&pen.x()[2]!=&pen[0][2])

	{

		printf("FAILURE! " __FILE__ "@%d: On request for horizontal iterator\n", __LINE__);

		return 1;

	}

	if(&pen.y()[2]!=&pen[2][0])

	{

		printf("FAILURE! " __FILE__ "@%d: On request for vertical iterator\n", __LINE__);

		return 1;

	}

	pen.move(1,1);

	pen2=pen;

	if(pen!=pen2)

	{

		printf("FAILURE! " __FILE__ "@%d: On pen assignment or pen comparison\n", __LINE__);

		return 1;

	}

	pen2.move(w,h);

	generic_pen<float>::difference_type diff(pen2-pen);</float>

	if(diff.x!=w || diff.y!=h)

	{

		printf("FAILURE! " __FILE__ "@%d: pen difference inconsistency ([%d,%d]!=[%d,%d])\n", __LINE__, diff.x, diff.y, w, h);

		return 1;

	}

	if(pen.end_x()-pen.x()!=w-1)

	{

		printf("FAILURE! " __FILE__ "@%d: iterator_x inconsistency (%ld!=%d)\n", __LINE__, pen.end_x()-pen.x(), w);

		return 1;

	}

	if(pen.end_y()-pen.y()!=h-1)

	{

		printf("FAILURE! " __FILE__ "@%d: iterator_y inconsistency (%d!=%d)\n", __LINE__, pen.end_y()-pen.y(), h);

		return 1;

	}

	if(&pen.end_y()[-1]!=&pen.y()[(h-2)])

	{

		printf("FAILURE! " __FILE__ "@%d: iterator_y inconsistency\n", __LINE__);

		return 1;

	}

	if(&pen.end_x()[-1]!=&pen.x()[(w-2)])

	{

		printf("FAILURE! " __FILE__ "@%d: iterator_x inconsistency\n", __LINE__);

		return 1;

	}

	printf("PASSED\n");

	return 0;

}

int alpha_pen_test(void)

{

	printf("alpha_pen: ");

	printf("SKIPPED\n");

	return 0;

}

int bbox_pen_test(void)

{

	printf("bbox_pen: ");

	printf("SKIPPED\n");

	return 0;

}

int display_pen(generic_pen<float> pen, int w, int h)</float>

{

	int ret=0;

	int x, y;

	// print out the after pic

	for(y=0;y

	{

		printf("|");

		for(x=0;x

		{

			if(pen.get_value()>=2.0f)

				printf("#");

			else if(pen.get_value()>=1.0f)

				printf("@");

			else if(pen.get_value()>=0.8f)

				printf("%%");

			else if(pen.get_value()>=0.6f)

				printf("O");

			else if(pen.get_value()>=0.4f)

				printf(":");

			else if(pen.get_value()>=0.2f)

				printf(".");

			else if(pen.get_value()>=-0.0001f)

				printf(" ");

			else

				printf("X"),ret++;

		}

		pen.dec_x(x);

		printf("|\n");

	}

	pen.dec_y(y);

	return ret;

}

int display_pen(generic_pen<double> pen, int w, int h)</double>

{

	int ret=0;

	int x, y;

	// print out the after pic

	for(y=0;y

	{

		printf("|");

		for(x=0;x

		{

			if(pen.get_value()>=2.0f)

				printf("#");

			else if(pen.get_value()>=1.0f)

				printf("@");

			else if(pen.get_value()>=0.8f)

				printf("%%");

			else if(pen.get_value()>=0.6f)

				printf("O");

			else if(pen.get_value()>=0.4f)

				printf(":");

			else if(pen.get_value()>=0.2f)

				printf(".");

			else if(pen.get_value()>=-0.0001f)

				printf(" ");

			else

				printf("X"),ret++;

		}

		pen.dec_x(x);

		printf("|\n");

	}

	pen.dec_y(y);

	return ret;

}

void emptyfunction(int v)

{

	static int stupid = 0;

	stupid = v;

	if (stupid == 0) return; // disable unused warning

	//printf("Called... %d\n",v);

}

int box_blur_test(void)

{

	typedef float boxblur_float;

	printf("box_blur: ");

	int w=25,h=25;

	//unique_ptr<boxblur_float> data(new boxblur_float[w*h]);</boxblur_float>

	//unique_ptr<boxblur_float> data2(new boxblur_float[w*h]);</boxblur_float>

	auto_ptr<boxblur_float> data(new boxblur_float[w*h]);</boxblur_float>

	auto_ptr<boxblur_float> data2(new boxblur_float[w*h]);</boxblur_float>

	if(!data.get())

	{

		printf("Um..... malloc failure on line %d of " __FILE__ "...\n", __LINE__);

		abort();

	}

	generic_pen<boxblur_float> pen(data.get(),w,h);</boxblur_float>

	generic_pen<boxblur_float> pen2;</boxblur_float>

	generic_pen<boxblur_float> pen3(data2.get(),w,h);</boxblur_float>

	int x,y;

	for(y=0;y

	{

		for(x=0;x

		{

			if( (x-y<=1 && y-x<=1) || y==h/2 || x==w/2)

				pen.put_value(2);

			else

				pen.put_value(0);

		}

		pen.dec_x(x);

	}

	pen.dec_y(y);

	int bad_values=0;

	printf("\nBEFORE BOX BLUR:\n");

	// print out the before pic

	display_pen(pen,w,h);

	// Pen 2 will be the end

	pen2=pen;

	pen2.move(w,h);

	//temporary

	vbox_blur(pen,pen2,2,pen3);

	printf("\n VBLUR ONLY:\n");

	display_pen(pen3,w,h);

//	box_blur(pen,w,h,4);

	hbox_blur(pen,pen2,2,pen3);

	printf("\n HBLUR ONLY:\n");

	display_pen(pen3,w,h);

	pen2=pen3;

	pen2.move(w,h);

	vbox_blur(pen3,pen2,2,pen);

	printf("\nAFTER BOX BLUR:\n");

	// print out the after pic

	bad_values=display_pen(pen,w,h);

	if(bad_values)

	{

		printf("FAILURE! " __FILE__ "@%d: blur result contained %d bad values\n", __LINE__, bad_values);

		return 1;

	}

	boxblur_float max=0;

	printf("CHECK BOXBLUR RESULTS %d,%d:\n",pen.diff_begin().x, pen.diff_begin().y);

	for(y=0;y

	{

		for(x=0;x

		{

			boxblur_float f = 0;

			for(int oy=-2; oy <= 2; ++oy)

			{

				int iy = y+oy;

				if(iy < 0) iy = 0;

				if(iy >= h) iy = h-1;

				for(int ox=-2; ox <= 2; ++ox)

				{

					int ix = x+ox;

					if(ix < 0) ix = 0;

					if(ix >= w) ix = w-1;

					if( (ix-iy<=1 && iy-ix<=1) || iy==h/2 || ix==w/2)

						f += 2;

				}

			}

			//print out if the relative error is high

			/*f /= 25;

			float rf = pen.get_value() - f/25;

			if(f && rf > 0.3)

			{

				printf("pixel (%d,%d) off by %f\n",x,y,rf);

			}*/

			boxblur_float diff = fabs(pen.get_value() - f/25);

			if(diff > max) max = diff;

			pen.put_value(f/25); //if length = 2 then dim = 5.. area = 25

		}

		pen.dec_x(x);

	}

	pen.dec_y(y);

	/*if(max)

	{

		for(y=0;y

		{

			for(x=0;x

			{

				pen.put_value(pen.get_value()/max);

			}

			pen.dec_x(x);

		}

		pen.dec_y(y);

	}*/

	//printf("\nHBOXBLUR ERROR (max = %e):\n",max);

	printf("\nCorrect results:\n");

	display_pen(pen,w,h);

	printf("PASSED\n");

	return 0;

}

/*

float:

|@@%O.     :::::          |

|@@@%:.    :::::          |

|%@@%O:.   :::::          |

|O%%@%O:.  :::::          |

|.:O%@%O:. :::::          |

| .:O%@%O:.:::::          |

|  .:O%@%O:O::::          |

|   .:O%@%O%O:::          |

|    .:O%@%@%O::          |

|     .:O%@@@%::          |

|::.:::O%@@@@@%O::::::::::|

|::.::::O%@@@@@%::::::::::|

|::.:::::OO@@@@@%O::::::::|

|::.:::::::%@@@@@%O:::::::|

|::.::::::.O%@@@@@%O::::::|

|          ::%@@@%O:.     |

|          ::O%@%@%O:.    |

|          :.:O%O%@%O:.   |

|          :.::O:O%@%O:.  |

|          :.:.:.:O%@%O:. |

|          :.:.: .:O%@%O:.|

|          :.:.:  .:O%@%%O|

|          :.:.:   .:O%@@%|

|          :.:.:    .:%@@@|

|          :.:.:     .O%@@|

double:

|@@%O.     .....          |

|@@@O:.    .....          |

|%@@%O:.   .....          |

|OO%@%O:.  .....          |

|.:O%@%O:. .....          |

| .:O%@%O:.:....          |

|  .:O%@%O:O:...          |

|   .:O%@%O%O:..          |

|    .:O%@%@%O:.          |

|     .:O%@@@O:.          |

|.....:O%@@@@@%O..........|

|......:O%@@@@@%::........|

|.......:OO@@@@@OO:.......|

|........::%@@@@@%O:......|

|..........O%@@@@@%O:.....|

|          .:O@@@%O:.     |

|          .:O%@%@%O:.    |

|          ..:O%O%@%O:.   |

|          ...:O:O%@%O:.  |

|          ....:.:O%@%O:. |

|          ..... .:O%@%O:.|

|          .....  .:O%@%OO|

|          .....   .:O%@@%|

|          .....    .:O@@@|

|          .....     .O%@@|

*/

int gaussian_blur_test(void)

{

	printf("gaussian_blur: ");

#if 0

	int w=25,h=25;

	int bad_values=0;

	auto_ptr<float> data(new float[w*h]);</float>

	if(!data.get())

	{

		printf("Um..... malloc failure on line %d of " __FILE__ "...\n",__LINE__);

		abort();

	}

	generic_pen<float> pen(data.get(),w,h);</float>

	generic_pen<float> pen2;</float>

	int x,y;

	for(y=0;y

	{

		for(x=0;x

		{

			if((x-y<=1 && y-x<=1) || y==h/2)

				pen.put_value(2);

			else

				pen.put_value(0);

		}

		pen.dec_x(x);

	}

	pen.dec_y(y);

	printf("\nBEFORE GAUSSIAN BLUR:\n");

	// print out the before pic

	for(y=0;y

	{

		printf("|");

		for(x=0;x

		{

			if(pen.get_value()>=2.0f)

				printf("#");

			else if(pen.get_value()>=1.0f)

				printf("@");

			else if(pen.get_value()>=0.8f)

				printf("%%");

			else if(pen.get_value()>=0.6f)

				printf("O");

			else if(pen.get_value()>=0.4f)

				printf(":");

			else if(pen.get_value()>=0.2f)

				printf(".");

			else if(pen.get_value()>=0.0f)

				printf(" ");

			else

				printf("X"),bad_values++;

		}

		pen.dec_x(x);

		printf("|\n");

	}

	pen.dec_y(y);

	// Pen 2 will be the end

	pen2=pen;

	pen2.move(w,h);

#if 0

	gaussian_blur_5x5(pen,pen2);

	gaussian_blur_5x5(pen,pen2);

	gaussian_blur_5x5(pen,pen2);

#endif

#if 0

	gaussian_blur_3x3(pen,pen2);

	gaussian_blur_3x3(pen,pen2);

	gaussian_blur_3x3(pen,pen2);

	gaussian_blur_3x3(pen,pen2);

	gaussian_blur_3x3(pen,pen2);

#endif

//	gaussian_blur(pen,pen2,15);

	gaussian_blur(pen,pen2,10,10);

	printf("\nAFTER GAUSSIAN BLUR:\n");

	// print out the after pic

	for(y=0;y

	{

		printf("|");

		for(x=0;x

		{

			if(pen.get_value()>=2.0f)

				printf("#");

			else if(pen.get_value()>=1.0f)

				printf("@");

			else if(pen.get_value()>=0.8f)

				printf("%%");

			else if(pen.get_value()>=0.6f)

				printf("O");

			else if(pen.get_value()>=0.4f)

				printf(":");

			else if(pen.get_value()>=0.2f)

				printf(".");

			else if(pen.get_value()>=0.0f)

				printf(" ");

			else

				printf("X"),bad_values++;

		}

		pen.dec_x(x);

		printf("|\n");

	}

	pen.dec_y(y);

	if(bad_values)

	{

		printf("FAILURE! " __FILE__ "@%d: blur result contained bad values\n",__LINE__);

		return 1;

	}

#endif

	printf("PASSED\n");

	return 0;

}

/* === E N T R Y P O I N T ================================================= */

int main()

{

	int error=0;

	error+=generic_pen_test(40,40);

	error+=generic_pen_test(10,40);

	error+=generic_pen_test(40,10);

    if(error)return error;

	error+=alpha_pen_test();

	error+=bbox_pen_test();

	error+=box_blur_test();

    if(error)return error;

	error+=gaussian_blur_test();

	return error;

}