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')
