Blame IlmImfTest/testYca.cpp

Packit 0d464f
//////////////////////////////////////////////////////////////////////////////
Packit 0d464f
//
Packit 0d464f
// Copyright (c) 2004-2012, Industrial Light & Magic, a division of Lucasfilm
Packit 0d464f
// Entertainment Company Ltd.  Portions contributed and copyright held by
Packit 0d464f
// others as indicated.  All rights reserved.
Packit 0d464f
//
Packit 0d464f
// Redistribution and use in source and binary forms, with or without
Packit 0d464f
// modification, are permitted provided that the following conditions are
Packit 0d464f
// met:
Packit 0d464f
//
Packit 0d464f
//     * Redistributions of source code must retain the above
Packit 0d464f
//       copyright notice, this list of conditions and the following
Packit 0d464f
//       disclaimer.
Packit 0d464f
//
Packit 0d464f
//     * Redistributions in binary form must reproduce the above
Packit 0d464f
//       copyright notice, this list of conditions and the following
Packit 0d464f
//       disclaimer in the documentation and/or other materials provided with
Packit 0d464f
//       the distribution.
Packit 0d464f
//
Packit 0d464f
//     * Neither the name of Industrial Light & Magic nor the names of
Packit 0d464f
//       any other contributors to this software may be used to endorse or
Packit 0d464f
//       promote products derived from this software without specific prior
Packit 0d464f
//       written permission.
Packit 0d464f
//
Packit 0d464f
// THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS
Packit 0d464f
// IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO,
Packit 0d464f
// THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
Packit 0d464f
// PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR
Packit 0d464f
// CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
Packit 0d464f
// EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
Packit 0d464f
// PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
Packit 0d464f
// PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF
Packit 0d464f
// LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
Packit 0d464f
// NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
Packit 0d464f
// SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
Packit 0d464f
//
Packit 0d464f
//////////////////////////////////////////////////////////////////////////////
Packit 0d464f
Packit 0d464f
#include <ImfRgbaFile.h>
Packit 0d464f
#include <ImfArray.h>
Packit 0d464f
#include <ImfThreading.h>
Packit 0d464f
#include "IlmThread.h"
Packit 0d464f
#include "ImathMath.h"
Packit 0d464f
#include <stdio.h>
Packit 0d464f
#include <assert.h>
Packit 0d464f
#include <algorithm>
Packit 0d464f
Packit 0d464f
Packit 0d464f
using namespace OPENEXR_IMF_NAMESPACE;
Packit 0d464f
using namespace std;
Packit 0d464f
using namespace IMATH_NAMESPACE;
Packit 0d464f
Packit 0d464f
namespace {
Packit 0d464f
Packit 0d464f
void
Packit 0d464f
fillPixelsColor (Array2D <Rgba> &pixels, int w, int h)
Packit 0d464f
{
Packit 0d464f
    for (int y = 0; y < h; ++y)
Packit 0d464f
    {
Packit 0d464f
	for (int x = 0; x < w; ++x)
Packit 0d464f
	{
Packit 0d464f
	    Rgba &p = pixels[y][x];
Packit 0d464f
Packit 0d464f
	    p.r = 0.8 + 0.5 * sin (x * 0.05);
Packit 0d464f
	    p.g = 0.8 + 0.5 * sin (x * 0.02 + y * 0.02);
Packit 0d464f
	    p.b = 0.8 + 0.5 * sin (y * 0.03);
Packit 0d464f
Packit 0d464f
	    float t = 0.8 + 0.5 * sin (x * 0.05 - y * 0.05);
Packit 0d464f
Packit 0d464f
	    p.r *= t;
Packit 0d464f
	    p.g *= t;
Packit 0d464f
	    p.b *= t;
Packit 0d464f
	    p.a  = t;
Packit 0d464f
	}
Packit 0d464f
    }
Packit 0d464f
}
Packit 0d464f
Packit 0d464f
Packit 0d464f
void
Packit 0d464f
fillPixelsGray (Array2D <Rgba> &pixels, int w, int h)
Packit 0d464f
{
Packit 0d464f
    for (int y = 0; y < h; ++y)
Packit 0d464f
    {
Packit 0d464f
	for (int x = 0; x < w; ++x)
Packit 0d464f
	{
Packit 0d464f
	    Rgba &p = pixels[y][x];
Packit 0d464f
Packit 0d464f
	    p.r = 0.8 + 0.5 * sin (x * 0.05 - y * 0.05);
Packit 0d464f
	    p.g = p.r;
Packit 0d464f
	    p.b = p.r;
Packit 0d464f
	    p.a = 0.5 + 0.5 * cos (x * 0.05 - y * 0.05);
Packit 0d464f
	}
Packit 0d464f
    }
Packit 0d464f
}
Packit 0d464f
Packit 0d464f
Packit 0d464f
void
Packit 0d464f
writeReadYca (const char fileName[],
Packit 0d464f
	      const Box2i &dw,
Packit 0d464f
	      RgbaChannels channels,
Packit 0d464f
	      LineOrder writeOrder,
Packit 0d464f
	      LineOrder readOrder,
Packit 0d464f
	      void (* fillPixels) (Array2D <Rgba> &pixels, int w, int h))
Packit 0d464f
{
Packit 0d464f
    int w = dw.max.x - dw.min.x + 1;
Packit 0d464f
    int h = dw.max.y - dw.min.y + 1;
Packit 0d464f
    Array2D <Rgba> pixels1 (h, w);
Packit 0d464f
    Array2D <Rgba> pixels2 (h, w);
Packit 0d464f
Packit 0d464f
    cout << w << " by " << h << " pixels, "
Packit 0d464f
	    "channels " << channels << ", "
Packit 0d464f
	    "write order " << writeOrder << ", "
Packit 0d464f
	    "read order " << readOrder <<
Packit 0d464f
	    endl;
Packit 0d464f
Packit 0d464f
    fillPixels (pixels1, w, h);
Packit 0d464f
Packit 0d464f
    cout << "writing " << flush;
Packit 0d464f
Packit 0d464f
    {
Packit 0d464f
	RgbaOutputFile out (fileName,
Packit 0d464f
			    dw, dw,	// display window, data window
Packit 0d464f
			    channels,
Packit 0d464f
			    1,		// pixelAspectRatio
Packit 0d464f
			    V2f (0, 0),	// screenWindowCenter
Packit 0d464f
			    1,		// screenWindowWidth
Packit 0d464f
			    writeOrder);
Packit 0d464f
Packit 0d464f
	out.setYCRounding (9, 9);
Packit 0d464f
	out.setFrameBuffer (&pixels1[-dw.min.y][-dw.min.x], 1, w);
Packit 0d464f
	out.writePixels (h);
Packit 0d464f
    }
Packit 0d464f
Packit 0d464f
    cout << "reading " << flush;
Packit 0d464f
Packit 0d464f
    {
Packit 0d464f
	RgbaInputFile in (fileName);
Packit 0d464f
Packit 0d464f
	in.setFrameBuffer (&pixels2[-dw.min.y][-dw.min.x], 1, w);
Packit 0d464f
Packit 0d464f
	switch (readOrder)
Packit 0d464f
	{
Packit 0d464f
	  case INCREASING_Y:
Packit 0d464f
Packit 0d464f
	    for (int y = dw.min.y; y <= dw.max.y; ++y)
Packit 0d464f
		in.readPixels (y);
Packit 0d464f
Packit 0d464f
	    break;
Packit 0d464f
Packit 0d464f
	  case DECREASING_Y:
Packit 0d464f
Packit 0d464f
	    for (int y = dw.max.y; y >= dw.min.y; --y)
Packit 0d464f
		in.readPixels (y);
Packit 0d464f
Packit 0d464f
	    break;
Packit 0d464f
Packit 0d464f
	  case RANDOM_Y:
Packit 0d464f
Packit 0d464f
	    assert (h % 5 != 0);
Packit 0d464f
Packit 0d464f
	    for (int i = 0; i < h; ++i)
Packit 0d464f
	    {
Packit 0d464f
		int y = dw.min.y + (i * 5) % h;
Packit 0d464f
		in.readPixels (y);
Packit 0d464f
	    }
Packit 0d464f
Packit 0d464f
	    break;
Packit 0d464f
	}
Packit 0d464f
    }
Packit 0d464f
Packit 0d464f
    cout << "comparing" << endl;
Packit 0d464f
Packit 0d464f
    for (int y = 0; y < h; ++y)
Packit 0d464f
    {
Packit 0d464f
	for (int x = 0; x < w; ++x)
Packit 0d464f
	{
Packit 0d464f
	    const Rgba &p1 = pixels1[y][x];
Packit 0d464f
	    const Rgba &p2 = pixels2[y][x];
Packit 0d464f
Packit 0d464f
	    if (channels & WRITE_C)
Packit 0d464f
	    {
Packit 0d464f
		float p1Max = max (p1.r, max (p1.g, p1.b));
Packit 0d464f
		float p2Max = max (p2.r, max (p2.g, p2.b));
Packit 0d464f
Packit 0d464f
		assert (equalWithAbsError (p1Max, p2Max, 0.03f));
Packit 0d464f
	    }
Packit 0d464f
	    else
Packit 0d464f
	    {
Packit 0d464f
		assert (p1.g == p2.g);
Packit 0d464f
	    }
Packit 0d464f
Packit 0d464f
	    if (channels & WRITE_A)
Packit 0d464f
	    {
Packit 0d464f
		assert (p1.a == p2.a);
Packit 0d464f
	    }
Packit 0d464f
	}
Packit 0d464f
    }
Packit 0d464f
Packit 0d464f
    remove (fileName);
Packit 0d464f
}
Packit 0d464f
Packit 0d464f
} // namespace
Packit 0d464f
Packit 0d464f
Packit 0d464f
void
Packit 0d464f
testYca (const std::string &tempDir)
Packit 0d464f
{
Packit 0d464f
    try
Packit 0d464f
    {
Packit 0d464f
        cout << "Testing luminance/chroma input and output" << endl;
Packit 0d464f
Packit 0d464f
	std::string fileName = tempDir + "imf_test_yca.exr";
Packit 0d464f
Packit 0d464f
	Box2i dataWindow[6];
Packit 0d464f
	dataWindow[0] = Box2i (V2i (0, 0), V2i (1, 17));
Packit 0d464f
	dataWindow[1] = Box2i (V2i (0, 0), V2i (5, 17));
Packit 0d464f
	dataWindow[2] = Box2i (V2i (0, 0), V2i (17, 1));
Packit 0d464f
	dataWindow[3] = Box2i (V2i (0, 0), V2i (17, 5));
Packit 0d464f
	dataWindow[4] = Box2i (V2i (0, 0), V2i (1, 1));
Packit 0d464f
	dataWindow[5] = Box2i (V2i (-18, -28), V2i (247, 255));
Packit 0d464f
Packit 0d464f
	int maxThreads = ILMTHREAD_NAMESPACE::supportsThreads()? 3: 0;
Packit 0d464f
Packit 0d464f
	for (int n = 0; n <= maxThreads; ++n)
Packit 0d464f
	{
Packit 0d464f
	    if (ILMTHREAD_NAMESPACE::supportsThreads())
Packit 0d464f
	    {
Packit 0d464f
		setGlobalThreadCount (n);
Packit 0d464f
		cout << "\nnumber of threads: " << globalThreadCount() << endl;
Packit 0d464f
	    }
Packit 0d464f
Packit 0d464f
	    for (int i = 0; i < 6; ++i)
Packit 0d464f
	    {
Packit 0d464f
		for (int writeOrder = INCREASING_Y;
Packit 0d464f
		     writeOrder <= DECREASING_Y;
Packit 0d464f
		     ++writeOrder)
Packit 0d464f
		{
Packit 0d464f
		    for (int readOrder = INCREASING_Y;
Packit 0d464f
			 readOrder <= RANDOM_Y;
Packit 0d464f
			 ++readOrder)
Packit 0d464f
		    {
Packit 0d464f
			writeReadYca (fileName.c_str(), dataWindow[i],
Packit 0d464f
				      WRITE_YCA,
Packit 0d464f
				      LineOrder (writeOrder),
Packit 0d464f
				      LineOrder (readOrder),
Packit 0d464f
				      fillPixelsColor);
Packit 0d464f
Packit 0d464f
			writeReadYca (fileName.c_str(), dataWindow[i],
Packit 0d464f
				      WRITE_YC,
Packit 0d464f
				      LineOrder (writeOrder),
Packit 0d464f
				      LineOrder (readOrder),
Packit 0d464f
				      fillPixelsColor);
Packit 0d464f
Packit 0d464f
			writeReadYca (fileName.c_str(), dataWindow[i],
Packit 0d464f
				      WRITE_YA,
Packit 0d464f
				      LineOrder (writeOrder),
Packit 0d464f
				      LineOrder (readOrder),
Packit 0d464f
				      fillPixelsGray);
Packit 0d464f
Packit 0d464f
			writeReadYca (fileName.c_str(), dataWindow[i],
Packit 0d464f
				      WRITE_Y,
Packit 0d464f
				      LineOrder (writeOrder),
Packit 0d464f
				      LineOrder (readOrder),
Packit 0d464f
				      fillPixelsGray);
Packit 0d464f
		    }
Packit 0d464f
		}
Packit 0d464f
	    }
Packit 0d464f
	}
Packit 0d464f
Packit 0d464f
	cout << "ok\n" << endl;
Packit 0d464f
    }
Packit 0d464f
    catch (const std::exception &e)
Packit 0d464f
    {
Packit 0d464f
	cerr << "ERROR -- caught exception: " << e.what() << endl;
Packit 0d464f
	assert (false);
Packit 0d464f
    }
Packit 0d464f
}