|
Packit Service |
6754ca |
///////////////////////////////////////////////////////////////////////////
|
|
Packit Service |
6754ca |
//
|
|
Packit Service |
6754ca |
// Copyright (c) 2009-2014 DreamWorks Animation LLC.
|
|
Packit Service |
6754ca |
//
|
|
Packit Service |
6754ca |
// All rights reserved.
|
|
Packit Service |
6754ca |
//
|
|
Packit Service |
6754ca |
// Redistribution and use in source and binary forms, with or without
|
|
Packit Service |
6754ca |
// modification, are permitted provided that the following conditions are
|
|
Packit Service |
6754ca |
// met:
|
|
Packit Service |
6754ca |
// * Redistributions of source code must retain the above copyright
|
|
Packit Service |
6754ca |
// notice, this list of conditions and the following disclaimer.
|
|
Packit Service |
6754ca |
// * Redistributions in binary form must reproduce the above
|
|
Packit Service |
6754ca |
// copyright notice, this list of conditions and the following disclaimer
|
|
Packit Service |
6754ca |
// in the documentation and/or other materials provided with the
|
|
Packit Service |
6754ca |
// distribution.
|
|
Packit Service |
6754ca |
// * Neither the name of DreamWorks Animation nor the names of
|
|
Packit Service |
6754ca |
// its contributors may be used to endorse or promote products derived
|
|
Packit Service |
6754ca |
// from this software without specific prior written permission.
|
|
Packit Service |
6754ca |
//
|
|
Packit Service |
6754ca |
// THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
|
|
Packit Service |
6754ca |
// "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
|
|
Packit Service |
6754ca |
// LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR
|
|
Packit Service |
6754ca |
// A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT
|
|
Packit Service |
6754ca |
// OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
|
|
Packit Service |
6754ca |
// SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT
|
|
Packit Service |
6754ca |
// LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE,
|
|
Packit Service |
6754ca |
// DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY
|
|
Packit Service |
6754ca |
// THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
|
|
Packit Service |
6754ca |
// (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
|
|
Packit Service |
6754ca |
// OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
|
|
Packit Service |
6754ca |
//
|
|
Packit Service |
6754ca |
///////////////////////////////////////////////////////////////////////////
|
|
Packit Service |
6754ca |
|
|
Packit Service |
6754ca |
#include <math.h>
|
|
Packit Service |
6754ca |
#include <assert.h>
|
|
Packit Service |
6754ca |
|
|
Packit Service |
6754ca |
#include "half.h"
|
|
Packit Service |
6754ca |
#include "compareDwa.h"
|
|
Packit Service |
6754ca |
|
|
Packit Service |
6754ca |
using namespace OPENEXR_IMF_NAMESPACE;
|
|
Packit Service |
6754ca |
using namespace std;
|
|
Packit Service |
6754ca |
|
|
Packit Service |
6754ca |
//
|
|
Packit Service |
6754ca |
// Convert from linear to a nonlinear representation,
|
|
Packit Service |
6754ca |
//
|
|
Packit Service |
6754ca |
|
|
Packit Service |
6754ca |
half
|
|
Packit Service |
6754ca |
toNonlinear(half linear)
|
|
Packit Service |
6754ca |
{
|
|
Packit Service |
6754ca |
float sign = 1;
|
|
Packit Service |
6754ca |
float logBase = pow(2.7182818, 2.2);
|
|
Packit Service |
6754ca |
|
|
Packit Service |
6754ca |
if ((float)linear < 0) {
|
|
Packit Service |
6754ca |
sign = -1;
|
|
Packit Service |
6754ca |
}
|
|
Packit Service |
6754ca |
|
|
Packit Service |
6754ca |
if ((linear.bits() & 0x7c00) == 0x7c00) {
|
|
Packit Service |
6754ca |
return (half)0.0;
|
|
Packit Service |
6754ca |
}
|
|
Packit Service |
6754ca |
|
|
Packit Service |
6754ca |
if ( fabs( (float)linear ) <= 1.0f) {
|
|
Packit Service |
6754ca |
return (half)(sign * pow(fabs((float)linear), 1.f/2.2f));
|
|
Packit Service |
6754ca |
} else {
|
|
Packit Service |
6754ca |
return (half)(sign * ( log(fabs((float)linear)) / log(logBase) + 1.0f) );
|
|
Packit Service |
6754ca |
}
|
|
Packit Service |
6754ca |
|
|
Packit Service |
6754ca |
return (half)0.0f;
|
|
Packit Service |
6754ca |
}
|
|
Packit Service |
6754ca |
|
|
Packit Service |
6754ca |
|
|
Packit Service |
6754ca |
void
|
|
Packit Service |
6754ca |
compareDwa(int width,
|
|
Packit Service |
6754ca |
int height,
|
|
Packit Service |
6754ca |
const Array2D<Rgba> &src,
|
|
Packit Service |
6754ca |
const Array2D<Rgba> &test,
|
|
Packit Service |
6754ca |
RgbaChannels channels)
|
|
Packit Service |
6754ca |
{
|
|
Packit Service |
6754ca |
half srcNonlin, testNonlin;
|
|
Packit Service |
6754ca |
float relError;
|
|
Packit Service |
6754ca |
|
|
Packit Service |
6754ca |
for (int y=0; y
|
|
Packit Service |
6754ca |
for (int x=0; x
|
|
Packit Service |
6754ca |
|
|
Packit Service |
6754ca |
for (int comp=0; comp<3; ++comp) {
|
|
Packit Service |
6754ca |
switch (comp) {
|
|
Packit Service |
6754ca |
case 0:
|
|
Packit Service |
6754ca |
if (!(channels & WRITE_R)) continue;
|
|
Packit Service |
6754ca |
|
|
Packit Service |
6754ca |
srcNonlin = toNonlinear(src[y][x].r);
|
|
Packit Service |
6754ca |
testNonlin = toNonlinear(test[y][x].r);
|
|
Packit Service |
6754ca |
break;
|
|
Packit Service |
6754ca |
case 1:
|
|
Packit Service |
6754ca |
if (!(channels & WRITE_G)) continue;
|
|
Packit Service |
6754ca |
|
|
Packit Service |
6754ca |
srcNonlin = toNonlinear(src[y][x].g);
|
|
Packit Service |
6754ca |
testNonlin = toNonlinear(test[y][x].g);
|
|
Packit Service |
6754ca |
break;
|
|
Packit Service |
6754ca |
case 2:
|
|
Packit Service |
6754ca |
if (!(channels & WRITE_B)) continue;
|
|
Packit Service |
6754ca |
|
|
Packit Service |
6754ca |
srcNonlin = toNonlinear(src[y][x].b);
|
|
Packit Service |
6754ca |
testNonlin = toNonlinear(test[y][x].b);
|
|
Packit Service |
6754ca |
break;
|
|
Packit Service |
6754ca |
}
|
|
Packit Service |
6754ca |
|
|
Packit Service |
6754ca |
//
|
|
Packit Service |
6754ca |
// Try to compare with relative error. This breaks down
|
|
Packit Service |
6754ca |
// for small numbers, which could be quantiezed to 0
|
|
Packit Service |
6754ca |
// giving 100% error.
|
|
Packit Service |
6754ca |
//
|
|
Packit Service |
6754ca |
if (srcNonlin.bits() != 0x00) {
|
|
Packit Service |
6754ca |
|
|
Packit Service |
6754ca |
relError = fabs( (float)srcNonlin - (float)testNonlin ) /
|
|
Packit Service |
6754ca |
fabs((float)srcNonlin);
|
|
Packit Service |
6754ca |
|
|
Packit Service |
6754ca |
|
|
Packit Service |
6754ca |
if (fabs(srcNonlin) < .1) continue;
|
|
Packit Service |
6754ca |
|
|
Packit Service |
6754ca |
if (fabs(srcNonlin) < .25) {
|
|
Packit Service |
6754ca |
assert( relError < .25);
|
|
Packit Service |
6754ca |
} else {
|
|
Packit Service |
6754ca |
assert( relError < .1);
|
|
Packit Service |
6754ca |
}
|
|
Packit Service |
6754ca |
|
|
Packit Service |
6754ca |
} else {
|
|
Packit Service |
6754ca |
assert( srcNonlin != testNonlin );
|
|
Packit Service |
6754ca |
}
|
|
Packit Service |
6754ca |
}
|
|
Packit Service |
6754ca |
|
|
Packit Service |
6754ca |
}
|
|
Packit Service |
6754ca |
}
|
|
Packit Service |
6754ca |
|
|
Packit Service |
6754ca |
|
|
Packit Service |
6754ca |
//
|
|
Packit Service |
6754ca |
// Test alpha, if necessary
|
|
Packit Service |
6754ca |
//
|
|
Packit Service |
6754ca |
if (channels & WRITE_A) {
|
|
Packit Service |
6754ca |
for (int y=0; y
|
|
Packit Service |
6754ca |
for (int x=0; x
|
|
Packit Service |
6754ca |
assert( src[y][x].a == test[y][x].a );
|
|
Packit Service |
6754ca |
}
|
|
Packit Service |
6754ca |
}
|
|
Packit Service |
6754ca |
}
|
|
Packit Service |
6754ca |
}
|