58   PetscScalar Rp, UU, VV;  in SolKxSolution()  local
88   Rp = sqrt((r + a) / 2.0);  in SolKxSolution()
90   UU = Rp - B;  in SolKxSolution()
91   VV = Rp + B;  in SolKxSolution()
122   t2   = B - Rp;  in SolKxSolution()
123   t4   = Rp + B;  in SolKxSolution()
125   t9   = exp(t6 - 0.4e1 * Rp);  in SolKxSolution()
128   t18  = Rp * Rp;  in SolKxSolution()
130   t20  = t15 * Rp;  in SolKxSolution()
131   t22  = t1 * Rp;  in SolKxSolution()
133   t32  = 0.8e1 * t15 * BB * kn * Rp;  in SolKxSolution()
136   t37  = Rm * Rp;  in SolKxSolution()
138   t63  = exp(t6 - 0.2e1 * Rp);  in SolKxSolution()
144   t75  = 0.2e1 * B * Rp;  in SolKxSolution()
155   t108 = exp(t6 - 0.3e1 * Rp - B);  in SolKxSolution()
162   t140 = exp(t6 - Rp - B);  in SolKxSolution()
163 …Rp * (-B * t13 + 0.3e1 * t15 * B - t19 - 0.2e1 * t20 - 0.2e1 * t22 - t24) * AA - t32) * t35 + (0.2…  in SolKxSolution()
165   t1   = Rp + B;  in SolKxSolution()
169   t5   = B * Rp;  in SolKxSolution()
172   t9   = Rp * Rp;  in SolKxSolution()
177   t23  = exp(t20 - 0.4e1 * Rp);  in SolKxSolution()
178   t25  = Rm * Rp;  in SolKxSolution()
183   t44  = t3 * Rp;  in SolKxSolution()
186   t54  = kn * Rp;  in SolKxSolution()
188   t69  = exp(t20 - 0.2e1 * Rp);  in SolKxSolution()
189   t71  = t9 * Rp;  in SolKxSolution()
192   t75  = 0.5e1 * t73 * Rp;  in SolKxSolution()
198   t82  = t79 * Rp;  in SolKxSolution()
206   t92  = t48 * Rp;  in SolKxSolution()
208   t96  = 0.3e1 * Rp * t41;  in SolKxSolution()
220   t132 = exp(t20 - 0.3e1 * Rp - B);  in SolKxSolution()
225   t162 = exp(t20 - Rp - B);  in SolKxSolution()
226 …) * AA + 0.8e1 * B * BB * t32 * Rp) * t38 + (0.2e1 * Rp * (-B * t8 + 0.3e1 * t41 - t9 * B - 0.2e1 …  in SolKxSolution()
229   t2   = B - Rp;  in SolKxSolution()
230   t4   = Rp + B;  in SolKxSolution()
234   t12  = Rp * Rp;  in SolKxSolution()
237   t15  = t14 * Rp;  in SolKxSolution()
239   t21  = t1 * Rp;  in SolKxSolution()
240   t30  = 0.8e1 * t14 * BB * kn * Rp;  in SolKxSolution()
243   t35  = Rm * Rp;  in SolKxSolution()
245   t61  = exp(-t6 - 0.2e1 * Rp);  in SolKxSolution()
251   t73  = 0.2e1 * B * Rp;  in SolKxSolution()
262   t106 = exp(-t6 - 0.3e1 * Rp - B);  in SolKxSolution()
269   t138 = exp(-t6 - Rp - B);  in SolKxSolution()
270 …Rp * (-t11 - t13 + 0.2e1 * t15 + 0.3e1 * t14 * B - B * t19 + 0.2e1 * t21) * AA + t30) * t33 + (-0.…  in SolKxSolution()
272   t1   = B - Rp;  in SolKxSolution()
276   t5   = B * Rp;  in SolKxSolution()
279   t9   = Rp * Rp;  in SolKxSolution()
285   t23  = Rm * Rp;  in SolKxSolution()
290   t40  = t3 * Rp;  in SolKxSolution()
293   t52  = kn * Rp;  in SolKxSolution()
295   t67  = exp(-t20 - 0.2e1 * Rp);  in SolKxSolution()
306   t104 = exp(-t20 - 0.3e1 * Rp - B);  in SolKxSolution()
307   t106 = t9 * Rp;  in SolKxSolution()
310   t110 = 0.5e1 * t108 * Rp;  in SolKxSolution()
316   t117 = t114 * Rp;  in SolKxSolution()
320   t123 = t38 * Rp;  in SolKxSolution()
322   t127 = 0.3e1 * Rp * t42;  in SolKxSolution()
332   t160 = exp(-t20 - Rp - B);  in SolKxSolution()
333 …AA - 0.8e1 * B * BB * t30 * Rp) * t36 + (-0.2e1 * Rp * (-t38 - t9 * B + 0.2e1 * t40 + 0.3e1 * t42 …  in SolKxSolution()
336   t2   = Rp * Rp;  in SolKxSolution()
340   t9   = exp(-0.4e1 * Rp);  in SolKxSolution()
342   t22  = exp(-0.2e1 * Rp);  in SolKxSolution()