myprojects-hg-reborn

Fork

#! /usr/bin/env python


import numpy as np
import scipy as s


def break_segment(n):
        x = np.arange(1, n + 1, 1)
        d = 1. / n * (1. / (n - x + 1)).cumsum()
        return d


# def break_segment(n, r_to_n):
#         x = np.arange(1, n + 1, 1)
#         d = 1. / n * (1. / (n - x + 1)).cumsum()
#         return d



number_frag=100


nvec=s.zeros(0)
frag_vec=s.zeros(0)
kvec=s.zeros(0)

# fragments=break_segment(n)

# print(fragments)

for j in xrange(number_frag):
        n=j+1
        temp= break_segment(n)
        temp2=s.tile(n, len(temp))
        temp3=s.arange(1, len(temp)+1,1)
        temp3=s.flipud(temp3)
        
        frag_vec=s.hstack((frag_vec, temp))
        nvec=s.hstack((nvec, temp2))
        kvec=s.hstack((kvec, temp3))


print(frag_vec)

s.savetxt("fragmentation.dat", frag_vec)
s.savetxt("fragmentation-number-cuts.dat", nvec,fmt='%d')
s.savetxt("fragmentation-number-of-the-piece.dat", kvec,fmt='%d')

res=s.transpose(s.vstack((frag_vec, nvec, kvec)))
s.savetxt("results-frag.dat", res,"%f, %d, %d ")



print ("So far so good")