from coef_setting import declare_coef #1
import main_code as main #2
#import main_code_con as main2
import picture_2d as pic #3
#from picture_2d_con import pic_2d
#from picture_3d import pic_3d
import numpy
from timeit import default_timer as timer
import time
import matplotlib.pyplot as plt
from mpmath.functions.rszeta import coef
#declare coefficients & initial settings
coef, set, sol = declare_coef() #1
#to plot interactively
plt.ion()
while set['t'] <= set['T'] and set['k'] < set['Nt']:
if set['con'] == True:
sol = main.continuous_iter(coef, set, sol)
'''PRINTING RESULT AS CONTROL'''
print 'at Time', set['t']
if not coef['Si'] == 0:
print 'NILAI C, b, n MAX', sol['c'].max(),',', sol['b'].max(),',', sol['n'].max()
print 'NILAI C, b, n MIN', sol['c'].min(),',', sol['b'].min(),',', sol['n'].min()
else:
sol = main.boolean_1_iter(coef, set, sol) #2
if sol['stop_iter'] >=100000:
from timeit import default_timer as timer
import time
import matplotlib.pyplot as plt
from matplotlib import cm
from matplotlib.ticker import LinearLocator, FormatStrFormatter
from mpl_toolkits.mplot3d import Axes3D
#from mpmath.functions.rszeta import coef
#import discrete_run as disc
from matplotlib import cm
from matplotlib.ticker import LinearLocator, FormatStrFormatter
from mpl_toolkits.mplot3d import Axes3D
#declare coefficients & initial settings
coef, set, sol = declare_coef()
#to plot interactively
plt.ion()
#hybrid part
while set['t'] <= set['T'] and set['k'] < set['Nt']:
start1 = timer()
sol = disc.boolean_1_iter(coef, set, sol)
start2 = timer()
if sol['stop_iter'] >= 100000:
set['k'] = sol['stop_iter']
'''PRINTING RESULT AS CONTROL'''
print 'at Time', set['t']
print 'Total Tips:', len(sol['matrix_tip'])
'''Import all files. We run all the files through run_main_code.py'''
from coef_setting import declare_coef #Ref.1
from initial_conditions import initial_prof #Ref.2
import main_code as main #Ref.4
import picture_2d as pic #Ref.5
'''Import additional functions to support calculation and graphics'''
import numpy
from timeit import default_timer as timer
import time
import matplotlib.pyplot as plt
plt.ion() #to plot interactively
'''declare coefficients & initial settings'''
coef, set, sol = declare_coef() #Ref.1
'''Initialize functions for k=0 (first iteration), t=0 (initial time)'''
sol = initial_prof(coef, set, sol) #Ref.2
pic.pic_2d(coef, set,
sol) #to generate the picture of initial condition #Ref.5
print '''All Coefficients:'''
print coef
print '''All Set'''
print set
'''PRINTING RESULT AS CONTROL'''
print 'at Time', set['t']
print 'Total Tips:', len(sol['matrix_tip'])
print 'Total Stop Tips:', len(sol['sp_stop'])
print 'Tip Stop:', sol['sp_stop']
print 'Tip Cell Pos:', sol['tip_cell']
print 'Max Value of c, f', sol['c'].max(), ',', sol['f'].max()
print 'Min Value of c, f', sol['c'].min(), ',', sol['f'].min()
'''Main part of the Hybrid calculation'''