swh:1:snp:d2bcff616bbf538fe8ce2a9c384200307730292a
Tip revision: 2f07c209fd34606453c6d5275423fb5e5cb3e91d authored by Joachim van der Herten on 16 August 2017, 13:12:14 UTC
Merge branch 'master' into nested_models_recompilation
Merge branch 'master' into nested_models_recompilation
Tip revision: 2f07c20
test_triang.py
import unittest
from GPflow.tf_wraps import vec_to_tri
import tensorflow as tf
import numpy as np
class TestVecToTri(unittest.TestCase):
def setUp(self):
tf.reset_default_graph()
def referenceInverse(self, matrices):
#this is the inverse operation of the vec_to_tri
#op being tested.
D, N, ignored = matrices.shape
M = (N * (N + 1)) // 2
tril_indices = np.tril_indices(N )
output = np.zeros( (D, M) )
for vector_index in range(D):
matrix = matrices[vector_index,:]
output[vector_index,:] = matrix[tril_indices]
return output
def getExampleMatrices(self, D, N ):
rng = np.random.RandomState(1)
random_matrices = rng.randn( D,N,N )
for matrix_index in range(D):
for row_index in range(N):
for col_index in range(N):
if col_index>row_index:
random_matrices[matrix_index,row_index,col_index] = 0.
return random_matrices
def testBasicFunctionality(self):
N = 3
D = 3
reference_matrices = self.getExampleMatrices(D,N)
input_vector_tensor = tf.constant(self.referenceInverse(reference_matrices))
test_matrices_tensor = vec_to_tri(input_vector_tensor,N)
test_matrices = tf.Session().run(test_matrices_tensor)
np.testing.assert_array_almost_equal( reference_matrices, test_matrices)
def testDifferentiable(self):
N = 3
D = 3
reference_matrices = self.getExampleMatrices(D,N)
input_vector_array = self.referenceInverse(reference_matrices)
input_vector_tensor = tf.constant(input_vector_array)
test_matrices_tensor = vec_to_tri(input_vector_tensor,N)
reduced_sum = tf.reduce_sum( test_matrices_tensor )
gradient = tf.gradients( reduced_sum, input_vector_tensor )[0]
reference_gradient = np.ones_like( input_vector_array )
test_gradient = tf.Session().run(gradient)
np.testing.assert_array_almost_equal( reference_gradient, test_gradient)
if __name__ == "__main__":
unittest.main()