# Authors:
#     Loic Gouarin <loic.gouarin@cmap.polytechnique.fr>
#     Nicole Spillane <nicole.spillane@cmap.polytechnique.fr>
#
# License: BSD 3 clause
# cython: boundscheck = False
# cython: wraparound = False
# cython: cdivision = True
import numpy as np
def getMatElemElasticity_2d_nonconst(double hx, double hy, double[:] lamb, double[:] mu):
	output = np.empty((lamb.size, 8, 8))
	cdef:
		int ie
		int size = lamb.size
		double[:, :, ::1] vout = output
	for ie in range(size):
		vout[ie, 0, 0] = hx*mu[ie]/(3*hy) + hy*lamb[ie]/(3*hx) + 2*hy*mu[ie]/(3*hx)
		vout[ie, 0, 1] = lamb[ie]/4 + mu[ie]/4
		vout[ie, 0, 2] = hx*mu[ie]/(6*hy) - hy*lamb[ie]/(3*hx) - 2*hy*mu[ie]/(3*hx)
		vout[ie, 0, 3] = lamb[ie]/4 - mu[ie]/4
		vout[ie, 0, 4] = -hx*mu[ie]/(6*hy) - hy*lamb[ie]/(6*hx) - hy*mu[ie]/(3*hx)
		vout[ie, 0, 5] = -lamb[ie]/4 - mu[ie]/4
		vout[ie, 0, 6] = -hx*mu[ie]/(3*hy) + hy*lamb[ie]/(6*hx) + hy*mu[ie]/(3*hx)
		vout[ie, 0, 7] = -lamb[ie]/4 + mu[ie]/4
		vout[ie, 1, 0] = lamb[ie]/4 + mu[ie]/4
		vout[ie, 1, 1] = hx*lamb[ie]/(3*hy) + 2*hx*mu[ie]/(3*hy) + hy*mu[ie]/(3*hx)
		vout[ie, 1, 2] = -lamb[ie]/4 + mu[ie]/4
		vout[ie, 1, 3] = hx*lamb[ie]/(6*hy) + hx*mu[ie]/(3*hy) - hy*mu[ie]/(3*hx)
		vout[ie, 1, 4] = -lamb[ie]/4 - mu[ie]/4
		vout[ie, 1, 5] = -hx*lamb[ie]/(6*hy) - hx*mu[ie]/(3*hy) - hy*mu[ie]/(6*hx)
		vout[ie, 1, 6] = lamb[ie]/4 - mu[ie]/4
		vout[ie, 1, 7] = -hx*lamb[ie]/(3*hy) - 2*hx*mu[ie]/(3*hy) + hy*mu[ie]/(6*hx)
		vout[ie, 2, 0] = hx*mu[ie]/(6*hy) - hy*lamb[ie]/(3*hx) - 2*hy*mu[ie]/(3*hx)
		vout[ie, 2, 1] = -lamb[ie]/4 + mu[ie]/4
		vout[ie, 2, 2] = hx*mu[ie]/(3*hy) + hy*lamb[ie]/(3*hx) + 2*hy*mu[ie]/(3*hx)
		vout[ie, 2, 3] = -lamb[ie]/4 - mu[ie]/4
		vout[ie, 2, 4] = -hx*mu[ie]/(3*hy) + hy*lamb[ie]/(6*hx) + hy*mu[ie]/(3*hx)
		vout[ie, 2, 5] = lamb[ie]/4 - mu[ie]/4
		vout[ie, 2, 6] = -hx*mu[ie]/(6*hy) - hy*lamb[ie]/(6*hx) - hy*mu[ie]/(3*hx)
		vout[ie, 2, 7] = lamb[ie]/4 + mu[ie]/4
		vout[ie, 3, 0] = lamb[ie]/4 - mu[ie]/4
		vout[ie, 3, 1] = hx*lamb[ie]/(6*hy) + hx*mu[ie]/(3*hy) - hy*mu[ie]/(3*hx)
		vout[ie, 3, 2] = -lamb[ie]/4 - mu[ie]/4
		vout[ie, 3, 3] = hx*lamb[ie]/(3*hy) + 2*hx*mu[ie]/(3*hy) + hy*mu[ie]/(3*hx)
		vout[ie, 3, 4] = -lamb[ie]/4 + mu[ie]/4
		vout[ie, 3, 5] = -hx*lamb[ie]/(3*hy) - 2*hx*mu[ie]/(3*hy) + hy*mu[ie]/(6*hx)
		vout[ie, 3, 6] = lamb[ie]/4 + mu[ie]/4
		vout[ie, 3, 7] = -hx*lamb[ie]/(6*hy) - hx*mu[ie]/(3*hy) - hy*mu[ie]/(6*hx)
		vout[ie, 4, 0] = -hx*mu[ie]/(6*hy) - hy*lamb[ie]/(6*hx) - hy*mu[ie]/(3*hx)
		vout[ie, 4, 1] = -lamb[ie]/4 - mu[ie]/4
		vout[ie, 4, 2] = -hx*mu[ie]/(3*hy) + hy*lamb[ie]/(6*hx) + hy*mu[ie]/(3*hx)
		vout[ie, 4, 3] = -lamb[ie]/4 + mu[ie]/4
		vout[ie, 4, 4] = hx*mu[ie]/(3*hy) + hy*lamb[ie]/(3*hx) + 2*hy*mu[ie]/(3*hx)
		vout[ie, 4, 5] = lamb[ie]/4 + mu[ie]/4
		vout[ie, 4, 6] = hx*mu[ie]/(6*hy) - hy*lamb[ie]/(3*hx) - 2*hy*mu[ie]/(3*hx)
		vout[ie, 4, 7] = lamb[ie]/4 - mu[ie]/4
		vout[ie, 5, 0] = -lamb[ie]/4 - mu[ie]/4
		vout[ie, 5, 1] = -hx*lamb[ie]/(6*hy) - hx*mu[ie]/(3*hy) - hy*mu[ie]/(6*hx)
		vout[ie, 5, 2] = lamb[ie]/4 - mu[ie]/4
		vout[ie, 5, 3] = -hx*lamb[ie]/(3*hy) - 2*hx*mu[ie]/(3*hy) + hy*mu[ie]/(6*hx)
		vout[ie, 5, 4] = lamb[ie]/4 + mu[ie]/4
		vout[ie, 5, 5] = hx*lamb[ie]/(3*hy) + 2*hx*mu[ie]/(3*hy) + hy*mu[ie]/(3*hx)
		vout[ie, 5, 6] = -lamb[ie]/4 + mu[ie]/4
		vout[ie, 5, 7] = hx*lamb[ie]/(6*hy) + hx*mu[ie]/(3*hy) - hy*mu[ie]/(3*hx)
		vout[ie, 6, 0] = -hx*mu[ie]/(3*hy) + hy*lamb[ie]/(6*hx) + hy*mu[ie]/(3*hx)
		vout[ie, 6, 1] = lamb[ie]/4 - mu[ie]/4
		vout[ie, 6, 2] = -hx*mu[ie]/(6*hy) - hy*lamb[ie]/(6*hx) - hy*mu[ie]/(3*hx)
		vout[ie, 6, 3] = lamb[ie]/4 + mu[ie]/4
		vout[ie, 6, 4] = hx*mu[ie]/(6*hy) - hy*lamb[ie]/(3*hx) - 2*hy*mu[ie]/(3*hx)
		vout[ie, 6, 5] = -lamb[ie]/4 + mu[ie]/4
		vout[ie, 6, 6] = hx*mu[ie]/(3*hy) + hy*lamb[ie]/(3*hx) + 2*hy*mu[ie]/(3*hx)
		vout[ie, 6, 7] = -lamb[ie]/4 - mu[ie]/4
		vout[ie, 7, 0] = -lamb[ie]/4 + mu[ie]/4
		vout[ie, 7, 1] = -hx*lamb[ie]/(3*hy) - 2*hx*mu[ie]/(3*hy) + hy*mu[ie]/(6*hx)
		vout[ie, 7, 2] = lamb[ie]/4 + mu[ie]/4
		vout[ie, 7, 3] = -hx*lamb[ie]/(6*hy) - hx*mu[ie]/(3*hy) - hy*mu[ie]/(6*hx)
		vout[ie, 7, 4] = lamb[ie]/4 - mu[ie]/4
		vout[ie, 7, 5] = hx*lamb[ie]/(6*hy) + hx*mu[ie]/(3*hy) - hy*mu[ie]/(3*hx)
		vout[ie, 7, 6] = -lamb[ie]/4 - mu[ie]/4
		vout[ie, 7, 7] = hx*lamb[ie]/(3*hy) + 2*hx*mu[ie]/(3*hy) + hy*mu[ie]/(3*hx)
	return output

def getMatElemElasticity_2d_const(double hx, double hy, double lamb, double mu):
	output = np.empty((8, 8))
	cdef:
		double[:, ::1] vout = output
	vout[0, 0] = hx*mu/(3*hy) + hy*lamb/(3*hx) + 2*hy*mu/(3*hx)
	vout[0, 1] = lamb/4 + mu/4
	vout[0, 2] = hx*mu/(6*hy) - hy*lamb/(3*hx) - 2*hy*mu/(3*hx)
	vout[0, 3] = lamb/4 - mu/4
	vout[0, 4] = -hx*mu/(6*hy) - hy*lamb/(6*hx) - hy*mu/(3*hx)
	vout[0, 5] = -lamb/4 - mu/4
	vout[0, 6] = -hx*mu/(3*hy) + hy*lamb/(6*hx) + hy*mu/(3*hx)
	vout[0, 7] = -lamb/4 + mu/4
	vout[1, 0] = lamb/4 + mu/4
	vout[1, 1] = hx*lamb/(3*hy) + 2*hx*mu/(3*hy) + hy*mu/(3*hx)
	vout[1, 2] = -lamb/4 + mu/4
	vout[1, 3] = hx*lamb/(6*hy) + hx*mu/(3*hy) - hy*mu/(3*hx)
	vout[1, 4] = -lamb/4 - mu/4
	vout[1, 5] = -hx*lamb/(6*hy) - hx*mu/(3*hy) - hy*mu/(6*hx)
	vout[1, 6] = lamb/4 - mu/4
	vout[1, 7] = -hx*lamb/(3*hy) - 2*hx*mu/(3*hy) + hy*mu/(6*hx)
	vout[2, 0] = hx*mu/(6*hy) - hy*lamb/(3*hx) - 2*hy*mu/(3*hx)
	vout[2, 1] = -lamb/4 + mu/4
	vout[2, 2] = hx*mu/(3*hy) + hy*lamb/(3*hx) + 2*hy*mu/(3*hx)
	vout[2, 3] = -lamb/4 - mu/4
	vout[2, 4] = -hx*mu/(3*hy) + hy*lamb/(6*hx) + hy*mu/(3*hx)
	vout[2, 5] = lamb/4 - mu/4
	vout[2, 6] = -hx*mu/(6*hy) - hy*lamb/(6*hx) - hy*mu/(3*hx)
	vout[2, 7] = lamb/4 + mu/4
	vout[3, 0] = lamb/4 - mu/4
	vout[3, 1] = hx*lamb/(6*hy) + hx*mu/(3*hy) - hy*mu/(3*hx)
	vout[3, 2] = -lamb/4 - mu/4
	vout[3, 3] = hx*lamb/(3*hy) + 2*hx*mu/(3*hy) + hy*mu/(3*hx)
	vout[3, 4] = -lamb/4 + mu/4
	vout[3, 5] = -hx*lamb/(3*hy) - 2*hx*mu/(3*hy) + hy*mu/(6*hx)
	vout[3, 6] = lamb/4 + mu/4
	vout[3, 7] = -hx*lamb/(6*hy) - hx*mu/(3*hy) - hy*mu/(6*hx)
	vout[4, 0] = -hx*mu/(6*hy) - hy*lamb/(6*hx) - hy*mu/(3*hx)
	vout[4, 1] = -lamb/4 - mu/4
	vout[4, 2] = -hx*mu/(3*hy) + hy*lamb/(6*hx) + hy*mu/(3*hx)
	vout[4, 3] = -lamb/4 + mu/4
	vout[4, 4] = hx*mu/(3*hy) + hy*lamb/(3*hx) + 2*hy*mu/(3*hx)
	vout[4, 5] = lamb/4 + mu/4
	vout[4, 6] = hx*mu/(6*hy) - hy*lamb/(3*hx) - 2*hy*mu/(3*hx)
	vout[4, 7] = lamb/4 - mu/4
	vout[5, 0] = -lamb/4 - mu/4
	vout[5, 1] = -hx*lamb/(6*hy) - hx*mu/(3*hy) - hy*mu/(6*hx)
	vout[5, 2] = lamb/4 - mu/4
	vout[5, 3] = -hx*lamb/(3*hy) - 2*hx*mu/(3*hy) + hy*mu/(6*hx)
	vout[5, 4] = lamb/4 + mu/4
	vout[5, 5] = hx*lamb/(3*hy) + 2*hx*mu/(3*hy) + hy*mu/(3*hx)
	vout[5, 6] = -lamb/4 + mu/4
	vout[5, 7] = hx*lamb/(6*hy) + hx*mu/(3*hy) - hy*mu/(3*hx)
	vout[6, 0] = -hx*mu/(3*hy) + hy*lamb/(6*hx) + hy*mu/(3*hx)
	vout[6, 1] = lamb/4 - mu/4
	vout[6, 2] = -hx*mu/(6*hy) - hy*lamb/(6*hx) - hy*mu/(3*hx)
	vout[6, 3] = lamb/4 + mu/4
	vout[6, 4] = hx*mu/(6*hy) - hy*lamb/(3*hx) - 2*hy*mu/(3*hx)
	vout[6, 5] = -lamb/4 + mu/4
	vout[6, 6] = hx*mu/(3*hy) + hy*lamb/(3*hx) + 2*hy*mu/(3*hx)
	vout[6, 7] = -lamb/4 - mu/4
	vout[7, 0] = -lamb/4 + mu/4
	vout[7, 1] = -hx*lamb/(3*hy) - 2*hx*mu/(3*hy) + hy*mu/(6*hx)
	vout[7, 2] = lamb/4 + mu/4
	vout[7, 3] = -hx*lamb/(6*hy) - hx*mu/(3*hy) - hy*mu/(6*hx)
	vout[7, 4] = lamb/4 - mu/4
	vout[7, 5] = hx*lamb/(6*hy) + hx*mu/(3*hy) - hy*mu/(3*hx)
	vout[7, 6] = -lamb/4 - mu/4
	vout[7, 7] = hx*lamb/(3*hy) + 2*hx*mu/(3*hy) + hy*mu/(3*hx)
	return output

