#-------------------------------------------------------------------------
#
# Local Imports
#
#-------------------------------------------------------------------------
import const
import lib.utils.utils as utils
DIR = const.DATA_DIR + os.sep + 'sandbox_domainsplit'
#create outputdir if not existing, delete if present
if not os.path.isdir(DIR):
os.mkdir(DIR)
else:
shutil.rmtree(DIR)
os.mkdir(DIR)
utils.set_outputdir(DIR)
#-------------------------------------------------------------------------
#
#Fipy Imports
#
#-------------------------------------------------------------------------
from fipy import *
#-------------------------------------------------------------------------
#
# DomainslitTestModel class
#
#-------------------------------------------------------------------------
class DomainsplitTestModel(object):
"""
pylab.ylim(0., 0.8)
pylab.axis()
pylab.show()
raw_input("check the figure")
#set up a yarn computation
from fiber.config import FiberConfigManager
from yarn.config import YarnConfigManager
from lib.utils.utils import set_outputdir
cfgf1 = FiberConfigManager.get_instance('tmpfiber1.ini', realdatastr=ini_fiber1)
cfgf2 = FiberConfigManager.get_instance('tmpfiber2.ini', realdatastr=ini_fiber2)
#cfgf3 = FiberConfigManager.get_instance('tmpfiber3.ini', realdatastr=ini_fiber3)
cfg = YarnConfigManager.get_instance('tmpyarn.ini', realdatastr=ini_yarn)
#create outputdir if not existing
if not os.path.isdir('temp'):
os.mkdir('temp')
set_outputdir('temp')
#create 10 2D grids for statistics
from yarn2d.yarn2dgrid import Yarn2dGrid
grid = Yarn2dGrid(cfg)
ouroptions = {
'x_central' : grid.x_central,
'y_central' : grid.y_central,
'number_fiber' : grid.number_fiber,
'number_fiber_blend' : grid.number_fiber_blend,
'radius_fiber' : grid.radius_fiber,
'radius_yarn' : grid.radius_yarn,
'theta_value' : grid.theta_value,
'beta_value' : grid.beta_value,
'mean_deviation': grid.mean_deviation,
'prob_area': grid.prob_area,
'radius_first_center': cfg.get(
