def test_1(): bg = 0 # background val = 2 # values space = 10 thick = 1 NWorldFine = np.array([128, 128]) CoefClass = buildcoef2d.Coefficient2d(NWorldFine, bg=bg, val=val, length=1, thick=thick, space=space, probfactor=1, right=1, down=0, diagr1=0, diagr2=0, diagl1=0, diagl2=0, LenSwitch=None, thickSwitch=None, equidistant=True, ChannelHorizontal=None, ChannelVertical=True, BoundarySpace=True) A = CoefClass.BuildCoefficient().flatten() plt.figure("Coefficient_test_1") drawCoefficient(NWorldFine, A)
def test_2(): bg = 0 # background val = 2 # values space = 10 thick = 1 fine = 64 NWorldFine = np.array([fine, fine]) ChoosingShapes = np.array([ # shape , len, thick, space [1, 1, 1, 1], [1, 1, 1, 1], [1, 1, 1, 1], [1, 1, 1, 1], [1, 1, 1, 1], [1, 1, 1, 1], [1, 1, 1, 1], [1, 1, 1, 1], [1, 1, 1, 1], [1, 1, 1, 1], [1, 1, 1, 1], [1, 1, 1, 1], [4, fine - 2, 1, 5], [4, fine - 2, 1, 5], [4, fine - 2, 1, 5] ]) CoefClass = buildcoef2d.Coefficient2d(NWorldFine, bg=bg, val=val, right=1, LenSwitch=None, thickSwitch=None, equidistant=True, ChannelHorizontal=None, ChannelVertical=None, BoundarySpace=True, probfactor=1, ChoosingShapes=ChoosingShapes) A = CoefClass.BuildCoefficient().flatten() Correct = CoefClass.SpecificVanish(Number=[10, 11]) plt.figure("Coefficient_test_1") drawCoefficient(NWorldFine, Correct.flatten()) Number = [10, 11] random.seed(32) lis = np.zeros(80) lis[0] = 1 for i in range(np.shape(CoefClass.ShapeRemember)[0]): Number.append(i * random.sample(list(lis), 1)[0]) plt.figure("Perturbation") Perturbed = CoefClass.SpecificVanish(Number=Number, Original=True).flatten() drawCoefficient(NWorldFine, Perturbed)
world = World(NWorldCoarse, NCoarseElement, boundaryConditions) #righthandside f = np.ones(NpCoarse) #coefficient CoefClass = buildcoef2d.Coefficient2d(NWorldFine, bg = 0.05, val = 1, length = 8, thick = 8, space = 8, probfactor = 1, right = 1, down = 0, diagr1 = 0, diagr2 = 0, diagl1 = 0, diagl2 = 0, LenSwitch = None, thickSwitch = None, equidistant = True, ChannelHorizontal = None, ChannelVertical = None, BoundarySpace = True) A = CoefClass.BuildCoefficient() ABase = A.flatten() plt.figure("OriginalCoefficient") drawCoefficient(NWorldFine, ABase)
val = 1 #values NWorldFine = np.array([11, 11]) NWorldFine2 = np.array([38, 42]) ######################## # Print all shapes ####################### CoefClass = buildcoef2d.Coefficient2d(NWorldFine, bg=0.05, val=1, length=7, thick=1, space=2, probfactor=1, right=1, down=0, diagr1=0, diagr2=0, diagl1=0, diagl2=0, BoundarySpace=True) CoefClass1 = buildcoef2d.Coefficient2d(NWorldFine, bg=0.05, val=1, length=7, thick=1, space=2, probfactor=1, right=0,
world = World(NWorldCoarse, NCoarseElement, boundaryConditions) #righthandside f = np.ones(NpCoarse) #Coefficient 4 CoefClass = buildcoef2d.Coefficient2d(NWorldFine, bg=bg, val=val, thick=1, space=0, probfactor=1, right=0, down=0, diagr1=1, diagr2=0, diagl1=1, diagl2=0, LenSwitch=[4, 5, 6, 7, 8], thickSwitch=None, equidistant=None, ChannelHorizontal=None, ChannelVertical=None, BoundarySpace=None) A = CoefClass.BuildCoefficient() ABase = A.flatten() ROOT = '../../../test_data/Coef4/' #safe NworldFine
bg = 0.05 #background val = 1 #values NWorldFine = np.array([256, 256]) #Coefficient 1 CoefClass1 = buildcoef2d.Coefficient2d(NWorldFine, bg = bg, val = val, length = 2, thick = 2, space = 2, probfactor = 1, right = 1, down = 0, diagr1 = 0, diagr2 = 0, diagl1 = 0, diagl2 = 0, LenSwitch = None, thickSwitch = None, equidistant = True, ChannelHorizontal = None, ChannelVertical = None, BoundarySpace = True) #Coefficient 2 CoefClass2 = buildcoef2d.Coefficient2d(NWorldFine, bg = bg, val = val, length = 1,
bg = 0.05 #background val = 1 #values NWorldFine = np.array([42, 42]) CoefClass = buildcoef2d.Coefficient2d( NWorldFine, bg=bg, # background val=val, # values length=2, # length thick=2, # thickness space=2, # space between values probfactor=1, # probability of an value right=1, # shape 1 down=0, # shape 2 diagr1=0, # shape 3 diagr2=0, # shape 4 diagl1=0, # shape 5 diagl2=0, # shape 6 LenSwitch=None, # various length thickSwitch=None, # various thickness ChannelHorizontal=None, # horizontal Channels ChannelVertical=None, # vertical Channels BoundarySpace=True # additional space on the boundary ) A = CoefClass.BuildCoefficient() # coefficient in a numpy array numbers = [13, 20, 27, 44, 73] #What entries will be changed # Change in Value
xtFine_pert = psi.evaluate(xtFine) xtFine_ref = psi.inverse_evaluate(xtFine) bg = 0.01 #background val = 1 #values CoefClass = buildcoef2d.Coefficient2d(NCoeff, bg=bg, val=val, length=4, thick=4, space=4, probfactor=1, right=1, down=0, diagr1=0, diagr2=0, diagl1=0, diagl2=0, LenSwitch=None, thickSwitch=None, equidistant=True, ChannelHorizontal=None, ChannelVertical=None, BoundarySpace=True) # Set reference coefficient aCoarse_ref_shaped = CoefClass.BuildCoefficient() aCoarse_ref = aCoarse_ref_shaped.flatten() aFine_ref = func.evaluateDQ0(NCoeff, aCoarse_ref, xtFine_ref)
#righthandside f = np.ones(NpCoarse) #coefficient WormsUltCoefClass = buildcoef2d.Coefficient2d(NWorldFine, bg=bg, val=val, length=1, thick=1, space=1, probfactor=1, right=1, down=1, diagr1=1, diagr2=1, diagl1=1, diagl2=1, LenSwitch=[1, 2, 3, 4, 5, 6], thickSwitch=[1, 2, 3], equidistant=None, ChannelHorizontal=None, ChannelVertical=None, BoundarySpace=True, Boxes2n=None, Channels2n=None) A = WormsUltCoefClass.BuildCoefficient() ABase = A.flatten() plt.figure("OriginalCoefficient")
NList = [fine, 2, 4, 8, 16, 32, 64, 128, 256, 512] boundaryConditions = np.array([[0, 0], [0, 0]]) #coefficient WormsUltCoefClass = buildcoef2d.Coefficient2d(NWorldFine, bg=bg, val=val, length=1, thick=65, space=65, probfactor=1, right=1, down=0, diagr1=0, diagr2=0, diagl1=0, diagl2=0, LenSwitch=None, thickSwitch=None, equidistant=True, ChannelHorizontal=None, ChannelVertical=True, BoundarySpace=None, Boxes2n=None, Channels2n=None, TestExample=True) A = WormsUltCoefClass.BuildCoefficient() ABase = A.flatten() plt.figure("OriginalCoefficient") drawCoefficient(NWorldFine, ABase)