contracted_tensor_check.py
No OneTemporary
Actions

Subscribers

None

File Metadata


	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')