///////////////////////////////////////////////////////////////////////////

//

// Copyright (c) 2004-2012, Industrial Light & Magic, a division of Lucas

// Digital Ltd. LLC

// 

// All rights reserved.

// 

// Redistribution and use in source and binary forms, with or without

// modification, are permitted provided that the following conditions are

// met:

// *       Redistributions of source code must retain the above copyright

// notice, this list of conditions and the following disclaimer.

// *       Redistributions in binary form must reproduce the above

// copyright notice, this list of conditions and the following disclaimer

// in the documentation and/or other materials provided with the

// distribution.

// *       Neither the name of Industrial Light & Magic nor the names of

// its contributors may be used to endorse or promote products derived

// from this software without specific prior written permission. 

// 

// THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS

// "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT

// LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR

// A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT

// OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,

// SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT

// LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE,

// DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY

// THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT

// (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE

// OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.

//

///////////////////////////////////////////////////////////////////////////

#include "compareB44.h"

#include <ImfOutputFile.h>

#include <ImfInputFile.h>

#include <ImfChannelList.h>

#include <ImfArray.h>

#include <ImathRandom.h>

#include <half.h>

#include "compareFloat.h"

#include <stdio.h>

#include <assert.h>

namespace IMF = OPENEXR_IMF_NAMESPACE;

using namespace IMF;

using namespace std;

using namespace IMATH_NAMESPACE;