def getMatElemElasticity_3d_nonconst(double hx, double hy, double hz, double[:] lamb, double[:] mu):
	output = np.empty((lamb.size, 24, 24))
	cdef:
		int ie
		int size = lamb.size
		double[:, :, ::1] vout = output
	for ie in range(size):
		vout[ie, 0, 0] = hx*hy*mu[ie]/(9*hz) + hx*hz*mu[ie]/(9*hy) + hy*hz*lamb[ie]/(9*hx) + 2*hy*hz*mu[ie]/(9*hx)
		vout[ie, 0, 1] = hz*lamb[ie]/12 + hz*mu[ie]/12
		vout[ie, 0, 2] = hy*lamb[ie]/12 + hy*mu[ie]/12
		vout[ie, 0, 3] = hx*hy*mu[ie]/(18*hz) + hx*hz*mu[ie]/(18*hy) - hy*hz*lamb[ie]/(9*hx) - 2*hy*hz*mu[ie]/(9*hx)
		vout[ie, 0, 4] = hz*lamb[ie]/12 - hz*mu[ie]/12
		vout[ie, 0, 5] = hy*lamb[ie]/12 - hy*mu[ie]/12
		vout[ie, 0, 6] = hx*hy*mu[ie]/(36*hz) - hx*hz*mu[ie]/(18*hy) - hy*hz*lamb[ie]/(18*hx) - hy*hz*mu[ie]/(9*hx)
		vout[ie, 0, 7] = -hz*lamb[ie]/12 - hz*mu[ie]/12
		vout[ie, 0, 8] = hy*lamb[ie]/24 - hy*mu[ie]/24
		vout[ie, 0, 9] = hx*hy*mu[ie]/(18*hz) - hx*hz*mu[ie]/(9*hy) + hy*hz*lamb[ie]/(18*hx) + hy*hz*mu[ie]/(9*hx)
		vout[ie, 0, 10] = -hz*lamb[ie]/12 + hz*mu[ie]/12
		vout[ie, 0, 11] = hy*lamb[ie]/24 + hy*mu[ie]/24
		vout[ie, 0, 12] = -hx*hy*mu[ie]/(9*hz) + hx*hz*mu[ie]/(18*hy) + hy*hz*lamb[ie]/(18*hx) + hy*hz*mu[ie]/(9*hx)
		vout[ie, 0, 13] = hz*lamb[ie]/24 + hz*mu[ie]/24
		vout[ie, 0, 14] = -hy*lamb[ie]/12 + hy*mu[ie]/12
		vout[ie, 0, 15] = -hx*hy*mu[ie]/(18*hz) + hx*hz*mu[ie]/(36*hy) - hy*hz*lamb[ie]/(18*hx) - hy*hz*mu[ie]/(9*hx)
		vout[ie, 0, 16] = hz*lamb[ie]/24 - hz*mu[ie]/24
		vout[ie, 0, 17] = -hy*lamb[ie]/12 - hy*mu[ie]/12
		vout[ie, 0, 18] = -hx*hy*mu[ie]/(36*hz) - hx*hz*mu[ie]/(36*hy) - hy*hz*lamb[ie]/(36*hx) - hy*hz*mu[ie]/(18*hx)
		vout[ie, 0, 19] = -hz*lamb[ie]/24 - hz*mu[ie]/24
		vout[ie, 0, 20] = -hy*lamb[ie]/24 - hy*mu[ie]/24
		vout[ie, 0, 21] = -hx*hy*mu[ie]/(18*hz) - hx*hz*mu[ie]/(18*hy) + hy*hz*lamb[ie]/(36*hx) + hy*hz*mu[ie]/(18*hx)
		vout[ie, 0, 22] = -hz*lamb[ie]/24 + hz*mu[ie]/24
		vout[ie, 0, 23] = -hy*lamb[ie]/24 + hy*mu[ie]/24
		vout[ie, 1, 0] = hz*lamb[ie]/12 + hz*mu[ie]/12
		vout[ie, 1, 1] = hx*hy*mu[ie]/(9*hz) + hx*hz*lamb[ie]/(9*hy) + 2*hx*hz*mu[ie]/(9*hy) + hy*hz*mu[ie]/(9*hx)
		vout[ie, 1, 2] = hx*lamb[ie]/12 + hx*mu[ie]/12
		vout[ie, 1, 3] = -hz*lamb[ie]/12 + hz*mu[ie]/12
		vout[ie, 1, 4] = hx*hy*mu[ie]/(18*hz) + hx*hz*lamb[ie]/(18*hy) + hx*hz*mu[ie]/(9*hy) - hy*hz*mu[ie]/(9*hx)
		vout[ie, 1, 5] = hx*lamb[ie]/24 + hx*mu[ie]/24
		vout[ie, 1, 6] = -hz*lamb[ie]/12 - hz*mu[ie]/12
		vout[ie, 1, 7] = hx*hy*mu[ie]/(36*hz) - hx*hz*lamb[ie]/(18*hy) - hx*hz*mu[ie]/(9*hy) - hy*hz*mu[ie]/(18*hx)
		vout[ie, 1, 8] = hx*lamb[ie]/24 - hx*mu[ie]/24
		vout[ie, 1, 9] = hz*lamb[ie]/12 - hz*mu[ie]/12
		vout[ie, 1, 10] = hx*hy*mu[ie]/(18*hz) - hx*hz*lamb[ie]/(9*hy) - 2*hx*hz*mu[ie]/(9*hy) + hy*hz*mu[ie]/(18*hx)
		vout[ie, 1, 11] = hx*lamb[ie]/12 - hx*mu[ie]/12
		vout[ie, 1, 12] = hz*lamb[ie]/24 + hz*mu[ie]/24
		vout[ie, 1, 13] = -hx*hy*mu[ie]/(9*hz) + hx*hz*lamb[ie]/(18*hy) + hx*hz*mu[ie]/(9*hy) + hy*hz*mu[ie]/(18*hx)
		vout[ie, 1, 14] = -hx*lamb[ie]/12 + hx*mu[ie]/12
		vout[ie, 1, 15] = -hz*lamb[ie]/24 + hz*mu[ie]/24
		vout[ie, 1, 16] = -hx*hy*mu[ie]/(18*hz) + hx*hz*lamb[ie]/(36*hy) + hx*hz*mu[ie]/(18*hy) - hy*hz*mu[ie]/(18*hx)
		vout[ie, 1, 17] = -hx*lamb[ie]/24 + hx*mu[ie]/24
		vout[ie, 1, 18] = -hz*lamb[ie]/24 - hz*mu[ie]/24
		vout[ie, 1, 19] = -hx*hy*mu[ie]/(36*hz) - hx*hz*lamb[ie]/(36*hy) - hx*hz*mu[ie]/(18*hy) - hy*hz*mu[ie]/(36*hx)
		vout[ie, 1, 20] = -hx*lamb[ie]/24 - hx*mu[ie]/24
		vout[ie, 1, 21] = hz*lamb[ie]/24 - hz*mu[ie]/24
		vout[ie, 1, 22] = -hx*hy*mu[ie]/(18*hz) - hx*hz*lamb[ie]/(18*hy) - hx*hz*mu[ie]/(9*hy) + hy*hz*mu[ie]/(36*hx)
		vout[ie, 1, 23] = -hx*lamb[ie]/12 - hx*mu[ie]/12
		vout[ie, 2, 0] = hy*lamb[ie]/12 + hy*mu[ie]/12
		vout[ie, 2, 1] = hx*lamb[ie]/12 + hx*mu[ie]/12
		vout[ie, 2, 2] = hx*hy*lamb[ie]/(9*hz) + 2*hx*hy*mu[ie]/(9*hz) + hx*hz*mu[ie]/(9*hy) + hy*hz*mu[ie]/(9*hx)
		vout[ie, 2, 3] = -hy*lamb[ie]/12 + hy*mu[ie]/12
		vout[ie, 2, 4] = hx*lamb[ie]/24 + hx*mu[ie]/24
		vout[ie, 2, 5] = hx*hy*lamb[ie]/(18*hz) + hx*hy*mu[ie]/(9*hz) + hx*hz*mu[ie]/(18*hy) - hy*hz*mu[ie]/(9*hx)
		vout[ie, 2, 6] = -hy*lamb[ie]/24 + hy*mu[ie]/24
		vout[ie, 2, 7] = -hx*lamb[ie]/24 + hx*mu[ie]/24
		vout[ie, 2, 8] = hx*hy*lamb[ie]/(36*hz) + hx*hy*mu[ie]/(18*hz) - hx*hz*mu[ie]/(18*hy) - hy*hz*mu[ie]/(18*hx)
		vout[ie, 2, 9] = hy*lamb[ie]/24 + hy*mu[ie]/24
		vout[ie, 2, 10] = -hx*lamb[ie]/12 + hx*mu[ie]/12
		vout[ie, 2, 11] = hx*hy*lamb[ie]/(18*hz) + hx*hy*mu[ie]/(9*hz) - hx*hz*mu[ie]/(9*hy) + hy*hz*mu[ie]/(18*hx)
		vout[ie, 2, 12] = hy*lamb[ie]/12 - hy*mu[ie]/12
		vout[ie, 2, 13] = hx*lamb[ie]/12 - hx*mu[ie]/12
		vout[ie, 2, 14] = -hx*hy*lamb[ie]/(9*hz) - 2*hx*hy*mu[ie]/(9*hz) + hx*hz*mu[ie]/(18*hy) + hy*hz*mu[ie]/(18*hx)
		vout[ie, 2, 15] = -hy*lamb[ie]/12 - hy*mu[ie]/12
		vout[ie, 2, 16] = hx*lamb[ie]/24 - hx*mu[ie]/24
		vout[ie, 2, 17] = -hx*hy*lamb[ie]/(18*hz) - hx*hy*mu[ie]/(9*hz) + hx*hz*mu[ie]/(36*hy) - hy*hz*mu[ie]/(18*hx)
		vout[ie, 2, 18] = -hy*lamb[ie]/24 - hy*mu[ie]/24
		vout[ie, 2, 19] = -hx*lamb[ie]/24 - hx*mu[ie]/24
		vout[ie, 2, 20] = -hx*hy*lamb[ie]/(36*hz) - hx*hy*mu[ie]/(18*hz) - hx*hz*mu[ie]/(36*hy) - hy*hz*mu[ie]/(36*hx)
		vout[ie, 2, 21] = hy*lamb[ie]/24 - hy*mu[ie]/24
		vout[ie, 2, 22] = -hx*lamb[ie]/12 - hx*mu[ie]/12
		vout[ie, 2, 23] = -hx*hy*lamb[ie]/(18*hz) - hx*hy*mu[ie]/(9*hz) - hx*hz*mu[ie]/(18*hy) + hy*hz*mu[ie]/(36*hx)
		vout[ie, 3, 0] = hx*hy*mu[ie]/(18*hz) + hx*hz*mu[ie]/(18*hy) - hy*hz*lamb[ie]/(9*hx) - 2*hy*hz*mu[ie]/(9*hx)
		vout[ie, 3, 1] = -hz*lamb[ie]/12 + hz*mu[ie]/12
		vout[ie, 3, 2] = -hy*lamb[ie]/12 + hy*mu[ie]/12
		vout[ie, 3, 3] = hx*hy*mu[ie]/(9*hz) + hx*hz*mu[ie]/(9*hy) + hy*hz*lamb[ie]/(9*hx) + 2*hy*hz*mu[ie]/(9*hx)
		vout[ie, 3, 4] = -hz*lamb[ie]/12 - hz*mu[ie]/12
		vout[ie, 3, 5] = -hy*lamb[ie]/12 - hy*mu[ie]/12
		vout[ie, 3, 6] = hx*hy*mu[ie]/(18*hz) - hx*hz*mu[ie]/(9*hy) + hy*hz*lamb[ie]/(18*hx) + hy*hz*mu[ie]/(9*hx)
		vout[ie, 3, 7] = hz*lamb[ie]/12 - hz*mu[ie]/12
		vout[ie, 3, 8] = -hy*lamb[ie]/24 - hy*mu[ie]/24
		vout[ie, 3, 9] = hx*hy*mu[ie]/(36*hz) - hx*hz*mu[ie]/(18*hy) - hy*hz*lamb[ie]/(18*hx) - hy*hz*mu[ie]/(9*hx)
		vout[ie, 3, 10] = hz*lamb[ie]/12 + hz*mu[ie]/12
		vout[ie, 3, 11] = -hy*lamb[ie]/24 + hy*mu[ie]/24
		vout[ie, 3, 12] = -hx*hy*mu[ie]/(18*hz) + hx*hz*mu[ie]/(36*hy) - hy*hz*lamb[ie]/(18*hx) - hy*hz*mu[ie]/(9*hx)
		vout[ie, 3, 13] = -hz*lamb[ie]/24 + hz*mu[ie]/24
		vout[ie, 3, 14] = hy*lamb[ie]/12 + hy*mu[ie]/12
		vout[ie, 3, 15] = -hx*hy*mu[ie]/(9*hz) + hx*hz*mu[ie]/(18*hy) + hy*hz*lamb[ie]/(18*hx) + hy*hz*mu[ie]/(9*hx)
		vout[ie, 3, 16] = -hz*lamb[ie]/24 - hz*mu[ie]/24
		vout[ie, 3, 17] = hy*lamb[ie]/12 - hy*mu[ie]/12
		vout[ie, 3, 18] = -hx*hy*mu[ie]/(18*hz) - hx*hz*mu[ie]/(18*hy) + hy*hz*lamb[ie]/(36*hx) + hy*hz*mu[ie]/(18*hx)
		vout[ie, 3, 19] = hz*lamb[ie]/24 - hz*mu[ie]/24
		vout[ie, 3, 20] = hy*lamb[ie]/24 - hy*mu[ie]/24
		vout[ie, 3, 21] = -hx*hy*mu[ie]/(36*hz) - hx*hz*mu[ie]/(36*hy) - hy*hz*lamb[ie]/(36*hx) - hy*hz*mu[ie]/(18*hx)
		vout[ie, 3, 22] = hz*lamb[ie]/24 + hz*mu[ie]/24
		vout[ie, 3, 23] = hy*lamb[ie]/24 + hy*mu[ie]/24
		vout[ie, 4, 0] = hz*lamb[ie]/12 - hz*mu[ie]/12
		vout[ie, 4, 1] = hx*hy*mu[ie]/(18*hz) + hx*hz*lamb[ie]/(18*hy) + hx*hz*mu[ie]/(9*hy) - hy*hz*mu[ie]/(9*hx)
		vout[ie, 4, 2] = hx*lamb[ie]/24 + hx*mu[ie]/24
		vout[ie, 4, 3] = -hz*lamb[ie]/12 - hz*mu[ie]/12
		vout[ie, 4, 4] = hx*hy*mu[ie]/(9*hz) + hx*hz*lamb[ie]/(9*hy) + 2*hx*hz*mu[ie]/(9*hy) + hy*hz*mu[ie]/(9*hx)
		vout[ie, 4, 5] = hx*lamb[ie]/12 + hx*mu[ie]/12
		vout[ie, 4, 6] = -hz*lamb[ie]/12 + hz*mu[ie]/12
		vout[ie, 4, 7] = hx*hy*mu[ie]/(18*hz) - hx*hz*lamb[ie]/(9*hy) - 2*hx*hz*mu[ie]/(9*hy) + hy*hz*mu[ie]/(18*hx)
		vout[ie, 4, 8] = hx*lamb[ie]/12 - hx*mu[ie]/12
		vout[ie, 4, 9] = hz*lamb[ie]/12 + hz*mu[ie]/12
		vout[ie, 4, 10] = hx*hy*mu[ie]/(36*hz) - hx*hz*lamb[ie]/(18*hy) - hx*hz*mu[ie]/(9*hy) - hy*hz*mu[ie]/(18*hx)
		vout[ie, 4, 11] = hx*lamb[ie]/24 - hx*mu[ie]/24
		vout[ie, 4, 12] = hz*lamb[ie]/24 - hz*mu[ie]/24
		vout[ie, 4, 13] = -hx*hy*mu[ie]/(18*hz) + hx*hz*lamb[ie]/(36*hy) + hx*hz*mu[ie]/(18*hy) - hy*hz*mu[ie]/(18*hx)
		vout[ie, 4, 14] = -hx*lamb[ie]/24 + hx*mu[ie]/24
		vout[ie, 4, 15] = -hz*lamb[ie]/24 - hz*mu[ie]/24
		vout[ie, 4, 16] = -hx*hy*mu[ie]/(9*hz) + hx*hz*lamb[ie]/(18*hy) + hx*hz*mu[ie]/(9*hy) + hy*hz*mu[ie]/(18*hx)
		vout[ie, 4, 17] = -hx*lamb[ie]/12 + hx*mu[ie]/12
		vout[ie, 4, 18] = -hz*lamb[ie]/24 + hz*mu[ie]/24
		vout[ie, 4, 19] = -hx*hy*mu[ie]/(18*hz) - hx*hz*lamb[ie]/(18*hy) - hx*hz*mu[ie]/(9*hy) + hy*hz*mu[ie]/(36*hx)
		vout[ie, 4, 20] = -hx*lamb[ie]/12 - hx*mu[ie]/12
		vout[ie, 4, 21] = hz*lamb[ie]/24 + hz*mu[ie]/24
		vout[ie, 4, 22] = -hx*hy*mu[ie]/(36*hz) - hx*hz*lamb[ie]/(36*hy) - hx*hz*mu[ie]/(18*hy) - hy*hz*mu[ie]/(36*hx)
		vout[ie, 4, 23] = -hx*lamb[ie]/24 - hx*mu[ie]/24
		vout[ie, 5, 0] = hy*lamb[ie]/12 - hy*mu[ie]/12
		vout[ie, 5, 1] = hx*lamb[ie]/24 + hx*mu[ie]/24
		vout[ie, 5, 2] = hx*hy*lamb[ie]/(18*hz) + hx*hy*mu[ie]/(9*hz) + hx*hz*mu[ie]/(18*hy) - hy*hz*mu[ie]/(9*hx)
		vout[ie, 5, 3] = -hy*lamb[ie]/12 - hy*mu[ie]/12
		vout[ie, 5, 4] = hx*lamb[ie]/12 + hx*mu[ie]/12
		vout[ie, 5, 5] = hx*hy*lamb[ie]/(9*hz) + 2*hx*hy*mu[ie]/(9*hz) + hx*hz*mu[ie]/(9*hy) + hy*hz*mu[ie]/(9*hx)
		vout[ie, 5, 6] = -hy*lamb[ie]/24 - hy*mu[ie]/24
		vout[ie, 5, 7] = -hx*lamb[ie]/12 + hx*mu[ie]/12
		vout[ie, 5, 8] = hx*hy*lamb[ie]/(18*hz) + hx*hy*mu[ie]/(9*hz) - hx*hz*mu[ie]/(9*hy) + hy*hz*mu[ie]/(18*hx)
		vout[ie, 5, 9] = hy*lamb[ie]/24 - hy*mu[ie]/24
		vout[ie, 5, 10] = -hx*lamb[ie]/24 + hx*mu[ie]/24
		vout[ie, 5, 11] = hx*hy*lamb[ie]/(36*hz) + hx*hy*mu[ie]/(18*hz) - hx*hz*mu[ie]/(18*hy) - hy*hz*mu[ie]/(18*hx)
		vout[ie, 5, 12] = hy*lamb[ie]/12 + hy*mu[ie]/12
		vout[ie, 5, 13] = hx*lamb[ie]/24 - hx*mu[ie]/24
		vout[ie, 5, 14] = -hx*hy*lamb[ie]/(18*hz) - hx*hy*mu[ie]/(9*hz) + hx*hz*mu[ie]/(36*hy) - hy*hz*mu[ie]/(18*hx)
		vout[ie, 5, 15] = -hy*lamb[ie]/12 + hy*mu[ie]/12
		vout[ie, 5, 16] = hx*lamb[ie]/12 - hx*mu[ie]/12
		vout[ie, 5, 17] = -hx*hy*lamb[ie]/(9*hz) - 2*hx*hy*mu[ie]/(9*hz) + hx*hz*mu[ie]/(18*hy) + hy*hz*mu[ie]/(18*hx)
		vout[ie, 5, 18] = -hy*lamb[ie]/24 + hy*mu[ie]/24
		vout[ie, 5, 19] = -hx*lamb[ie]/12 - hx*mu[ie]/12
		vout[ie, 5, 20] = -hx*hy*lamb[ie]/(18*hz) - hx*hy*mu[ie]/(9*hz) - hx*hz*mu[ie]/(18*hy) + hy*hz*mu[ie]/(36*hx)
		vout[ie, 5, 21] = hy*lamb[ie]/24 + hy*mu[ie]/24
		vout[ie, 5, 22] = -hx*lamb[ie]/24 - hx*mu[ie]/24
		vout[ie, 5, 23] = -hx*hy*lamb[ie]/(36*hz) - hx*hy*mu[ie]/(18*hz) - hx*hz*mu[ie]/(36*hy) - hy*hz*mu[ie]/(36*hx)
		vout[ie, 6, 0] = hx*hy*mu[ie]/(36*hz) - hx*hz*mu[ie]/(18*hy) - hy*hz*lamb[ie]/(18*hx) - hy*hz*mu[ie]/(9*hx)
		vout[ie, 6, 1] = -hz*lamb[ie]/12 - hz*mu[ie]/12
		vout[ie, 6, 2] = -hy*lamb[ie]/24 + hy*mu[ie]/24
		vout[ie, 6, 3] = hx*hy*mu[ie]/(18*hz) - hx*hz*mu[ie]/(9*hy) + hy*hz*lamb[ie]/(18*hx) + hy*hz*mu[ie]/(9*hx)
		vout[ie, 6, 4] = -hz*lamb[ie]/12 + hz*mu[ie]/12
		vout[ie, 6, 5] = -hy*lamb[ie]/24 - hy*mu[ie]/24
		vout[ie, 6, 6] = hx*hy*mu[ie]/(9*hz) + hx*hz*mu[ie]/(9*hy) + hy*hz*lamb[ie]/(9*hx) + 2*hy*hz*mu[ie]/(9*hx)
		vout[ie, 6, 7] = hz*lamb[ie]/12 + hz*mu[ie]/12
		vout[ie, 6, 8] = -hy*lamb[ie]/12 - hy*mu[ie]/12
		vout[ie, 6, 9] = hx*hy*mu[ie]/(18*hz) + hx*hz*mu[ie]/(18*hy) - hy*hz*lamb[ie]/(9*hx) - 2*hy*hz*mu[ie]/(9*hx)
		vout[ie, 6, 10] = hz*lamb[ie]/12 - hz*mu[ie]/12
		vout[ie, 6, 11] = -hy*lamb[ie]/12 + hy*mu[ie]/12
		vout[ie, 6, 12] = -hx*hy*mu[ie]/(36*hz) - hx*hz*mu[ie]/(36*hy) - hy*hz*lamb[ie]/(36*hx) - hy*hz*mu[ie]/(18*hx)
		vout[ie, 6, 13] = -hz*lamb[ie]/24 - hz*mu[ie]/24
		vout[ie, 6, 14] = hy*lamb[ie]/24 + hy*mu[ie]/24
		vout[ie, 6, 15] = -hx*hy*mu[ie]/(18*hz) - hx*hz*mu[ie]/(18*hy) + hy*hz*lamb[ie]/(36*hx) + hy*hz*mu[ie]/(18*hx)
		vout[ie, 6, 16] = -hz*lamb[ie]/24 + hz*mu[ie]/24
		vout[ie, 6, 17] = hy*lamb[ie]/24 - hy*mu[ie]/24
		vout[ie, 6, 18] = -hx*hy*mu[ie]/(9*hz) + hx*hz*mu[ie]/(18*hy) + hy*hz*lamb[ie]/(18*hx) + hy*hz*mu[ie]/(9*hx)
		vout[ie, 6, 19] = hz*lamb[ie]/24 + hz*mu[ie]/24
		vout[ie, 6, 20] = hy*lamb[ie]/12 - hy*mu[ie]/12
		vout[ie, 6, 21] = -hx*hy*mu[ie]/(18*hz) + hx*hz*mu[ie]/(36*hy) - hy*hz*lamb[ie]/(18*hx) - hy*hz*mu[ie]/(9*hx)
		vout[ie, 6, 22] = hz*lamb[ie]/24 - hz*mu[ie]/24
		vout[ie, 6, 23] = hy*lamb[ie]/12 + hy*mu[ie]/12
		vout[ie, 7, 0] = -hz*lamb[ie]/12 - hz*mu[ie]/12
		vout[ie, 7, 1] = hx*hy*mu[ie]/(36*hz) - hx*hz*lamb[ie]/(18*hy) - hx*hz*mu[ie]/(9*hy) - hy*hz*mu[ie]/(18*hx)
		vout[ie, 7, 2] = -hx*lamb[ie]/24 + hx*mu[ie]/24
		vout[ie, 7, 3] = hz*lamb[ie]/12 - hz*mu[ie]/12
		vout[ie, 7, 4] = hx*hy*mu[ie]/(18*hz) - hx*hz*lamb[ie]/(9*hy) - 2*hx*hz*mu[ie]/(9*hy) + hy*hz*mu[ie]/(18*hx)
		vout[ie, 7, 5] = -hx*lamb[ie]/12 + hx*mu[ie]/12
		vout[ie, 7, 6] = hz*lamb[ie]/12 + hz*mu[ie]/12
		vout[ie, 7, 7] = hx*hy*mu[ie]/(9*hz) + hx*hz*lamb[ie]/(9*hy) + 2*hx*hz*mu[ie]/(9*hy) + hy*hz*mu[ie]/(9*hx)
		vout[ie, 7, 8] = -hx*lamb[ie]/12 - hx*mu[ie]/12
		vout[ie, 7, 9] = -hz*lamb[ie]/12 + hz*mu[ie]/12
		vout[ie, 7, 10] = hx*hy*mu[ie]/(18*hz) + hx*hz*lamb[ie]/(18*hy) + hx*hz*mu[ie]/(9*hy) - hy*hz*mu[ie]/(9*hx)
		vout[ie, 7, 11] = -hx*lamb[ie]/24 - hx*mu[ie]/24
		vout[ie, 7, 12] = -hz*lamb[ie]/24 - hz*mu[ie]/24
		vout[ie, 7, 13] = -hx*hy*mu[ie]/(36*hz) - hx*hz*lamb[ie]/(36*hy) - hx*hz*mu[ie]/(18*hy) - hy*hz*mu[ie]/(36*hx)
		vout[ie, 7, 14] = hx*lamb[ie]/24 + hx*mu[ie]/24
		vout[ie, 7, 15] = hz*lamb[ie]/24 - hz*mu[ie]/24
		vout[ie, 7, 16] = -hx*hy*mu[ie]/(18*hz) - hx*hz*lamb[ie]/(18*hy) - hx*hz*mu[ie]/(9*hy) + hy*hz*mu[ie]/(36*hx)
		vout[ie, 7, 17] = hx*lamb[ie]/12 + hx*mu[ie]/12
		vout[ie, 7, 18] = hz*lamb[ie]/24 + hz*mu[ie]/24
		vout[ie, 7, 19] = -hx*hy*mu[ie]/(9*hz) + hx*hz*lamb[ie]/(18*hy) + hx*hz*mu[ie]/(9*hy) + hy*hz*mu[ie]/(18*hx)
		vout[ie, 7, 20] = hx*lamb[ie]/12 - hx*mu[ie]/12
		vout[ie, 7, 21] = -hz*lamb[ie]/24 + hz*mu[ie]/24
		vout[ie, 7, 22] = -hx*hy*mu[ie]/(18*hz) + hx*hz*lamb[ie]/(36*hy) + hx*hz*mu[ie]/(18*hy) - hy*hz*mu[ie]/(18*hx)
		vout[ie, 7, 23] = hx*lamb[ie]/24 - hx*mu[ie]/24
		vout[ie, 8, 0] = hy*lamb[ie]/24 - hy*mu[ie]/24
		vout[ie, 8, 1] = hx*lamb[ie]/24 - hx*mu[ie]/24
		vout[ie, 8, 2] = hx*hy*lamb[ie]/(36*hz) + hx*hy*mu[ie]/(18*hz) - hx*hz*mu[ie]/(18*hy) - hy*hz*mu[ie]/(18*hx)
		vout[ie, 8, 3] = -hy*lamb[ie]/24 - hy*mu[ie]/24
		vout[ie, 8, 4] = hx*lamb[ie]/12 - hx*mu[ie]/12
		vout[ie, 8, 5] = hx*hy*lamb[ie]/(18*hz) + hx*hy*mu[ie]/(9*hz) - hx*hz*mu[ie]/(9*hy) + hy*hz*mu[ie]/(18*hx)
		vout[ie, 8, 6] = -hy*lamb[ie]/12 - hy*mu[ie]/12
		vout[ie, 8, 7] = -hx*lamb[ie]/12 - hx*mu[ie]/12
		vout[ie, 8, 8] = hx*hy*lamb[ie]/(9*hz) + 2*hx*hy*mu[ie]/(9*hz) + hx*hz*mu[ie]/(9*hy) + hy*hz*mu[ie]/(9*hx)
		vout[ie, 8, 9] = hy*lamb[ie]/12 - hy*mu[ie]/12
		vout[ie, 8, 10] = -hx*lamb[ie]/24 - hx*mu[ie]/24
		vout[ie, 8, 11] = hx*hy*lamb[ie]/(18*hz) + hx*hy*mu[ie]/(9*hz) + hx*hz*mu[ie]/(18*hy) - hy*hz*mu[ie]/(9*hx)
		vout[ie, 8, 12] = hy*lamb[ie]/24 + hy*mu[ie]/24
		vout[ie, 8, 13] = hx*lamb[ie]/24 + hx*mu[ie]/24
		vout[ie, 8, 14] = -hx*hy*lamb[ie]/(36*hz) - hx*hy*mu[ie]/(18*hz) - hx*hz*mu[ie]/(36*hy) - hy*hz*mu[ie]/(36*hx)
		vout[ie, 8, 15] = -hy*lamb[ie]/24 + hy*mu[ie]/24
		vout[ie, 8, 16] = hx*lamb[ie]/12 + hx*mu[ie]/12
		vout[ie, 8, 17] = -hx*hy*lamb[ie]/(18*hz) - hx*hy*mu[ie]/(9*hz) - hx*hz*mu[ie]/(18*hy) + hy*hz*mu[ie]/(36*hx)
		vout[ie, 8, 18] = -hy*lamb[ie]/12 + hy*mu[ie]/12
		vout[ie, 8, 19] = -hx*lamb[ie]/12 + hx*mu[ie]/12
		vout[ie, 8, 20] = -hx*hy*lamb[ie]/(9*hz) - 2*hx*hy*mu[ie]/(9*hz) + hx*hz*mu[ie]/(18*hy) + hy*hz*mu[ie]/(18*hx)
		vout[ie, 8, 21] = hy*lamb[ie]/12 + hy*mu[ie]/12
		vout[ie, 8, 22] = -hx*lamb[ie]/24 + hx*mu[ie]/24
		vout[ie, 8, 23] = -hx*hy*lamb[ie]/(18*hz) - hx*hy*mu[ie]/(9*hz) + hx*hz*mu[ie]/(36*hy) - hy*hz*mu[ie]/(18*hx)
		vout[ie, 9, 0] = hx*hy*mu[ie]/(18*hz) - hx*hz*mu[ie]/(9*hy) + hy*hz*lamb[ie]/(18*hx) + hy*hz*mu[ie]/(9*hx)
		vout[ie, 9, 1] = hz*lamb[ie]/12 - hz*mu[ie]/12
		vout[ie, 9, 2] = hy*lamb[ie]/24 + hy*mu[ie]/24
		vout[ie, 9, 3] = hx*hy*mu[ie]/(36*hz) - hx*hz*mu[ie]/(18*hy) - hy*hz*lamb[ie]/(18*hx) - hy*hz*mu[ie]/(9*hx)
		vout[ie, 9, 4] = hz*lamb[ie]/12 + hz*mu[ie]/12
		vout[ie, 9, 5] = hy*lamb[ie]/24 - hy*mu[ie]/24
		vout[ie, 9, 6] = hx*hy*mu[ie]/(18*hz) + hx*hz*mu[ie]/(18*hy) - hy*hz*lamb[ie]/(9*hx) - 2*hy*hz*mu[ie]/(9*hx)
		vout[ie, 9, 7] = -hz*lamb[ie]/12 + hz*mu[ie]/12
		vout[ie, 9, 8] = hy*lamb[ie]/12 - hy*mu[ie]/12
		vout[ie, 9, 9] = hx*hy*mu[ie]/(9*hz) + hx*hz*mu[ie]/(9*hy) + hy*hz*lamb[ie]/(9*hx) + 2*hy*hz*mu[ie]/(9*hx)
		vout[ie, 9, 10] = -hz*lamb[ie]/12 - hz*mu[ie]/12
		vout[ie, 9, 11] = hy*lamb[ie]/12 + hy*mu[ie]/12
		vout[ie, 9, 12] = -hx*hy*mu[ie]/(18*hz) - hx*hz*mu[ie]/(18*hy) + hy*hz*lamb[ie]/(36*hx) + hy*hz*mu[ie]/(18*hx)
		vout[ie, 9, 13] = hz*lamb[ie]/24 - hz*mu[ie]/24
		vout[ie, 9, 14] = -hy*lamb[ie]/24 + hy*mu[ie]/24
		vout[ie, 9, 15] = -hx*hy*mu[ie]/(36*hz) - hx*hz*mu[ie]/(36*hy) - hy*hz*lamb[ie]/(36*hx) - hy*hz*mu[ie]/(18*hx)
		vout[ie, 9, 16] = hz*lamb[ie]/24 + hz*mu[ie]/24
		vout[ie, 9, 17] = -hy*lamb[ie]/24 - hy*mu[ie]/24
		vout[ie, 9, 18] = -hx*hy*mu[ie]/(18*hz) + hx*hz*mu[ie]/(36*hy) - hy*hz*lamb[ie]/(18*hx) - hy*hz*mu[ie]/(9*hx)
		vout[ie, 9, 19] = -hz*lamb[ie]/24 + hz*mu[ie]/24
		vout[ie, 9, 20] = -hy*lamb[ie]/12 - hy*mu[ie]/12
		vout[ie, 9, 21] = -hx*hy*mu[ie]/(9*hz) + hx*hz*mu[ie]/(18*hy) + hy*hz*lamb[ie]/(18*hx) + hy*hz*mu[ie]/(9*hx)
		vout[ie, 9, 22] = -hz*lamb[ie]/24 - hz*mu[ie]/24
		vout[ie, 9, 23] = -hy*lamb[ie]/12 + hy*mu[ie]/12
		vout[ie, 10, 0] = -hz*lamb[ie]/12 + hz*mu[ie]/12
		vout[ie, 10, 1] = hx*hy*mu[ie]/(18*hz) - hx*hz*lamb[ie]/(9*hy) - 2*hx*hz*mu[ie]/(9*hy) + hy*hz*mu[ie]/(18*hx)
		vout[ie, 10, 2] = -hx*lamb[ie]/12 + hx*mu[ie]/12
		vout[ie, 10, 3] = hz*lamb[ie]/12 + hz*mu[ie]/12
		vout[ie, 10, 4] = hx*hy*mu[ie]/(36*hz) - hx*hz*lamb[ie]/(18*hy) - hx*hz*mu[ie]/(9*hy) - hy*hz*mu[ie]/(18*hx)
		vout[ie, 10, 5] = -hx*lamb[ie]/24 + hx*mu[ie]/24
		vout[ie, 10, 6] = hz*lamb[ie]/12 - hz*mu[ie]/12
		vout[ie, 10, 7] = hx*hy*mu[ie]/(18*hz) + hx*hz*lamb[ie]/(18*hy) + hx*hz*mu[ie]/(9*hy) - hy*hz*mu[ie]/(9*hx)
		vout[ie, 10, 8] = -hx*lamb[ie]/24 - hx*mu[ie]/24
		vout[ie, 10, 9] = -hz*lamb[ie]/12 - hz*mu[ie]/12
		vout[ie, 10, 10] = hx*hy*mu[ie]/(9*hz) + hx*hz*lamb[ie]/(9*hy) + 2*hx*hz*mu[ie]/(9*hy) + hy*hz*mu[ie]/(9*hx)
		vout[ie, 10, 11] = -hx*lamb[ie]/12 - hx*mu[ie]/12
		vout[ie, 10, 12] = -hz*lamb[ie]/24 + hz*mu[ie]/24
		vout[ie, 10, 13] = -hx*hy*mu[ie]/(18*hz) - hx*hz*lamb[ie]/(18*hy) - hx*hz*mu[ie]/(9*hy) + hy*hz*mu[ie]/(36*hx)
		vout[ie, 10, 14] = hx*lamb[ie]/12 + hx*mu[ie]/12
		vout[ie, 10, 15] = hz*lamb[ie]/24 + hz*mu[ie]/24
		vout[ie, 10, 16] = -hx*hy*mu[ie]/(36*hz) - hx*hz*lamb[ie]/(36*hy) - hx*hz*mu[ie]/(18*hy) - hy*hz*mu[ie]/(36*hx)
		vout[ie, 10, 17] = hx*lamb[ie]/24 + hx*mu[ie]/24
		vout[ie, 10, 18] = hz*lamb[ie]/24 - hz*mu[ie]/24
		vout[ie, 10, 19] = -hx*hy*mu[ie]/(18*hz) + hx*hz*lamb[ie]/(36*hy) + hx*hz*mu[ie]/(18*hy) - hy*hz*mu[ie]/(18*hx)
		vout[ie, 10, 20] = hx*lamb[ie]/24 - hx*mu[ie]/24
		vout[ie, 10, 21] = -hz*lamb[ie]/24 - hz*mu[ie]/24
		vout[ie, 10, 22] = -hx*hy*mu[ie]/(9*hz) + hx*hz*lamb[ie]/(18*hy) + hx*hz*mu[ie]/(9*hy) + hy*hz*mu[ie]/(18*hx)
		vout[ie, 10, 23] = hx*lamb[ie]/12 - hx*mu[ie]/12
		vout[ie, 11, 0] = hy*lamb[ie]/24 + hy*mu[ie]/24
		vout[ie, 11, 1] = hx*lamb[ie]/12 - hx*mu[ie]/12
		vout[ie, 11, 2] = hx*hy*lamb[ie]/(18*hz) + hx*hy*mu[ie]/(9*hz) - hx*hz*mu[ie]/(9*hy) + hy*hz*mu[ie]/(18*hx)
		vout[ie, 11, 3] = -hy*lamb[ie]/24 + hy*mu[ie]/24
		vout[ie, 11, 4] = hx*lamb[ie]/24 - hx*mu[ie]/24
		vout[ie, 11, 5] = hx*hy*lamb[ie]/(36*hz) + hx*hy*mu[ie]/(18*hz) - hx*hz*mu[ie]/(18*hy) - hy*hz*mu[ie]/(18*hx)
		vout[ie, 11, 6] = -hy*lamb[ie]/12 + hy*mu[ie]/12
		vout[ie, 11, 7] = -hx*lamb[ie]/24 - hx*mu[ie]/24
		vout[ie, 11, 8] = hx*hy*lamb[ie]/(18*hz) + hx*hy*mu[ie]/(9*hz) + hx*hz*mu[ie]/(18*hy) - hy*hz*mu[ie]/(9*hx)
		vout[ie, 11, 9] = hy*lamb[ie]/12 + hy*mu[ie]/12
		vout[ie, 11, 10] = -hx*lamb[ie]/12 - hx*mu[ie]/12
		vout[ie, 11, 11] = hx*hy*lamb[ie]/(9*hz) + 2*hx*hy*mu[ie]/(9*hz) + hx*hz*mu[ie]/(9*hy) + hy*hz*mu[ie]/(9*hx)
		vout[ie, 11, 12] = hy*lamb[ie]/24 - hy*mu[ie]/24
		vout[ie, 11, 13] = hx*lamb[ie]/12 + hx*mu[ie]/12
		vout[ie, 11, 14] = -hx*hy*lamb[ie]/(18*hz) - hx*hy*mu[ie]/(9*hz) - hx*hz*mu[ie]/(18*hy) + hy*hz*mu[ie]/(36*hx)
		vout[ie, 11, 15] = -hy*lamb[ie]/24 - hy*mu[ie]/24
		vout[ie, 11, 16] = hx*lamb[ie]/24 + hx*mu[ie]/24
		vout[ie, 11, 17] = -hx*hy*lamb[ie]/(36*hz) - hx*hy*mu[ie]/(18*hz) - hx*hz*mu[ie]/(36*hy) - hy*hz*mu[ie]/(36*hx)
		vout[ie, 11, 18] = -hy*lamb[ie]/12 - hy*mu[ie]/12
		vout[ie, 11, 19] = -hx*lamb[ie]/24 + hx*mu[ie]/24
		vout[ie, 11, 20] = -hx*hy*lamb[ie]/(18*hz) - hx*hy*mu[ie]/(9*hz) + hx*hz*mu[ie]/(36*hy) - hy*hz*mu[ie]/(18*hx)
		vout[ie, 11, 21] = hy*lamb[ie]/12 - hy*mu[ie]/12
		vout[ie, 11, 22] = -hx*lamb[ie]/12 + hx*mu[ie]/12
		vout[ie, 11, 23] = -hx*hy*lamb[ie]/(9*hz) - 2*hx*hy*mu[ie]/(9*hz) + hx*hz*mu[ie]/(18*hy) + hy*hz*mu[ie]/(18*hx)
		vout[ie, 12, 0] = -hx*hy*mu[ie]/(9*hz) + hx*hz*mu[ie]/(18*hy) + hy*hz*lamb[ie]/(18*hx) + hy*hz*mu[ie]/(9*hx)
		vout[ie, 12, 1] = hz*lamb[ie]/24 + hz*mu[ie]/24
		vout[ie, 12, 2] = hy*lamb[ie]/12 - hy*mu[ie]/12
		vout[ie, 12, 3] = -hx*hy*mu[ie]/(18*hz) + hx*hz*mu[ie]/(36*hy) - hy*hz*lamb[ie]/(18*hx) - hy*hz*mu[ie]/(9*hx)
		vout[ie, 12, 4] = hz*lamb[ie]/24 - hz*mu[ie]/24
		vout[ie, 12, 5] = hy*lamb[ie]/12 + hy*mu[ie]/12
		vout[ie, 12, 6] = -hx*hy*mu[ie]/(36*hz) - hx*hz*mu[ie]/(36*hy) - hy*hz*lamb[ie]/(36*hx) - hy*hz*mu[ie]/(18*hx)
		vout[ie, 12, 7] = -hz*lamb[ie]/24 - hz*mu[ie]/24
		vout[ie, 12, 8] = hy*lamb[ie]/24 + hy*mu[ie]/24
		vout[ie, 12, 9] = -hx*hy*mu[ie]/(18*hz) - hx*hz*mu[ie]/(18*hy) + hy*hz*lamb[ie]/(36*hx) + hy*hz*mu[ie]/(18*hx)
		vout[ie, 12, 10] = -hz*lamb[ie]/24 + hz*mu[ie]/24
		vout[ie, 12, 11] = hy*lamb[ie]/24 - hy*mu[ie]/24
		vout[ie, 12, 12] = hx*hy*mu[ie]/(9*hz) + hx*hz*mu[ie]/(9*hy) + hy*hz*lamb[ie]/(9*hx) + 2*hy*hz*mu[ie]/(9*hx)
		vout[ie, 12, 13] = hz*lamb[ie]/12 + hz*mu[ie]/12
		vout[ie, 12, 14] = -hy*lamb[ie]/12 - hy*mu[ie]/12
		vout[ie, 12, 15] = hx*hy*mu[ie]/(18*hz) + hx*hz*mu[ie]/(18*hy) - hy*hz*lamb[ie]/(9*hx) - 2*hy*hz*mu[ie]/(9*hx)
		vout[ie, 12, 16] = hz*lamb[ie]/12 - hz*mu[ie]/12
		vout[ie, 12, 17] = -hy*lamb[ie]/12 + hy*mu[ie]/12
		vout[ie, 12, 18] = hx*hy*mu[ie]/(36*hz) - hx*hz*mu[ie]/(18*hy) - hy*hz*lamb[ie]/(18*hx) - hy*hz*mu[ie]/(9*hx)
		vout[ie, 12, 19] = -hz*lamb[ie]/12 - hz*mu[ie]/12
		vout[ie, 12, 20] = -hy*lamb[ie]/24 + hy*mu[ie]/24
		vout[ie, 12, 21] = hx*hy*mu[ie]/(18*hz) - hx*hz*mu[ie]/(9*hy) + hy*hz*lamb[ie]/(18*hx) + hy*hz*mu[ie]/(9*hx)
		vout[ie, 12, 22] = -hz*lamb[ie]/12 + hz*mu[ie]/12
		vout[ie, 12, 23] = -hy*lamb[ie]/24 - hy*mu[ie]/24
		vout[ie, 13, 0] = hz*lamb[ie]/24 + hz*mu[ie]/24
		vout[ie, 13, 1] = -hx*hy*mu[ie]/(9*hz) + hx*hz*lamb[ie]/(18*hy) + hx*hz*mu[ie]/(9*hy) + hy*hz*mu[ie]/(18*hx)
		vout[ie, 13, 2] = hx*lamb[ie]/12 - hx*mu[ie]/12
		vout[ie, 13, 3] = -hz*lamb[ie]/24 + hz*mu[ie]/24
		vout[ie, 13, 4] = -hx*hy*mu[ie]/(18*hz) + hx*hz*lamb[ie]/(36*hy) + hx*hz*mu[ie]/(18*hy) - hy*hz*mu[ie]/(18*hx)
		vout[ie, 13, 5] = hx*lamb[ie]/24 - hx*mu[ie]/24
		vout[ie, 13, 6] = -hz*lamb[ie]/24 - hz*mu[ie]/24
		vout[ie, 13, 7] = -hx*hy*mu[ie]/(36*hz) - hx*hz*lamb[ie]/(36*hy) - hx*hz*mu[ie]/(18*hy) - hy*hz*mu[ie]/(36*hx)
		vout[ie, 13, 8] = hx*lamb[ie]/24 + hx*mu[ie]/24
		vout[ie, 13, 9] = hz*lamb[ie]/24 - hz*mu[ie]/24
		vout[ie, 13, 10] = -hx*hy*mu[ie]/(18*hz) - hx*hz*lamb[ie]/(18*hy) - hx*hz*mu[ie]/(9*hy) + hy*hz*mu[ie]/(36*hx)
		vout[ie, 13, 11] = hx*lamb[ie]/12 + hx*mu[ie]/12
		vout[ie, 13, 12] = hz*lamb[ie]/12 + hz*mu[ie]/12
		vout[ie, 13, 13] = hx*hy*mu[ie]/(9*hz) + hx*hz*lamb[ie]/(9*hy) + 2*hx*hz*mu[ie]/(9*hy) + hy*hz*mu[ie]/(9*hx)
		vout[ie, 13, 14] = -hx*lamb[ie]/12 - hx*mu[ie]/12
		vout[ie, 13, 15] = -hz*lamb[ie]/12 + hz*mu[ie]/12
		vout[ie, 13, 16] = hx*hy*mu[ie]/(18*hz) + hx*hz*lamb[ie]/(18*hy) + hx*hz*mu[ie]/(9*hy) - hy*hz*mu[ie]/(9*hx)
		vout[ie, 13, 17] = -hx*lamb[ie]/24 - hx*mu[ie]/24
		vout[ie, 13, 18] = -hz*lamb[ie]/12 - hz*mu[ie]/12
		vout[ie, 13, 19] = hx*hy*mu[ie]/(36*hz) - hx*hz*lamb[ie]/(18*hy) - hx*hz*mu[ie]/(9*hy) - hy*hz*mu[ie]/(18*hx)
		vout[ie, 13, 20] = -hx*lamb[ie]/24 + hx*mu[ie]/24
		vout[ie, 13, 21] = hz*lamb[ie]/12 - hz*mu[ie]/12
		vout[ie, 13, 22] = hx*hy*mu[ie]/(18*hz) - hx*hz*lamb[ie]/(9*hy) - 2*hx*hz*mu[ie]/(9*hy) + hy*hz*mu[ie]/(18*hx)
		vout[ie, 13, 23] = -hx*lamb[ie]/12 + hx*mu[ie]/12
		vout[ie, 14, 0] = -hy*lamb[ie]/12 + hy*mu[ie]/12
		vout[ie, 14, 1] = -hx*lamb[ie]/12 + hx*mu[ie]/12
		vout[ie, 14, 2] = -hx*hy*lamb[ie]/(9*hz) - 2*hx*hy*mu[ie]/(9*hz) + hx*hz*mu[ie]/(18*hy) + hy*hz*mu[ie]/(18*hx)
		vout[ie, 14, 3] = hy*lamb[ie]/12 + hy*mu[ie]/12
		vout[ie, 14, 4] = -hx*lamb[ie]/24 + hx*mu[ie]/24
		vout[ie, 14, 5] = -hx*hy*lamb[ie]/(18*hz) - hx*hy*mu[ie]/(9*hz) + hx*hz*mu[ie]/(36*hy) - hy*hz*mu[ie]/(18*hx)
		vout[ie, 14, 6] = hy*lamb[ie]/24 + hy*mu[ie]/24
		vout[ie, 14, 7] = hx*lamb[ie]/24 + hx*mu[ie]/24
		vout[ie, 14, 8] = -hx*hy*lamb[ie]/(36*hz) - hx*hy*mu[ie]/(18*hz) - hx*hz*mu[ie]/(36*hy) - hy*hz*mu[ie]/(36*hx)
		vout[ie, 14, 9] = -hy*lamb[ie]/24 + hy*mu[ie]/24
		vout[ie, 14, 10] = hx*lamb[ie]/12 + hx*mu[ie]/12
		vout[ie, 14, 11] = -hx*hy*lamb[ie]/(18*hz) - hx*hy*mu[ie]/(9*hz) - hx*hz*mu[ie]/(18*hy) + hy*hz*mu[ie]/(36*hx)
		vout[ie, 14, 12] = -hy*lamb[ie]/12 - hy*mu[ie]/12
		vout[ie, 14, 13] = -hx*lamb[ie]/12 - hx*mu[ie]/12
		vout[ie, 14, 14] = hx*hy*lamb[ie]/(9*hz) + 2*hx*hy*mu[ie]/(9*hz) + hx*hz*mu[ie]/(9*hy) + hy*hz*mu[ie]/(9*hx)
		vout[ie, 14, 15] = hy*lamb[ie]/12 - hy*mu[ie]/12
		vout[ie, 14, 16] = -hx*lamb[ie]/24 - hx*mu[ie]/24
		vout[ie, 14, 17] = hx*hy*lamb[ie]/(18*hz) + hx*hy*mu[ie]/(9*hz) + hx*hz*mu[ie]/(18*hy) - hy*hz*mu[ie]/(9*hx)
		vout[ie, 14, 18] = hy*lamb[ie]/24 - hy*mu[ie]/24
		vout[ie, 14, 19] = hx*lamb[ie]/24 - hx*mu[ie]/24
		vout[ie, 14, 20] = hx*hy*lamb[ie]/(36*hz) + hx*hy*mu[ie]/(18*hz) - hx*hz*mu[ie]/(18*hy) - hy*hz*mu[ie]/(18*hx)
		vout[ie, 14, 21] = -hy*lamb[ie]/24 - hy*mu[ie]/24
		vout[ie, 14, 22] = hx*lamb[ie]/12 - hx*mu[ie]/12
		vout[ie, 14, 23] = hx*hy*lamb[ie]/(18*hz) + hx*hy*mu[ie]/(9*hz) - hx*hz*mu[ie]/(9*hy) + hy*hz*mu[ie]/(18*hx)
		vout[ie, 15, 0] = -hx*hy*mu[ie]/(18*hz) + hx*hz*mu[ie]/(36*hy) - hy*hz*lamb[ie]/(18*hx) - hy*hz*mu[ie]/(9*hx)
		vout[ie, 15, 1] = -hz*lamb[ie]/24 + hz*mu[ie]/24
		vout[ie, 15, 2] = -hy*lamb[ie]/12 - hy*mu[ie]/12
		vout[ie, 15, 3] = -hx*hy*mu[ie]/(9*hz) + hx*hz*mu[ie]/(18*hy) + hy*hz*lamb[ie]/(18*hx) + hy*hz*mu[ie]/(9*hx)
		vout[ie, 15, 4] = -hz*lamb[ie]/24 - hz*mu[ie]/24
		vout[ie, 15, 5] = -hy*lamb[ie]/12 + hy*mu[ie]/12
		vout[ie, 15, 6] = -hx*hy*mu[ie]/(18*hz) - hx*hz*mu[ie]/(18*hy) + hy*hz*lamb[ie]/(36*hx) + hy*hz*mu[ie]/(18*hx)
		vout[ie, 15, 7] = hz*lamb[ie]/24 - hz*mu[ie]/24
		vout[ie, 15, 8] = -hy*lamb[ie]/24 + hy*mu[ie]/24
		vout[ie, 15, 9] = -hx*hy*mu[ie]/(36*hz) - hx*hz*mu[ie]/(36*hy) - hy*hz*lamb[ie]/(36*hx) - hy*hz*mu[ie]/(18*hx)
		vout[ie, 15, 10] = hz*lamb[ie]/24 + hz*mu[ie]/24
		vout[ie, 15, 11] = -hy*lamb[ie]/24 - hy*mu[ie]/24
		vout[ie, 15, 12] = hx*hy*mu[ie]/(18*hz) + hx*hz*mu[ie]/(18*hy) - hy*hz*lamb[ie]/(9*hx) - 2*hy*hz*mu[ie]/(9*hx)
		vout[ie, 15, 13] = -hz*lamb[ie]/12 + hz*mu[ie]/12
		vout[ie, 15, 14] = hy*lamb[ie]/12 - hy*mu[ie]/12
		vout[ie, 15, 15] = hx*hy*mu[ie]/(9*hz) + hx*hz*mu[ie]/(9*hy) + hy*hz*lamb[ie]/(9*hx) + 2*hy*hz*mu[ie]/(9*hx)
		vout[ie, 15, 16] = -hz*lamb[ie]/12 - hz*mu[ie]/12
		vout[ie, 15, 17] = hy*lamb[ie]/12 + hy*mu[ie]/12
		vout[ie, 15, 18] = hx*hy*mu[ie]/(18*hz) - hx*hz*mu[ie]/(9*hy) + hy*hz*lamb[ie]/(18*hx) + hy*hz*mu[ie]/(9*hx)
		vout[ie, 15, 19] = hz*lamb[ie]/12 - hz*mu[ie]/12
		vout[ie, 15, 20] = hy*lamb[ie]/24 + hy*mu[ie]/24
		vout[ie, 15, 21] = hx*hy*mu[ie]/(36*hz) - hx*hz*mu[ie]/(18*hy) - hy*hz*lamb[ie]/(18*hx) - hy*hz*mu[ie]/(9*hx)
		vout[ie, 15, 22] = hz*lamb[ie]/12 + hz*mu[ie]/12
		vout[ie, 15, 23] = hy*lamb[ie]/24 - hy*mu[ie]/24
		vout[ie, 16, 0] = hz*lamb[ie]/24 - hz*mu[ie]/24
		vout[ie, 16, 1] = -hx*hy*mu[ie]/(18*hz) + hx*hz*lamb[ie]/(36*hy) + hx*hz*mu[ie]/(18*hy) - hy*hz*mu[ie]/(18*hx)
		vout[ie, 16, 2] = hx*lamb[ie]/24 - hx*mu[ie]/24
		vout[ie, 16, 3] = -hz*lamb[ie]/24 - hz*mu[ie]/24
		vout[ie, 16, 4] = -hx*hy*mu[ie]/(9*hz) + hx*hz*lamb[ie]/(18*hy) + hx*hz*mu[ie]/(9*hy) + hy*hz*mu[ie]/(18*hx)
		vout[ie, 16, 5] = hx*lamb[ie]/12 - hx*mu[ie]/12
		vout[ie, 16, 6] = -hz*lamb[ie]/24 + hz*mu[ie]/24
		vout[ie, 16, 7] = -hx*hy*mu[ie]/(18*hz) - hx*hz*lamb[ie]/(18*hy) - hx*hz*mu[ie]/(9*hy) + hy*hz*mu[ie]/(36*hx)
		vout[ie, 16, 8] = hx*lamb[ie]/12 + hx*mu[ie]/12
		vout[ie, 16, 9] = hz*lamb[ie]/24 + hz*mu[ie]/24
		vout[ie, 16, 10] = -hx*hy*mu[ie]/(36*hz) - hx*hz*lamb[ie]/(36*hy) - hx*hz*mu[ie]/(18*hy) - hy*hz*mu[ie]/(36*hx)
		vout[ie, 16, 11] = hx*lamb[ie]/24 + hx*mu[ie]/24
		vout[ie, 16, 12] = hz*lamb[ie]/12 - hz*mu[ie]/12
		vout[ie, 16, 13] = hx*hy*mu[ie]/(18*hz) + hx*hz*lamb[ie]/(18*hy) + hx*hz*mu[ie]/(9*hy) - hy*hz*mu[ie]/(9*hx)
		vout[ie, 16, 14] = -hx*lamb[ie]/24 - hx*mu[ie]/24
		vout[ie, 16, 15] = -hz*lamb[ie]/12 - hz*mu[ie]/12
		vout[ie, 16, 16] = hx*hy*mu[ie]/(9*hz) + hx*hz*lamb[ie]/(9*hy) + 2*hx*hz*mu[ie]/(9*hy) + hy*hz*mu[ie]/(9*hx)
		vout[ie, 16, 17] = -hx*lamb[ie]/12 - hx*mu[ie]/12
		vout[ie, 16, 18] = -hz*lamb[ie]/12 + hz*mu[ie]/12
		vout[ie, 16, 19] = hx*hy*mu[ie]/(18*hz) - hx*hz*lamb[ie]/(9*hy) - 2*hx*hz*mu[ie]/(9*hy) + hy*hz*mu[ie]/(18*hx)
		vout[ie, 16, 20] = -hx*lamb[ie]/12 + hx*mu[ie]/12
		vout[ie, 16, 21] = hz*lamb[ie]/12 + hz*mu[ie]/12
		vout[ie, 16, 22] = hx*hy*mu[ie]/(36*hz) - hx*hz*lamb[ie]/(18*hy) - hx*hz*mu[ie]/(9*hy) - hy*hz*mu[ie]/(18*hx)
		vout[ie, 16, 23] = -hx*lamb[ie]/24 + hx*mu[ie]/24
		vout[ie, 17, 0] = -hy*lamb[ie]/12 - hy*mu[ie]/12
		vout[ie, 17, 1] = -hx*lamb[ie]/24 + hx*mu[ie]/24
		vout[ie, 17, 2] = -hx*hy*lamb[ie]/(18*hz) - hx*hy*mu[ie]/(9*hz) + hx*hz*mu[ie]/(36*hy) - hy*hz*mu[ie]/(18*hx)
		vout[ie, 17, 3] = hy*lamb[ie]/12 - hy*mu[ie]/12
		vout[ie, 17, 4] = -hx*lamb[ie]/12 + hx*mu[ie]/12
		vout[ie, 17, 5] = -hx*hy*lamb[ie]/(9*hz) - 2*hx*hy*mu[ie]/(9*hz) + hx*hz*mu[ie]/(18*hy) + hy*hz*mu[ie]/(18*hx)
		vout[ie, 17, 6] = hy*lamb[ie]/24 - hy*mu[ie]/24
		vout[ie, 17, 7] = hx*lamb[ie]/12 + hx*mu[ie]/12
		vout[ie, 17, 8] = -hx*hy*lamb[ie]/(18*hz) - hx*hy*mu[ie]/(9*hz) - hx*hz*mu[ie]/(18*hy) + hy*hz*mu[ie]/(36*hx)
		vout[ie, 17, 9] = -hy*lamb[ie]/24 - hy*mu[ie]/24
		vout[ie, 17, 10] = hx*lamb[ie]/24 + hx*mu[ie]/24
		vout[ie, 17, 11] = -hx*hy*lamb[ie]/(36*hz) - hx*hy*mu[ie]/(18*hz) - hx*hz*mu[ie]/(36*hy) - hy*hz*mu[ie]/(36*hx)
		vout[ie, 17, 12] = -hy*lamb[ie]/12 + hy*mu[ie]/12
		vout[ie, 17, 13] = -hx*lamb[ie]/24 - hx*mu[ie]/24
		vout[ie, 17, 14] = hx*hy*lamb[ie]/(18*hz) + hx*hy*mu[ie]/(9*hz) + hx*hz*mu[ie]/(18*hy) - hy*hz*mu[ie]/(9*hx)
		vout[ie, 17, 15] = hy*lamb[ie]/12 + hy*mu[ie]/12
		vout[ie, 17, 16] = -hx*lamb[ie]/12 - hx*mu[ie]/12
		vout[ie, 17, 17] = hx*hy*lamb[ie]/(9*hz) + 2*hx*hy*mu[ie]/(9*hz) + hx*hz*mu[ie]/(9*hy) + hy*hz*mu[ie]/(9*hx)
		vout[ie, 17, 18] = hy*lamb[ie]/24 + hy*mu[ie]/24
		vout[ie, 17, 19] = hx*lamb[ie]/12 - hx*mu[ie]/12
		vout[ie, 17, 20] = hx*hy*lamb[ie]/(18*hz) + hx*hy*mu[ie]/(9*hz) - hx*hz*mu[ie]/(9*hy) + hy*hz*mu[ie]/(18*hx)
		vout[ie, 17, 21] = -hy*lamb[ie]/24 + hy*mu[ie]/24
		vout[ie, 17, 22] = hx*lamb[ie]/24 - hx*mu[ie]/24
		vout[ie, 17, 23] = hx*hy*lamb[ie]/(36*hz) + hx*hy*mu[ie]/(18*hz) - hx*hz*mu[ie]/(18*hy) - hy*hz*mu[ie]/(18*hx)
		vout[ie, 18, 0] = -hx*hy*mu[ie]/(36*hz) - hx*hz*mu[ie]/(36*hy) - hy*hz*lamb[ie]/(36*hx) - hy*hz*mu[ie]/(18*hx)
		vout[ie, 18, 1] = -hz*lamb[ie]/24 - hz*mu[ie]/24
		vout[ie, 18, 2] = -hy*lamb[ie]/24 - hy*mu[ie]/24
		vout[ie, 18, 3] = -hx*hy*mu[ie]/(18*hz) - hx*hz*mu[ie]/(18*hy) + hy*hz*lamb[ie]/(36*hx) + hy*hz*mu[ie]/(18*hx)
		vout[ie, 18, 4] = -hz*lamb[ie]/24 + hz*mu[ie]/24
		vout[ie, 18, 5] = -hy*lamb[ie]/24 + hy*mu[ie]/24
		vout[ie, 18, 6] = -hx*hy*mu[ie]/(9*hz) + hx*hz*mu[ie]/(18*hy) + hy*hz*lamb[ie]/(18*hx) + hy*hz*mu[ie]/(9*hx)
		vout[ie, 18, 7] = hz*lamb[ie]/24 + hz*mu[ie]/24
		vout[ie, 18, 8] = -hy*lamb[ie]/12 + hy*mu[ie]/12
		vout[ie, 18, 9] = -hx*hy*mu[ie]/(18*hz) + hx*hz*mu[ie]/(36*hy) - hy*hz*lamb[ie]/(18*hx) - hy*hz*mu[ie]/(9*hx)
		vout[ie, 18, 10] = hz*lamb[ie]/24 - hz*mu[ie]/24
		vout[ie, 18, 11] = -hy*lamb[ie]/12 - hy*mu[ie]/12
		vout[ie, 18, 12] = hx*hy*mu[ie]/(36*hz) - hx*hz*mu[ie]/(18*hy) - hy*hz*lamb[ie]/(18*hx) - hy*hz*mu[ie]/(9*hx)
		vout[ie, 18, 13] = -hz*lamb[ie]/12 - hz*mu[ie]/12
		vout[ie, 18, 14] = hy*lamb[ie]/24 - hy*mu[ie]/24
		vout[ie, 18, 15] = hx*hy*mu[ie]/(18*hz) - hx*hz*mu[ie]/(9*hy) + hy*hz*lamb[ie]/(18*hx) + hy*hz*mu[ie]/(9*hx)
		vout[ie, 18, 16] = -hz*lamb[ie]/12 + hz*mu[ie]/12
		vout[ie, 18, 17] = hy*lamb[ie]/24 + hy*mu[ie]/24
		vout[ie, 18, 18] = hx*hy*mu[ie]/(9*hz) + hx*hz*mu[ie]/(9*hy) + hy*hz*lamb[ie]/(9*hx) + 2*hy*hz*mu[ie]/(9*hx)
		vout[ie, 18, 19] = hz*lamb[ie]/12 + hz*mu[ie]/12
		vout[ie, 18, 20] = hy*lamb[ie]/12 + hy*mu[ie]/12
		vout[ie, 18, 21] = hx*hy*mu[ie]/(18*hz) + hx*hz*mu[ie]/(18*hy) - hy*hz*lamb[ie]/(9*hx) - 2*hy*hz*mu[ie]/(9*hx)
		vout[ie, 18, 22] = hz*lamb[ie]/12 - hz*mu[ie]/12
		vout[ie, 18, 23] = hy*lamb[ie]/12 - hy*mu[ie]/12
		vout[ie, 19, 0] = -hz*lamb[ie]/24 - hz*mu[ie]/24
		vout[ie, 19, 1] = -hx*hy*mu[ie]/(36*hz) - hx*hz*lamb[ie]/(36*hy) - hx*hz*mu[ie]/(18*hy) - hy*hz*mu[ie]/(36*hx)
		vout[ie, 19, 2] = -hx*lamb[ie]/24 - hx*mu[ie]/24
		vout[ie, 19, 3] = hz*lamb[ie]/24 - hz*mu[ie]/24
		vout[ie, 19, 4] = -hx*hy*mu[ie]/(18*hz) - hx*hz*lamb[ie]/(18*hy) - hx*hz*mu[ie]/(9*hy) + hy*hz*mu[ie]/(36*hx)
		vout[ie, 19, 5] = -hx*lamb[ie]/12 - hx*mu[ie]/12
		vout[ie, 19, 6] = hz*lamb[ie]/24 + hz*mu[ie]/24
		vout[ie, 19, 7] = -hx*hy*mu[ie]/(9*hz) + hx*hz*lamb[ie]/(18*hy) + hx*hz*mu[ie]/(9*hy) + hy*hz*mu[ie]/(18*hx)
		vout[ie, 19, 8] = -hx*lamb[ie]/12 + hx*mu[ie]/12
		vout[ie, 19, 9] = -hz*lamb[ie]/24 + hz*mu[ie]/24
		vout[ie, 19, 10] = -hx*hy*mu[ie]/(18*hz) + hx*hz*lamb[ie]/(36*hy) + hx*hz*mu[ie]/(18*hy) - hy*hz*mu[ie]/(18*hx)
		vout[ie, 19, 11] = -hx*lamb[ie]/24 + hx*mu[ie]/24
		vout[ie, 19, 12] = -hz*lamb[ie]/12 - hz*mu[ie]/12
		vout[ie, 19, 13] = hx*hy*mu[ie]/(36*hz) - hx*hz*lamb[ie]/(18*hy) - hx*hz*mu[ie]/(9*hy) - hy*hz*mu[ie]/(18*hx)
		vout[ie, 19, 14] = hx*lamb[ie]/24 - hx*mu[ie]/24
		vout[ie, 19, 15] = hz*lamb[ie]/12 - hz*mu[ie]/12
		vout[ie, 19, 16] = hx*hy*mu[ie]/(18*hz) - hx*hz*lamb[ie]/(9*hy) - 2*hx*hz*mu[ie]/(9*hy) + hy*hz*mu[ie]/(18*hx)
		vout[ie, 19, 17] = hx*lamb[ie]/12 - hx*mu[ie]/12
		vout[ie, 19, 18] = hz*lamb[ie]/12 + hz*mu[ie]/12
		vout[ie, 19, 19] = hx*hy*mu[ie]/(9*hz) + hx*hz*lamb[ie]/(9*hy) + 2*hx*hz*mu[ie]/(9*hy) + hy*hz*mu[ie]/(9*hx)
		vout[ie, 19, 20] = hx*lamb[ie]/12 + hx*mu[ie]/12
		vout[ie, 19, 21] = -hz*lamb[ie]/12 + hz*mu[ie]/12
		vout[ie, 19, 22] = hx*hy*mu[ie]/(18*hz) + hx*hz*lamb[ie]/(18*hy) + hx*hz*mu[ie]/(9*hy) - hy*hz*mu[ie]/(9*hx)
		vout[ie, 19, 23] = hx*lamb[ie]/24 + hx*mu[ie]/24
		vout[ie, 20, 0] = -hy*lamb[ie]/24 - hy*mu[ie]/24
		vout[ie, 20, 1] = -hx*lamb[ie]/24 - hx*mu[ie]/24
		vout[ie, 20, 2] = -hx*hy*lamb[ie]/(36*hz) - hx*hy*mu[ie]/(18*hz) - hx*hz*mu[ie]/(36*hy) - hy*hz*mu[ie]/(36*hx)
		vout[ie, 20, 3] = hy*lamb[ie]/24 - hy*mu[ie]/24
		vout[ie, 20, 4] = -hx*lamb[ie]/12 - hx*mu[ie]/12
		vout[ie, 20, 5] = -hx*hy*lamb[ie]/(18*hz) - hx*hy*mu[ie]/(9*hz) - hx*hz*mu[ie]/(18*hy) + hy*hz*mu[ie]/(36*hx)
		vout[ie, 20, 6] = hy*lamb[ie]/12 - hy*mu[ie]/12
		vout[ie, 20, 7] = hx*lamb[ie]/12 - hx*mu[ie]/12
		vout[ie, 20, 8] = -hx*hy*lamb[ie]/(9*hz) - 2*hx*hy*mu[ie]/(9*hz) + hx*hz*mu[ie]/(18*hy) + hy*hz*mu[ie]/(18*hx)
		vout[ie, 20, 9] = -hy*lamb[ie]/12 - hy*mu[ie]/12
		vout[ie, 20, 10] = hx*lamb[ie]/24 - hx*mu[ie]/24
		vout[ie, 20, 11] = -hx*hy*lamb[ie]/(18*hz) - hx*hy*mu[ie]/(9*hz) + hx*hz*mu[ie]/(36*hy) - hy*hz*mu[ie]/(18*hx)
		vout[ie, 20, 12] = -hy*lamb[ie]/24 + hy*mu[ie]/24
		vout[ie, 20, 13] = -hx*lamb[ie]/24 + hx*mu[ie]/24
		vout[ie, 20, 14] = hx*hy*lamb[ie]/(36*hz) + hx*hy*mu[ie]/(18*hz) - hx*hz*mu[ie]/(18*hy) - hy*hz*mu[ie]/(18*hx)
		vout[ie, 20, 15] = hy*lamb[ie]/24 + hy*mu[ie]/24
		vout[ie, 20, 16] = -hx*lamb[ie]/12 + hx*mu[ie]/12
		vout[ie, 20, 17] = hx*hy*lamb[ie]/(18*hz) + hx*hy*mu[ie]/(9*hz) - hx*hz*mu[ie]/(9*hy) + hy*hz*mu[ie]/(18*hx)
		vout[ie, 20, 18] = hy*lamb[ie]/12 + hy*mu[ie]/12
		vout[ie, 20, 19] = hx*lamb[ie]/12 + hx*mu[ie]/12
		vout[ie, 20, 20] = hx*hy*lamb[ie]/(9*hz) + 2*hx*hy*mu[ie]/(9*hz) + hx*hz*mu[ie]/(9*hy) + hy*hz*mu[ie]/(9*hx)
		vout[ie, 20, 21] = -hy*lamb[ie]/12 + hy*mu[ie]/12
		vout[ie, 20, 22] = hx*lamb[ie]/24 + hx*mu[ie]/24
		vout[ie, 20, 23] = hx*hy*lamb[ie]/(18*hz) + hx*hy*mu[ie]/(9*hz) + hx*hz*mu[ie]/(18*hy) - hy*hz*mu[ie]/(9*hx)
		vout[ie, 21, 0] = -hx*hy*mu[ie]/(18*hz) - hx*hz*mu[ie]/(18*hy) + hy*hz*lamb[ie]/(36*hx) + hy*hz*mu[ie]/(18*hx)
		vout[ie, 21, 1] = hz*lamb[ie]/24 - hz*mu[ie]/24
		vout[ie, 21, 2] = hy*lamb[ie]/24 - hy*mu[ie]/24
		vout[ie, 21, 3] = -hx*hy*mu[ie]/(36*hz) - hx*hz*mu[ie]/(36*hy) - hy*hz*lamb[ie]/(36*hx) - hy*hz*mu[ie]/(18*hx)
		vout[ie, 21, 4] = hz*lamb[ie]/24 + hz*mu[ie]/24
		vout[ie, 21, 5] = hy*lamb[ie]/24 + hy*mu[ie]/24
		vout[ie, 21, 6] = -hx*hy*mu[ie]/(18*hz) + hx*hz*mu[ie]/(36*hy) - hy*hz*lamb[ie]/(18*hx) - hy*hz*mu[ie]/(9*hx)
		vout[ie, 21, 7] = -hz*lamb[ie]/24 + hz*mu[ie]/24
		vout[ie, 21, 8] = hy*lamb[ie]/12 + hy*mu[ie]/12
		vout[ie, 21, 9] = -hx*hy*mu[ie]/(9*hz) + hx*hz*mu[ie]/(18*hy) + hy*hz*lamb[ie]/(18*hx) + hy*hz*mu[ie]/(9*hx)
		vout[ie, 21, 10] = -hz*lamb[ie]/24 - hz*mu[ie]/24
		vout[ie, 21, 11] = hy*lamb[ie]/12 - hy*mu[ie]/12
		vout[ie, 21, 12] = hx*hy*mu[ie]/(18*hz) - hx*hz*mu[ie]/(9*hy) + hy*hz*lamb[ie]/(18*hx) + hy*hz*mu[ie]/(9*hx)
		vout[ie, 21, 13] = hz*lamb[ie]/12 - hz*mu[ie]/12
		vout[ie, 21, 14] = -hy*lamb[ie]/24 - hy*mu[ie]/24
		vout[ie, 21, 15] = hx*hy*mu[ie]/(36*hz) - hx*hz*mu[ie]/(18*hy) - hy*hz*lamb[ie]/(18*hx) - hy*hz*mu[ie]/(9*hx)
		vout[ie, 21, 16] = hz*lamb[ie]/12 + hz*mu[ie]/12
		vout[ie, 21, 17] = -hy*lamb[ie]/24 + hy*mu[ie]/24
		vout[ie, 21, 18] = hx*hy*mu[ie]/(18*hz) + hx*hz*mu[ie]/(18*hy) - hy*hz*lamb[ie]/(9*hx) - 2*hy*hz*mu[ie]/(9*hx)
		vout[ie, 21, 19] = -hz*lamb[ie]/12 + hz*mu[ie]/12
		vout[ie, 21, 20] = -hy*lamb[ie]/12 + hy*mu[ie]/12
		vout[ie, 21, 21] = hx*hy*mu[ie]/(9*hz) + hx*hz*mu[ie]/(9*hy) + hy*hz*lamb[ie]/(9*hx) + 2*hy*hz*mu[ie]/(9*hx)
		vout[ie, 21, 22] = -hz*lamb[ie]/12 - hz*mu[ie]/12
		vout[ie, 21, 23] = -hy*lamb[ie]/12 - hy*mu[ie]/12
		vout[ie, 22, 0] = -hz*lamb[ie]/24 + hz*mu[ie]/24
		vout[ie, 22, 1] = -hx*hy*mu[ie]/(18*hz) - hx*hz*lamb[ie]/(18*hy) - hx*hz*mu[ie]/(9*hy) + hy*hz*mu[ie]/(36*hx)
		vout[ie, 22, 2] = -hx*lamb[ie]/12 - hx*mu[ie]/12
		vout[ie, 22, 3] = hz*lamb[ie]/24 + hz*mu[ie]/24
		vout[ie, 22, 4] = -hx*hy*mu[ie]/(36*hz) - hx*hz*lamb[ie]/(36*hy) - hx*hz*mu[ie]/(18*hy) - hy*hz*mu[ie]/(36*hx)
		vout[ie, 22, 5] = -hx*lamb[ie]/24 - hx*mu[ie]/24
		vout[ie, 22, 6] = hz*lamb[ie]/24 - hz*mu[ie]/24
		vout[ie, 22, 7] = -hx*hy*mu[ie]/(18*hz) + hx*hz*lamb[ie]/(36*hy) + hx*hz*mu[ie]/(18*hy) - hy*hz*mu[ie]/(18*hx)
		vout[ie, 22, 8] = -hx*lamb[ie]/24 + hx*mu[ie]/24
		vout[ie, 22, 9] = -hz*lamb[ie]/24 - hz*mu[ie]/24
		vout[ie, 22, 10] = -hx*hy*mu[ie]/(9*hz) + hx*hz*lamb[ie]/(18*hy) + hx*hz*mu[ie]/(9*hy) + hy*hz*mu[ie]/(18*hx)
		vout[ie, 22, 11] = -hx*lamb[ie]/12 + hx*mu[ie]/12
		vout[ie, 22, 12] = -hz*lamb[ie]/12 + hz*mu[ie]/12
		vout[ie, 22, 13] = hx*hy*mu[ie]/(18*hz) - hx*hz*lamb[ie]/(9*hy) - 2*hx*hz*mu[ie]/(9*hy) + hy*hz*mu[ie]/(18*hx)
		vout[ie, 22, 14] = hx*lamb[ie]/12 - hx*mu[ie]/12
		vout[ie, 22, 15] = hz*lamb[ie]/12 + hz*mu[ie]/12
		vout[ie, 22, 16] = hx*hy*mu[ie]/(36*hz) - hx*hz*lamb[ie]/(18*hy) - hx*hz*mu[ie]/(9*hy) - hy*hz*mu[ie]/(18*hx)
		vout[ie, 22, 17] = hx*lamb[ie]/24 - hx*mu[ie]/24
		vout[ie, 22, 18] = hz*lamb[ie]/12 - hz*mu[ie]/12
		vout[ie, 22, 19] = hx*hy*mu[ie]/(18*hz) + hx*hz*lamb[ie]/(18*hy) + hx*hz*mu[ie]/(9*hy) - hy*hz*mu[ie]/(9*hx)
		vout[ie, 22, 20] = hx*lamb[ie]/24 + hx*mu[ie]/24
		vout[ie, 22, 21] = -hz*lamb[ie]/12 - hz*mu[ie]/12
		vout[ie, 22, 22] = hx*hy*mu[ie]/(9*hz) + hx*hz*lamb[ie]/(9*hy) + 2*hx*hz*mu[ie]/(9*hy) + hy*hz*mu[ie]/(9*hx)
		vout[ie, 22, 23] = hx*lamb[ie]/12 + hx*mu[ie]/12
		vout[ie, 23, 0] = -hy*lamb[ie]/24 + hy*mu[ie]/24
		vout[ie, 23, 1] = -hx*lamb[ie]/12 - hx*mu[ie]/12
		vout[ie, 23, 2] = -hx*hy*lamb[ie]/(18*hz) - hx*hy*mu[ie]/(9*hz) - hx*hz*mu[ie]/(18*hy) + hy*hz*mu[ie]/(36*hx)
		vout[ie, 23, 3] = hy*lamb[ie]/24 + hy*mu[ie]/24
		vout[ie, 23, 4] = -hx*lamb[ie]/24 - hx*mu[ie]/24
		vout[ie, 23, 5] = -hx*hy*lamb[ie]/(36*hz) - hx*hy*mu[ie]/(18*hz) - hx*hz*mu[ie]/(36*hy) - hy*hz*mu[ie]/(36*hx)
		vout[ie, 23, 6] = hy*lamb[ie]/12 + hy*mu[ie]/12
		vout[ie, 23, 7] = hx*lamb[ie]/24 - hx*mu[ie]/24
		vout[ie, 23, 8] = -hx*hy*lamb[ie]/(18*hz) - hx*hy*mu[ie]/(9*hz) + hx*hz*mu[ie]/(36*hy) - hy*hz*mu[ie]/(18*hx)
		vout[ie, 23, 9] = -hy*lamb[ie]/12 + hy*mu[ie]/12
		vout[ie, 23, 10] = hx*lamb[ie]/12 - hx*mu[ie]/12
		vout[ie, 23, 11] = -hx*hy*lamb[ie]/(9*hz) - 2*hx*hy*mu[ie]/(9*hz) + hx*hz*mu[ie]/(18*hy) + hy*hz*mu[ie]/(18*hx)
		vout[ie, 23, 12] = -hy*lamb[ie]/24 - hy*mu[ie]/24
		vout[ie, 23, 13] = -hx*lamb[ie]/12 + hx*mu[ie]/12
		vout[ie, 23, 14] = hx*hy*lamb[ie]/(18*hz) + hx*hy*mu[ie]/(9*hz) - hx*hz*mu[ie]/(9*hy) + hy*hz*mu[ie]/(18*hx)
		vout[ie, 23, 15] = hy*lamb[ie]/24 - hy*mu[ie]/24
		vout[ie, 23, 16] = -hx*lamb[ie]/24 + hx*mu[ie]/24
		vout[ie, 23, 17] = hx*hy*lamb[ie]/(36*hz) + hx*hy*mu[ie]/(18*hz) - hx*hz*mu[ie]/(18*hy) - hy*hz*mu[ie]/(18*hx)
		vout[ie, 23, 18] = hy*lamb[ie]/12 - hy*mu[ie]/12
		vout[ie, 23, 19] = hx*lamb[ie]/24 + hx*mu[ie]/24
		vout[ie, 23, 20] = hx*hy*lamb[ie]/(18*hz) + hx*hy*mu[ie]/(9*hz) + hx*hz*mu[ie]/(18*hy) - hy*hz*mu[ie]/(9*hx)
		vout[ie, 23, 21] = -hy*lamb[ie]/12 - hy*mu[ie]/12
		vout[ie, 23, 22] = hx*lamb[ie]/12 + hx*mu[ie]/12
		vout[ie, 23, 23] = hx*hy*lamb[ie]/(9*hz) + 2*hx*hy*mu[ie]/(9*hz) + hx*hz*mu[ie]/(9*hy) + hy*hz*mu[ie]/(9*hx)
	return output

