import BPV
import pattern_manipulation as pm
import numpy as np
import matplotlib.pyplot as plt
import pandas as pd
import ipdb
n = 10
N = 5
W = 0.025
data = BPV.Data()
data.read_csv("pixel.dist.csv", False)
data.df.sort_index(by="p", inplace=True, ascending=False)
data.df.set_index(pd.Index([j for j in range(len(data.df))]), inplace=True)
dftmp = data.data_head(n).df
p = np.zeros(2 * n)
for i in range(2 * n):
p[i] = dftmp['p'][i % n]
p /= p.sum()
df = BPV.Data(pd.DataFrame(p, columns=['p']))
df.df['plog1onp'] = df.df['p'] * np.log(1 / df.df['p'])
df.df.sort_index(by='p', inplace=True, ascending=False)
df.df.set_index(pd.Index([j for j in range(len(df.df))]), inplace=True)
prbl = BPV.BPV("decgraphV", df, N, W, time_solver=False)
prbl.solve()
#cProfile.run('prbl.solve()',sort=1)
prbl.pprint_solution()
import BPV
import pattern_manipulation as pm
import numpy as np
import matplotlib.pyplot as plt
import pandas as pd
import ipdb
n = 10
N = 5
W = 0.025
epsilon = 1000
df = BPV.Data()
df.read_csv("pixel.dist.csv", False)
df.df.sort_index(by="p", inplace=True, ascending=False)
df.df.set_index(pd.Index([j for j in range(len(df.df))]), inplace=True)
df = df.data_head(n)
#p = np.zeros(2*n)
#for i in range(2*n):
# p[i] = df.df['p'][i%n]
#p /= p.sum()
#df = BPV.Data(pd.DataFrame(p, columns=['p']))
#df.df['plog1onp'] = df.df['p']*np.log(1/df.df['p'])
#df.df.sort_index(by='p',inplace=True,ascending=False)
#df.df.set_index(pd.Index([j for j in range(len(df.df))]), inplace=True)
min_entropy = df.df['plog1onp'].min()
c = 2 * N / (epsilon * min_entropy)
print("\n c = ", c, "\n\n")
scaler = lambda x: (1 / c) * (int(c * x) + 1)
import BPV
import pattern_manipulation as pm
import numpy as np
import matplotlib.pyplot as plt
import pandas as pd
import cProfile
n = 100
N=20
W=0.03 #rate
data = BPV.Data()
data.read_csv("pixel.dist.csv",False)
data.df.sort_index(by="p",inplace=True,ascending=True)
data.df.set_index(pd.Index([j for j in range(len(data.df))]), inplace=True)
data_head = data.data_head(n)
#data.df.sort_index(by="p")
#prbl_pulp = BPV.BPV("pulp",data_head,N,W,time_solver=False)
#prbl_pulp.solve()
#prbl_pulp.pprint_solution()
#
prbl_decH = BPV.BPV("decgraphH",data_head,N,W,time_solver=True)
prbl_decH.solve()
#cProfile.run('prbl_decH.solve()',sort=1)
prbl_decH.pprint_solution()
#prbl_decW = BPV.BPV("decgraphW",data_head,N,W,time_solver=False)
#prbl_decW.solve()
##cProfile.run('prbl_decW.solve()',sort=1)
