contracted_tensor_check.py
55115a530536master

contracted_tensor_check.py
View Options


			from os import listdir, linesep
			from os.path import isfile, join

			CTcsv = '/home/sernst/tensor_contractor/bin/contracted_tensors.csv'
			CTcsv_ref = '/home/sernst/tensor_contractor/bin/tests/contracted_tensors_mathematica.csv'

			refs = {};

			with open(CTcsv_ref,'r') as file:
			for line in file.readlines():
			splited = line.split(',')
			p = splited[0]
			u = splited[1]
			pp = splited[2]
			val = float(splited[3])
			refs[(p,u,pp)] = val;


			max_error = 0;
			count = 0;
			with open(CTcsv,'r') as file:
			for line in file.readlines():
			splited = line.split(',')
			p = splited[0]
			u = splited[1]
			pp = splited[2]
			val = float(splited[3])
			error = abs(val - refs[(p,u,pp)]);
			count += 1
			max_error = max(max_error, error)
			#print('Error: %e' % error)

			print('Max error:', max_error);
			print(len(refs), 'in CTcsv_ref')
			print(count, 'in CTcsv')