namespace {

void

fillPixels1 (Array2D<unsigned int> &pi,

	     Array2D<half> &ph,

	     Array2D<float> &pf,

	     int width,

	     int height)

{

    cout << "only zeroes" << endl;

    for (int y = 0; y < height; ++y)

	for (int x = 0; x < width; ++x)

	{

	    pi[y][x] = 0;

	    ph[y][x] = 0;

	    pf[y][x] = 0;

	}

}

void

fillPixels2 (Array2D<unsigned int> &pi,

	     Array2D<half> &ph,

	     Array2D<float> &pf,

	     int width,

	     int height)

{

    cout << "pattern 1" << endl;

    for (int y = 0; y < height; ++y)

	for (int x = 0; x < width; ++x)

	{

	    pi[y][x] = (x + y) & 1;

	    ph[y][x] = pi[y][x];

	    pf[y][x] = pi[y][x];

	}

}

void

fillPixels3 (Array2D<unsigned int> &pi,

	     Array2D<half> &ph,

	     Array2D<float> &pf,

	     int width,

	     int height)

{

    cout << "pattern 2" << endl;

    for (int y = 0; y < height; ++y)

	for (int x = 0; x < width; ++x)

	{

	    pi[y][x] = x % 100 + 100 * (y % 100);

	    ph[y][x] = sin (double (x)) + sin (y * 0.5);

	    pf[y][x] = sin (double (y)) + sin (x * 0.5);

	}

}

void

fillPixels4 (Array2D<unsigned int> &pi,

	     Array2D<half> &ph,

	     Array2D<float> &pf,

	     int width,

	     int height)

{

    cout << "random bits" << endl;

    //

    // Use of a union to extract the bit pattern from a float, as is

    // done below, works only if int and float have the same size.

    //

    assert (sizeof (int) == sizeof (float));

    Rand48 rand;

    for (int y = 0; y < height; ++y)

	for (int x = 0; x < width; ++x)

	{

	    pi[y][x] = rand.nexti();

	    ph[y][x].setBits (rand.nexti());

	    union {int i; float f;} u;

	    u.i = rand.nexti();

	    pf[y][x] = u.f;

	}

}

void

writeRead (const Array2D<unsigned int> &pi1,

	   const Array2D<half> &ph1,

	   const Array2D<float> &pf1,

	   const char fileName[],

	   int width,

	   int height,

	   int xOffset,

	   int yOffset,

	   Compression comp,

	   int xs,

	   int ys)

{

    //

    // Write the pixel data in pi1, ph1 and ph2 to an

    // image file using the specified compression type

    // and subsampling rates.  Read the pixel data back

    // from the file and verify that the data did not

    // change.

    //

    cout << "compression " << comp <<

	    ", x sampling " << xs <<

	    ", y sampling " << ys <<

	    ":" << flush;

    Header hdr ((Box2i (V2i (0, 0),			// display window

		        V2i (width - 1, height -1))),

		(Box2i (V2i (xOffset, yOffset),		// data window

		        V2i (xOffset + width - 1, yOffset + height - 1))));

    hdr.compression() = comp;

    hdr.channels().insert ("I",			// name

			   Channel (IMF::UINT,	// type

				    xs,		// xSampling

				    ys)		// ySampling

			  );

    hdr.channels().insert ("H",			// name

			   Channel (IMF::HALF,	// type

				    xs,		// xSampling

				    ys)		// ySampling

			  );

    hdr.channels().insert ("F",			// name

			   Channel (IMF::FLOAT,	// type

				    xs,		// xSampling

				    ys)		// ySampling

			  );

    {

	FrameBuffer fb; 

	fb.insert ("I",						// name

		   Slice (IMF::UINT,				// type

			  (char *) &pi1[-yOffset / ys][-xOffset / xs], // base

			  sizeof (pi1[0][0]), 			// xStride

			  sizeof (pi1[0][0]) * (width / xs),	// yStride

			  xs,					// xSampling

			  ys)					// ySampling

		  );

	fb.insert ("H",						// name

		   Slice (IMF::HALF,				// type

			  (char *) &ph1[-yOffset / ys][-xOffset / xs], // base

			  sizeof (ph1[0][0]), 			// xStride

			  sizeof (ph1[0][0]) * (width / xs),	// yStride

			  xs,					// xSampling

			  ys)					// ySampling

		  );

	fb.insert ("F",						// name

		   Slice (IMF::FLOAT,				// type

			  (char *) &pf1[-yOffset / ys][-xOffset / xs], // base

			  sizeof (pf1[0][0]), 			// xStride

			  sizeof (pf1[0][0]) * (width / xs),	// yStride

			  xs,					// xSampling

			  ys)					// ySampling

		  );

	cout << " writing" << flush;

	remove (fileName);

	OutputFile out (fileName, hdr);

	out.setFrameBuffer (fb);

	out.writePixels (height);

    }

    {

	cout << " reading" << flush;

	InputFile in (fileName);

	const Box2i &dw = in.header().dataWindow();

	int w = dw.max.x - dw.min.x + 1;

	int h = dw.max.y - dw.min.y + 1;

	int dx = dw.min.x;

	int dy = dw.min.y;

	Array2D<unsigned int> pi2 (h / ys, w / xs);

	Array2D<half>         ph2 (h / ys, w / xs);

	Array2D<float>        pf2 (h / ys, w / xs);

	FrameBuffer fb;

	{

	    int xs = in.header().channels()["I"].xSampling;

	    int ys = in.header().channels()["I"].ySampling;

	    fb.insert ("I",					// name

		       Slice (IMF::UINT,			// type

			      (char *) &pi2[-dy / ys][-dx / xs], // base

			      sizeof (pi2[0][0]), 		// xStride

			      sizeof (pi2[0][0]) * (w / xs),	// yStride

			      xs,				// xSampling

			      ys)				// ySampling

		      );

	}

	{

	    int xs = in.header().channels()["H"].xSampling;

	    int ys = in.header().channels()["H"].ySampling;

	    fb.insert ("H",					// name

		       Slice (IMF::HALF,			// type

			      (char *) &ph2[-dy / ys][-dx / xs], // base

			      sizeof (ph2[0][0]), 		// xStride

			      sizeof (ph2[0][0]) * (w / xs),	// yStride

			      xs,				// xSampling

			      ys)				// ySampling

		      );

	}

	{

	    int xs = in.header().channels()["F"].xSampling;

	    int ys = in.header().channels()["F"].ySampling;

	    fb.insert ("F",					// name

		       Slice (IMF::FLOAT,			// type

			      (char *) &pf2[-dy / ys][-dx / xs], // base

			      sizeof (pf2[0][0]), 		// xStride

			      sizeof (pf2[0][0]) * (w / xs),	// yStride

			      xs,				// xSampling

			      ys)				// ySampling

		      );

	}

	in.setFrameBuffer (fb);

	in.readPixels (dw.min.y, dw.max.y);

	cout << " comparing" << flush;

	assert (in.header().displayWindow() == hdr.displayWindow());

	assert (in.header().dataWindow() == hdr.dataWindow());

	assert (in.header().pixelAspectRatio() == hdr.pixelAspectRatio());

	assert (in.header().screenWindowCenter() == hdr.screenWindowCenter());

	assert (in.header().screenWindowWidth() == hdr.screenWindowWidth());

	assert (in.header().lineOrder() == hdr.lineOrder());

	assert (in.header().compression() == hdr.compression());

	ChannelList::ConstIterator hi = hdr.channels().begin();

	ChannelList::ConstIterator ii = in.header().channels().begin();

	while (hi != hdr.channels().end())

	{

	    assert (!strcmp (hi.name(), ii.name()));

	    assert (hi.channel().type == ii.channel().type);

	    assert (hi.channel().xSampling == ii.channel().xSampling);

	    assert (hi.channel().ySampling == ii.channel().ySampling);

	    ++hi;

	    ++ii;

	}

	assert (ii == in.header().channels().end());

	for (int y = 0; y < h / ys; ++y)

	{

	    for (int x = 0; x < w / xs; ++x)

	    {

		assert (pi1[y][x] == pi2[y][x]);

		assert (equivalent (pf1[y][x], pf2[y][x], comp));

		if (comp != B44_COMPRESSION &&

                    comp != B44A_COMPRESSION)

                {

		    assert (ph1[y][x].bits() == ph2[y][x].bits());

                }

	    }

	}

	if (comp == B44_COMPRESSION ||

            comp == B44A_COMPRESSION)

	{

	    Array2D<half> ph3 (h / ys, w / xs);

	    for (int y = 0; y < h / ys; ++y)

		for (int x = 0; x < w / xs; ++x)

		    ph3[y][x] = ph1[y][x];

	    compareB44 (w / xs, h / ys, ph3, ph2);

	}

    }

    remove (fileName);

    cout << endl;

}

void

writeRead (const std::string &tempDir,

           const Array2D<unsigned int> &pi,

	   const Array2D<half> &ph,

	   const Array2D<float> &pf,

	   int w,

	   int h,

	   int dx,

	   int dy)

{

    std::string filename = tempDir + "imf_test_comp.exr";

    for (int xs = 1; xs <= 2; ++xs)

    {

	for (int ys = 1; ys <= 2; ++ys)

	{

	    for (int comp = 0; comp < NUM_COMPRESSION_METHODS; ++comp)

	    {

		writeRead (pi, ph, pf,

			   filename.c_str(),

			   w  * xs, h  * ys,

			   dx * xs, dy * ys,

			   Compression (comp),

			   xs, ys);

	    }

	}

    }

}

} // namespace

