def check_uv_primes_inverse_1(self, xyz0, verbose):
''' Check that uv_primes() and uv_primes_inverse() are inverses. '''
up0, vp0 = colormodels.uv_primes(xyz0)
xyz1 = colormodels.uv_primes_inverse(up0, vp0, xyz0[1])
# Check error.
dxyz = (xyz1 - xyz0)
error = math.sqrt(numpy.dot(dxyz, dxyz))
tolerance = 1.0e-13
self.assertLessEqual(error, tolerance)
msg = 'xyz0: %s up: %g vp: %g xyz(up,vp): %s Error: %g' % (
str(xyz0), up0, vp0, str(xyz1), error)
if verbose:
print(msg)
def check_uv_primes_inverse_1(self, xyz0, verbose):
''' Check that uv_primes() and uv_primes_inverse() are inverses. '''
up0, vp0 = colormodels.uv_primes (xyz0)
xyz1 = colormodels.uv_primes_inverse (up0, vp0, xyz0[1])
# Check error.
dxyz = (xyz1 - xyz0)
error = math.sqrt (numpy.dot (dxyz, dxyz))
tolerance = 1.0e-13
self.assertLessEqual(error, tolerance)
msg = 'xyz0: %s up: %g vp: %g xyz(up,vp): %s Error: %g' % (
str (xyz0), up0, vp0, str(xyz1), error)
if verbose:
print (msg)
def check_uv_primes_inverse_2(self, up0, vp0, y0, verbose):
''' Check that uv_primes() and uv_primes_inverse() are inverses. '''
xyz0 = colormodels.uv_primes_inverse (up0, vp0, y0)
up1, vp1 = colormodels.uv_primes (xyz0)
# Check error.
error_up = up1 - up0
error_vp = vp1 - vp0
error = numpy.hypot (error_up, error_vp)
tolerance = 1.0e-13
self.assertLessEqual(error, tolerance)
msg = 'up0, vp0, y0: %g, %g, %g xyz(up0,vp0,y0): %s (up,vp)(xyz): %g, %g Error: %g' % (
up0, vp0, y0, str (xyz0), up1, vp1, error)
if verbose:
print (msg)
def check_uv_primes_inverse_2(self, up0, vp0, y0, verbose):
''' Check that uv_primes() and uv_primes_inverse() are inverses. '''
xyz0 = colormodels.uv_primes_inverse(up0, vp0, y0)
up1, vp1 = colormodels.uv_primes(xyz0)
# Check error.
error_up = up1 - up0
error_vp = vp1 - vp0
error = numpy.hypot(error_up, error_vp)
tolerance = 1.0e-13
self.assertLessEqual(error, tolerance)
msg = 'up0, vp0, y0: %g, %g, %g xyz(up0,vp0,y0): %s (up,vp)(xyz): %g, %g Error: %g' % (
up0, vp0, y0, str(xyz0), up1, vp1, error)
if verbose:
print(msg)
def test_A(xyz0, tolerance=0.0, verbose=1):
(up0, vp0) = colormodels.uv_primes(xyz0)
xyz1 = colormodels.uv_primes_inverse(up0, vp0, xyz0[1])
# check error
dxyz = (xyz1 - xyz0)
error = math.sqrt(numpy.dot(dxyz, dxyz))
passed = (error <= tolerance)
if passed:
status = 'pass'
else:
status = 'FAILED'
msg = 'test_uv_primes.test_A() : xyz0 = %s, (up,vp) = (%g,%g), xyz(up,vp) = %s, error = %g, %s' % (
str(xyz0), up0, vp0, str(xyz1), error, status)
if verbose >= 1:
print(msg)
if not passed:
pass
raise ValueError(msg)
return passed
def test_A (xyz0, tolerance=0.0, verbose=1):
(up0, vp0) = colormodels.uv_primes (xyz0)
xyz1 = colormodels.uv_primes_inverse (up0, vp0, xyz0[1])
# check error
dxyz = (xyz1 - xyz0)
error = math.sqrt (numpy.dot (dxyz, dxyz))
passed = (error <= tolerance)
if passed:
status = 'pass'
else:
status = 'FAILED'
msg = 'test_uv_primes.test_A() : xyz0 = %s, (up,vp) = (%g,%g), xyz(up,vp) = %s, error = %g, %s' % (
str (xyz0), up0, vp0, str(xyz1), error, status)
if verbose >= 1:
print msg
if not passed:
pass
raise ValueError, msg
return passed
def test_B(up0, vp0, y0, tolerance=0.0, verbose=1):
xyz0 = colormodels.uv_primes_inverse(up0, vp0, y0)
(up1, vp1) = colormodels.uv_primes(xyz0)
# check error
error_up = up1 - up0
error_vp = vp1 - vp0
error = numpy.hypot(error_up, error_vp)
passed = (error <= tolerance)
if passed:
status = 'pass'
else:
status = 'FAILED'
msg = 'test_uv_primes.test_B() : (up0,vp0,y0) = (%g,%g,%g), xyz (up0,vp0,y0) = %s, (up,vp)(xyz) = (%g,%g), error = %g, %s' % (
up0, vp0, y0, str(xyz0), up1, vp1, error, status)
if verbose >= 1:
print(msg)
if not passed:
pass
raise ValueError(msg)
return passed
def test_B (up0, vp0, y0, tolerance=0.0, verbose=1):
xyz0 = colormodels.uv_primes_inverse (up0, vp0, y0)
(up1, vp1) = colormodels.uv_primes (xyz0)
# check error
error_up = up1 - up0
error_vp = vp1 - vp0
error = numpy.hypot (error_up, error_vp)
passed = (error <= tolerance)
if passed:
status = 'pass'
else:
status = 'FAILED'
msg = 'test_uv_primes.test_B() : (up0,vp0,y0) = (%g,%g,%g), xyz (up0,vp0,y0) = %s, (up,vp)(xyz) = (%g,%g), error = %g, %s' % (
up0, vp0, y0, str (xyz0), up1, vp1, error, status)
if verbose >= 1:
print msg
if not passed:
pass
raise ValueError, msg
return passed
