示例#1
0
def force_dr(F0,f0,t):
	f = np.zeros(len(t))
	for ii in range(len(t)):
		f[ii] = forces.driven(F0,f0,t[ii])

	return f
示例#2
0
def force_ts(option_mask,z,v,t,f0,k,Q,Rt,a0,E,H,mu,m,D,F0,zc,eta,sigmas,sigmat,landadeb,epsilonr,mutip,musample):
	f = np.zeros(len(z))
	for ii in range(len(z)):
		f[ii] = forces.sumatory(option_mask,z[ii],v[ii],t[ii],f0,k,Q,Rt,a0,E,H,mu,m,D,F0,zc,eta,sigmas,sigmat,landadeb,epsilonr,mutip,musample) - forces.driven(F0,f0,t[ii])

	return f