def thomsen_voigt_2d(dummy, keys, vals):
input = Struct(zip(keys, vals))
output = Struct()
vp = input.vp
vs = input.vs
rho = input.rho
epslion = input.epsilon
delta = input.delta
output.c11 = rho * vp**2. * (1. + 2. * epsilon)
output.c12 = rho * vp**2. * (1. + 2. * epsilon) - 2. * rho * vs**2. * (
1. + 2. * epsilon)
output.c13 = rho * vp**2. * (1. + 2. * delta) - 2. * rho * vs**2.
output.c33 = rho * vp**2.
output.c55 = rho * vs**2.
return output
def tti_voight_2d(dummy, keys, vals):
input = Struct(zip(keys, vals))
output = Struct()
vp = input.vp
vs = input.vs
rho = input.rho
epsilon = input.epsilon
delta = input.delta
theta = input.theta
c11 = rho * vp**2. * (1. + 2. * epsilon)
c12 = rho * vp**2. * (1. + 2. * epsilon) - 2. * rho * vs**2. * (
1. + 2. * epsilon)
c13 = rho * vp**2. * (1. + 2. * delta) - 2. * rho * vs**2.
c33 = rho * vp**2.
c55 = rho * vs**2.
sint = sin(PI / 180. * theta)
cost = cos(PI / 180. * theta)
sin2t = sin(2. * PI / 180. * theta)
cos2t = cos(2. * PI / 180. * theta)
output.c11 = c11 * cost**4 + c33 * sint**4 + 2 * c13 * sint**2 * cost**2 + c55 * sin2t**2
output.c12 = 1.e-6
output.c13 = c13 * (cost**4 + sint**4) + (
c11 + c33) * sint**2 * cost**2 - c55 * sin2t**2
output.c15 = ((c11 - c13) * cost**2 +
(c13 - c33) * sint**2) * sint * cost - c55 * sin2t * cos2t
output.c23 = 1.e-6
output.c25 = 0.0
output.c33 = c11 * sint**4 + c33 * cost**4 + 2 * c13 * sint**2 * cost**2 + c55 * sin2t**2
output.c35 = ((c11 - c13) * sint**2 +
(c13 - c33) * cost**2) * sint * cost + c55 * sin2t * cos2t
output.c55 = (c11 - 2. * c13 + c33) * sint**2 * cost**2 + c55 * cos2t**2
return output