Add files via upload
This commit is contained in:
Silviu Marian Udrescu
GitHub
parent
be88c45a8a
commit
f1f13db1b9
2 changed files with 266 additions and 355 deletions
|
|
@ -2,212 +2,141 @@ import numpy as np
|
||||||
import os
|
import os
|
||||||
from S_run_bf_polyfit import run_bf_polyfit
|
from S_run_bf_polyfit import run_bf_polyfit
|
||||||
|
|
||||||
def get_acos(pathdir,pathdir_write_to,filename,BF_try_time,BF_ops_file_type, PA, polyfit_deg=4):
|
def get_acos(pathdir,pathdir_write_to,filename,BF_try_time,BF_ops_file_type, PA, polyfit_deg=3):
|
||||||
try:
|
try:
|
||||||
os.mkdir(pathdir_write_to)
|
os.mkdir(pathdir_write_to)
|
||||||
except:
|
except:
|
||||||
pass
|
pass
|
||||||
|
data = np.loadtxt(pathdir+filename)
|
||||||
try:
|
try:
|
||||||
n_variables = np.loadtxt(pathdir+"%s" %filename, dtype='str').shape[1]-1
|
data[:,-1] = np.arccos(data[:,-1])
|
||||||
variables = np.loadtxt(pathdir+"%s" %filename, usecols=(0,))
|
np.savetxt(pathdir_write_to+filename,data)
|
||||||
for j in range(1,n_variables):
|
|
||||||
v = np.loadtxt(pathdir+"%s" %filename, usecols=(j,))
|
|
||||||
variables = np.column_stack((variables,v))
|
|
||||||
f_dependent = np.loadtxt(pathdir+"%s" %filename, usecols=(n_variables,))
|
|
||||||
|
|
||||||
dt = np.column_stack((variables,np.arccos(f_dependent)))
|
|
||||||
np.savetxt(pathdir_write_to+filename,dt)
|
|
||||||
|
|
||||||
PA = run_bf_polyfit(pathdir,pathdir_write_to,filename,BF_try_time,BF_ops_file_type, PA, polyfit_deg, "acos")
|
PA = run_bf_polyfit(pathdir,pathdir_write_to,filename,BF_try_time,BF_ops_file_type, PA, polyfit_deg, "acos")
|
||||||
|
|
||||||
except:
|
except:
|
||||||
return PA
|
return PA
|
||||||
|
|
||||||
return PA
|
return PA
|
||||||
|
|
||||||
def get_asin(pathdir,pathdir_write_to,filename,BF_try_time,BF_ops_file_type, PA, polyfit_deg=4):
|
def get_asin(pathdir,pathdir_write_to,filename,BF_try_time,BF_ops_file_type, PA, polyfit_deg=3):
|
||||||
try:
|
try:
|
||||||
os.mkdir(pathdir_write_to)
|
os.mkdir(pathdir_write_to)
|
||||||
except:
|
except:
|
||||||
pass
|
pass
|
||||||
|
data = np.loadtxt(pathdir+filename)
|
||||||
try:
|
try:
|
||||||
n_variables = np.loadtxt(pathdir+"%s" %filename, dtype='str').shape[1]-1
|
data[:,-1] = np.arcsin(data[:,-1])
|
||||||
variables = np.loadtxt(pathdir+"%s" %filename, usecols=(0,))
|
np.savetxt(pathdir_write_to+filename,data)
|
||||||
for j in range(1,n_variables):
|
|
||||||
v = np.loadtxt(pathdir+"%s" %filename, usecols=(j,))
|
|
||||||
variables = np.column_stack((variables,v))
|
|
||||||
f_dependent = np.loadtxt(pathdir+"%s" %filename, usecols=(n_variables,))
|
|
||||||
|
|
||||||
dt = np.column_stack((variables,np.arcsin(f_dependent)))
|
|
||||||
np.savetxt(pathdir_write_to+filename,dt)
|
|
||||||
|
|
||||||
PA = run_bf_polyfit(pathdir,pathdir_write_to,filename,BF_try_time,BF_ops_file_type, PA, polyfit_deg, "asin")
|
PA = run_bf_polyfit(pathdir,pathdir_write_to,filename,BF_try_time,BF_ops_file_type, PA, polyfit_deg, "asin")
|
||||||
|
|
||||||
except:
|
except:
|
||||||
return PA
|
return PA
|
||||||
|
|
||||||
return PA
|
return PA
|
||||||
|
|
||||||
def get_atan(pathdir,pathdir_write_to,filename,BF_try_time,BF_ops_file_type, PA, polyfit_deg=4):
|
def get_atan(pathdir,pathdir_write_to,filename,BF_try_time,BF_ops_file_type, PA, polyfit_deg=3):
|
||||||
try:
|
try:
|
||||||
os.mkdir(pathdir_write_to)
|
os.mkdir(pathdir_write_to)
|
||||||
except:
|
except:
|
||||||
pass
|
pass
|
||||||
|
data = np.loadtxt(pathdir+filename)
|
||||||
try:
|
try:
|
||||||
n_variables = np.loadtxt(pathdir+"%s" %filename, dtype='str').shape[1]-1
|
data[:,-1] = np.arctan(data[:,-1])
|
||||||
variables = np.loadtxt(pathdir+"%s" %filename, usecols=(0,))
|
np.savetxt(pathdir_write_to+filename,data)
|
||||||
for j in range(1,n_variables):
|
|
||||||
v = np.loadtxt(pathdir+"%s" %filename, usecols=(j,))
|
|
||||||
variables = np.column_stack((variables,v))
|
|
||||||
f_dependent = np.loadtxt(pathdir+"%s" %filename, usecols=(n_variables,))
|
|
||||||
|
|
||||||
dt = np.column_stack((variables,np.arctan(f_dependent)))
|
|
||||||
np.savetxt(pathdir_write_to+filename,dt)
|
|
||||||
|
|
||||||
PA = run_bf_polyfit(pathdir,pathdir_write_to,filename,BF_try_time,BF_ops_file_type, PA, polyfit_deg, "atan")
|
PA = run_bf_polyfit(pathdir,pathdir_write_to,filename,BF_try_time,BF_ops_file_type, PA, polyfit_deg, "atan")
|
||||||
|
|
||||||
except:
|
except:
|
||||||
return PA
|
return PA
|
||||||
|
|
||||||
return PA
|
return PA
|
||||||
|
|
||||||
|
|
||||||
def get_cos(pathdir,pathdir_write_to,filename,BF_try_time,BF_ops_file_type, PA, polyfit_deg=4):
|
def get_cos(pathdir,pathdir_write_to,filename,BF_try_time,BF_ops_file_type, PA, polyfit_deg=3):
|
||||||
try:
|
try:
|
||||||
os.mkdir(pathdir_write_to)
|
os.mkdir(pathdir_write_to)
|
||||||
except:
|
except:
|
||||||
pass
|
pass
|
||||||
|
data = np.loadtxt(pathdir+filename)
|
||||||
try:
|
try:
|
||||||
n_variables = np.loadtxt(pathdir+"%s" %filename, dtype='str').shape[1]-1
|
data[:,-1] = np.cos(data[:,-1])
|
||||||
variables = np.loadtxt(pathdir+"%s" %filename, usecols=(0,))
|
np.savetxt(pathdir_write_to+filename,data)
|
||||||
for j in range(1,n_variables):
|
|
||||||
v = np.loadtxt(pathdir+"%s" %filename, usecols=(j,))
|
|
||||||
variables = np.column_stack((variables,v))
|
|
||||||
f_dependent = np.loadtxt(pathdir+"%s" %filename, usecols=(n_variables,))
|
|
||||||
|
|
||||||
dt = np.column_stack((variables,np.cos(f_dependent)))
|
|
||||||
np.savetxt(pathdir_write_to+filename,dt)
|
|
||||||
|
|
||||||
PA = run_bf_polyfit(pathdir,pathdir_write_to,filename,BF_try_time,BF_ops_file_type, PA, polyfit_deg, "cos")
|
PA = run_bf_polyfit(pathdir,pathdir_write_to,filename,BF_try_time,BF_ops_file_type, PA, polyfit_deg, "cos")
|
||||||
|
|
||||||
except:
|
except:
|
||||||
return PA
|
return PA
|
||||||
|
|
||||||
return PA
|
return PA
|
||||||
|
|
||||||
|
|
||||||
def get_exp(pathdir,pathdir_write_to,filename,BF_try_time,BF_ops_file_type, PA, polyfit_deg=4):
|
def get_exp(pathdir,pathdir_write_to,filename,BF_try_time,BF_ops_file_type, PA, polyfit_deg=3):
|
||||||
try:
|
try:
|
||||||
os.mkdir(pathdir_write_to)
|
os.mkdir(pathdir_write_to)
|
||||||
except:
|
except:
|
||||||
pass
|
pass
|
||||||
|
data = np.loadtxt(pathdir+filename)
|
||||||
try:
|
try:
|
||||||
n_variables = np.loadtxt(pathdir+"%s" %filename, dtype='str').shape[1]-1
|
data[:,-1] = np.exp(data[:,-1])
|
||||||
variables = np.loadtxt(pathdir+"%s" %filename, usecols=(0,))
|
np.savetxt(pathdir_write_to+filename,data)
|
||||||
for j in range(1,n_variables):
|
PA = run_bf_polyfit(pathdir,pathdir_write_to,filename,BF_try_time,BF_ops_file_type, PA, polyfit_deg, "exp")
|
||||||
v = np.loadtxt(pathdir+"%s" %filename, usecols=(j,))
|
|
||||||
variables = np.column_stack((variables,v))
|
|
||||||
f_dependent = np.loadtxt(pathdir+"%s" %filename, usecols=(n_variables,))
|
|
||||||
|
|
||||||
dt = np.column_stack((variables,np.exp(f_dependent)))
|
|
||||||
np.savetxt(pathdir_write_to+filename,dt)
|
|
||||||
|
|
||||||
PA = run_bf_polyfit(pathdir,pathdir_write_to,filename,BF_try_time,BF_ops_file_type, PA, polyfit_deg, "exp")
|
|
||||||
|
|
||||||
except:
|
except:
|
||||||
return PA
|
return PA
|
||||||
|
|
||||||
return PA
|
return PA
|
||||||
|
|
||||||
|
|
||||||
def get_inverse(pathdir,pathdir_write_to,filename,BF_try_time,BF_ops_file_type, PA, polyfit_deg=4):
|
def get_inverse(pathdir,pathdir_write_to,filename,BF_try_time,BF_ops_file_type, PA, polyfit_deg=3):
|
||||||
try:
|
try:
|
||||||
os.mkdir(pathdir_write_to)
|
os.mkdir(pathdir_write_to)
|
||||||
except:
|
except:
|
||||||
pass
|
pass
|
||||||
|
data = np.loadtxt(pathdir+filename)
|
||||||
try:
|
try:
|
||||||
n_variables = np.loadtxt(pathdir+"%s" %filename, dtype='str').shape[1]-1
|
data[:,-1] = 1/data[:,-1]
|
||||||
variables = np.loadtxt(pathdir+"%s" %filename, usecols=(0,))
|
np.savetxt(pathdir_write_to+filename,data)
|
||||||
for j in range(1,n_variables):
|
PA = run_bf_polyfit(pathdir,pathdir_write_to,filename,BF_try_time,BF_ops_file_type, PA, polyfit_deg, "inverse")
|
||||||
v = np.loadtxt(pathdir+"%s" %filename, usecols=(j,))
|
|
||||||
variables = np.column_stack((variables,v))
|
|
||||||
f_dependent = np.loadtxt(pathdir+"%s" %filename, usecols=(n_variables,))
|
|
||||||
|
|
||||||
dt = np.column_stack((variables,1/f_dependent))
|
|
||||||
np.savetxt(pathdir_write_to+filename,dt)
|
|
||||||
|
|
||||||
PA = run_bf_polyfit(pathdir,pathdir_write_to,filename,BF_try_time,BF_ops_file_type, PA, polyfit_deg, "inverse")
|
|
||||||
|
|
||||||
except:
|
except:
|
||||||
return PA
|
return PA
|
||||||
|
|
||||||
return PA
|
return PA
|
||||||
|
|
||||||
|
|
||||||
def get_log(pathdir,pathdir_write_to,filename,BF_try_time,BF_ops_file_type, PA, polyfit_deg=4):
|
def get_log(pathdir,pathdir_write_to,filename,BF_try_time,BF_ops_file_type, PA, polyfit_deg=3):
|
||||||
try:
|
try:
|
||||||
os.mkdir(pathdir_write_to)
|
os.mkdir(pathdir_write_to)
|
||||||
except:
|
except:
|
||||||
pass
|
pass
|
||||||
|
data = np.loadtxt(pathdir+filename)
|
||||||
try:
|
try:
|
||||||
n_variables = np.loadtxt(pathdir+"%s" %filename, dtype='str').shape[1]-1
|
data[:,-1] = np.log(data[:,-1])
|
||||||
variables = np.loadtxt(pathdir+"%s" %filename, usecols=(0,))
|
np.savetxt(pathdir_write_to+filename,data)
|
||||||
for j in range(1,n_variables):
|
PA = run_bf_polyfit(pathdir,pathdir_write_to,filename,BF_try_time,BF_ops_file_type, PA, polyfit_deg, "log")
|
||||||
v = np.loadtxt(pathdir+"%s" %filename, usecols=(j,))
|
|
||||||
variables = np.column_stack((variables,v))
|
|
||||||
f_dependent = np.loadtxt(pathdir+"%s" %filename, usecols=(n_variables,))
|
|
||||||
|
|
||||||
dt = np.column_stack((variables,np.log(f_dependent)))
|
|
||||||
np.savetxt(pathdir_write_to+filename,dt)
|
|
||||||
|
|
||||||
PA = run_bf_polyfit(pathdir,pathdir_write_to,filename,BF_try_time,BF_ops_file_type, PA, polyfit_deg, "log")
|
|
||||||
|
|
||||||
except:
|
except:
|
||||||
return PA
|
return PA
|
||||||
|
|
||||||
return PA
|
return PA
|
||||||
|
|
||||||
|
|
||||||
def get_sin(pathdir,pathdir_write_to,filename,BF_try_time,BF_ops_file_type, PA, polyfit_deg=4):
|
def get_sin(pathdir,pathdir_write_to,filename,BF_try_time,BF_ops_file_type, PA, polyfit_deg=3):
|
||||||
try:
|
try:
|
||||||
os.mkdir(pathdir_write_to)
|
os.mkdir(pathdir_write_to)
|
||||||
except:
|
except:
|
||||||
pass
|
pass
|
||||||
|
data = np.loadtxt(pathdir+filename)
|
||||||
try:
|
try:
|
||||||
n_variables = np.loadtxt(pathdir+"%s" %filename, dtype='str').shape[1]-1
|
data[:,-1] = np.sin(data[:,-1])
|
||||||
variables = np.loadtxt(pathdir+"%s" %filename, usecols=(0,))
|
np.savetxt(pathdir_write_to+filename,data)
|
||||||
for j in range(1,n_variables):
|
PA = run_bf_polyfit(pathdir,pathdir_write_to,filename,BF_try_time,BF_ops_file_type, PA, polyfit_deg, "sin")
|
||||||
v = np.loadtxt(pathdir+"%s" %filename, usecols=(j,))
|
|
||||||
variables = np.column_stack((variables,v))
|
|
||||||
f_dependent = np.loadtxt(pathdir+"%s" %filename, usecols=(n_variables,))
|
|
||||||
|
|
||||||
dt = np.column_stack((variables,np.sin(f_dependent)))
|
|
||||||
np.savetxt(pathdir_write_to+filename,dt)
|
|
||||||
|
|
||||||
PA = run_bf_polyfit(pathdir,pathdir_write_to,filename,BF_try_time,BF_ops_file_type, PA, polyfit_deg, "sin")
|
|
||||||
|
|
||||||
except:
|
except:
|
||||||
return PA
|
return PA
|
||||||
|
|
||||||
return PA
|
return PA
|
||||||
|
|
||||||
|
|
||||||
def get_sqrt(pathdir,pathdir_write_to,filename,BF_try_time,BF_ops_file_type, PA, polyfit_deg=4):
|
def get_sqrt(pathdir,pathdir_write_to,filename,BF_try_time,BF_ops_file_type, PA, polyfit_deg=3):
|
||||||
try:
|
try:
|
||||||
os.mkdir(pathdir_write_to)
|
os.mkdir(pathdir_write_to)
|
||||||
except:
|
except:
|
||||||
pass
|
pass
|
||||||
|
data = np.loadtxt(pathdir+filename)
|
||||||
try:
|
try:
|
||||||
n_variables = np.loadtxt(pathdir+"%s" %filename, dtype='str').shape[1]-1
|
data[:,-1] = np.sqrt(data[:,-1])
|
||||||
variables = np.loadtxt(pathdir+"%s" %filename, usecols=(0,))
|
np.savetxt(pathdir_write_to+filename,data)
|
||||||
for j in range(1,n_variables):
|
|
||||||
v = np.loadtxt(pathdir+"%s" %filename, usecols=(j,))
|
|
||||||
variables = np.column_stack((variables,v))
|
|
||||||
f_dependent = np.loadtxt(pathdir+"%s" %filename, usecols=(n_variables,))
|
|
||||||
|
|
||||||
dt = np.column_stack((variables,np.sqrt(f_dependent)))
|
|
||||||
np.savetxt(pathdir_write_to+filename,dt)
|
|
||||||
|
|
||||||
PA = run_bf_polyfit(pathdir,pathdir_write_to,filename,BF_try_time,BF_ops_file_type, PA, polyfit_deg, "sqrt")
|
PA = run_bf_polyfit(pathdir,pathdir_write_to,filename,BF_try_time,BF_ops_file_type, PA, polyfit_deg, "sqrt")
|
||||||
|
|
||||||
except:
|
except:
|
||||||
|
|
@ -216,22 +145,15 @@ def get_sqrt(pathdir,pathdir_write_to,filename,BF_try_time,BF_ops_file_type, PA,
|
||||||
return PA
|
return PA
|
||||||
|
|
||||||
|
|
||||||
def get_squared(pathdir,pathdir_write_to,filename,BF_try_time,BF_ops_file_type, PA, polyfit_deg=4):
|
def get_squared(pathdir,pathdir_write_to,filename,BF_try_time,BF_ops_file_type, PA, polyfit_deg=3):
|
||||||
try:
|
try:
|
||||||
os.mkdir(pathdir_write_to)
|
os.mkdir(pathdir_write_to)
|
||||||
except:
|
except:
|
||||||
pass
|
pass
|
||||||
|
data = np.loadtxt(pathdir+filename)
|
||||||
try:
|
try:
|
||||||
n_variables = np.loadtxt(pathdir+"%s" %filename, dtype='str').shape[1]-1
|
data[:,-1] = data[:,-1]**2
|
||||||
variables = np.loadtxt(pathdir+"%s" %filename, usecols=(0,))
|
np.savetxt(pathdir_write_to+filename,data)
|
||||||
for j in range(1,n_variables):
|
|
||||||
v = np.loadtxt(pathdir+"%s" %filename, usecols=(j,))
|
|
||||||
variables = np.column_stack((variables,v))
|
|
||||||
f_dependent = np.loadtxt(pathdir+"%s" %filename, usecols=(n_variables,))
|
|
||||||
|
|
||||||
dt = np.column_stack((variables,f_dependent**2))
|
|
||||||
np.savetxt(pathdir_write_to+filename,dt)
|
|
||||||
|
|
||||||
PA = run_bf_polyfit(pathdir,pathdir_write_to,filename,BF_try_time,BF_ops_file_type, PA, polyfit_deg, "squared")
|
PA = run_bf_polyfit(pathdir,pathdir_write_to,filename,BF_try_time,BF_ops_file_type, PA, polyfit_deg, "squared")
|
||||||
|
|
||||||
except:
|
except:
|
||||||
|
|
@ -240,22 +162,15 @@ def get_squared(pathdir,pathdir_write_to,filename,BF_try_time,BF_ops_file_type,
|
||||||
return PA
|
return PA
|
||||||
|
|
||||||
|
|
||||||
def get_tan(pathdir,pathdir_write_to,filename,BF_try_time,BF_ops_file_type, PA, polyfit_deg=4):
|
def get_tan(pathdir,pathdir_write_to,filename,BF_try_time,BF_ops_file_type, PA, polyfit_deg=3):
|
||||||
try:
|
try:
|
||||||
os.mkdir(pathdir_write_to)
|
os.mkdir(pathdir_write_to)
|
||||||
except:
|
except:
|
||||||
pass
|
pass
|
||||||
|
data = np.loadtxt(pathdir+filename)
|
||||||
try:
|
try:
|
||||||
n_variables = np.loadtxt(pathdir+"%s" %filename, dtype='str').shape[1]-1
|
data[:,-1] = np.tan(data[:,-1])
|
||||||
variables = np.loadtxt(pathdir+"%s" %filename, usecols=(0,))
|
np.savetxt(pathdir_write_to+filename,data)
|
||||||
for j in range(1,n_variables):
|
|
||||||
v = np.loadtxt(pathdir+"%s" %filename, usecols=(j,))
|
|
||||||
variables = np.column_stack((variables,v))
|
|
||||||
f_dependent = np.loadtxt(pathdir+"%s" %filename, usecols=(n_variables,))
|
|
||||||
|
|
||||||
dt = np.column_stack((variables,np.tan(f_dependent)))
|
|
||||||
np.savetxt(pathdir_write_to+filename,dt)
|
|
||||||
|
|
||||||
PA = run_bf_polyfit(pathdir,pathdir_write_to,filename,BF_try_time,BF_ops_file_type, PA, polyfit_deg, "tan")
|
PA = run_bf_polyfit(pathdir,pathdir_write_to,filename,BF_try_time,BF_ops_file_type, PA, polyfit_deg, "tan")
|
||||||
|
|
||||||
except:
|
except:
|
||||||
|
|
@ -265,9 +180,3 @@ def get_tan(pathdir,pathdir_write_to,filename,BF_try_time,BF_ops_file_type, PA,
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
|
|
||||||
|
|
@ -20,225 +20,227 @@ from os import path
|
||||||
def run_bf_polyfit(pathdir,pathdir_transformed,filename,BF_try_time,BF_ops_file_type, PA, polyfit_deg=3, output_type=""):
|
def run_bf_polyfit(pathdir,pathdir_transformed,filename,BF_try_time,BF_ops_file_type, PA, polyfit_deg=3, output_type=""):
|
||||||
input_data = np.loadtxt(pathdir_transformed+filename)
|
input_data = np.loadtxt(pathdir_transformed+filename)
|
||||||
#############################################################################################################################
|
#############################################################################################################################
|
||||||
|
if np.isnan(input_data).any()==False:
|
||||||
# run BF on the data (+)
|
# run BF on the data (+)
|
||||||
print("Checking for brute force + \n")
|
print("Checking for brute force + \n")
|
||||||
brute_force(pathdir_transformed,filename,BF_try_time,BF_ops_file_type,"+")
|
brute_force(pathdir_transformed,filename,BF_try_time,BF_ops_file_type,"+")
|
||||||
|
|
||||||
try:
|
|
||||||
# load the BF output data
|
|
||||||
bf_all_output = np.loadtxt("results.dat", dtype="str")
|
|
||||||
express = bf_all_output[:,2]
|
|
||||||
prefactors = bf_all_output[:,1]
|
|
||||||
prefactors = [str(i) for i in prefactors]
|
|
||||||
|
|
||||||
# Calculate the complexity of the bf expression the same way as for gradient descent case
|
|
||||||
complexity = []
|
|
||||||
errors = []
|
|
||||||
eqns = []
|
|
||||||
for i in range(len(prefactors)):
|
|
||||||
try:
|
|
||||||
if output_type=="":
|
|
||||||
eqn = prefactors[i] + "+" + RPN_to_eq(express[i])
|
|
||||||
elif output_type=="acos":
|
|
||||||
eqn = "cos(" + prefactors[i] + "+" + RPN_to_eq(express[i]) + ")"
|
|
||||||
elif output_type=="asin":
|
|
||||||
eqn = "sin(" + prefactors[i] + "+" + RPN_to_eq(express[i]) + ")"
|
|
||||||
elif output_type=="atan":
|
|
||||||
eqn = "tan(" + prefactors[i] + "+" + RPN_to_eq(express[i]) + ")"
|
|
||||||
elif output_type=="cos":
|
|
||||||
eqn = "acos(" + prefactors[i] + "+" + RPN_to_eq(express[i]) + ")"
|
|
||||||
elif output_type=="exp":
|
|
||||||
eqn = "log(" + prefactors[i] + "+" + RPN_to_eq(express[i]) + ")"
|
|
||||||
elif output_type=="inverse":
|
|
||||||
eqn = "1/(" + prefactors[i] + "+" + RPN_to_eq(express[i]) + ")"
|
|
||||||
elif output_type=="log":
|
|
||||||
eqn = "exp(" + prefactors[i] + "+" + RPN_to_eq(express[i]) + ")"
|
|
||||||
elif output_type=="sin":
|
|
||||||
eqn = "asin(" + prefactors[i] + "+" + RPN_to_eq(express[i]) + ")"
|
|
||||||
elif output_type=="sqrt":
|
|
||||||
eqn = "(" + prefactors[i] + "+" + RPN_to_eq(express[i]) + ")**2"
|
|
||||||
elif output_type=="squared":
|
|
||||||
eqn = "sqrt(" + prefactors[i] + "+" + RPN_to_eq(express[i]) + ")"
|
|
||||||
elif output_type=="tan":
|
|
||||||
eqn = "atan(" + prefactors[i] + "+" + RPN_to_eq(express[i]) + ")"
|
|
||||||
|
|
||||||
eqns = eqns + [eqn]
|
|
||||||
errors = errors + [get_symbolic_expr_error(input_data,eqn)]
|
|
||||||
expr = parse_expr(eqn)
|
|
||||||
is_atomic_number = lambda expr: expr.is_Atom and expr.is_number
|
|
||||||
numbers_expr = [subexpression for subexpression in preorder_traversal(expr) if is_atomic_number(subexpression)]
|
|
||||||
compl = 0
|
|
||||||
for j in numbers_expr:
|
|
||||||
try:
|
|
||||||
compl = compl + get_number_DL_snapped(float(j))
|
|
||||||
except:
|
|
||||||
compl = compl + 1000000
|
|
||||||
|
|
||||||
# Add the complexity due to symbols
|
|
||||||
n_variables = len(expr.free_symbols)
|
|
||||||
n_operations = len(count_ops(expr,visual=True).free_symbols)
|
|
||||||
if n_operations!=0 or n_variables!=0:
|
|
||||||
compl = compl + (n_variables+n_operations)*np.log2((n_variables+n_operations))
|
|
||||||
|
|
||||||
complexity = complexity + [compl]
|
|
||||||
except:
|
|
||||||
continue
|
|
||||||
|
|
||||||
for i in range(len(complexity)):
|
|
||||||
PA.add(Point(x=complexity[i], y=errors[i], data=eqns[i]))
|
|
||||||
|
|
||||||
# run gradient descent of BF output parameters and add the results to the Pareto plot
|
|
||||||
for i in range(len(express)):
|
|
||||||
try:
|
|
||||||
bf_gd_update = RPN_to_pytorch(input_data,eqns[i])
|
|
||||||
PA.add(Point(x=bf_gd_update[1],y=bf_gd_update[0],data=bf_gd_update[2]))
|
|
||||||
except:
|
|
||||||
continue
|
|
||||||
except:
|
|
||||||
pass
|
|
||||||
|
|
||||||
#############################################################################################################################
|
|
||||||
# run BF on the data (*)
|
|
||||||
print("Checking for brute force * \n")
|
|
||||||
brute_force(pathdir_transformed,filename,BF_try_time,BF_ops_file_type,"*")
|
|
||||||
|
|
||||||
try:
|
|
||||||
# load the BF output data
|
|
||||||
bf_all_output = np.loadtxt("results.dat", dtype="str")
|
|
||||||
express = bf_all_output[:,2]
|
|
||||||
prefactors = bf_all_output[:,1]
|
|
||||||
prefactors = [str(i) for i in prefactors]
|
|
||||||
|
|
||||||
# Calculate the complexity of the bf expression the same way as for gradient descent case
|
|
||||||
complexity = []
|
|
||||||
errors = []
|
|
||||||
eqns = []
|
|
||||||
for i in range(len(prefactors)):
|
|
||||||
try:
|
|
||||||
if output_type=="":
|
|
||||||
eqn = prefactors[i] + "*" + RPN_to_eq(express[i])
|
|
||||||
elif output_type=="acos":
|
|
||||||
eqn = "cos(" + prefactors[i] + "*" + RPN_to_eq(express[i]) + ")"
|
|
||||||
elif output_type=="asin":
|
|
||||||
eqn = "sin(" + prefactors[i] + "*" + RPN_to_eq(express[i]) + ")"
|
|
||||||
elif output_type=="atan":
|
|
||||||
eqn = "tan(" + prefactors[i] + "*" + RPN_to_eq(express[i]) + ")"
|
|
||||||
elif output_type=="cos":
|
|
||||||
eqn = "acos(" + prefactors[i] + "*" + RPN_to_eq(express[i]) + ")"
|
|
||||||
elif output_type=="exp":
|
|
||||||
eqn = "log(" + prefactors[i] + "*" + RPN_to_eq(express[i]) + ")"
|
|
||||||
elif output_type=="inverse":
|
|
||||||
eqn = "1/(" + prefactors[i] + "*" + RPN_to_eq(express[i]) + ")"
|
|
||||||
elif output_type=="log":
|
|
||||||
eqn = "exp(" + prefactors[i] + "*" + RPN_to_eq(express[i]) + ")"
|
|
||||||
elif output_type=="sin":
|
|
||||||
eqn = "asin(" + prefactors[i] + "*" + RPN_to_eq(express[i]) + ")"
|
|
||||||
elif output_type=="sqrt":
|
|
||||||
eqn = "(" + prefactors[i] + "*" + RPN_to_eq(express[i]) + ")**2"
|
|
||||||
elif output_type=="squared":
|
|
||||||
eqn = "sqrt(" + prefactors[i] + "*" + RPN_to_eq(express[i]) + ")"
|
|
||||||
elif output_type=="tan":
|
|
||||||
eqn = "atan(" + prefactors[i] + "*" + RPN_to_eq(express[i]) + ")"
|
|
||||||
|
|
||||||
eqns = eqns + [eqn]
|
|
||||||
errors = errors + [get_symbolic_expr_error(input_data,eqn)]
|
|
||||||
expr = parse_expr(eqn)
|
|
||||||
is_atomic_number = lambda expr: expr.is_Atom and expr.is_number
|
|
||||||
numbers_expr = [subexpression for subexpression in preorder_traversal(expr) if is_atomic_number(subexpression)]
|
|
||||||
compl = 0
|
|
||||||
for j in numbers_expr:
|
|
||||||
try:
|
|
||||||
compl = compl + get_number_DL_snapped(float(j))
|
|
||||||
except:
|
|
||||||
compl = compl + 1000000
|
|
||||||
|
|
||||||
# Add the complexity due to symbols
|
|
||||||
n_variables = len(expr.free_symbols)
|
|
||||||
n_operations = len(count_ops(expr,visual=True).free_symbols)
|
|
||||||
if n_operations!=0 or n_variables!=0:
|
|
||||||
compl = compl + (n_variables+n_operations)*np.log2((n_variables+n_operations))
|
|
||||||
|
|
||||||
complexity = complexity + [compl]
|
|
||||||
except:
|
|
||||||
continue
|
|
||||||
|
|
||||||
# add the BF output to the Pareto plot
|
|
||||||
for i in range(len(complexity)):
|
|
||||||
PA.add(Point(x=complexity[i], y=errors[i], data=eqns[i]))
|
|
||||||
|
|
||||||
# run gradient descent of BF output parameters and add the results to the Pareto plot
|
|
||||||
for i in range(len(express)):
|
|
||||||
try:
|
|
||||||
bf_gd_update = RPN_to_pytorch(input_data,eqns[i])
|
|
||||||
PA.add(Point(x=bf_gd_update[1],y=bf_gd_update[0],data=bf_gd_update[2]))
|
|
||||||
except:
|
|
||||||
continue
|
|
||||||
except:
|
|
||||||
pass
|
|
||||||
|
|
||||||
#############################################################################################################################
|
|
||||||
# run polyfit on the data
|
|
||||||
print("Checking polyfit \n")
|
|
||||||
try:
|
|
||||||
polyfit_result = polyfit(polyfit_deg, pathdir_transformed+filename)
|
|
||||||
eqn = str(polyfit_result[0])
|
|
||||||
|
|
||||||
# Calculate the complexity of the polyfit expression the same way as for gradient descent case
|
|
||||||
if output_type=="":
|
|
||||||
eqn = eqn
|
|
||||||
elif output_type=="acos":
|
|
||||||
eqn = "cos(" + eqn + ")"
|
|
||||||
elif output_type=="asin":
|
|
||||||
eqn = "sin(" + eqn + ")"
|
|
||||||
elif output_type=="atan":
|
|
||||||
eqn = "tan(" + eqn + ")"
|
|
||||||
elif output_type=="cos":
|
|
||||||
eqn = "acos(" + eqn + ")"
|
|
||||||
elif output_type=="exp":
|
|
||||||
eqn = "log(" + eqn + ")"
|
|
||||||
elif output_type=="inverse":
|
|
||||||
eqn = "1/(" + eqn + ")"
|
|
||||||
elif output_type=="log":
|
|
||||||
eqn = "exp(" + eqn + ")"
|
|
||||||
elif output_type=="sin":
|
|
||||||
eqn = "asin(" + eqn + ")"
|
|
||||||
elif output_type=="sqrt":
|
|
||||||
eqn = "(" + eqn + ")**2"
|
|
||||||
elif output_type=="squared":
|
|
||||||
eqn = "sqrt(" + eqn + ")"
|
|
||||||
elif output_type=="tan":
|
|
||||||
eqn = "atan(" + eqn + ")"
|
|
||||||
|
|
||||||
polyfit_err = get_symbolic_expr_error(input_data,eqn)
|
|
||||||
expr = parse_expr(eqn)
|
|
||||||
is_atomic_number = lambda expr: expr.is_Atom and expr.is_number
|
|
||||||
numbers_expr = [subexpression for subexpression in preorder_traversal(expr) if is_atomic_number(subexpression)]
|
|
||||||
complexity = 0
|
|
||||||
for j in numbers_expr:
|
|
||||||
complexity = complexity + get_number_DL_snapped(float(j))
|
|
||||||
try:
|
try:
|
||||||
# Add the complexity due to symbols
|
# load the BF output data
|
||||||
n_variables = len(polyfit_result[0].free_symbols)
|
bf_all_output = np.loadtxt("results.dat", dtype="str")
|
||||||
n_operations = len(count_ops(polyfit_result[0],visual=True).free_symbols)
|
express = bf_all_output[:,2]
|
||||||
if n_operations!=0 or n_variables!=0:
|
prefactors = bf_all_output[:,1]
|
||||||
complexity = complexity + (n_variables+n_operations)*np.log2((n_variables+n_operations))
|
prefactors = [str(i) for i in prefactors]
|
||||||
|
|
||||||
|
# Calculate the complexity of the bf expression the same way as for gradient descent case
|
||||||
|
complexity = []
|
||||||
|
errors = []
|
||||||
|
eqns = []
|
||||||
|
for i in range(len(prefactors)):
|
||||||
|
try:
|
||||||
|
if output_type=="":
|
||||||
|
eqn = prefactors[i] + "+" + RPN_to_eq(express[i])
|
||||||
|
elif output_type=="acos":
|
||||||
|
eqn = "cos(" + prefactors[i] + "+" + RPN_to_eq(express[i]) + ")"
|
||||||
|
elif output_type=="asin":
|
||||||
|
eqn = "sin(" + prefactors[i] + "+" + RPN_to_eq(express[i]) + ")"
|
||||||
|
elif output_type=="atan":
|
||||||
|
eqn = "tan(" + prefactors[i] + "+" + RPN_to_eq(express[i]) + ")"
|
||||||
|
elif output_type=="cos":
|
||||||
|
eqn = "acos(" + prefactors[i] + "+" + RPN_to_eq(express[i]) + ")"
|
||||||
|
elif output_type=="exp":
|
||||||
|
eqn = "log(" + prefactors[i] + "+" + RPN_to_eq(express[i]) + ")"
|
||||||
|
elif output_type=="inverse":
|
||||||
|
eqn = "1/(" + prefactors[i] + "+" + RPN_to_eq(express[i]) + ")"
|
||||||
|
elif output_type=="log":
|
||||||
|
eqn = "exp(" + prefactors[i] + "+" + RPN_to_eq(express[i]) + ")"
|
||||||
|
elif output_type=="sin":
|
||||||
|
eqn = "asin(" + prefactors[i] + "+" + RPN_to_eq(express[i]) + ")"
|
||||||
|
elif output_type=="sqrt":
|
||||||
|
eqn = "(" + prefactors[i] + "+" + RPN_to_eq(express[i]) + ")**2"
|
||||||
|
elif output_type=="squared":
|
||||||
|
eqn = "sqrt(" + prefactors[i] + "+" + RPN_to_eq(express[i]) + ")"
|
||||||
|
elif output_type=="tan":
|
||||||
|
eqn = "atan(" + prefactors[i] + "+" + RPN_to_eq(express[i]) + ")"
|
||||||
|
|
||||||
|
eqns = eqns + [eqn]
|
||||||
|
errors = errors + [get_symbolic_expr_error(input_data,eqn)]
|
||||||
|
expr = parse_expr(eqn)
|
||||||
|
is_atomic_number = lambda expr: expr.is_Atom and expr.is_number
|
||||||
|
numbers_expr = [subexpression for subexpression in preorder_traversal(expr) if is_atomic_number(subexpression)]
|
||||||
|
compl = 0
|
||||||
|
for j in numbers_expr:
|
||||||
|
try:
|
||||||
|
compl = compl + get_number_DL_snapped(float(j))
|
||||||
|
except:
|
||||||
|
compl = compl + 1000000
|
||||||
|
|
||||||
|
# Add the complexity due to symbols
|
||||||
|
n_variables = len(expr.free_symbols)
|
||||||
|
n_operations = len(count_ops(expr,visual=True).free_symbols)
|
||||||
|
if n_operations!=0 or n_variables!=0:
|
||||||
|
compl = compl + (n_variables+n_operations)*np.log2((n_variables+n_operations))
|
||||||
|
|
||||||
|
complexity = complexity + [compl]
|
||||||
|
except:
|
||||||
|
continue
|
||||||
|
|
||||||
|
for i in range(len(complexity)):
|
||||||
|
PA.add(Point(x=complexity[i], y=errors[i], data=eqns[i]))
|
||||||
|
|
||||||
|
# run gradient descent of BF output parameters and add the results to the Pareto plot
|
||||||
|
for i in range(len(express)):
|
||||||
|
try:
|
||||||
|
bf_gd_update = RPN_to_pytorch(input_data,eqns[i])
|
||||||
|
PA.add(Point(x=bf_gd_update[1],y=bf_gd_update[0],data=bf_gd_update[2]))
|
||||||
|
except:
|
||||||
|
continue
|
||||||
except:
|
except:
|
||||||
pass
|
pass
|
||||||
|
|
||||||
#run zero snap on polyfit output
|
#############################################################################################################################
|
||||||
PA_poly = ParetoSet()
|
# run BF on the data (*)
|
||||||
PA_poly.add(Point(x=complexity, y=polyfit_err, data=str(eqn)))
|
print("Checking for brute force * \n")
|
||||||
PA_poly = add_snap_expr_on_pareto_polyfit(pathdir, filename, str(eqn), PA_poly)
|
brute_force(pathdir_transformed,filename,BF_try_time,BF_ops_file_type,"*")
|
||||||
|
|
||||||
for l in range(len(PA_poly.get_pareto_points())):
|
try:
|
||||||
PA.add(Point(PA_poly.get_pareto_points()[l][0],PA_poly.get_pareto_points()[l][1],PA_poly.get_pareto_points()[l][2]))
|
# load the BF output data
|
||||||
|
bf_all_output = np.loadtxt("results.dat", dtype="str")
|
||||||
except:
|
express = bf_all_output[:,2]
|
||||||
pass
|
prefactors = bf_all_output[:,1]
|
||||||
|
prefactors = [str(i) for i in prefactors]
|
||||||
print("Complexity RMSE Expression")
|
|
||||||
for pareto_i in range(len(PA.get_pareto_points())):
|
# Calculate the complexity of the bf expression the same way as for gradient descent case
|
||||||
print(PA.get_pareto_points()[pareto_i])
|
complexity = []
|
||||||
|
errors = []
|
||||||
return PA
|
eqns = []
|
||||||
|
for i in range(len(prefactors)):
|
||||||
|
try:
|
||||||
|
if output_type=="":
|
||||||
|
eqn = prefactors[i] + "*" + RPN_to_eq(express[i])
|
||||||
|
elif output_type=="acos":
|
||||||
|
eqn = "cos(" + prefactors[i] + "*" + RPN_to_eq(express[i]) + ")"
|
||||||
|
elif output_type=="asin":
|
||||||
|
eqn = "sin(" + prefactors[i] + "*" + RPN_to_eq(express[i]) + ")"
|
||||||
|
elif output_type=="atan":
|
||||||
|
eqn = "tan(" + prefactors[i] + "*" + RPN_to_eq(express[i]) + ")"
|
||||||
|
elif output_type=="cos":
|
||||||
|
eqn = "acos(" + prefactors[i] + "*" + RPN_to_eq(express[i]) + ")"
|
||||||
|
elif output_type=="exp":
|
||||||
|
eqn = "log(" + prefactors[i] + "*" + RPN_to_eq(express[i]) + ")"
|
||||||
|
elif output_type=="inverse":
|
||||||
|
eqn = "1/(" + prefactors[i] + "*" + RPN_to_eq(express[i]) + ")"
|
||||||
|
elif output_type=="log":
|
||||||
|
eqn = "exp(" + prefactors[i] + "*" + RPN_to_eq(express[i]) + ")"
|
||||||
|
elif output_type=="sin":
|
||||||
|
eqn = "asin(" + prefactors[i] + "*" + RPN_to_eq(express[i]) + ")"
|
||||||
|
elif output_type=="sqrt":
|
||||||
|
eqn = "(" + prefactors[i] + "*" + RPN_to_eq(express[i]) + ")**2"
|
||||||
|
elif output_type=="squared":
|
||||||
|
eqn = "sqrt(" + prefactors[i] + "*" + RPN_to_eq(express[i]) + ")"
|
||||||
|
elif output_type=="tan":
|
||||||
|
eqn = "atan(" + prefactors[i] + "*" + RPN_to_eq(express[i]) + ")"
|
||||||
|
|
||||||
|
eqns = eqns + [eqn]
|
||||||
|
errors = errors + [get_symbolic_expr_error(input_data,eqn)]
|
||||||
|
expr = parse_expr(eqn)
|
||||||
|
is_atomic_number = lambda expr: expr.is_Atom and expr.is_number
|
||||||
|
numbers_expr = [subexpression for subexpression in preorder_traversal(expr) if is_atomic_number(subexpression)]
|
||||||
|
compl = 0
|
||||||
|
for j in numbers_expr:
|
||||||
|
try:
|
||||||
|
compl = compl + get_number_DL_snapped(float(j))
|
||||||
|
except:
|
||||||
|
compl = compl + 1000000
|
||||||
|
|
||||||
|
# Add the complexity due to symbols
|
||||||
|
n_variables = len(expr.free_symbols)
|
||||||
|
n_operations = len(count_ops(expr,visual=True).free_symbols)
|
||||||
|
if n_operations!=0 or n_variables!=0:
|
||||||
|
compl = compl + (n_variables+n_operations)*np.log2((n_variables+n_operations))
|
||||||
|
|
||||||
|
complexity = complexity + [compl]
|
||||||
|
except:
|
||||||
|
continue
|
||||||
|
|
||||||
|
# add the BF output to the Pareto plot
|
||||||
|
for i in range(len(complexity)):
|
||||||
|
PA.add(Point(x=complexity[i], y=errors[i], data=eqns[i]))
|
||||||
|
|
||||||
|
# run gradient descent of BF output parameters and add the results to the Pareto plot
|
||||||
|
for i in range(len(express)):
|
||||||
|
try:
|
||||||
|
bf_gd_update = RPN_to_pytorch(input_data,eqns[i])
|
||||||
|
PA.add(Point(x=bf_gd_update[1],y=bf_gd_update[0],data=bf_gd_update[2]))
|
||||||
|
except:
|
||||||
|
continue
|
||||||
|
except:
|
||||||
|
pass
|
||||||
|
|
||||||
|
#############################################################################################################################
|
||||||
|
# run polyfit on the data
|
||||||
|
print("Checking polyfit \n")
|
||||||
|
try:
|
||||||
|
polyfit_result = polyfit(polyfit_deg, pathdir_transformed+filename)
|
||||||
|
eqn = str(polyfit_result[0])
|
||||||
|
|
||||||
|
# Calculate the complexity of the polyfit expression the same way as for gradient descent case
|
||||||
|
if output_type=="":
|
||||||
|
eqn = eqn
|
||||||
|
elif output_type=="acos":
|
||||||
|
eqn = "cos(" + eqn + ")"
|
||||||
|
elif output_type=="asin":
|
||||||
|
eqn = "sin(" + eqn + ")"
|
||||||
|
elif output_type=="atan":
|
||||||
|
eqn = "tan(" + eqn + ")"
|
||||||
|
elif output_type=="cos":
|
||||||
|
eqn = "acos(" + eqn + ")"
|
||||||
|
elif output_type=="exp":
|
||||||
|
eqn = "log(" + eqn + ")"
|
||||||
|
elif output_type=="inverse":
|
||||||
|
eqn = "1/(" + eqn + ")"
|
||||||
|
elif output_type=="log":
|
||||||
|
eqn = "exp(" + eqn + ")"
|
||||||
|
elif output_type=="sin":
|
||||||
|
eqn = "asin(" + eqn + ")"
|
||||||
|
elif output_type=="sqrt":
|
||||||
|
eqn = "(" + eqn + ")**2"
|
||||||
|
elif output_type=="squared":
|
||||||
|
eqn = "sqrt(" + eqn + ")"
|
||||||
|
elif output_type=="tan":
|
||||||
|
eqn = "atan(" + eqn + ")"
|
||||||
|
|
||||||
|
polyfit_err = get_symbolic_expr_error(input_data,eqn)
|
||||||
|
expr = parse_expr(eqn)
|
||||||
|
is_atomic_number = lambda expr: expr.is_Atom and expr.is_number
|
||||||
|
numbers_expr = [subexpression for subexpression in preorder_traversal(expr) if is_atomic_number(subexpression)]
|
||||||
|
complexity = 0
|
||||||
|
for j in numbers_expr:
|
||||||
|
complexity = complexity + get_number_DL_snapped(float(j))
|
||||||
|
try:
|
||||||
|
# Add the complexity due to symbols
|
||||||
|
n_variables = len(polyfit_result[0].free_symbols)
|
||||||
|
n_operations = len(count_ops(polyfit_result[0],visual=True).free_symbols)
|
||||||
|
if n_operations!=0 or n_variables!=0:
|
||||||
|
complexity = complexity + (n_variables+n_operations)*np.log2((n_variables+n_operations))
|
||||||
|
except:
|
||||||
|
pass
|
||||||
|
|
||||||
|
#run zero snap on polyfit output
|
||||||
|
PA_poly = ParetoSet()
|
||||||
|
PA_poly.add(Point(x=complexity, y=polyfit_err, data=str(eqn)))
|
||||||
|
PA_poly = add_snap_expr_on_pareto_polyfit(pathdir, filename, str(eqn), PA_poly)
|
||||||
|
|
||||||
|
for l in range(len(PA_poly.get_pareto_points())):
|
||||||
|
PA.add(Point(PA_poly.get_pareto_points()[l][0],PA_poly.get_pareto_points()[l][1],PA_poly.get_pareto_points()[l][2]))
|
||||||
|
|
||||||
|
except:
|
||||||
|
pass
|
||||||
|
|
||||||
|
print("Complexity RMSE Expression")
|
||||||
|
for pareto_i in range(len(PA.get_pareto_points())):
|
||||||
|
print(PA.get_pareto_points()[pareto_i])
|
||||||
|
|
||||||
|
return PA
|
||||||
|
else:
|
||||||
|
return PA
|
||||||
|
|
|
||||||
Loading…
Add table
Add a link
Reference in a new issue