void

testCompression (const std::string &tempDir)

{

    try

    {

	cout << "Testing pixel data types, "

		"subsampling and "

		"compression schemes" << endl;

	const int W = 1371;

	const int H = 159;

	const int DX = 17;

	const int DY = 29;

	Array2D<unsigned int> pi (H, W);

	Array2D<half> ph (H, W);

	Array2D<float> pf (H, W);

	//

	// If the following assertion fails, new pixel types have

	// been added to the Imf library; testing code for the new

	// pixel types should be added to this file.

	//

	assert (NUM_PIXELTYPES == 3);

	fillPixels1 (pi, ph, pf, W, H);

	writeRead (tempDir, pi, ph, pf, W, H, DX, DY);

	fillPixels2 (pi, ph, pf, W, H);

	writeRead (tempDir, pi, ph, pf, W, H, DX, DY);

	fillPixels3 (pi, ph, pf, W, H);

	writeRead (tempDir, pi, ph, pf, W, H, DX, DY);

	fillPixels4 (pi, ph, pf, W, H);

	writeRead (tempDir, pi, ph, pf, W, H, DX, DY);

	cout << "ok\n" << endl;

    }

    catch (const std::exception &e)

    {

	cerr << "ERROR -- caught exception: " << e.what() << endl;

	assert (false);

    }

}