marchwicka/signature_function
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

refactoring

Browse Source
This commit is contained in:
Maria Marchwicka 2020-09-08 17:28:56 +02:00
parent ac1a1d3475
commit 1824b075e9
2 changed files with 114 additions and 186 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

290
cable_signature.sage
View File

@ -4,10 +4,13 @@ import numpy as np
import itertools as it import itertools as it
SIGNATURE = 0
SIGMA = 1
class TorusCable(object): class TorusCable(object):
def __init__(self, knot_formula, k_vector=None, q_vector=None): def __init__(self, knot_formula, k_vector=None, q_vector=None):
# q_i = 2 * k_i + 1 # q_i = 2 * k_i + 1
if k_vector is None: if k_vector is None:
if q_vector is None: if q_vector is None:
# TBD docstring # TBD docstring
@ -27,8 +30,39 @@ class TorusCable(object):
self.__sigma_function = None self.__sigma_function = None
self.__signature_as_function_of_theta = None self.__signature_as_function_of_theta = None
def __get_sigma_function(self):
k_1, k_2, k_3, k_4 = [abs(k) for k in self.k_vector]
last_q = 2 * k_4 + 1
ksi = 1/last_q
sigma_q_1 = self.get_untwisted_signature_function(k_1)
sigma_q_2 = self.get_untwisted_signature_function(k_2)
sigma_q_3 = self.get_untwisted_signature_function(k_3)
def get_untwisted_signature_function(self, j, q=None): def sigma_function(theta_vector, print_results=False):
# "untwisted" part (Levine-Tristram signatures)
a_1, a_2, a_3, a_4 = theta_vector
untwisted_part = 2 * (sigma_q_2(ksi * a_1) -
sigma_q_2(ksi * a_2) +
sigma_q_3(ksi * a_3) -
sigma_q_3(ksi * a_4) +
sigma_q_1(ksi * a_1 * 2) -
sigma_q_1(ksi * a_4 * 2))
# "twisted" part
tp = [0, 0, 0, 0]
for i, a in enumerate(theta_vector):
if a:
tp[i] = -last_q + 2 * a - 2 * (a^2/last_q)
twisted_part = tp[0] - tp[1] + tp[2] - tp[3]
# if print_results:
# self.print_results_LT(theta_vector, untwisted_part)
# self.print_results_LT(theta_vector, twisted_part)
sigma_v = untwisted_part + twisted_part
return sigma_v
return sigma_function
@staticmethod
def get_untwisted_signature_function(j, q=None):
# return the signature function of the T_{2,2k+1} torus knot # return the signature function of the T_{2,2k+1} torus knot
k = abs(j) k = abs(j)
q = 2 * k + 1 q = 2 * k + 1
@ -37,7 +71,6 @@ class TorusCable(object):
for a in range(k + 1, 2 * k + 1)]) for a in range(k + 1, 2 * k + 1)])
return SignatureFunction(values=w) return SignatureFunction(values=w)
def get_knot_descrption(self): def get_knot_descrption(self):
description = "" description = ""
for knot in self.knot_sum: for knot in self.knot_sum:
@ -116,30 +149,16 @@ class TorusCable(object):
number_of_null_comb += 2^m number_of_null_comb += 2^m
return number_of_null_comb, list_of_good_vectors return number_of_null_comb, list_of_good_vectors
# searching for sigma > 5 + #(v_i != 0)
def eval_cable_for_large_signature(self, list_of_ranges, def __check_all_combinations_in_ranges(self, list_of_ranges,
print_results=False, sigma_or_sign,
verbose=False): print_results=False):
if self.__signature_as_function_of_theta is None:
self.__signature_as_function_of_theta= \
self.__get_signature_as_function_of_theta()
if print_results:
print()
print(self.knot_description, end="\t\t\t")
print()
f = self.__signature_as_function_of_theta
if self.s__check_all_combinations_in_ranges(list_of_ranges,
print_results=print_results):
return True
return False
def s__check_all_combinations_in_ranges(self, list_of_ranges,
print_results=False):
all_combinations_pass = True all_combinations_pass = True
all_bad_vectors = [] all_bad_vectors = []
number_of_all_good_v = 0 number_of_all_good_v = 0
for i, range_product in enumerate(list_of_ranges): for i, range_product in enumerate(list_of_ranges):
good_v, bad_v = self.s__check_combinations_in_range(range_product) good_v, bad_v = self.__check_combinations_in_range(range_product,
sigma_or_sign)
number_of_all_good_v += len(good_v) number_of_all_good_v += len(good_v)
all_bad_vectors = list(it.chain(all_bad_vectors, bad_v)) all_bad_vectors = list(it.chain(all_bad_vectors, bad_v))
if bad_v: if bad_v:
@ -157,53 +176,27 @@ class TorusCable(object):
" : " + str(len(all_bad_vectors))) " : " + str(len(all_bad_vectors)))
return all_combinations_pass return all_combinations_pass
def s__check_combinations_in_range(self, range_product): # searching for signature or sigma > 5 + #(v_i != 0)
large_sigma_for_all_combinations = True def __check_combinations_in_range(self, range_product, sigma_or_sign):
bad_vectors = [] bad_vectors = []
good_vectors = [] good_vectors = []
q_4 = self.q_vector[-1] last_q = self.q_vector[-1]
for vector in range_product: for vector in range_product:
a_1, a_2, a_3, a_4 = vector a_1, a_2, a_3, a_4 = vector
if (a_1^2 - a_2^2 + a_3^2 - a_4^2) % q_4: if (a_1^2 - a_2^2 + a_3^2 - a_4^2) % last_q:
continue continue
if all(a in [1, q_4 - 1] for a in vector): # if all(a in [1, last_q - 1] for a in vector):
pass # pass
else: # else:
continue # continue
if self.s__is_sigma_for_vector_class_big(vector): if self.__is_sigma_for_vector_class_big(vector, sigma_or_sign):
good_vectors.append(vector) good_vectors.append(vector)
else: else:
# print(vector) # print(vector)
bad_vectors.append(vector) bad_vectors.append(vector)
return good_vectors, bad_vectors return good_vectors, bad_vectors
def s__is_sigma_for_vector_class_big(self, theta_vector):
[a_1, a_2, a_3, a_4] = theta_vector
q_4 = self.q_vector[-1]
k_4 = self.k_vector[-1]
max_sigma = 0
print(theta_vector, end="\t")
for shift in range(1, k_4 + 1):
shifted_theta = [(shift * a) % q_4 for a in
[a_1, a_2, a_3, a_4]]
sf = self.__signature_as_function_of_theta(shifted_theta)
sigma_v = sf.is_big()
print(sigma_v, end=" ")
if abs(sigma_v) > abs(max_sigma):
max_sigma = sigma_v
if abs(sigma_v) > 5 + np.count_nonzero(shifted_theta):
print("\tok " + str(sigma_v))
return True
print("\tbad class " + str(max_sigma))
return False
# searching for signature == 0 # searching for signature == 0
def eval_cable_for_null_signature(self, print_results=False, verbose=False): def eval_cable_for_null_signature(self, print_results=False, verbose=False):
# search for zero combinations # search for zero combinations
@ -224,17 +217,35 @@ class TorusCable(object):
return number_of_null_comb, number_of_all_comb return number_of_null_comb, number_of_all_comb
return None return None
# check sigma for all v = s * [a_1, a_2, a_3, a_4] for s in [1, q_4 - 1] # check sigma for all v = s * [a_1, a_2, a_3, a_4] for s in [1, last_q - 1]
def __is_sigma_for_vector_class_big(self, theta_vector): def __is_sigma_for_vector_class_big(self, theta_vector, sigma_or_sign):
[a_1, a_2, a_3, a_4] = theta_vector [a_1, a_2, a_3, a_4] = theta_vector
q_4 = self.q_vector[-1] q_4 = self.q_vector[-1]
k_4 = self.k_vector[-1] k_4 = self.k_vector[-1]
max_sigma = 0
if sigma_or_sign == SIGNATURE:
f = self.__signature_as_function_of_theta
else:
f = self.__sigma_function
print(theta_vector, end="\t")
for shift in range(1, k_4 + 1): for shift in range(1, k_4 + 1):
shifted_theta = [(shift * a) % q_4 for a in shifted_theta = [(shift * a) % q_4 for a in
[a_1, a_2, a_3, a_4]] [a_1, a_2, a_3, a_4]]
sigma_v = self.__sigma_function(shifted_theta) if sigma_or_sign == SIGNATURE:
sf = f(shifted_theta)
sigma_v = sf.is_big()
else:
sigma_v = f(shifted_theta)
print(sigma_v, end=" ")
if abs(sigma_v) > abs(max_sigma):
max_sigma = sigma_v
if abs(sigma_v) > 5 + np.count_nonzero(shifted_theta): if abs(sigma_v) > 5 + np.count_nonzero(shifted_theta):
print("\tok " + str(sigma_v))
return True return True
print("\tbad class " + str(max_sigma))
return False return False
def __tmp_print_all_sigma_for_vector_class(self, theta_vector): def __tmp_print_all_sigma_for_vector_class(self, theta_vector):
@ -242,9 +253,9 @@ class TorusCable(object):
print(self.knot_description) print(self.knot_description)
print("vector = " + str(theta_vector)) print("vector = " + str(theta_vector))
[a_1, a_2, a_3, a_4] = theta_vector [a_1, a_2, a_3, a_4] = theta_vector
q_4 = self.q_vector[3] last_q = self.q_vector[-1]
for shift in range(1, q_4): for shift in range(1, last_q):
shifted_theta = [(shift * a) % q_4 for a in shifted_theta = [(shift * a) % last_q for a in
[a_1, a_2, a_3, a_4]] [a_1, a_2, a_3, a_4]]
print(str(shifted_theta) + "\t\t" + \ print(str(shifted_theta) + "\t\t" + \
str(self.__sigma_function(shifted_theta))) str(self.__sigma_function(shifted_theta)))
@ -253,9 +264,9 @@ class TorusCable(object):
def __tmp_get_max_sigma_for_vector_class(self, theta_vector): def __tmp_get_max_sigma_for_vector_class(self, theta_vector):
max_sigma = (theta_vector, 0) max_sigma = (theta_vector, 0)
[a_1, a_2, a_3, a_4] = theta_vector [a_1, a_2, a_3, a_4] = theta_vector
q_4 = self.q_vector[3] last_q = self.q_vector[-1]
for shift in range(1, q_4): for shift in range(1, last_q):
shifted_theta = [(shift * a) % q_4 for a in shifted_theta = [(shift * a) % last_q for a in
[a_1, a_2, a_3, a_4]] [a_1, a_2, a_3, a_4]]
sigma = self.__sigma_function(shifted_theta) sigma = self.__sigma_function(shifted_theta)
if abs(sigma) > abs(max_sigma[1]): if abs(sigma) > abs(max_sigma[1]):
@ -265,45 +276,14 @@ class TorusCable(object):
def is_sigma_for_vector_class_big(self, theta_vector): def is_sigma_for_vector_class_big(self, theta_vector):
if self.__sigma_function is None: if self.__sigma_function is None:
self.__sigma_function = self.__get_sigma_function() self.__sigma_function = self.__get_sigma_function()
return self.__is_sigma_for_vector_class_big(theta_vector) return self.__is_sigma_for_vector_class_big(theta_vector, SIGMA)
def __get_sigma_function(self):
k_1, k_2, k_3, k_4 = [abs(k) for k in self.k_vector]
q_4 = 2 * k_4 + 1
ksi = 1/q_4
sigma_q_1 = self.get_untwisted_signature_function(k_1)
sigma_q_2 = self.get_untwisted_signature_function(k_2)
sigma_q_3 = self.get_untwisted_signature_function(k_3)
def sigma_function(theta_vector, print_results=False):
# "untwisted" part (Levine-Tristram signatures)
a_1, a_2, a_3, a_4 = theta_vector
untwisted_part = 2 * (sigma_q_2(ksi * a_1) -
sigma_q_2(ksi * a_2) +
sigma_q_3(ksi * a_3) -
sigma_q_3(ksi * a_4) +
sigma_q_1(ksi * a_1 * 2) -
sigma_q_1(ksi * a_4 * 2))
# "twisted" part
tp = [0, 0, 0, 0]
for i, a in enumerate(theta_vector):
if a:
tp[i] = -q_4 + 2 * a - 2 * (a^2/q_4)
twisted_part = tp[0] - tp[1] + tp[2] - tp[3]
# if print_results:
# self.print_results_LT(theta_vector, untwisted_part)
# self.print_results_LT(theta_vector, twisted_part)
sigma_v = untwisted_part + twisted_part
return sigma_v
return sigma_function
def print_results_LT(self, theta_vector, untwisted_part): def print_results_LT(self, theta_vector, untwisted_part):
knot_description = self.knot_description knot_description = self.knot_description
k_1, k_2, k_3, k_4 = [abs(k) for k in self.k_vector] k_1, k_2, k_3, k_4 = [abs(k) for k in self.k_vector]
a_1, a_2, a_3, a_4 = theta_vector a_1, a_2, a_3, a_4 = theta_vector
q_4 = 2 * k_4 + 1 last_q = 2 * k_4 + 1
ksi = 1/q_4 ksi = 1/last_q
sigma_q_1 = self.get_untwisted_signature_function(k_1) sigma_q_1 = self.get_untwisted_signature_function(k_1)
sigma_q_2 = self.get_untwisted_signature_function(k_2) sigma_q_2 = self.get_untwisted_signature_function(k_2)
sigma_q_3 = self.get_untwisted_signature_function(k_3) sigma_q_3 = self.get_untwisted_signature_function(k_3)
@ -378,7 +358,7 @@ class TorusCable(object):
def print_results_sigma(self, theta_vector, twisted_part): def print_results_sigma(self, theta_vector, twisted_part):
a_1, a_2, a_3, a_4 = theta_vector a_1, a_2, a_3, a_4 = theta_vector
knot_description = self.knot_description knot_description = self.knot_description
q_4 = self.q_vector[-1] last_q = self.q_vector[-1]
print("\n\nSigma values for the cable sum: ") print("\n\nSigma values for the cable sum: ")
print(knot_description) print(knot_description)
print("and characters: " + str(v_theta)) print("and characters: " + str(v_theta))
@ -388,29 +368,29 @@ class TorusCable(object):
" = 0 if a == 0.") " = 0 if a == 0.")
print("\nsigma(T_{2, q_4}, chi_a_1) = ", end="") print("\nsigma(T_{2, q_4}, chi_a_1) = ", end="")
if a_1: if a_1:
print("- (" + str(q_4) + ") + 2 * " + str(a_1) + " + " +\ print("- (" + str(last_q) + ") + 2 * " + str(a_1) + " + " +\
"- 2 * " + str(a_1^2) + "/" + str(q_4) + \ "- 2 * " + str(a_1^2) + "/" + str(last_q) + \
" = " + str(tp[0])) " = " + str(tp[0]))
else: else:
print("0") print("0")
print("\nsigma(T_{2, q_4}, chi_a_2) = ", end ="") print("\nsigma(T_{2, q_4}, chi_a_2) = ", end ="")
if a_2: if a_2:
print("- (" + str(q_4) + ") + 2 * " + str(a_2) + " + " +\ print("- (" + str(last_q) + ") + 2 * " + str(a_2) + " + " +\
"- 2 * " + str(a_2^2) + "/" + str(q_4) + \ "- 2 * " + str(a_2^2) + "/" + str(last_q) + \
" = " + str(tp[1])) " = " + str(tp[1]))
else: else:
print("0", end="") print("0", end="")
print("\nsigma(T_{2, q_4}, chi_a_3) = ", end="") print("\nsigma(T_{2, q_4}, chi_a_3) = ", end="")
if a_3: if a_3:
print("- (" + str(q_4) + ") + 2 * " + str(a_3) + " + " +\ print("- (" + str(last_q) + ") + 2 * " + str(a_3) + " + " +\
"- 2 * " + str(a_3^2) + "/" + str(q_4) + \ "- 2 * " + str(a_3^2) + "/" + str(last_q) + \
" = " + str(tp[2])) " = " + str(tp[2]))
else: else:
print("0", end="") print("0", end="")
print("\nsigma(T_{2, q_4}, chi_a_4) = ", end="") print("\nsigma(T_{2, q_4}, chi_a_4) = ", end="")
if a_4: if a_4:
print("- (" + str(q_4) + ") + 2 * " + str(a_4) + " + " +\ print("- (" + str(last_q) + ") + 2 * " + str(a_4) + " + " +\
"- 2 * " + str(a_4^2) + "/" + str(q_4) + \ "- 2 * " + str(a_4^2) + "/" + str(last_q) + \
" = " + str(tp[3])) " = " + str(tp[3]))
else: else:
print("0") print("0")
@ -431,68 +411,34 @@ class TorusCable(object):
self.__sigma_function = self.__get_sigma_function() self.__sigma_function = self.__get_sigma_function()
return self.__sigma_function(theta_vector) return self.__sigma_function(theta_vector)
# searching for sigma > 5 + #(v_i != 0)
def __check_combinations_in_range(self, range_product):
bad_vectors = []
good_vectors = []
q_4 = self.q_vector[-1]
for vector in range_product:
a_1, a_2, a_3, a_4 = vector
if (a_1^2 - a_2^2 + a_3^2 - a_4^2) % q_4:
continue
# if all(a in [1, q_4 - 1] for a in vector):
# continue
if self.__is_sigma_for_vector_class_big(vector):
good_vectors.append(vector)
else:
bad_vectors.append(vector)
return good_vectors, bad_vectors
# searching for sigma > 5 + #(v_i != 0) # searching for sigma > 5 + #(v_i != 0)
def check_combinations_in_range(self, range_product): def check_combinations_in_range(self, range_product):
if self.__sigma_function is None: if self.__sigma_function is None:
self.__sigma_function = self.__get_sigma_function() self.__sigma_function = self.__get_sigma_function()
return self.__check_combinations_in_range(range_product) return self.__check_combinations_in_range(range_product, SIGMA)
# searching for sigma > 5 + #(v_i != 0) # searching for sigma > 5 + #(v_i != 0)
def __check_all_combinations_in_ranges(self, list_of_ranges, def eval_cable_for_large_values(self, list_of_ranges,
print_results=False): sigma_or_sign,
all_combinations_pass = True print_results=False,
all_bad_vectors = [] verbose=False):
number_of_all_good_v = 0
for i, range_product in enumerate(list_of_ranges):
good_v, bad_v = self.__check_combinations_in_range(range_product)
number_of_all_good_v += len(good_v)
all_bad_vectors = list(it.chain(all_bad_vectors, bad_v))
if bad_v:
all_combinations_pass = False
# if print_results:
# print("good : bad:\t " + str(len(good_v)) +\
# " : " + str(len(bad_v)))
# if i in [0, 4,]:
# print()
# if bad_v:
# print(bad_v)
if print_results:
print("good : bad:\t " + str(number_of_all_good_v) +\
" : " + str(len(all_bad_vectors)))
return all_combinations_pass
# searching for sigma > 5 + #(v_i != 0)
def eval_cable_for_large_sigma(self, list_of_ranges,
print_results=False, verbose=False):
if self.__sigma_function is None:
self.__sigma_function = self.__get_sigma_function()
if print_results: if print_results:
print(self.knot_description, end="\t\t\t") print(self.knot_description, end="\t\t\t")
if sigma_or_sign == SIGMA:
if self.__sigma_function is None:
self.__sigma_function = self.__get_sigma_function()
else:
if self.__signature_as_function_of_theta is None:
self.__signature_as_function_of_theta= \
self.__get_signature_as_function_of_theta()
if self.__check_all_combinations_in_ranges(list_of_ranges, if self.__check_all_combinations_in_ranges(list_of_ranges,
sigma_or_sign,
print_results=print_results): print_results=print_results):
return True return True
return False return False
class SignatureFunction(object): class SignatureFunction(object):
def __init__(self, values=None, counter=None): def __init__(self, values=None, counter=None):
@ -586,7 +532,8 @@ class SignatureFunction(object):
max = 0 max = 0
items = self.cnt_signature_jumps.items() items = self.cnt_signature_jumps.items()
for arg, _ in items: for arg, _ in items:
current = sum([jump for jump_arg, jump in items if jump_arg <= arg]) # current = sum([jump for jump_arg, jump in items if jump_arg <= arg])
current = self(arg)
if abs(current) > abs(max): if abs(current) > abs(max):
max = current max = current
if abs(max) > 9: if abs(max) > 9:
@ -594,22 +541,10 @@ class SignatureFunction(object):
return max return max
def mod_one(n): def mod_one(n):
return n - floor(n) return n - floor(n)
def get_untwisted_signature_function(j):
# return the signature function of the T_{2,2k+1} torus knot
k = abs(j)
w = ([((2 * a + 1)/(4 * k + 2), -1 * sgn(j)) for a in range(k)] +
[((2 * a + 1)/(4 * k + 2), 1 * sgn(j))
for a in range(k + 1, 2 * k + 1)])
return SignatureFunction(values=w)
def get_summand_signature_as_theta_function(*arg): def get_summand_signature_as_theta_function(*arg):
def get_summand_signture_function(theta): def get_summand_signture_function(theta):
# TBD: another formula (for t^2) description # TBD: another formula (for t^2) description
@ -628,7 +563,7 @@ def get_summand_signature_as_theta_function(*arg):
for i, k in enumerate(arg[:-1][::-1]): for i, k in enumerate(arg[:-1][::-1]):
ksi = 1/(2 * k_n + 1) ksi = 1/(2 * k_n + 1)
power = 2^i power = 2^i
a = get_untwisted_signature_function(k) a = TorusCable.get_untwisted_signature_function(k)
shift = theta * ksi * power shift = theta * ksi * power
b = a >> shift b = a >> shift
c = a << shift c = a << shift
@ -709,15 +644,6 @@ def get_summand_signature_as_theta_function(*arg):
return get_summand_signture_function return get_summand_signture_function
def get_untwisted_signature_function(j):
# return the signature function of the T_{2,2k+1} torus knot
k = abs(j)
w = ([((2 * a + 1)/(4 * k + 2), -1 * sgn(j)) for a in range(k)] +
[((2 * a + 1)/(4 * k + 2), 1 * sgn(j))
for a in range(k + 1, 2 * k + 1)])
return SignatureFunction(values=w)
TorusCable.get_number_of_combinations_of_theta.__doc__ = \ TorusCable.get_number_of_combinations_of_theta.__doc__ = \
""" """
Arguments: Arguments:

