def test_bot_vs_time_top_vs_time_drn_1_double_loads(self): assert_allclose(dim1sin_E_Igamv_the_BC_deltaf_linear( **{'drn': 1, 'm': self.m, 'eigs':self.eigs, 'tvals':self.tvals,'Igamv':self.Igamv, 'zvals': self.zvals, 'pseudo_k':self.pseudo_k, 'top_vs_time': [self.top_vs_time,self.top_vs_time], 'bot_vs_time':[self.bot_vs_time,self.bot_vs_time]}), np.array([[ 175.40726974, 801.94431061], [ 343.82160516, 1571.91763224], [ 498.52885818, 2279.22355834]]))
def test_bot_vs_time_top_vs_time_drn_1(self): assert_allclose(dim1sin_E_Igamv_the_BC_deltaf_linear( **{'drn': 1, 'm': self.m, 'eigs':self.eigs, 'tvals':self.tvals,'Igamv':self.Igamv, 'zvals': self.zvals, 'pseudo_k':self.pseudo_k, 'top_vs_time': [self.top_vs_time], 'bot_vs_time':[self.bot_vs_time]}), np.array([[ 87.70363487, 400.97215531], [ 171.91080258, 785.95881612], [ 249.26442909, 1139.61177917]]))
def test_bot_vs_time_drn_0(self): assert_allclose(dim1sin_E_Igamv_the_BC_deltaf_linear( **{'drn': 0, 'm': self.m, 'eigs':self.eigs, 'tvals':self.tvals,'Igamv':self.Igamv, 'zvals': self.zvals, 'pseudo_k':self.pseudo_k, 'top_vs_time': None, 'bot_vs_time':[self.bot_vs_time]}), np.array([[ 14.61727248, 66.82869255], [ 28.65180043, 130.99313602], [ 41.54407151, 189.93529653]]))
def test_top_vs_time_drn_1(self): assert_allclose(dim1sin_E_Igamv_the_BC_deltaf_linear( **{'drn': 1, 'm': self.m, 'eigs':self.eigs, 'tvals':self.tvals,'Igamv':self.Igamv, 'zvals': self.zvals, 'pseudo_k':self.pseudo_k, 'top_vs_time': [self.top_vs_time], 'bot_vs_time':None}), np.array([[ 73.08636239, 334.14346275], [ 143.25900215, 654.9656801 ], [ 207.72035757, 949.67648264]]))
def test_top_vs_time_drn_0(self): assert_allclose(dim1sin_E_Igamv_the_BC_deltaf_linear( **{'drn': 0, 'm': self.m, 'eigs':self.eigs, 'tvals':self.tvals,'Igamv':self.Igamv, 'zvals': self.zvals, 'pseudo_k':self.pseudo_k, 'top_vs_time': [self.top_vs_time], 'bot_vs_time':None}), np.array([[ 58.46908991, 267.3147702 ], [ 114.60720172, 523.97254408], [ 166.17628606, 759.74118611]]))
def test_no_BC(self): assert_allclose(dim1sin_E_Igamv_the_BC_deltaf_linear( **{'drn': 0, 'm': self.m, 'eigs':self.eigs, 'tvals':self.tvals,'Igamv':self.Igamv, 'zvals': self.zvals, 'pseudo_k':self.pseudo_k, 'top_vs_time': None, 'bot_vs_time':None}), np.array([[ 0., 0.], [ 0., 0.], [ 0., 0.]]))
def test_bot_vs_time_drn_0_omega_phase(self): assert_allclose(dim1sin_E_Igamv_the_BC_deltaf_linear( **{'drn': 0, 'm': self.m, 'eigs':self.eigs, 'tvals':self.tvals,'Igamv':self.Igamv, 'zvals': self.zvals, 'pseudo_k':self.pseudo_k, 'top_vs_time': None, 'bot_vs_time':[self.bot_vs_time], 'bot_omega_phase': [self.omega_phase]}), np.array([[-13.02239055, -20.40017837], [-25.52561949, -39.98706601], [-37.01122253, -57.97979552]]))
def test_top_vs_time_drn_0_omega_phase(self): assert_allclose(dim1sin_E_Igamv_the_BC_deltaf_linear( **{'drn': 0, 'm': self.m, 'eigs':self.eigs, 'tvals':self.tvals,'Igamv':self.Igamv, 'zvals': self.zvals, 'pseudo_k':self.pseudo_k, 'top_vs_time': [self.top_vs_time], 'bot_vs_time':None, 'top_omega_phase': [self.omega_phase]}), np.array([[ -52.08956221, -81.6007135 ], [-102.10247796, -159.94826406], [-148.04489011, -231.91918207]]))
def gen_dim1sin_E_Igamv_the_BC_deltaf_linear(): """Test case generation for dim1sin_E_Igamv_the_BC_deltaf_linear. 2014-03-22""" #dim1sin_E_Igamv_the_BC_deltaf_linear(drn, m, eigs, tvals, Igamv, zvals, pseudo_k, top_vs_time, bot_vs_time, top_omega_phase=None, bot_omega_phase=None, dT=1.0): outz = np.array([[0.2, 0.4], [0.4, 0.6]]) z1 = outz[:, 0] z2 = outz[:, 1] m = np.array([1.0,2.0, 3.0]) v_E_Igamv_the = np.ones((3,2), dtype=float) tvals = np.array([1.0, 3]) top_vs_time = PolyLine([0,2,4],[0,2,2]) bot_vs_time = PolyLine([0,2,4],[0,2,2]) omega_phase = (1,2) a = PolyLine([0, 1], [1, 2])# y = 1 + z b = PolyLine([0, 1], [1, 2])# y = 1 + z g = np.array([1.0,2.0])# this is interpolated from top_vs_time at t = 1, 3 Igamv = np.identity(3) eigs = np.ones(3) zvals=[0.2] pseudo_k=[1000] fn=OrderedDict() fn['no_bc'] = {'drn': 0, 'm': m, 'eigs':eigs, 'tvals':tvals,'Igamv':Igamv, 'zvals':zvals, 'pseudo_k': pseudo_k, 'top_vs_time': None, 'bot_vs_time':None} fn['top_vs_time_drn_0']={'drn': 0, 'm': m, 'eigs':eigs, 'tvals':tvals,'Igamv':Igamv, 'zvals':zvals, 'pseudo_k': pseudo_k, 'top_vs_time': [top_vs_time], 'bot_vs_time':None} fn['top_vs_time_drn_1']={'drn': 1, 'm': m, 'eigs':eigs, 'tvals':tvals,'Igamv':Igamv, 'zvals':zvals, 'pseudo_k': pseudo_k, 'top_vs_time': [top_vs_time], 'bot_vs_time':None} fn['bot_vs_time_drn_0']={'drn': 0, 'm': m, 'eigs':eigs, 'tvals':tvals,'Igamv':Igamv, 'zvals':zvals, 'pseudo_k': pseudo_k, 'top_vs_time': None,'bot_vs_time':[bot_vs_time]} fn['bot_vs_time_top_vs_time_drn_1']={'drn': 1, 'm': m, 'eigs':eigs, 'tvals':tvals,'Igamv':Igamv, 'zvals':zvals, 'pseudo_k': pseudo_k, 'top_vs_time': [top_vs_time],'bot_vs_time':[bot_vs_time]} fn['test_top_vs_time_drn_0_omega_phase']={'drn': 0, 'm': m, 'eigs':eigs, 'tvals':tvals,'Igamv':Igamv, 'zvals':zvals, 'pseudo_k': pseudo_k, 'top_vs_time': [top_vs_time], 'bot_vs_time':None, 'top_omega_phase': [omega_phase]} fn['test_bot_vs_time_drn_0_omega_phase']={'drn': 0, 'm': m, 'eigs':eigs, 'tvals':tvals,'Igamv':Igamv, 'zvals':zvals, 'pseudo_k': pseudo_k, 'top_vs_time': None, 'bot_vs_time':[bot_vs_time], 'bot_omega_phase': [omega_phase]} fn['bot_vs_time_top_vs_time_drn_1_double_loads']={'drn': 1, 'm': m, 'eigs':eigs, 'tvals':tvals,'Igamv':Igamv, 'zvals':zvals, 'pseudo_k': pseudo_k, 'top_vs_time': [top_vs_time, top_vs_time],'bot_vs_time':[bot_vs_time, bot_vs_time]} # print(dim1sin_E_Igamv_the_BC_abf_linear(0,m,eigs,tvals,Igamv, a,b, top_vs_time=None, bot_vs_time=None)) for k, v in fn.items(): print(k) print (dim1sin_E_Igamv_the_BC_deltaf_linear(**v)) print('#'*10+'\n')
def gen_dim1sin_E_Igamv_the_BC_deltaf_linear(): """Test case generation for dim1sin_E_Igamv_the_BC_deltaf_linear. 2014-03-22""" #dim1sin_E_Igamv_the_BC_deltaf_linear(drn, m, eigs, tvals, Igamv, zvals, pseudo_k, top_vs_time, bot_vs_time, top_omega_phase=None, bot_omega_phase=None, dT=1.0): outz = np.array([[0.2, 0.4], [0.4, 0.6]]) z1 = outz[:, 0] z2 = outz[:, 1] m = np.array([1.0, 2.0, 3.0]) v_E_Igamv_the = np.ones((3, 2), dtype=float) tvals = np.array([1.0, 3]) top_vs_time = PolyLine([0, 2, 4], [0, 2, 2]) bot_vs_time = PolyLine([0, 2, 4], [0, 2, 2]) omega_phase = (1, 2) a = PolyLine([0, 1], [1, 2]) # y = 1 + z b = PolyLine([0, 1], [1, 2]) # y = 1 + z g = np.array([1.0, 2.0]) # this is interpolated from top_vs_time at t = 1, 3 Igamv = np.identity(3) eigs = np.ones(3) zvals = [0.2] pseudo_k = [1000] fn = OrderedDict() fn['no_bc'] = { 'drn': 0, 'm': m, 'eigs': eigs, 'tvals': tvals, 'Igamv': Igamv, 'zvals': zvals, 'pseudo_k': pseudo_k, 'top_vs_time': None, 'bot_vs_time': None } fn['top_vs_time_drn_0'] = { 'drn': 0, 'm': m, 'eigs': eigs, 'tvals': tvals, 'Igamv': Igamv, 'zvals': zvals, 'pseudo_k': pseudo_k, 'top_vs_time': [top_vs_time], 'bot_vs_time': None } fn['top_vs_time_drn_1'] = { 'drn': 1, 'm': m, 'eigs': eigs, 'tvals': tvals, 'Igamv': Igamv, 'zvals': zvals, 'pseudo_k': pseudo_k, 'top_vs_time': [top_vs_time], 'bot_vs_time': None } fn['bot_vs_time_drn_0'] = { 'drn': 0, 'm': m, 'eigs': eigs, 'tvals': tvals, 'Igamv': Igamv, 'zvals': zvals, 'pseudo_k': pseudo_k, 'top_vs_time': None, 'bot_vs_time': [bot_vs_time] } fn['bot_vs_time_top_vs_time_drn_1'] = { 'drn': 1, 'm': m, 'eigs': eigs, 'tvals': tvals, 'Igamv': Igamv, 'zvals': zvals, 'pseudo_k': pseudo_k, 'top_vs_time': [top_vs_time], 'bot_vs_time': [bot_vs_time] } fn['test_top_vs_time_drn_0_omega_phase'] = { 'drn': 0, 'm': m, 'eigs': eigs, 'tvals': tvals, 'Igamv': Igamv, 'zvals': zvals, 'pseudo_k': pseudo_k, 'top_vs_time': [top_vs_time], 'bot_vs_time': None, 'top_omega_phase': [omega_phase] } fn['test_bot_vs_time_drn_0_omega_phase'] = { 'drn': 0, 'm': m, 'eigs': eigs, 'tvals': tvals, 'Igamv': Igamv, 'zvals': zvals, 'pseudo_k': pseudo_k, 'top_vs_time': None, 'bot_vs_time': [bot_vs_time], 'bot_omega_phase': [omega_phase] } fn['bot_vs_time_top_vs_time_drn_1_double_loads'] = { 'drn': 1, 'm': m, 'eigs': eigs, 'tvals': tvals, 'Igamv': Igamv, 'zvals': zvals, 'pseudo_k': pseudo_k, 'top_vs_time': [top_vs_time, top_vs_time], 'bot_vs_time': [bot_vs_time, bot_vs_time] } # print(dim1sin_E_Igamv_the_BC_abf_linear(0,m,eigs,tvals,Igamv, a,b, top_vs_time=None, bot_vs_time=None)) for k, v in fn.items(): print(k) print(dim1sin_E_Igamv_the_BC_deltaf_linear(**v)) print('#' * 10 + '\n')