200   PetscReal kn, km, x, z;  in SolKxSolution()  local
227   km = m * PETSC_PI; /* solution valid for km not zero -- should get trivial solution if km=0 */  in SolKxSolution()
231   a  = B * B + km * km;  in SolKxSolution()
232   b  = 2.0 * km * B;  in SolKxSolution()
252   t4  = km * km;  in SolKxSolution()
255   t8  = pow(km + kn, 0.2e1);  in SolKxSolution()
257   t16 = pow(km - kn, 0.2e1);  in SolKxSolution()
260   t2  = km * km;  in SolKxSolution()
263   t6  = km - kn;  in SolKxSolution()
264   t7  = km + kn;  in SolKxSolution()
508   u1  = -km * (_PC1 * t2 + _PC2 * t4 + _PC3 * t2 + _PC4 * t4 + t10 * AA * t13 + t10 * BB * t16);  in SolKxSolution()
522   t8   = km * km;  in SolKxSolution()
545 …- 0.4e1 * t83 * BB * kn + 0.3e1 * t8 * BB * kn - sigma + BB * t87 - 0.8e1 * t76 * AA) * t106) / km;  in SolKxSolution()
548   t4  = km * km;  in SolKxSolution()
577   t8   = km * km;  in SolKxSolution()
600 …t76 * B - 0.4e1 * t79 * BB * kn + t8 * BB * kn - sigma + BB * t83 - 0.4e1 * t72 * AA) * t101) / km;  in SolKxSolution()
604   t8  = km * km;  in SolKxSolution()
626 …84 + (-0.4e1 * AA * t71 * B - t8 * BB * kn - 0.4e1 * t74 * BB * kn - sigma + BB * t80) * t96) / km;  in SolKxSolution()
632   sum5 += u5 * PetscCosReal(km * z); /* pressure */  in SolKxSolution()
633   sum6 += u6 * PetscCosReal(km * z); /* zz total stress */  in SolKxSolution()
635   u1 *= PetscCosReal(km * z); /* x velocity */  in SolKxSolution()
637   u2 *= PetscSinReal(km * z); /* z velocity */  in SolKxSolution()
640   u3 *= PetscCosReal(km * z); /* xx total stress */  in SolKxSolution()
642   u4 *= PetscSinReal(km * z); /* zx stress */  in SolKxSolution()