10
my_signature.sage
View File

@ -7,11 +7,11 @@ import sys
import itertools as it import itertools as it
import re import re
#
# if not os.path.isfile('cable_signature.py'): # if not os.path.isfile('cable_signature.py'):
# os.system('sage --preparse cable_signature.sage') # os.system('sage --preparse cable_signature.sage')
# os.system('mv cable_signature.sage.py cable_signature.py') # os.system('mv cable_signature.sage.py cable_signature.py')
# from cable_signature import SignatureFunction, TorusCable # from cable_signature import SignatureFunction, TorusCable, SIGNATURE, SIGMA
class Config(object): class Config(object):
@ -132,7 +132,8 @@ def search_for_large_sigma_value(knot_formula=None, limit=None,
continue continue
cable = TorusCable(knot_formula=knot_formula, q_vector=q) cable = TorusCable(knot_formula=knot_formula, q_vector=q)
list_of_ranges = config.get_list_of_ranges(cable.k_vector[-1] + 1) list_of_ranges = config.get_list_of_ranges(cable.k_vector[-1] + 1)
if cable.eval_cable_for_large_sigma(list_of_ranges, verbose=verbose, if cable.eval_cable_for_large_values(list_of_ranges, SIGMA,
verbose=verbose,
print_results=print_results): print_results=print_results):
good_knots.append(cable.knot_description) good_knots.append(cable.knot_description)
return good_knots return good_knots
@ -187,7 +188,8 @@ def search_for_large_signature_value(knot_formula=None, limit=None,
continue continue
cable = TorusCable(knot_formula=knot_formula, q_vector=q) cable = TorusCable(knot_formula=knot_formula, q_vector=q)
list_of_ranges = config.get_list_of_ranges(cable.q_vector[-1]) list_of_ranges = config.get_list_of_ranges(cable.q_vector[-1])
if cable.eval_cable_for_large_signature(list_of_ranges, verbose=verbose, if cable.eval_cable_for_large_values(list_of_ranges, SIGNATURE,
verbose=verbose,
print_results=print_results): print_results=print_results):
good_knots.append(cable.knot_description) good_knots.append(cable.knot_description)