def getMatElemElasticity_3d_const(double hx, double hy, double hz, double lamb, double mu):
	output = np.empty((24, 24))
	cdef:
		double[:, ::1] vout = output
	vout[0, 0] = hx*hy*mu/(9*hz) + hx*hz*mu/(9*hy) + hy*hz*lamb/(9*hx) + 2*hy*hz*mu/(9*hx)
	vout[0, 1] = hz*lamb/12 + hz*mu/12
	vout[0, 2] = hy*lamb/12 + hy*mu/12
	vout[0, 3] = hx*hy*mu/(18*hz) + hx*hz*mu/(18*hy) - hy*hz*lamb/(9*hx) - 2*hy*hz*mu/(9*hx)
	vout[0, 4] = hz*lamb/12 - hz*mu/12
	vout[0, 5] = hy*lamb/12 - hy*mu/12
	vout[0, 6] = hx*hy*mu/(36*hz) - hx*hz*mu/(18*hy) - hy*hz*lamb/(18*hx) - hy*hz*mu/(9*hx)
	vout[0, 7] = -hz*lamb/12 - hz*mu/12
	vout[0, 8] = hy*lamb/24 - hy*mu/24
	vout[0, 9] = hx*hy*mu/(18*hz) - hx*hz*mu/(9*hy) + hy*hz*lamb/(18*hx) + hy*hz*mu/(9*hx)
	vout[0, 10] = -hz*lamb/12 + hz*mu/12
	vout[0, 11] = hy*lamb/24 + hy*mu/24
	vout[0, 12] = -hx*hy*mu/(9*hz) + hx*hz*mu/(18*hy) + hy*hz*lamb/(18*hx) + hy*hz*mu/(9*hx)
	vout[0, 13] = hz*lamb/24 + hz*mu/24
	vout[0, 14] = -hy*lamb/12 + hy*mu/12
	vout[0, 15] = -hx*hy*mu/(18*hz) + hx*hz*mu/(36*hy) - hy*hz*lamb/(18*hx) - hy*hz*mu/(9*hx)
	vout[0, 16] = hz*lamb/24 - hz*mu/24
	vout[0, 17] = -hy*lamb/12 - hy*mu/12
	vout[0, 18] = -hx*hy*mu/(36*hz) - hx*hz*mu/(36*hy) - hy*hz*lamb/(36*hx) - hy*hz*mu/(18*hx)
	vout[0, 19] = -hz*lamb/24 - hz*mu/24
	vout[0, 20] = -hy*lamb/24 - hy*mu/24
	vout[0, 21] = -hx*hy*mu/(18*hz) - hx*hz*mu/(18*hy) + hy*hz*lamb/(36*hx) + hy*hz*mu/(18*hx)
	vout[0, 22] = -hz*lamb/24 + hz*mu/24
	vout[0, 23] = -hy*lamb/24 + hy*mu/24
	vout[1, 0] = hz*lamb/12 + hz*mu/12
	vout[1, 1] = hx*hy*mu/(9*hz) + hx*hz*lamb/(9*hy) + 2*hx*hz*mu/(9*hy) + hy*hz*mu/(9*hx)
	vout[1, 2] = hx*lamb/12 + hx*mu/12
	vout[1, 3] = -hz*lamb/12 + hz*mu/12
	vout[1, 4] = hx*hy*mu/(18*hz) + hx*hz*lamb/(18*hy) + hx*hz*mu/(9*hy) - hy*hz*mu/(9*hx)
	vout[1, 5] = hx*lamb/24 + hx*mu/24
	vout[1, 6] = -hz*lamb/12 - hz*mu/12
	vout[1, 7] = hx*hy*mu/(36*hz) - hx*hz*lamb/(18*hy) - hx*hz*mu/(9*hy) - hy*hz*mu/(18*hx)
	vout[1, 8] = hx*lamb/24 - hx*mu/24
	vout[1, 9] = hz*lamb/12 - hz*mu/12
	vout[1, 10] = hx*hy*mu/(18*hz) - hx*hz*lamb/(9*hy) - 2*hx*hz*mu/(9*hy) + hy*hz*mu/(18*hx)
	vout[1, 11] = hx*lamb/12 - hx*mu/12
	vout[1, 12] = hz*lamb/24 + hz*mu/24
	vout[1, 13] = -hx*hy*mu/(9*hz) + hx*hz*lamb/(18*hy) + hx*hz*mu/(9*hy) + hy*hz*mu/(18*hx)
	vout[1, 14] = -hx*lamb/12 + hx*mu/12
	vout[1, 15] = -hz*lamb/24 + hz*mu/24
	vout[1, 16] = -hx*hy*mu/(18*hz) + hx*hz*lamb/(36*hy) + hx*hz*mu/(18*hy) - hy*hz*mu/(18*hx)
	vout[1, 17] = -hx*lamb/24 + hx*mu/24
	vout[1, 18] = -hz*lamb/24 - hz*mu/24
	vout[1, 19] = -hx*hy*mu/(36*hz) - hx*hz*lamb/(36*hy) - hx*hz*mu/(18*hy) - hy*hz*mu/(36*hx)
	vout[1, 20] = -hx*lamb/24 - hx*mu/24
	vout[1, 21] = hz*lamb/24 - hz*mu/24
	vout[1, 22] = -hx*hy*mu/(18*hz) - hx*hz*lamb/(18*hy) - hx*hz*mu/(9*hy) + hy*hz*mu/(36*hx)
	vout[1, 23] = -hx*lamb/12 - hx*mu/12
	vout[2, 0] = hy*lamb/12 + hy*mu/12
	vout[2, 1] = hx*lamb/12 + hx*mu/12
	vout[2, 2] = hx*hy*lamb/(9*hz) + 2*hx*hy*mu/(9*hz) + hx*hz*mu/(9*hy) + hy*hz*mu/(9*hx)
	vout[2, 3] = -hy*lamb/12 + hy*mu/12
	vout[2, 4] = hx*lamb/24 + hx*mu/24
	vout[2, 5] = hx*hy*lamb/(18*hz) + hx*hy*mu/(9*hz) + hx*hz*mu/(18*hy) - hy*hz*mu/(9*hx)
	vout[2, 6] = -hy*lamb/24 + hy*mu/24
	vout[2, 7] = -hx*lamb/24 + hx*mu/24
	vout[2, 8] = hx*hy*lamb/(36*hz) + hx*hy*mu/(18*hz) - hx*hz*mu/(18*hy) - hy*hz*mu/(18*hx)
	vout[2, 9] = hy*lamb/24 + hy*mu/24
	vout[2, 10] = -hx*lamb/12 + hx*mu/12
	vout[2, 11] = hx*hy*lamb/(18*hz) + hx*hy*mu/(9*hz) - hx*hz*mu/(9*hy) + hy*hz*mu/(18*hx)
	vout[2, 12] = hy*lamb/12 - hy*mu/12
	vout[2, 13] = hx*lamb/12 - hx*mu/12
	vout[2, 14] = -hx*hy*lamb/(9*hz) - 2*hx*hy*mu/(9*hz) + hx*hz*mu/(18*hy) + hy*hz*mu/(18*hx)
	vout[2, 15] = -hy*lamb/12 - hy*mu/12
	vout[2, 16] = hx*lamb/24 - hx*mu/24
	vout[2, 17] = -hx*hy*lamb/(18*hz) - hx*hy*mu/(9*hz) + hx*hz*mu/(36*hy) - hy*hz*mu/(18*hx)
	vout[2, 18] = -hy*lamb/24 - hy*mu/24
	vout[2, 19] = -hx*lamb/24 - hx*mu/24
	vout[2, 20] = -hx*hy*lamb/(36*hz) - hx*hy*mu/(18*hz) - hx*hz*mu/(36*hy) - hy*hz*mu/(36*hx)
	vout[2, 21] = hy*lamb/24 - hy*mu/24
	vout[2, 22] = -hx*lamb/12 - hx*mu/12
	vout[2, 23] = -hx*hy*lamb/(18*hz) - hx*hy*mu/(9*hz) - hx*hz*mu/(18*hy) + hy*hz*mu/(36*hx)
	vout[3, 0] = hx*hy*mu/(18*hz) + hx*hz*mu/(18*hy) - hy*hz*lamb/(9*hx) - 2*hy*hz*mu/(9*hx)
	vout[3, 1] = -hz*lamb/12 + hz*mu/12
	vout[3, 2] = -hy*lamb/12 + hy*mu/12
	vout[3, 3] = hx*hy*mu/(9*hz) + hx*hz*mu/(9*hy) + hy*hz*lamb/(9*hx) + 2*hy*hz*mu/(9*hx)
	vout[3, 4] = -hz*lamb/12 - hz*mu/12
	vout[3, 5] = -hy*lamb/12 - hy*mu/12
	vout[3, 6] = hx*hy*mu/(18*hz) - hx*hz*mu/(9*hy) + hy*hz*lamb/(18*hx) + hy*hz*mu/(9*hx)
	vout[3, 7] = hz*lamb/12 - hz*mu/12
	vout[3, 8] = -hy*lamb/24 - hy*mu/24
	vout[3, 9] = hx*hy*mu/(36*hz) - hx*hz*mu/(18*hy) - hy*hz*lamb/(18*hx) - hy*hz*mu/(9*hx)
	vout[3, 10] = hz*lamb/12 + hz*mu/12
	vout[3, 11] = -hy*lamb/24 + hy*mu/24
	vout[3, 12] = -hx*hy*mu/(18*hz) + hx*hz*mu/(36*hy) - hy*hz*lamb/(18*hx) - hy*hz*mu/(9*hx)
	vout[3, 13] = -hz*lamb/24 + hz*mu/24
	vout[3, 14] = hy*lamb/12 + hy*mu/12
	vout[3, 15] = -hx*hy*mu/(9*hz) + hx*hz*mu/(18*hy) + hy*hz*lamb/(18*hx) + hy*hz*mu/(9*hx)
	vout[3, 16] = -hz*lamb/24 - hz*mu/24
	vout[3, 17] = hy*lamb/12 - hy*mu/12
	vout[3, 18] = -hx*hy*mu/(18*hz) - hx*hz*mu/(18*hy) + hy*hz*lamb/(36*hx) + hy*hz*mu/(18*hx)
	vout[3, 19] = hz*lamb/24 - hz*mu/24
	vout[3, 20] = hy*lamb/24 - hy*mu/24
	vout[3, 21] = -hx*hy*mu/(36*hz) - hx*hz*mu/(36*hy) - hy*hz*lamb/(36*hx) - hy*hz*mu/(18*hx)
	vout[3, 22] = hz*lamb/24 + hz*mu/24
	vout[3, 23] = hy*lamb/24 + hy*mu/24
	vout[4, 0] = hz*lamb/12 - hz*mu/12
	vout[4, 1] = hx*hy*mu/(18*hz) + hx*hz*lamb/(18*hy) + hx*hz*mu/(9*hy) - hy*hz*mu/(9*hx)
	vout[4, 2] = hx*lamb/24 + hx*mu/24
	vout[4, 3] = -hz*lamb/12 - hz*mu/12
	vout[4, 4] = hx*hy*mu/(9*hz) + hx*hz*lamb/(9*hy) + 2*hx*hz*mu/(9*hy) + hy*hz*mu/(9*hx)
	vout[4, 5] = hx*lamb/12 + hx*mu/12
	vout[4, 6] = -hz*lamb/12 + hz*mu/12
	vout[4, 7] = hx*hy*mu/(18*hz) - hx*hz*lamb/(9*hy) - 2*hx*hz*mu/(9*hy) + hy*hz*mu/(18*hx)
	vout[4, 8] = hx*lamb/12 - hx*mu/12
	vout[4, 9] = hz*lamb/12 + hz*mu/12
	vout[4, 10] = hx*hy*mu/(36*hz) - hx*hz*lamb/(18*hy) - hx*hz*mu/(9*hy) - hy*hz*mu/(18*hx)
	vout[4, 11] = hx*lamb/24 - hx*mu/24
	vout[4, 12] = hz*lamb/24 - hz*mu/24
	vout[4, 13] = -hx*hy*mu/(18*hz) + hx*hz*lamb/(36*hy) + hx*hz*mu/(18*hy) - hy*hz*mu/(18*hx)
	vout[4, 14] = -hx*lamb/24 + hx*mu/24
	vout[4, 15] = -hz*lamb/24 - hz*mu/24
	vout[4, 16] = -hx*hy*mu/(9*hz) + hx*hz*lamb/(18*hy) + hx*hz*mu/(9*hy) + hy*hz*mu/(18*hx)
	vout[4, 17] = -hx*lamb/12 + hx*mu/12
	vout[4, 18] = -hz*lamb/24 + hz*mu/24
	vout[4, 19] = -hx*hy*mu/(18*hz) - hx*hz*lamb/(18*hy) - hx*hz*mu/(9*hy) + hy*hz*mu/(36*hx)
	vout[4, 20] = -hx*lamb/12 - hx*mu/12
	vout[4, 21] = hz*lamb/24 + hz*mu/24
	vout[4, 22] = -hx*hy*mu/(36*hz) - hx*hz*lamb/(36*hy) - hx*hz*mu/(18*hy) - hy*hz*mu/(36*hx)
	vout[4, 23] = -hx*lamb/24 - hx*mu/24
	vout[5, 0] = hy*lamb/12 - hy*mu/12
	vout[5, 1] = hx*lamb/24 + hx*mu/24
	vout[5, 2] = hx*hy*lamb/(18*hz) + hx*hy*mu/(9*hz) + hx*hz*mu/(18*hy) - hy*hz*mu/(9*hx)
	vout[5, 3] = -hy*lamb/12 - hy*mu/12
	vout[5, 4] = hx*lamb/12 + hx*mu/12
	vout[5, 5] = hx*hy*lamb/(9*hz) + 2*hx*hy*mu/(9*hz) + hx*hz*mu/(9*hy) + hy*hz*mu/(9*hx)
	vout[5, 6] = -hy*lamb/24 - hy*mu/24
	vout[5, 7] = -hx*lamb/12 + hx*mu/12
	vout[5, 8] = hx*hy*lamb/(18*hz) + hx*hy*mu/(9*hz) - hx*hz*mu/(9*hy) + hy*hz*mu/(18*hx)
	vout[5, 9] = hy*lamb/24 - hy*mu/24
	vout[5, 10] = -hx*lamb/24 + hx*mu/24
	vout[5, 11] = hx*hy*lamb/(36*hz) + hx*hy*mu/(18*hz) - hx*hz*mu/(18*hy) - hy*hz*mu/(18*hx)
	vout[5, 12] = hy*lamb/12 + hy*mu/12
	vout[5, 13] = hx*lamb/24 - hx*mu/24
	vout[5, 14] = -hx*hy*lamb/(18*hz) - hx*hy*mu/(9*hz) + hx*hz*mu/(36*hy) - hy*hz*mu/(18*hx)
	vout[5, 15] = -hy*lamb/12 + hy*mu/12
	vout[5, 16] = hx*lamb/12 - hx*mu/12
	vout[5, 17] = -hx*hy*lamb/(9*hz) - 2*hx*hy*mu/(9*hz) + hx*hz*mu/(18*hy) + hy*hz*mu/(18*hx)
	vout[5, 18] = -hy*lamb/24 + hy*mu/24
	vout[5, 19] = -hx*lamb/12 - hx*mu/12
	vout[5, 20] = -hx*hy*lamb/(18*hz) - hx*hy*mu/(9*hz) - hx*hz*mu/(18*hy) + hy*hz*mu/(36*hx)
	vout[5, 21] = hy*lamb/24 + hy*mu/24
	vout[5, 22] = -hx*lamb/24 - hx*mu/24
	vout[5, 23] = -hx*hy*lamb/(36*hz) - hx*hy*mu/(18*hz) - hx*hz*mu/(36*hy) - hy*hz*mu/(36*hx)
	vout[6, 0] = hx*hy*mu/(36*hz) - hx*hz*mu/(18*hy) - hy*hz*lamb/(18*hx) - hy*hz*mu/(9*hx)
	vout[6, 1] = -hz*lamb/12 - hz*mu/12
	vout[6, 2] = -hy*lamb/24 + hy*mu/24
	vout[6, 3] = hx*hy*mu/(18*hz) - hx*hz*mu/(9*hy) + hy*hz*lamb/(18*hx) + hy*hz*mu/(9*hx)
	vout[6, 4] = -hz*lamb/12 + hz*mu/12
	vout[6, 5] = -hy*lamb/24 - hy*mu/24
	vout[6, 6] = hx*hy*mu/(9*hz) + hx*hz*mu/(9*hy) + hy*hz*lamb/(9*hx) + 2*hy*hz*mu/(9*hx)
	vout[6, 7] = hz*lamb/12 + hz*mu/12
	vout[6, 8] = -hy*lamb/12 - hy*mu/12
	vout[6, 9] = hx*hy*mu/(18*hz) + hx*hz*mu/(18*hy) - hy*hz*lamb/(9*hx) - 2*hy*hz*mu/(9*hx)
	vout[6, 10] = hz*lamb/12 - hz*mu/12
	vout[6, 11] = -hy*lamb/12 + hy*mu/12
	vout[6, 12] = -hx*hy*mu/(36*hz) - hx*hz*mu/(36*hy) - hy*hz*lamb/(36*hx) - hy*hz*mu/(18*hx)
	vout[6, 13] = -hz*lamb/24 - hz*mu/24
	vout[6, 14] = hy*lamb/24 + hy*mu/24
	vout[6, 15] = -hx*hy*mu/(18*hz) - hx*hz*mu/(18*hy) + hy*hz*lamb/(36*hx) + hy*hz*mu/(18*hx)
	vout[6, 16] = -hz*lamb/24 + hz*mu/24
	vout[6, 17] = hy*lamb/24 - hy*mu/24
	vout[6, 18] = -hx*hy*mu/(9*hz) + hx*hz*mu/(18*hy) + hy*hz*lamb/(18*hx) + hy*hz*mu/(9*hx)
	vout[6, 19] = hz*lamb/24 + hz*mu/24
	vout[6, 20] = hy*lamb/12 - hy*mu/12
	vout[6, 21] = -hx*hy*mu/(18*hz) + hx*hz*mu/(36*hy) - hy*hz*lamb/(18*hx) - hy*hz*mu/(9*hx)
	vout[6, 22] = hz*lamb/24 - hz*mu/24
	vout[6, 23] = hy*lamb/12 + hy*mu/12
	vout[7, 0] = -hz*lamb/12 - hz*mu/12
	vout[7, 1] = hx*hy*mu/(36*hz) - hx*hz*lamb/(18*hy) - hx*hz*mu/(9*hy) - hy*hz*mu/(18*hx)
	vout[7, 2] = -hx*lamb/24 + hx*mu/24
	vout[7, 3] = hz*lamb/12 - hz*mu/12
	vout[7, 4] = hx*hy*mu/(18*hz) - hx*hz*lamb/(9*hy) - 2*hx*hz*mu/(9*hy) + hy*hz*mu/(18*hx)
	vout[7, 5] = -hx*lamb/12 + hx*mu/12
	vout[7, 6] = hz*lamb/12 + hz*mu/12
	vout[7, 7] = hx*hy*mu/(9*hz) + hx*hz*lamb/(9*hy) + 2*hx*hz*mu/(9*hy) + hy*hz*mu/(9*hx)
	vout[7, 8] = -hx*lamb/12 - hx*mu/12
	vout[7, 9] = -hz*lamb/12 + hz*mu/12
	vout[7, 10] = hx*hy*mu/(18*hz) + hx*hz*lamb/(18*hy) + hx*hz*mu/(9*hy) - hy*hz*mu/(9*hx)
	vout[7, 11] = -hx*lamb/24 - hx*mu/24
	vout[7, 12] = -hz*lamb/24 - hz*mu/24
	vout[7, 13] = -hx*hy*mu/(36*hz) - hx*hz*lamb/(36*hy) - hx*hz*mu/(18*hy) - hy*hz*mu/(36*hx)
	vout[7, 14] = hx*lamb/24 + hx*mu/24
	vout[7, 15] = hz*lamb/24 - hz*mu/24
	vout[7, 16] = -hx*hy*mu/(18*hz) - hx*hz*lamb/(18*hy) - hx*hz*mu/(9*hy) + hy*hz*mu/(36*hx)
	vout[7, 17] = hx*lamb/12 + hx*mu/12
	vout[7, 18] = hz*lamb/24 + hz*mu/24
	vout[7, 19] = -hx*hy*mu/(9*hz) + hx*hz*lamb/(18*hy) + hx*hz*mu/(9*hy) + hy*hz*mu/(18*hx)
	vout[7, 20] = hx*lamb/12 - hx*mu/12
	vout[7, 21] = -hz*lamb/24 + hz*mu/24
	vout[7, 22] = -hx*hy*mu/(18*hz) + hx*hz*lamb/(36*hy) + hx*hz*mu/(18*hy) - hy*hz*mu/(18*hx)
	vout[7, 23] = hx*lamb/24 - hx*mu/24
	vout[8, 0] = hy*lamb/24 - hy*mu/24
	vout[8, 1] = hx*lamb/24 - hx*mu/24
	vout[8, 2] = hx*hy*lamb/(36*hz) + hx*hy*mu/(18*hz) - hx*hz*mu/(18*hy) - hy*hz*mu/(18*hx)
	vout[8, 3] = -hy*lamb/24 - hy*mu/24
	vout[8, 4] = hx*lamb/12 - hx*mu/12
	vout[8, 5] = hx*hy*lamb/(18*hz) + hx*hy*mu/(9*hz) - hx*hz*mu/(9*hy) + hy*hz*mu/(18*hx)
	vout[8, 6] = -hy*lamb/12 - hy*mu/12
	vout[8, 7] = -hx*lamb/12 - hx*mu/12
	vout[8, 8] = hx*hy*lamb/(9*hz) + 2*hx*hy*mu/(9*hz) + hx*hz*mu/(9*hy) + hy*hz*mu/(9*hx)
	vout[8, 9] = hy*lamb/12 - hy*mu/12
	vout[8, 10] = -hx*lamb/24 - hx*mu/24
	vout[8, 11] = hx*hy*lamb/(18*hz) + hx*hy*mu/(9*hz) + hx*hz*mu/(18*hy) - hy*hz*mu/(9*hx)
	vout[8, 12] = hy*lamb/24 + hy*mu/24
	vout[8, 13] = hx*lamb/24 + hx*mu/24
	vout[8, 14] = -hx*hy*lamb/(36*hz) - hx*hy*mu/(18*hz) - hx*hz*mu/(36*hy) - hy*hz*mu/(36*hx)
	vout[8, 15] = -hy*lamb/24 + hy*mu/24
	vout[8, 16] = hx*lamb/12 + hx*mu/12
	vout[8, 17] = -hx*hy*lamb/(18*hz) - hx*hy*mu/(9*hz) - hx*hz*mu/(18*hy) + hy*hz*mu/(36*hx)
	vout[8, 18] = -hy*lamb/12 + hy*mu/12
	vout[8, 19] = -hx*lamb/12 + hx*mu/12
	vout[8, 20] = -hx*hy*lamb/(9*hz) - 2*hx*hy*mu/(9*hz) + hx*hz*mu/(18*hy) + hy*hz*mu/(18*hx)
	vout[8, 21] = hy*lamb/12 + hy*mu/12
	vout[8, 22] = -hx*lamb/24 + hx*mu/24
	vout[8, 23] = -hx*hy*lamb/(18*hz) - hx*hy*mu/(9*hz) + hx*hz*mu/(36*hy) - hy*hz*mu/(18*hx)
	vout[9, 0] = hx*hy*mu/(18*hz) - hx*hz*mu/(9*hy) + hy*hz*lamb/(18*hx) + hy*hz*mu/(9*hx)
	vout[9, 1] = hz*lamb/12 - hz*mu/12
	vout[9, 2] = hy*lamb/24 + hy*mu/24
	vout[9, 3] = hx*hy*mu/(36*hz) - hx*hz*mu/(18*hy) - hy*hz*lamb/(18*hx) - hy*hz*mu/(9*hx)
	vout[9, 4] = hz*lamb/12 + hz*mu/12
	vout[9, 5] = hy*lamb/24 - hy*mu/24
	vout[9, 6] = hx*hy*mu/(18*hz) + hx*hz*mu/(18*hy) - hy*hz*lamb/(9*hx) - 2*hy*hz*mu/(9*hx)
	vout[9, 7] = -hz*lamb/12 + hz*mu/12
	vout[9, 8] = hy*lamb/12 - hy*mu/12
	vout[9, 9] = hx*hy*mu/(9*hz) + hx*hz*mu/(9*hy) + hy*hz*lamb/(9*hx) + 2*hy*hz*mu/(9*hx)
	vout[9, 10] = -hz*lamb/12 - hz*mu/12
	vout[9, 11] = hy*lamb/12 + hy*mu/12
	vout[9, 12] = -hx*hy*mu/(18*hz) - hx*hz*mu/(18*hy) + hy*hz*lamb/(36*hx) + hy*hz*mu/(18*hx)
	vout[9, 13] = hz*lamb/24 - hz*mu/24
	vout[9, 14] = -hy*lamb/24 + hy*mu/24
	vout[9, 15] = -hx*hy*mu/(36*hz) - hx*hz*mu/(36*hy) - hy*hz*lamb/(36*hx) - hy*hz*mu/(18*hx)
	vout[9, 16] = hz*lamb/24 + hz*mu/24
	vout[9, 17] = -hy*lamb/24 - hy*mu/24
	vout[9, 18] = -hx*hy*mu/(18*hz) + hx*hz*mu/(36*hy) - hy*hz*lamb/(18*hx) - hy*hz*mu/(9*hx)
	vout[9, 19] = -hz*lamb/24 + hz*mu/24
	vout[9, 20] = -hy*lamb/12 - hy*mu/12
	vout[9, 21] = -hx*hy*mu/(9*hz) + hx*hz*mu/(18*hy) + hy*hz*lamb/(18*hx) + hy*hz*mu/(9*hx)
	vout[9, 22] = -hz*lamb/24 - hz*mu/24
	vout[9, 23] = -hy*lamb/12 + hy*mu/12
	vout[10, 0] = -hz*lamb/12 + hz*mu/12
	vout[10, 1] = hx*hy*mu/(18*hz) - hx*hz*lamb/(9*hy) - 2*hx*hz*mu/(9*hy) + hy*hz*mu/(18*hx)
	vout[10, 2] = -hx*lamb/12 + hx*mu/12
	vout[10, 3] = hz*lamb/12 + hz*mu/12
	vout[10, 4] = hx*hy*mu/(36*hz) - hx*hz*lamb/(18*hy) - hx*hz*mu/(9*hy) - hy*hz*mu/(18*hx)
	vout[10, 5] = -hx*lamb/24 + hx*mu/24
	vout[10, 6] = hz*lamb/12 - hz*mu/12
	vout[10, 7] = hx*hy*mu/(18*hz) + hx*hz*lamb/(18*hy) + hx*hz*mu/(9*hy) - hy*hz*mu/(9*hx)
	vout[10, 8] = -hx*lamb/24 - hx*mu/24
	vout[10, 9] = -hz*lamb/12 - hz*mu/12
	vout[10, 10] = hx*hy*mu/(9*hz) + hx*hz*lamb/(9*hy) + 2*hx*hz*mu/(9*hy) + hy*hz*mu/(9*hx)
	vout[10, 11] = -hx*lamb/12 - hx*mu/12
	vout[10, 12] = -hz*lamb/24 + hz*mu/24
	vout[10, 13] = -hx*hy*mu/(18*hz) - hx*hz*lamb/(18*hy) - hx*hz*mu/(9*hy) + hy*hz*mu/(36*hx)
	vout[10, 14] = hx*lamb/12 + hx*mu/12
	vout[10, 15] = hz*lamb/24 + hz*mu/24
	vout[10, 16] = -hx*hy*mu/(36*hz) - hx*hz*lamb/(36*hy) - hx*hz*mu/(18*hy) - hy*hz*mu/(36*hx)
	vout[10, 17] = hx*lamb/24 + hx*mu/24
	vout[10, 18] = hz*lamb/24 - hz*mu/24
	vout[10, 19] = -hx*hy*mu/(18*hz) + hx*hz*lamb/(36*hy) + hx*hz*mu/(18*hy) - hy*hz*mu/(18*hx)
	vout[10, 20] = hx*lamb/24 - hx*mu/24
	vout[10, 21] = -hz*lamb/24 - hz*mu/24
	vout[10, 22] = -hx*hy*mu/(9*hz) + hx*hz*lamb/(18*hy) + hx*hz*mu/(9*hy) + hy*hz*mu/(18*hx)
	vout[10, 23] = hx*lamb/12 - hx*mu/12
	vout[11, 0] = hy*lamb/24 + hy*mu/24
	vout[11, 1] = hx*lamb/12 - hx*mu/12
	vout[11, 2] = hx*hy*lamb/(18*hz) + hx*hy*mu/(9*hz) - hx*hz*mu/(9*hy) + hy*hz*mu/(18*hx)
	vout[11, 3] = -hy*lamb/24 + hy*mu/24
	vout[11, 4] = hx*lamb/24 - hx*mu/24
	vout[11, 5] = hx*hy*lamb/(36*hz) + hx*hy*mu/(18*hz) - hx*hz*mu/(18*hy) - hy*hz*mu/(18*hx)
	vout[11, 6] = -hy*lamb/12 + hy*mu/12
	vout[11, 7] = -hx*lamb/24 - hx*mu/24
	vout[11, 8] = hx*hy*lamb/(18*hz) + hx*hy*mu/(9*hz) + hx*hz*mu/(18*hy) - hy*hz*mu/(9*hx)
	vout[11, 9] = hy*lamb/12 + hy*mu/12
	vout[11, 10] = -hx*lamb/12 - hx*mu/12
	vout[11, 11] = hx*hy*lamb/(9*hz) + 2*hx*hy*mu/(9*hz) + hx*hz*mu/(9*hy) + hy*hz*mu/(9*hx)
	vout[11, 12] = hy*lamb/24 - hy*mu/24
	vout[11, 13] = hx*lamb/12 + hx*mu/12
	vout[11, 14] = -hx*hy*lamb/(18*hz) - hx*hy*mu/(9*hz) - hx*hz*mu/(18*hy) + hy*hz*mu/(36*hx)
	vout[11, 15] = -hy*lamb/24 - hy*mu/24
	vout[11, 16] = hx*lamb/24 + hx*mu/24
	vout[11, 17] = -hx*hy*lamb/(36*hz) - hx*hy*mu/(18*hz) - hx*hz*mu/(36*hy) - hy*hz*mu/(36*hx)
	vout[11, 18] = -hy*lamb/12 - hy*mu/12
	vout[11, 19] = -hx*lamb/24 + hx*mu/24
	vout[11, 20] = -hx*hy*lamb/(18*hz) - hx*hy*mu/(9*hz) + hx*hz*mu/(36*hy) - hy*hz*mu/(18*hx)
	vout[11, 21] = hy*lamb/12 - hy*mu/12
	vout[11, 22] = -hx*lamb/12 + hx*mu/12
	vout[11, 23] = -hx*hy*lamb/(9*hz) - 2*hx*hy*mu/(9*hz) + hx*hz*mu/(18*hy) + hy*hz*mu/(18*hx)
	vout[12, 0] = -hx*hy*mu/(9*hz) + hx*hz*mu/(18*hy) + hy*hz*lamb/(18*hx) + hy*hz*mu/(9*hx)
	vout[12, 1] = hz*lamb/24 + hz*mu/24
	vout[12, 2] = hy*lamb/12 - hy*mu/12
	vout[12, 3] = -hx*hy*mu/(18*hz) + hx*hz*mu/(36*hy) - hy*hz*lamb/(18*hx) - hy*hz*mu/(9*hx)
	vout[12, 4] = hz*lamb/24 - hz*mu/24
	vout[12, 5] = hy*lamb/12 + hy*mu/12
	vout[12, 6] = -hx*hy*mu/(36*hz) - hx*hz*mu/(36*hy) - hy*hz*lamb/(36*hx) - hy*hz*mu/(18*hx)
	vout[12, 7] = -hz*lamb/24 - hz*mu/24
	vout[12, 8] = hy*lamb/24 + hy*mu/24
	vout[12, 9] = -hx*hy*mu/(18*hz) - hx*hz*mu/(18*hy) + hy*hz*lamb/(36*hx) + hy*hz*mu/(18*hx)
	vout[12, 10] = -hz*lamb/24 + hz*mu/24
	vout[12, 11] = hy*lamb/24 - hy*mu/24
	vout[12, 12] = hx*hy*mu/(9*hz) + hx*hz*mu/(9*hy) + hy*hz*lamb/(9*hx) + 2*hy*hz*mu/(9*hx)
	vout[12, 13] = hz*lamb/12 + hz*mu/12
	vout[12, 14] = -hy*lamb/12 - hy*mu/12
	vout[12, 15] = hx*hy*mu/(18*hz) + hx*hz*mu/(18*hy) - hy*hz*lamb/(9*hx) - 2*hy*hz*mu/(9*hx)
	vout[12, 16] = hz*lamb/12 - hz*mu/12
	vout[12, 17] = -hy*lamb/12 + hy*mu/12
	vout[12, 18] = hx*hy*mu/(36*hz) - hx*hz*mu/(18*hy) - hy*hz*lamb/(18*hx) - hy*hz*mu/(9*hx)
	vout[12, 19] = -hz*lamb/12 - hz*mu/12
	vout[12, 20] = -hy*lamb/24 + hy*mu/24
	vout[12, 21] = hx*hy*mu/(18*hz) - hx*hz*mu/(9*hy) + hy*hz*lamb/(18*hx) + hy*hz*mu/(9*hx)
	vout[12, 22] = -hz*lamb/12 + hz*mu/12
	vout[12, 23] = -hy*lamb/24 - hy*mu/24
	vout[13, 0] = hz*lamb/24 + hz*mu/24
	vout[13, 1] = -hx*hy*mu/(9*hz) + hx*hz*lamb/(18*hy) + hx*hz*mu/(9*hy) + hy*hz*mu/(18*hx)
	vout[13, 2] = hx*lamb/12 - hx*mu/12
	vout[13, 3] = -hz*lamb/24 + hz*mu/24
	vout[13, 4] = -hx*hy*mu/(18*hz) + hx*hz*lamb/(36*hy) + hx*hz*mu/(18*hy) - hy*hz*mu/(18*hx)
	vout[13, 5] = hx*lamb/24 - hx*mu/24
	vout[13, 6] = -hz*lamb/24 - hz*mu/24
	vout[13, 7] = -hx*hy*mu/(36*hz) - hx*hz*lamb/(36*hy) - hx*hz*mu/(18*hy) - hy*hz*mu/(36*hx)
	vout[13, 8] = hx*lamb/24 + hx*mu/24
	vout[13, 9] = hz*lamb/24 - hz*mu/24
	vout[13, 10] = -hx*hy*mu/(18*hz) - hx*hz*lamb/(18*hy) - hx*hz*mu/(9*hy) + hy*hz*mu/(36*hx)
	vout[13, 11] = hx*lamb/12 + hx*mu/12
	vout[13, 12] = hz*lamb/12 + hz*mu/12
	vout[13, 13] = hx*hy*mu/(9*hz) + hx*hz*lamb/(9*hy) + 2*hx*hz*mu/(9*hy) + hy*hz*mu/(9*hx)
	vout[13, 14] = -hx*lamb/12 - hx*mu/12
	vout[13, 15] = -hz*lamb/12 + hz*mu/12
	vout[13, 16] = hx*hy*mu/(18*hz) + hx*hz*lamb/(18*hy) + hx*hz*mu/(9*hy) - hy*hz*mu/(9*hx)
	vout[13, 17] = -hx*lamb/24 - hx*mu/24
	vout[13, 18] = -hz*lamb/12 - hz*mu/12
	vout[13, 19] = hx*hy*mu/(36*hz) - hx*hz*lamb/(18*hy) - hx*hz*mu/(9*hy) - hy*hz*mu/(18*hx)
	vout[13, 20] = -hx*lamb/24 + hx*mu/24
	vout[13, 21] = hz*lamb/12 - hz*mu/12
	vout[13, 22] = hx*hy*mu/(18*hz) - hx*hz*lamb/(9*hy) - 2*hx*hz*mu/(9*hy) + hy*hz*mu/(18*hx)
	vout[13, 23] = -hx*lamb/12 + hx*mu/12
	vout[14, 0] = -hy*lamb/12 + hy*mu/12
	vout[14, 1] = -hx*lamb/12 + hx*mu/12
	vout[14, 2] = -hx*hy*lamb/(9*hz) - 2*hx*hy*mu/(9*hz) + hx*hz*mu/(18*hy) + hy*hz*mu/(18*hx)
	vout[14, 3] = hy*lamb/12 + hy*mu/12
	vout[14, 4] = -hx*lamb/24 + hx*mu/24
	vout[14, 5] = -hx*hy*lamb/(18*hz) - hx*hy*mu/(9*hz) + hx*hz*mu/(36*hy) - hy*hz*mu/(18*hx)
	vout[14, 6] = hy*lamb/24 + hy*mu/24
	vout[14, 7] = hx*lamb/24 + hx*mu/24
	vout[14, 8] = -hx*hy*lamb/(36*hz) - hx*hy*mu/(18*hz) - hx*hz*mu/(36*hy) - hy*hz*mu/(36*hx)
	vout[14, 9] = -hy*lamb/24 + hy*mu/24
	vout[14, 10] = hx*lamb/12 + hx*mu/12
	vout[14, 11] = -hx*hy*lamb/(18*hz) - hx*hy*mu/(9*hz) - hx*hz*mu/(18*hy) + hy*hz*mu/(36*hx)
	vout[14, 12] = -hy*lamb/12 - hy*mu/12
	vout[14, 13] = -hx*lamb/12 - hx*mu/12
	vout[14, 14] = hx*hy*lamb/(9*hz) + 2*hx*hy*mu/(9*hz) + hx*hz*mu/(9*hy) + hy*hz*mu/(9*hx)
	vout[14, 15] = hy*lamb/12 - hy*mu/12
	vout[14, 16] = -hx*lamb/24 - hx*mu/24
	vout[14, 17] = hx*hy*lamb/(18*hz) + hx*hy*mu/(9*hz) + hx*hz*mu/(18*hy) - hy*hz*mu/(9*hx)
	vout[14, 18] = hy*lamb/24 - hy*mu/24
	vout[14, 19] = hx*lamb/24 - hx*mu/24
	vout[14, 20] = hx*hy*lamb/(36*hz) + hx*hy*mu/(18*hz) - hx*hz*mu/(18*hy) - hy*hz*mu/(18*hx)
	vout[14, 21] = -hy*lamb/24 - hy*mu/24
	vout[14, 22] = hx*lamb/12 - hx*mu/12
	vout[14, 23] = hx*hy*lamb/(18*hz) + hx*hy*mu/(9*hz) - hx*hz*mu/(9*hy) + hy*hz*mu/(18*hx)
	vout[15, 0] = -hx*hy*mu/(18*hz) + hx*hz*mu/(36*hy) - hy*hz*lamb/(18*hx) - hy*hz*mu/(9*hx)
	vout[15, 1] = -hz*lamb/24 + hz*mu/24
	vout[15, 2] = -hy*lamb/12 - hy*mu/12
	vout[15, 3] = -hx*hy*mu/(9*hz) + hx*hz*mu/(18*hy) + hy*hz*lamb/(18*hx) + hy*hz*mu/(9*hx)
	vout[15, 4] = -hz*lamb/24 - hz*mu/24
	vout[15, 5] = -hy*lamb/12 + hy*mu/12
	vout[15, 6] = -hx*hy*mu/(18*hz) - hx*hz*mu/(18*hy) + hy*hz*lamb/(36*hx) + hy*hz*mu/(18*hx)
	vout[15, 7] = hz*lamb/24 - hz*mu/24
	vout[15, 8] = -hy*lamb/24 + hy*mu/24
	vout[15, 9] = -hx*hy*mu/(36*hz) - hx*hz*mu/(36*hy) - hy*hz*lamb/(36*hx) - hy*hz*mu/(18*hx)
	vout[15, 10] = hz*lamb/24 + hz*mu/24
	vout[15, 11] = -hy*lamb/24 - hy*mu/24
	vout[15, 12] = hx*hy*mu/(18*hz) + hx*hz*mu/(18*hy) - hy*hz*lamb/(9*hx) - 2*hy*hz*mu/(9*hx)
	vout[15, 13] = -hz*lamb/12 + hz*mu/12
	vout[15, 14] = hy*lamb/12 - hy*mu/12
	vout[15, 15] = hx*hy*mu/(9*hz) + hx*hz*mu/(9*hy) + hy*hz*lamb/(9*hx) + 2*hy*hz*mu/(9*hx)
	vout[15, 16] = -hz*lamb/12 - hz*mu/12
	vout[15, 17] = hy*lamb/12 + hy*mu/12
	vout[15, 18] = hx*hy*mu/(18*hz) - hx*hz*mu/(9*hy) + hy*hz*lamb/(18*hx) + hy*hz*mu/(9*hx)
	vout[15, 19] = hz*lamb/12 - hz*mu/12
	vout[15, 20] = hy*lamb/24 + hy*mu/24
	vout[15, 21] = hx*hy*mu/(36*hz) - hx*hz*mu/(18*hy) - hy*hz*lamb/(18*hx) - hy*hz*mu/(9*hx)
	vout[15, 22] = hz*lamb/12 + hz*mu/12
	vout[15, 23] = hy*lamb/24 - hy*mu/24
	vout[16, 0] = hz*lamb/24 - hz*mu/24
	vout[16, 1] = -hx*hy*mu/(18*hz) + hx*hz*lamb/(36*hy) + hx*hz*mu/(18*hy) - hy*hz*mu/(18*hx)
	vout[16, 2] = hx*lamb/24 - hx*mu/24
	vout[16, 3] = -hz*lamb/24 - hz*mu/24
	vout[16, 4] = -hx*hy*mu/(9*hz) + hx*hz*lamb/(18*hy) + hx*hz*mu/(9*hy) + hy*hz*mu/(18*hx)
	vout[16, 5] = hx*lamb/12 - hx*mu/12
	vout[16, 6] = -hz*lamb/24 + hz*mu/24
	vout[16, 7] = -hx*hy*mu/(18*hz) - hx*hz*lamb/(18*hy) - hx*hz*mu/(9*hy) + hy*hz*mu/(36*hx)
	vout[16, 8] = hx*lamb/12 + hx*mu/12
	vout[16, 9] = hz*lamb/24 + hz*mu/24
	vout[16, 10] = -hx*hy*mu/(36*hz) - hx*hz*lamb/(36*hy) - hx*hz*mu/(18*hy) - hy*hz*mu/(36*hx)
	vout[16, 11] = hx*lamb/24 + hx*mu/24
	vout[16, 12] = hz*lamb/12 - hz*mu/12
	vout[16, 13] = hx*hy*mu/(18*hz) + hx*hz*lamb/(18*hy) + hx*hz*mu/(9*hy) - hy*hz*mu/(9*hx)
	vout[16, 14] = -hx*lamb/24 - hx*mu/24
	vout[16, 15] = -hz*lamb/12 - hz*mu/12
	vout[16, 16] = hx*hy*mu/(9*hz) + hx*hz*lamb/(9*hy) + 2*hx*hz*mu/(9*hy) + hy*hz*mu/(9*hx)
	vout[16, 17] = -hx*lamb/12 - hx*mu/12
	vout[16, 18] = -hz*lamb/12 + hz*mu/12
	vout[16, 19] = hx*hy*mu/(18*hz) - hx*hz*lamb/(9*hy) - 2*hx*hz*mu/(9*hy) + hy*hz*mu/(18*hx)
	vout[16, 20] = -hx*lamb/12 + hx*mu/12
	vout[16, 21] = hz*lamb/12 + hz*mu/12
	vout[16, 22] = hx*hy*mu/(36*hz) - hx*hz*lamb/(18*hy) - hx*hz*mu/(9*hy) - hy*hz*mu/(18*hx)
	vout[16, 23] = -hx*lamb/24 + hx*mu/24
	vout[17, 0] = -hy*lamb/12 - hy*mu/12
	vout[17, 1] = -hx*lamb/24 + hx*mu/24
	vout[17, 2] = -hx*hy*lamb/(18*hz) - hx*hy*mu/(9*hz) + hx*hz*mu/(36*hy) - hy*hz*mu/(18*hx)
	vout[17, 3] = hy*lamb/12 - hy*mu/12
	vout[17, 4] = -hx*lamb/12 + hx*mu/12
	vout[17, 5] = -hx*hy*lamb/(9*hz) - 2*hx*hy*mu/(9*hz) + hx*hz*mu/(18*hy) + hy*hz*mu/(18*hx)
	vout[17, 6] = hy*lamb/24 - hy*mu/24
	vout[17, 7] = hx*lamb/12 + hx*mu/12
	vout[17, 8] = -hx*hy*lamb/(18*hz) - hx*hy*mu/(9*hz) - hx*hz*mu/(18*hy) + hy*hz*mu/(36*hx)
	vout[17, 9] = -hy*lamb/24 - hy*mu/24
	vout[17, 10] = hx*lamb/24 + hx*mu/24
	vout[17, 11] = -hx*hy*lamb/(36*hz) - hx*hy*mu/(18*hz) - hx*hz*mu/(36*hy) - hy*hz*mu/(36*hx)
	vout[17, 12] = -hy*lamb/12 + hy*mu/12
	vout[17, 13] = -hx*lamb/24 - hx*mu/24
	vout[17, 14] = hx*hy*lamb/(18*hz) + hx*hy*mu/(9*hz) + hx*hz*mu/(18*hy) - hy*hz*mu/(9*hx)
	vout[17, 15] = hy*lamb/12 + hy*mu/12
	vout[17, 16] = -hx*lamb/12 - hx*mu/12
	vout[17, 17] = hx*hy*lamb/(9*hz) + 2*hx*hy*mu/(9*hz) + hx*hz*mu/(9*hy) + hy*hz*mu/(9*hx)
	vout[17, 18] = hy*lamb/24 + hy*mu/24
	vout[17, 19] = hx*lamb/12 - hx*mu/12
	vout[17, 20] = hx*hy*lamb/(18*hz) + hx*hy*mu/(9*hz) - hx*hz*mu/(9*hy) + hy*hz*mu/(18*hx)
	vout[17, 21] = -hy*lamb/24 + hy*mu/24
	vout[17, 22] = hx*lamb/24 - hx*mu/24
	vout[17, 23] = hx*hy*lamb/(36*hz) + hx*hy*mu/(18*hz) - hx*hz*mu/(18*hy) - hy*hz*mu/(18*hx)
	vout[18, 0] = -hx*hy*mu/(36*hz) - hx*hz*mu/(36*hy) - hy*hz*lamb/(36*hx) - hy*hz*mu/(18*hx)
	vout[18, 1] = -hz*lamb/24 - hz*mu/24
	vout[18, 2] = -hy*lamb/24 - hy*mu/24
	vout[18, 3] = -hx*hy*mu/(18*hz) - hx*hz*mu/(18*hy) + hy*hz*lamb/(36*hx) + hy*hz*mu/(18*hx)
	vout[18, 4] = -hz*lamb/24 + hz*mu/24
	vout[18, 5] = -hy*lamb/24 + hy*mu/24
	vout[18, 6] = -hx*hy*mu/(9*hz) + hx*hz*mu/(18*hy) + hy*hz*lamb/(18*hx) + hy*hz*mu/(9*hx)
	vout[18, 7] = hz*lamb/24 + hz*mu/24
	vout[18, 8] = -hy*lamb/12 + hy*mu/12
	vout[18, 9] = -hx*hy*mu/(18*hz) + hx*hz*mu/(36*hy) - hy*hz*lamb/(18*hx) - hy*hz*mu/(9*hx)
	vout[18, 10] = hz*lamb/24 - hz*mu/24
	vout[18, 11] = -hy*lamb/12 - hy*mu/12
	vout[18, 12] = hx*hy*mu/(36*hz) - hx*hz*mu/(18*hy) - hy*hz*lamb/(18*hx) - hy*hz*mu/(9*hx)
	vout[18, 13] = -hz*lamb/12 - hz*mu/12
	vout[18, 14] = hy*lamb/24 - hy*mu/24
	vout[18, 15] = hx*hy*mu/(18*hz) - hx*hz*mu/(9*hy) + hy*hz*lamb/(18*hx) + hy*hz*mu/(9*hx)
	vout[18, 16] = -hz*lamb/12 + hz*mu/12
	vout[18, 17] = hy*lamb/24 + hy*mu/24
	vout[18, 18] = hx*hy*mu/(9*hz) + hx*hz*mu/(9*hy) + hy*hz*lamb/(9*hx) + 2*hy*hz*mu/(9*hx)
	vout[18, 19] = hz*lamb/12 + hz*mu/12
	vout[18, 20] = hy*lamb/12 + hy*mu/12
	vout[18, 21] = hx*hy*mu/(18*hz) + hx*hz*mu/(18*hy) - hy*hz*lamb/(9*hx) - 2*hy*hz*mu/(9*hx)
	vout[18, 22] = hz*lamb/12 - hz*mu/12
	vout[18, 23] = hy*lamb/12 - hy*mu/12
	vout[19, 0] = -hz*lamb/24 - hz*mu/24
	vout[19, 1] = -hx*hy*mu/(36*hz) - hx*hz*lamb/(36*hy) - hx*hz*mu/(18*hy) - hy*hz*mu/(36*hx)
	vout[19, 2] = -hx*lamb/24 - hx*mu/24
	vout[19, 3] = hz*lamb/24 - hz*mu/24
	vout[19, 4] = -hx*hy*mu/(18*hz) - hx*hz*lamb/(18*hy) - hx*hz*mu/(9*hy) + hy*hz*mu/(36*hx)
	vout[19, 5] = -hx*lamb/12 - hx*mu/12
	vout[19, 6] = hz*lamb/24 + hz*mu/24
	vout[19, 7] = -hx*hy*mu/(9*hz) + hx*hz*lamb/(18*hy) + hx*hz*mu/(9*hy) + hy*hz*mu/(18*hx)
	vout[19, 8] = -hx*lamb/12 + hx*mu/12
	vout[19, 9] = -hz*lamb/24 + hz*mu/24
	vout[19, 10] = -hx*hy*mu/(18*hz) + hx*hz*lamb/(36*hy) + hx*hz*mu/(18*hy) - hy*hz*mu/(18*hx)
	vout[19, 11] = -hx*lamb/24 + hx*mu/24
	vout[19, 12] = -hz*lamb/12 - hz*mu/12
	vout[19, 13] = hx*hy*mu/(36*hz) - hx*hz*lamb/(18*hy) - hx*hz*mu/(9*hy) - hy*hz*mu/(18*hx)
	vout[19, 14] = hx*lamb/24 - hx*mu/24
	vout[19, 15] = hz*lamb/12 - hz*mu/12
	vout[19, 16] = hx*hy*mu/(18*hz) - hx*hz*lamb/(9*hy) - 2*hx*hz*mu/(9*hy) + hy*hz*mu/(18*hx)
	vout[19, 17] = hx*lamb/12 - hx*mu/12
	vout[19, 18] = hz*lamb/12 + hz*mu/12
	vout[19, 19] = hx*hy*mu/(9*hz) + hx*hz*lamb/(9*hy) + 2*hx*hz*mu/(9*hy) + hy*hz*mu/(9*hx)
	vout[19, 20] = hx*lamb/12 + hx*mu/12
	vout[19, 21] = -hz*lamb/12 + hz*mu/12
	vout[19, 22] = hx*hy*mu/(18*hz) + hx*hz*lamb/(18*hy) + hx*hz*mu/(9*hy) - hy*hz*mu/(9*hx)
	vout[19, 23] = hx*lamb/24 + hx*mu/24
	vout[20, 0] = -hy*lamb/24 - hy*mu/24
	vout[20, 1] = -hx*lamb/24 - hx*mu/24
	vout[20, 2] = -hx*hy*lamb/(36*hz) - hx*hy*mu/(18*hz) - hx*hz*mu/(36*hy) - hy*hz*mu/(36*hx)
	vout[20, 3] = hy*lamb/24 - hy*mu/24
	vout[20, 4] = -hx*lamb/12 - hx*mu/12
	vout[20, 5] = -hx*hy*lamb/(18*hz) - hx*hy*mu/(9*hz) - hx*hz*mu/(18*hy) + hy*hz*mu/(36*hx)
	vout[20, 6] = hy*lamb/12 - hy*mu/12
	vout[20, 7] = hx*lamb/12 - hx*mu/12
	vout[20, 8] = -hx*hy*lamb/(9*hz) - 2*hx*hy*mu/(9*hz) + hx*hz*mu/(18*hy) + hy*hz*mu/(18*hx)
	vout[20, 9] = -hy*lamb/12 - hy*mu/12
	vout[20, 10] = hx*lamb/24 - hx*mu/24
	vout[20, 11] = -hx*hy*lamb/(18*hz) - hx*hy*mu/(9*hz) + hx*hz*mu/(36*hy) - hy*hz*mu/(18*hx)
	vout[20, 12] = -hy*lamb/24 + hy*mu/24
	vout[20, 13] = -hx*lamb/24 + hx*mu/24
	vout[20, 14] = hx*hy*lamb/(36*hz) + hx*hy*mu/(18*hz) - hx*hz*mu/(18*hy) - hy*hz*mu/(18*hx)
	vout[20, 15] = hy*lamb/24 + hy*mu/24
	vout[20, 16] = -hx*lamb/12 + hx*mu/12
	vout[20, 17] = hx*hy*lamb/(18*hz) + hx*hy*mu/(9*hz) - hx*hz*mu/(9*hy) + hy*hz*mu/(18*hx)
	vout[20, 18] = hy*lamb/12 + hy*mu/12
	vout[20, 19] = hx*lamb/12 + hx*mu/12
	vout[20, 20] = hx*hy*lamb/(9*hz) + 2*hx*hy*mu/(9*hz) + hx*hz*mu/(9*hy) + hy*hz*mu/(9*hx)
	vout[20, 21] = -hy*lamb/12 + hy*mu/12
	vout[20, 22] = hx*lamb/24 + hx*mu/24
	vout[20, 23] = hx*hy*lamb/(18*hz) + hx*hy*mu/(9*hz) + hx*hz*mu/(18*hy) - hy*hz*mu/(9*hx)
	vout[21, 0] = -hx*hy*mu/(18*hz) - hx*hz*mu/(18*hy) + hy*hz*lamb/(36*hx) + hy*hz*mu/(18*hx)
	vout[21, 1] = hz*lamb/24 - hz*mu/24
	vout[21, 2] = hy*lamb/24 - hy*mu/24
	vout[21, 3] = -hx*hy*mu/(36*hz) - hx*hz*mu/(36*hy) - hy*hz*lamb/(36*hx) - hy*hz*mu/(18*hx)
	vout[21, 4] = hz*lamb/24 + hz*mu/24
	vout[21, 5] = hy*lamb/24 + hy*mu/24
	vout[21, 6] = -hx*hy*mu/(18*hz) + hx*hz*mu/(36*hy) - hy*hz*lamb/(18*hx) - hy*hz*mu/(9*hx)
	vout[21, 7] = -hz*lamb/24 + hz*mu/24
	vout[21, 8] = hy*lamb/12 + hy*mu/12
	vout[21, 9] = -hx*hy*mu/(9*hz) + hx*hz*mu/(18*hy) + hy*hz*lamb/(18*hx) + hy*hz*mu/(9*hx)
	vout[21, 10] = -hz*lamb/24 - hz*mu/24
	vout[21, 11] = hy*lamb/12 - hy*mu/12
	vout[21, 12] = hx*hy*mu/(18*hz) - hx*hz*mu/(9*hy) + hy*hz*lamb/(18*hx) + hy*hz*mu/(9*hx)
	vout[21, 13] = hz*lamb/12 - hz*mu/12
	vout[21, 14] = -hy*lamb/24 - hy*mu/24
	vout[21, 15] = hx*hy*mu/(36*hz) - hx*hz*mu/(18*hy) - hy*hz*lamb/(18*hx) - hy*hz*mu/(9*hx)
	vout[21, 16] = hz*lamb/12 + hz*mu/12
	vout[21, 17] = -hy*lamb/24 + hy*mu/24
	vout[21, 18] = hx*hy*mu/(18*hz) + hx*hz*mu/(18*hy) - hy*hz*lamb/(9*hx) - 2*hy*hz*mu/(9*hx)
	vout[21, 19] = -hz*lamb/12 + hz*mu/12
	vout[21, 20] = -hy*lamb/12 + hy*mu/12
	vout[21, 21] = hx*hy*mu/(9*hz) + hx*hz*mu/(9*hy) + hy*hz*lamb/(9*hx) + 2*hy*hz*mu/(9*hx)
	vout[21, 22] = -hz*lamb/12 - hz*mu/12
	vout[21, 23] = -hy*lamb/12 - hy*mu/12
	vout[22, 0] = -hz*lamb/24 + hz*mu/24
	vout[22, 1] = -hx*hy*mu/(18*hz) - hx*hz*lamb/(18*hy) - hx*hz*mu/(9*hy) + hy*hz*mu/(36*hx)
	vout[22, 2] = -hx*lamb/12 - hx*mu/12
	vout[22, 3] = hz*lamb/24 + hz*mu/24
	vout[22, 4] = -hx*hy*mu/(36*hz) - hx*hz*lamb/(36*hy) - hx*hz*mu/(18*hy) - hy*hz*mu/(36*hx)
	vout[22, 5] = -hx*lamb/24 - hx*mu/24
	vout[22, 6] = hz*lamb/24 - hz*mu/24
	vout[22, 7] = -hx*hy*mu/(18*hz) + hx*hz*lamb/(36*hy) + hx*hz*mu/(18*hy) - hy*hz*mu/(18*hx)
	vout[22, 8] = -hx*lamb/24 + hx*mu/24
	vout[22, 9] = -hz*lamb/24 - hz*mu/24
	vout[22, 10] = -hx*hy*mu/(9*hz) + hx*hz*lamb/(18*hy) + hx*hz*mu/(9*hy) + hy*hz*mu/(18*hx)
	vout[22, 11] = -hx*lamb/12 + hx*mu/12
	vout[22, 12] = -hz*lamb/12 + hz*mu/12
	vout[22, 13] = hx*hy*mu/(18*hz) - hx*hz*lamb/(9*hy) - 2*hx*hz*mu/(9*hy) + hy*hz*mu/(18*hx)
	vout[22, 14] = hx*lamb/12 - hx*mu/12
	vout[22, 15] = hz*lamb/12 + hz*mu/12
	vout[22, 16] = hx*hy*mu/(36*hz) - hx*hz*lamb/(18*hy) - hx*hz*mu/(9*hy) - hy*hz*mu/(18*hx)
	vout[22, 17] = hx*lamb/24 - hx*mu/24
	vout[22, 18] = hz*lamb/12 - hz*mu/12
	vout[22, 19] = hx*hy*mu/(18*hz) + hx*hz*lamb/(18*hy) + hx*hz*mu/(9*hy) - hy*hz*mu/(9*hx)
	vout[22, 20] = hx*lamb/24 + hx*mu/24
	vout[22, 21] = -hz*lamb/12 - hz*mu/12
	vout[22, 22] = hx*hy*mu/(9*hz) + hx*hz*lamb/(9*hy) + 2*hx*hz*mu/(9*hy) + hy*hz*mu/(9*hx)
	vout[22, 23] = hx*lamb/12 + hx*mu/12
	vout[23, 0] = -hy*lamb/24 + hy*mu/24
	vout[23, 1] = -hx*lamb/12 - hx*mu/12
	vout[23, 2] = -hx*hy*lamb/(18*hz) - hx*hy*mu/(9*hz) - hx*hz*mu/(18*hy) + hy*hz*mu/(36*hx)
	vout[23, 3] = hy*lamb/24 + hy*mu/24
	vout[23, 4] = -hx*lamb/24 - hx*mu/24
	vout[23, 5] = -hx*hy*lamb/(36*hz) - hx*hy*mu/(18*hz) - hx*hz*mu/(36*hy) - hy*hz*mu/(36*hx)
	vout[23, 6] = hy*lamb/12 + hy*mu/12
	vout[23, 7] = hx*lamb/24 - hx*mu/24
	vout[23, 8] = -hx*hy*lamb/(18*hz) - hx*hy*mu/(9*hz) + hx*hz*mu/(36*hy) - hy*hz*mu/(18*hx)
	vout[23, 9] = -hy*lamb/12 + hy*mu/12
	vout[23, 10] = hx*lamb/12 - hx*mu/12
	vout[23, 11] = -hx*hy*lamb/(9*hz) - 2*hx*hy*mu/(9*hz) + hx*hz*mu/(18*hy) + hy*hz*mu/(18*hx)
	vout[23, 12] = -hy*lamb/24 - hy*mu/24
	vout[23, 13] = -hx*lamb/12 + hx*mu/12
	vout[23, 14] = hx*hy*lamb/(18*hz) + hx*hy*mu/(9*hz) - hx*hz*mu/(9*hy) + hy*hz*mu/(18*hx)
	vout[23, 15] = hy*lamb/24 - hy*mu/24
	vout[23, 16] = -hx*lamb/24 + hx*mu/24
	vout[23, 17] = hx*hy*lamb/(36*hz) + hx*hy*mu/(18*hz) - hx*hz*mu/(18*hy) - hy*hz*mu/(18*hx)
	vout[23, 18] = hy*lamb/12 - hy*mu/12
	vout[23, 19] = hx*lamb/24 + hx*mu/24
	vout[23, 20] = hx*hy*lamb/(18*hz) + hx*hy*mu/(9*hz) + hx*hz*mu/(18*hy) - hy*hz*mu/(9*hx)
	vout[23, 21] = -hy*lamb/12 - hy*mu/12
	vout[23, 22] = hx*lamb/12 + hx*mu/12
	vout[23, 23] = hx*hy*lamb/(9*hz) + 2*hx*hy*mu/(9*hz) + hx*hz*mu/(9*hy) + hy*hz*mu/(9*hx)
	return output