From 83afb26d0bbe09df10fd3337336135d22d956814 Mon Sep 17 00:00:00 2001
From: Zizhou Huang <zizhouhuang@roblox.com>
Date: Tue, 21 Oct 2025 20:23:18 -0700
Subject: [PATCH 1/2] remove unreferenced code

---
 src/ipc/smooth_contact/distance/autogen.cpp | 958 --------------------
 1 file changed, 958 deletions(-)

diff --git a/src/ipc/smooth_contact/distance/autogen.cpp b/src/ipc/smooth_contact/distance/autogen.cpp
index 79a58bf39..8472e4023 100644
--- a/src/ipc/smooth_contact/distance/autogen.cpp
+++ b/src/ipc/smooth_contact/distance/autogen.cpp
@@ -1742,964 +1742,6 @@ void edge_edge_closest_point_hessian_b(
            + 4 * t16 * t27 * t551 * t94 - t269 * t426 - t550 - t551 * t559);
 }
 
-void face_normal_squared_norm_gradient(
-    double t0_x,
-    double t0_y,
-    double t0_z,
-    double t1_x,
-    double t1_y,
-    double t1_z,
-    double t2_x,
-    double t2_y,
-    double t2_z,
-    double grad[9])
-{
-    const auto t0 = -t2_y;
-    const auto t1 = t0 + t1_y;
-    const auto t2 = t0_x - t1_x;
-    const auto t3 = t0 + t0_y;
-    const auto t4 = -t2_x;
-    const auto t5 = t0_x + t4;
-    const auto t6 = t0_y - t1_y;
-    const auto t7 = t2 * t3 - t5 * t6;
-    const auto t8 = -t2_z;
-    const auto t9 = t1_z + t8;
-    const auto t10 = t0_z + t8;
-    const auto t11 = t0_z - t1_z;
-    const auto t12 = t10 * t2 - t11 * t5;
-    const auto t13 = t1_x + t4;
-    const auto t14 = t10 * t6 - t11 * t3;
-    const auto t15 = t10 * t2 - t11 * t5;
-    grad[0] = 2 * t1 * t7 + 2 * t12 * t9;
-    grad[1] = -2 * t13 * t7 + 2 * t14 * t9;
-    grad[2] = -2 * t1 * t14 - 2 * t13 * t15;
-    grad[3] = -2 * t10 * t15 - 2 * t3 * t7;
-    grad[4] = -2 * t10 * t14 + 2 * t5 * t7;
-    grad[5] = 2 * t12 * t5 + 2 * t14 * t3;
-    grad[6] = 2 * t11 * t12 + 2 * t6 * t7;
-    grad[7] = 2 * t11 * t14 - 2 * t2 * t7;
-    grad[8] = -2 * t14 * t6 - 2 * t15 * t2;
-}
-
-// hess is (81�1) flattened in column-major order
-void face_normal_squared_norm_hessian(
-    double t0_x,
-    double t0_y,
-    double t0_z,
-    double t1_x,
-    double t1_y,
-    double t1_z,
-    double t2_x,
-    double t2_y,
-    double t2_z,
-    double hess[81])
-{
-    const auto t0 = -t2_y;
-    const auto t1 = t0 + t1_y;
-    const auto t2 = std::pow(t1, 2);
-    const auto t3 = -t2_z;
-    const auto t4 = t1_z + t3;
-    const auto t5 = std::pow(t4, 2);
-    const auto t6 = -t2_x;
-    const auto t7 = t1_x + t6;
-    const auto t8 = 2 * t7;
-    const auto t9 = -t1 * t8;
-    const auto t10 = -t4 * t8;
-    const auto t11 = t0 + t0_y;
-    const auto t12 = t1 * t11;
-    const auto t13 = t0_z + t3;
-    const auto t14 = t13 * t4;
-    const auto t15 = -2 * t12 - 2 * t14;
-    const auto t16 = t0_x + t6;
-    const auto t17 = t0_y - t1_y;
-    const auto t18 = t16 * t17;
-    const auto t19 = -t18;
-    const auto t20 = t0_x - t1_x;
-    const auto t21 = t11 * t20;
-    const auto t22 = t19 + t21;
-    const auto t23 = 2 * t1 * t16 + 2 * t22;
-    const auto t24 = t0_z - t1_z;
-    const auto t25 = t16 * t24;
-    const auto t26 = -t25;
-    const auto t27 = t13 * t20;
-    const auto t28 = t26 + t27;
-    const auto t29 = 2 * t16 * t4 + 2 * t28;
-    const auto t30 = t1 * t17;
-    const auto t31 = t24 * t4;
-    const auto t32 = 2 * t30 + 2 * t31;
-    const auto t33 = -2 * t1 * t20 - 2 * t22;
-    const auto t34 = -2 * t20 * t4 - 2 * t28;
-    const auto t35 = std::pow(t7, 2);
-    const auto t36 = -2 * t1 * t4;
-    const auto t37 = t18 - t21;
-    const auto t38 = 2 * t11 * t7 + 2 * t37;
-    const auto t39 = t16 * t7;
-    const auto t40 = -2 * t14 - 2 * t39;
-    const auto t41 = t11 * t24;
-    const auto t42 = -t41;
-    const auto t43 = t13 * t17;
-    const auto t44 = t42 + t43;
-    const auto t45 = 2 * t11 * t4 + 2 * t44;
-    const auto t46 = -2 * t17 * t7 - 2 * t37;
-    const auto t47 = t20 * t7;
-    const auto t48 = 2 * t31 + 2 * t47;
-    const auto t49 = -2 * t17 * t4 - 2 * t44;
-    const auto t50 = t25 - t27;
-    const auto t51 = 2 * t13 * t7 + 2 * t50;
-    const auto t52 = t41 - t43;
-    const auto t53 = 2 * t1 * t13 + 2 * t52;
-    const auto t54 = -2 * t12 - 2 * t39;
-    const auto t55 = -2 * t24 * t7 - 2 * t50;
-    const auto t56 = -2 * t1 * t24 - 2 * t52;
-    const auto t57 = 2 * t30 + 2 * t47;
-    const auto t58 = std::pow(t11, 2);
-    const auto t59 = std::pow(t13, 2);
-    const auto t60 = 2 * t16;
-    const auto t61 = -t11 * t60;
-    const auto t62 = -t13 * t60;
-    const auto t63 = t11 * t17;
-    const auto t64 = t13 * t24;
-    const auto t65 = -2 * t63 - 2 * t64;
-    const auto t66 = 2 * t19 + 4 * t21;
-    const auto t67 = 2 * t26 + 4 * t27;
-    const auto t68 = std::pow(t16, 2);
-    const auto t69 = -2 * t11 * t13;
-    const auto t70 = 4 * t16 * t17 - 2 * t21;
-    const auto t71 = t16 * t20;
-    const auto t72 = -2 * t64 - 2 * t71;
-    const auto t73 = 2 * t42 + 4 * t43;
-    const auto t74 = 4 * t16 * t24 - 2 * t27;
-    const auto t75 = 4 * t11 * t24 - 2 * t43;
-    const auto t76 = -2 * t63 - 2 * t71;
-    const auto t77 = std::pow(t17, 2);
-    const auto t78 = std::pow(t24, 2);
-    const auto t79 = 2 * t20;
-    const auto t80 = -t17 * t79;
-    const auto t81 = -t24 * t79;
-    const auto t82 = std::pow(t20, 2);
-    const auto t83 = -2 * t17 * t24;
-    hess[0] = 2 * t2 + 2 * t5;
-    hess[1] = t9;
-    hess[2] = t10;
-    hess[3] = t15;
-    hess[4] = t23;
-    hess[5] = t29;
-    hess[6] = t32;
-    hess[7] = t33;
-    hess[8] = t34;
-    hess[9] = t9;
-    hess[10] = 2 * t35 + 2 * t5;
-    hess[11] = t36;
-    hess[12] = t38;
-    hess[13] = t40;
-    hess[14] = t45;
-    hess[15] = t46;
-    hess[16] = t48;
-    hess[17] = t49;
-    hess[18] = t10;
-    hess[19] = t36;
-    hess[20] = 2 * t2 + 2 * t35;
-    hess[21] = t51;
-    hess[22] = t53;
-    hess[23] = t54;
-    hess[24] = t55;
-    hess[25] = t56;
-    hess[26] = t57;
-    hess[27] = t15;
-    hess[28] = t38;
-    hess[29] = t51;
-    hess[30] = 2 * t58 + 2 * t59;
-    hess[31] = t61;
-    hess[32] = t62;
-    hess[33] = t65;
-    hess[34] = t66;
-    hess[35] = t67;
-    hess[36] = t23;
-    hess[37] = t40;
-    hess[38] = t53;
-    hess[39] = t61;
-    hess[40] = 2 * t59 + 2 * t68;
-    hess[41] = t69;
-    hess[42] = t70;
-    hess[43] = t72;
-    hess[44] = t73;
-    hess[45] = t29;
-    hess[46] = t45;
-    hess[47] = t54;
-    hess[48] = t62;
-    hess[49] = t69;
-    hess[50] = 2 * t58 + 2 * t68;
-    hess[51] = t74;
-    hess[52] = t75;
-    hess[53] = t76;
-    hess[54] = t32;
-    hess[55] = t46;
-    hess[56] = t55;
-    hess[57] = t65;
-    hess[58] = t70;
-    hess[59] = t74;
-    hess[60] = 2 * t77 + 2 * t78;
-    hess[61] = t80;
-    hess[62] = t81;
-    hess[63] = t33;
-    hess[64] = t48;
-    hess[65] = t56;
-    hess[66] = t66;
-    hess[67] = t72;
-    hess[68] = t75;
-    hess[69] = t80;
-    hess[70] = 2 * t78 + 2 * t82;
-    hess[71] = t83;
-    hess[72] = t34;
-    hess[73] = t49;
-    hess[74] = t57;
-    hess[75] = t67;
-    hess[76] = t73;
-    hess[77] = t76;
-    hess[78] = t81;
-    hess[79] = t83;
-    hess[80] = 2 * t77 + 2 * t82;
-}
-
-void face_term_aux_gradient(
-    double t0_x,
-    double t0_y,
-    double t0_z,
-    double t1_x,
-    double t1_y,
-    double t1_z,
-    double t2_x,
-    double t2_y,
-    double t2_z,
-    double p1_x,
-    double p1_y,
-    double p1_z,
-    double p2_x,
-    double p2_y,
-    double p2_z,
-    double grad[15])
-{
-    const auto t0 = p1_x - p2_x;
-    const auto t1 = p1_y - p2_y;
-    const auto t2 = p1_z - p2_z;
-    const auto t3 = std::pow(t0, 2) + std::pow(t1, 2) + std::pow(t2, 2);
-    const auto t4 = t0_x - t1_x;
-    const auto t5 = -t2_z;
-    const auto t6 = t0_z + t5;
-    const auto t7 = t4 * t6;
-    const auto t8 = -t2_x;
-    const auto t9 = t0_x + t8;
-    const auto t10 = t0_z - t1_z;
-    const auto t11 = t10 * t9;
-    const auto t12 = -t11 + t7;
-    const auto t13 = -t2_y;
-    const auto t14 = t0_y + t13;
-    const auto t15 = t14 * t4;
-    const auto t16 = t0_y - t1_y;
-    const auto t17 = t16 * t9;
-    const auto t18 = t15 - t17;
-    const auto t19 = -t10 * t14 + t16 * t6;
-    const auto t20 = std::pow(t18, 2) + std::pow(t19, 2);
-    const auto t21 = std::pow(t12, 2) + t20;
-    const auto t22 = std::pow(t21 * t3, -1.0 / 2.0);
-    const auto t23 = p2_y - t0_y;
-    const auto t24 = t1_z + t5;
-    const auto t25 = p2_z - t0_z;
-    const auto t26 = t13 + t1_y;
-    const auto t27 = p2_x - t0_x;
-    const auto t28 = t18 * t25 + t19 * t27;
-    const auto t29 = -t12 * t23 + t28;
-    const auto t30 = t29 / t21;
-    const auto t31 = t1_x + t8;
-    const auto t32 = t10 * t9;
-    const auto t33 = -t32 + t4 * t6;
-    const auto t34 = t20 + std::pow(t33, 2);
-    const auto t35 = std::pow(t3 * t34, -1.0 / 2.0);
-    const auto t36 = -t23 * t33 + t28;
-    const auto t37 = t36 / t34;
-    const auto t38 = 1.0 / t3;
-    const auto t39 = t35 * t36 * t38;
-    const auto t40 = t29 * t38;
-    grad[0] =
-        t22 * (-t19 - t23 * t24 + t25 * t26 - t30 * (t12 * t24 + t18 * t26));
-    grad[1] = t22
-        * (t24 * t27 - t25 * t31 - t30 * (-t18 * t31 + t19 * t24) - t32
-           + t4 * t6);
-    grad[2] = t22
-        * (-t15 + t17 + t23 * t31 - t26 * t27 + t30 * (t12 * t31 + t19 * t26));
-    grad[3] = t35 * (-t14 * t25 + t23 * t6 + t37 * (t12 * t6 + t14 * t18));
-    grad[4] = -t35 * (-t25 * t9 + t27 * t6 + t37 * (t18 * t9 - t19 * t6));
-    grad[5] = t35 * (t14 * t27 - t23 * t9 - t37 * (t14 * t19 + t33 * t9));
-    grad[6] = -t35 * (t10 * t23 - t16 * t25 + t37 * (t10 * t33 + t16 * t18));
-    grad[7] = t35 * (t10 * t27 - t25 * t4 + t37 * (-t10 * t19 + t18 * t4));
-    grad[8] = t35 * (-t16 * t27 + t23 * t4 + t37 * (t12 * t4 + t16 * t19));
-    grad[9] = -t0 * t39;
-    grad[10] = -t1 * t39;
-    grad[11] = -t2 * t39;
-    grad[12] = t22 * (t0 * t40 + t19);
-    grad[13] = t22 * (t1 * t40 + t11 - t7);
-    grad[14] = t22 * (t18 + t2 * t40);
-}
-
-// hess is (225�1) flattened in column-major order
-void face_term_aux_hessian(
-    double t0_x,
-    double t0_y,
-    double t0_z,
-    double t1_x,
-    double t1_y,
-    double t1_z,
-    double t2_x,
-    double t2_y,
-    double t2_z,
-    double p1_x,
-    double p1_y,
-    double p1_z,
-    double p2_x,
-    double p2_y,
-    double p2_z,
-    double hess[225])
-{
-    const auto t0 = -t2_y;
-    const auto t1 = t0 + t1_y;
-    const auto t2 = std::pow(t1, 2);
-    const auto t3 = -t2_z;
-    const auto t4 = t1_z + t3;
-    const auto t5 = std::pow(t4, 2);
-    const auto t6 = -t0_y;
-    const auto t7 = p2_y + t6;
-    const auto t8 = -t1_x;
-    const auto t9 = t0_x + t8;
-    const auto t10 = -t9;
-    const auto t11 = t0_z + t3;
-    const auto t12 = -t11;
-    const auto t13 = t10 * t12;
-    const auto t14 = -t2_x;
-    const auto t15 = t0_x + t14;
-    const auto t16 = -t15;
-    const auto t17 = -t1_z;
-    const auto t18 = t0_z + t17;
-    const auto t19 = -t18;
-    const auto t20 = t16 * t19;
-    const auto t21 = t13 - t20;
-    const auto t22 = t21 * t7;
-    const auto t23 = -t0_x;
-    const auto t24 = p2_x + t23;
-    const auto t25 = -t1_y;
-    const auto t26 = t0_y + t25;
-    const auto t27 = t11 * t26;
-    const auto t28 = t0 + t0_y;
-    const auto t29 = t18 * t28;
-    const auto t30 = -t29;
-    const auto t31 = t27 + t30;
-    const auto t32 = t24 * t31;
-    const auto t33 = -t0_z;
-    const auto t34 = p2_z + t33;
-    const auto t35 = t28 * t9;
-    const auto t36 = t15 * t26;
-    const auto t37 = -t36;
-    const auto t38 = t35 + t37;
-    const auto t39 = t34 * t38;
-    const auto t40 = t32 + t39;
-    const auto t41 = -t22 + t40;
-    const auto t42 = -t1 * t34 + t31 + t4 * t7;
-    const auto t43 = -t42;
-    const auto t44 = t1 * t38;
-    const auto t45 = t21 * t4 + t44;
-    const auto t46 = 2 * t45;
-    const auto t47 = std::pow(t31, 2) + std::pow(t38, 2);
-    const auto t48 = std::pow(t21, 2) + t47;
-    const auto t49 = 1.0 / t48;
-    const auto t50 = t41 * t49;
-    const auto t51 = -t28;
-    const auto t52 = -t26;
-    const auto t53 = t10 * t51 - t16 * t52;
-    const auto t54 = -t21;
-    const auto t55 = t1 * t53 - t4 * t54;
-    const auto t56 = t50 * t55;
-    const auto t57 = -p2_x;
-    const auto t58 = p1_x + t57;
-    const auto t59 = std::pow(t58, 2);
-    const auto t60 = -p2_y;
-    const auto t61 = p1_y + t60;
-    const auto t62 = std::pow(t61, 2);
-    const auto t63 = -p2_z;
-    const auto t64 = p1_z + t63;
-    const auto t65 = std::pow(t64, 2);
-    const auto t66 = t59 + t62 + t65;
-    const auto t67 = std::pow(t48 * t66, -1.0 / 2.0);
-    const auto t68 = t49 * t67;
-    const auto t69 = t14 + t1_x;
-    const auto t70 = -t31 * t4 + t38 * t69;
-    const auto t71 = -t70;
-    const auto t72 = 2 * t71;
-    const auto t73 = t11 * t9;
-    const auto t74 = t15 * t18;
-    const auto t75 = -t74;
-    const auto t76 = t73 + t75;
-    const auto t77 = t24 * t4 - t34 * t69 + t76;
-    const auto t78 = t41 * t69;
-    const auto t79 = t45 * t50;
-    const auto t80 = t1 * t78 - t43 * t71 - t45 * t77 + t71 * t79;
-    const auto t81 = t1 * t24;
-    const auto t82 = t69 * t7;
-    const auto t83 = t38 + t81 - t82;
-    const auto t84 = -t83;
-    const auto t85 = t45 * t84;
-    const auto t86 = t1 * t31 + t21 * t69;
-    const auto t87 = t79 * t86;
-    const auto t88 = 2 * t86;
-    const auto t89 = t11 * t7 - t28 * t34;
-    const auto t90 = t45 * t89;
-    const auto t91 = t11 * t21 + t28 * t38;
-    const auto t92 = t1 * t28;
-    const auto t93 = t11 * t4;
-    const auto t94 = t92 + t93;
-    const auto t95 = t50 * t91;
-    const auto t96 = t45 * t95;
-    const auto t97 = 2 * t91;
-    const auto t98 = -t11 * t31 + t15 * t38;
-    const auto t99 = std::pow(t48, -2);
-    const auto t100 = t11 * t24 - t15 * t34;
-    const auto t101 = -t100;
-    const auto t102 = t45 * t49;
-    const auto t103 = t43 * t49;
-    const auto t104 = t101 * t102 + t103 * t98 + t2_z - t41 * t45 * t98 * t99
-        + t50 * (t1 * t15 + t38) + t63;
-    const auto t105 = t15 * t21 + t28 * t31;
-    const auto t106 = t24 * t28;
-    const auto t107 = t15 * t7;
-    const auto t108 = t106 - t107;
-    const auto t109 = p2_y + t0;
-    const auto t110 = t102 * t108 + t103 * t105 - t105 * t41 * t45 * t99 + t109
-        + t50 * (t15 * t4 + t76);
-    const auto t111 = t18 * t21 + t26 * t38;
-    const auto t112 = t18 * t7 - t26 * t34;
-    const auto t113 = -t112;
-    const auto t114 = t1 * t26;
-    const auto t115 = t18 * t4;
-    const auto t116 =
-        -t111 * t41 * t45 * t49 + t111 * t43 + t113 * t45 + t41 * (t114 + t115);
-    const auto t117 = t18 * t31 - t38 * t9;
-    const auto t118 = t18 * t24;
-    const auto t119 = t34 * t9;
-    const auto t120 = t118 - t119;
-    const auto t121 = p2_z + t102 * t120 + t103 * t117 - t117 * t41 * t45 * t99
-        + t17 + t50 * (-t1 * t9 - t38);
-    const auto t122 = t21 * t9 + t26 * t31;
-    const auto t123 = 2 * t122;
-    const auto t124 = t41 * t99;
-    const auto t125 = t24 * t26;
-    const auto t126 = -t125 + t7 * t9;
-    const auto t127 = t102 * t126 + t122 * t124 * t45 - t122 * t43 * t49 + t1_y
-        + t50 * (-t4 * t9 - t76) + t60;
-    const auto t128 = t47 + std::pow(t76, 2);
-    const auto t129 = 1.0 / t128;
-    const auto t130 = t40 - t7 * t76;
-    const auto t131 = t129 * t130;
-    const auto t132 = t131 * (t4 * t76 + t44) + t42;
-    const auto t133 = std::pow(t128 * t66, -1.0 / 2.0);
-    const auto t134 = 1.0 / t66;
-    const auto t135 = t133 * t134;
-    const auto t136 = t132 * t135;
-    const auto t137 = t136 * t58;
-    const auto t138 = t134 * t61;
-    const auto t139 = t133 * t138;
-    const auto t140 = t132 * t139;
-    const auto t141 = t136 * t64;
-    const auto t142 = t134 * t79;
-    const auto t143 = t67 * (-t102 * t31 + t134 * t43 * t58 - t142 * t58);
-    const auto t144 = t17 + t2_z;
-    const auto t145 = t67 * (t102 * t21 + t134 * t43 * t61 - t138 * t79 + t144);
-    const auto t146 = t25 + t2_y;
-    const auto t147 =
-        t67 * (-t102 * t38 + t134 * t43 * t64 - t142 * t64 - t146);
-    const auto t148 = t12 * t52 - t19 * t51;
-    const auto t149 = t148 * t4 - t53 * t69;
-    const auto t150 = t149 * t50;
-    const auto t151 = std::pow(t69, 2);
-    const auto t152 = t71 * t84;
-    const auto t153 = t50 * t71;
-    const auto t154 = t153 * t86;
-    const auto t155 = t124 * t149;
-    const auto t156 = t49 * t71;
-    const auto t157 = -t35 + t36;
-    const auto t158 = t124 * t71;
-    const auto t159 = p2_z + t156 * t89 + t158 * t91 + t3 - t49 * t77 * t91
-        + t50 * (t157 + t28 * t69);
-    const auto t160 = 2 * t98;
-    const auto t161 = t15 * t69;
-    const auto t162 = t50 * t98;
-    const auto t163 = -t101 * t71 + t162 * t71 + t41 * (t161 + t93) - t77 * t98;
-    const auto t164 = t49 * t77;
-    const auto t165 = t105 * t164 - t105 * t41 * t71 * t99 + t108 * t156 + t2_x
-        + t50 * (t28 * t4 + t31) + t57;
-    const auto t166 = t111 * t164 - t111 * t41 * t71 * t99 + t113 * t156 + t1_z
-        + t50 * (-t157 - t26 * t69) + t63;
-    const auto t167 = t69 * t9;
-    const auto t168 =
-        -t117 * t41 * t49 * t71 + t117 * t77 + t120 * t71 + t41 * (t115 + t167);
-    const auto t169 = p2_x + t122 * t158 - t122 * t49 * t77 + t126 * t156
-        + t50 * (-t26 * t4 - t31) + t8;
-    const auto t170 = t131 * t70 + t77;
-    const auto t171 = t135 * t170;
-    const auto t172 = -t171 * t58;
-    const auto t173 = -t139 * t170;
-    const auto t174 = -t171 * t64;
-    const auto t175 = t134 * t153;
-    const auto t176 = t67 * (t134 * t58 * t77 - t144 - t156 * t31 - t175 * t58);
-    const auto t177 = t67 * (t134 * t61 * t77 - t138 * t153 + t156 * t21);
-    const auto t178 = t67 * (t134 * t64 * t77 - t156 * t38 - t175 * t64 - t69);
-    const auto t179 = -t1 * t148 + t54 * t69;
-    const auto t180 = t179 * t50;
-    const auto t181 = t49 * t86;
-    const auto t182 = t49 * t84;
-    const auto t183 = -t73 + t74;
-    const auto t184 = t124 * t86;
-    const auto t185 =
-        t109 + t181 * t89 + t182 * t91 + t184 * t91 - t50 * (t11 * t69 + t183);
-    const auto t186 = -t27 + t29;
-    const auto t187 = p2_x - t101 * t49 * t86 + t14 + t182 * t98 + t184 * t98
-        + t50 * (t1 * t11 + t186);
-    const auto t188 = t105 * t84;
-    const auto t189 = t161 + t92;
-    const auto t190 = 2 * t105;
-    const auto t191 = t105 * t50;
-    const auto t192 = t191 * t86;
-    const auto t193 = p2_y + t111 * t182 + t111 * t184 - t113 * t49 * t86 + t25
-        + t50 * (-t18 * t69 - t183);
-    const auto t194 = t117 * t182 + t117 * t184 - t120 * t49 * t86 + t1_x
-        + t50 * (-t1 * t18 - t186) + t57;
-    const auto t195 = t122 * t50;
-    const auto t196 =
-        t122 * t84 + t126 * t86 + t195 * t86 - t41 * (t114 + t167);
-    const auto t197 = -t131 * t86 - t157 + t81 - t82;
-    const auto t198 = t135 * t197;
-    const auto t199 = t198 * t58;
-    const auto t200 = t139 * t197;
-    const auto t201 = t198 * t64;
-    const auto t202 = t134 * t84;
-    const auto t203 = t134 * t58;
-    const auto t204 = t50 * t86;
-    const auto t205 = t67 * (t146 + t181 * t31 + t202 * t58 + t203 * t204);
-    const auto t206 = t133
-        * (t129 * t130 * t134 * t61 * t86 - t129 * t76 * t86 - t138 * t83 - t2_x
-           - t8);
-    const auto t207 = t134 * t64;
-    const auto t208 = t67 * (t181 * t38 + t202 * t64 + t204 * t207);
-    const auto t209 = t11 * t54 + t51 * t53;
-    const auto t210 = t209 * t50;
-    const auto t211 = t124 * t209;
-    const auto t212 = std::pow(t28, 2);
-    const auto t213 = std::pow(t11, 2);
-    const auto t214 = t89 * t98;
-    const auto t215 = t15 * t41;
-    const auto t216 = t95 * t98;
-    const auto t217 = t105 * t89;
-    const auto t218 = t105 * t95;
-    const auto t219 = t111 * t89;
-    const auto t220 = t26 * t28;
-    const auto t221 = t11 * t18;
-    const auto t222 = t220 + t221;
-    const auto t223 = 2 * t111;
-    const auto t224 = t111 * t95;
-    const auto t225 = t49 * t89;
-    const auto t226 = t124 * t91;
-    const auto t227 = t0_z + t117 * t225 + t117 * t226 - t120 * t49 * t91
-        + t50 * (2 * t35 + t37) + t63;
-    const auto t228 = t49 * t91;
-    const auto t229 = t0_y + t60;
-    const auto t230 =
-        t122 * t225 + t122 * t226 + t126 * t228 + t229 - t50 * (2 * t73 + t75);
-    const auto t231 = t89 + t95;
-    const auto t232 = t134 * t67;
-    const auto t233 = t231 * t232;
-    const auto t234 = -t233 * t58;
-    const auto t235 = t138 * t67;
-    const auto t236 = -t231 * t235;
-    const auto t237 = -t233 * t64;
-    const auto t238 = t134 * t89;
-    const auto t239 = t67 * (t203 * t95 + t228 * t31 + t238 * t58);
-    const auto t240 = t67 * (t11 + t138 * t89 + t138 * t95 - t21 * t228);
-    const auto t241 = t2_y + t6;
-    const auto t242 = t67 * (t207 * t95 + t228 * t38 + t238 * t64 + t241);
-    const auto t243 = t12 * t148 + t15 * t53;
-    const auto t244 = t243 * t50;
-    const auto t245 = t124 * t243;
-    const auto t246 = std::pow(t15, 2);
-    const auto t247 = -t101 * t105 + t105 * t162 - t108 * t98 + t11 * t28 * t41;
-    const auto t248 = t49 * t98;
-    const auto t249 = t101 * t111 * t49 - t111 * t41 * t98 * t99 + t113 * t248
-        + t34 + t50 * (2 * t15 * t26 - t35);
-    const auto t250 = 2 * t117;
-    const auto t251 = t15 * t9;
-    const auto t252 =
-        -t101 * t117 + t117 * t162 - t120 * t98 + t41 * (t221 + t251);
-    const auto t253 = t0_x - t101 * t122 * t49 + t122 * t124 * t98 + t126 * t248
-        + t50 * (2 * t27 + t30) + t57;
-    const auto t254 = t100 + t162;
-    const auto t255 = t232 * t254;
-    const auto t256 = t255 * t58;
-    const auto t257 = t235 * t254;
-    const auto t258 = t255 * t64;
-    const auto t259 =
-        t67 * (t101 * t134 * t58 - t11 - t162 * t203 - t248 * t31);
-    const auto t260 = t67 * (t101 * t138 - t138 * t162 + t21 * t248);
-    const auto t261 = t23 + t2_x;
-    const auto t262 =
-        t67 * (t101 * t134 * t64 - t162 * t207 - t248 * t38 - t261);
-    const auto t263 = t148 * t28 + t16 * t54;
-    const auto t264 = t263 * t50;
-    const auto t265 = t108 * t49;
-    const auto t266 = t105 * t49;
-    const auto t267 = -t105 * t111 * t41 * t99 + t111 * t265 + t113 * t266
-        + t229 + t50 * (2 * t15 * t18 - t73);
-    const auto t268 = -t105 * t117 * t41 * t99 + t117 * t265 + t120 * t266 + t24
-        + t50 * (2 * t18 * t28 - t27);
-    const auto t269 = t105 * t126;
-    const auto t270 = t220 + t251;
-    const auto t271 = t122 * t191;
-    const auto t272 = -t106 + t107 + t191;
-    const auto t273 = t232 * t272;
-    const auto t274 = t273 * t58;
-    const auto t275 = t235 * t272;
-    const auto t276 = t273 * t64;
-    const auto t277 =
-        t67 * (t108 * t134 * t58 - t191 * t203 - t241 - t266 * t31);
-    const auto t278 = t67 * (t108 * t138 - t138 * t191 + t21 * t266 + t261);
-    const auto t279 = t67 * (t108 * t134 * t64 - t191 * t207 - t266 * t38);
-    const auto t280 = t19 * t54 + t26 * t53;
-    const auto t281 = t280 * t50;
-    const auto t282 = std::pow(t26, 2);
-    const auto t283 = std::pow(t18, 2);
-    const auto t284 = t41 * t9;
-    const auto t285 = t111 * t50;
-    const auto t286 = -t111 * t120 - t113 * t117 + t117 * t285 + t26 * t284;
-    const auto t287 = t111 * t126;
-    const auto t288 = t122 * t285;
-    const auto t289 = t112 + t285;
-    const auto t290 = t232 * t289;
-    const auto t291 = t290 * t58;
-    const auto t292 = t235 * t289;
-    const auto t293 = t290 * t64;
-    const auto t294 = t111 * t49;
-    const auto t295 = t67 * (t113 * t134 * t58 - t203 * t285 - t294 * t31);
-    const auto t296 = t1_z + t33;
-    const auto t297 = t67 * (t113 * t138 - t138 * t285 + t21 * t294 + t296);
-    const auto t298 = t1_y + t6;
-    const auto t299 =
-        t67 * (t113 * t134 * t64 - t207 * t285 - t294 * t38 - t298);
-    const auto t300 = t10 * t53 + t148 * t18;
-    const auto t301 = t124 * t300;
-    const auto t302 = t300 * t50;
-    const auto t303 = std::pow(t9, 2);
-    const auto t304 = t117 * t126;
-    const auto t305 = t117 * t50;
-    const auto t306 = t122 * t305;
-    const auto t307 = -t118 + t119 + t305;
-    const auto t308 = t232 * t307;
-    const auto t309 = t308 * t58;
-    const auto t310 = t235 * t307;
-    const auto t311 = t308 * t64;
-    const auto t312 = t117 * t49;
-    const auto t313 =
-        t67 * (t120 * t134 * t58 - t203 * t305 - t296 - t31 * t312);
-    const auto t314 = t67 * (t120 * t138 - t138 * t305 + t21 * t312);
-    const auto t315 = t67 * (t120 * t134 * t64 - t207 * t305 - t312 * t38 - t9);
-    const auto t316 = t148 * t52 + t54 * t9;
-    const auto t317 = t316 * t50;
-    const auto t318 = -t125 + t195 + t7 * t9;
-    const auto t319 = t232 * t318;
-    const auto t320 = -t319 * t58;
-    const auto t321 = -t235 * t318;
-    const auto t322 = -t319 * t64;
-    const auto t323 = t126 * t134;
-    const auto t324 = t122 * t49;
-    const auto t325 = t67 * (t195 * t203 + t298 + t31 * t324 + t323 * t58);
-    const auto t326 = t67 * (t126 * t138 + t138 * t195 - t21 * t324 + t9);
-    const auto t327 = t67 * (t195 * t207 + t323 * t64 + t324 * t38);
-    const auto t328 = 3 * t134;
-    const auto t329 = t328 * t59;
-    const auto t330 = t232 * t41;
-    const auto t331 = std::pow(t66, -2);
-    const auto t332 = 3 * t41;
-    const auto t333 = t331 * t332 * t58 * t67;
-    const auto t334 = t333 * t61;
-    const auto t335 = t333 * t64;
-    const auto t336 = t31 * t58;
-    const auto t337 = t22 - t32 - t39;
-    const auto t338 = t329 * t41 + t337;
-    const auto t339 = t232 * (-t336 - t338);
-    const auto t340 = t138 * t332;
-    const auto t341 = t13 - t20 - t340;
-    const auto t342 = t232 * t58;
-    const auto t343 = t341 * t342;
-    const auto t344 = t328 * t41;
-    const auto t345 = t344 * t64;
-    const auto t346 = t345 + t38;
-    const auto t347 = -t342 * t346;
-    const auto t348 = t328 * t62;
-    const auto t349 = t331 * t332 * t61 * t64 * t67;
-    const auto t350 = t31 + t344 * t58;
-    const auto t351 = -t235 * t350;
-    const auto t352 = t21 * t61;
-    const auto t353 = -t348 * t41 + t41;
-    const auto t354 = -t235 * t346;
-    const auto t355 = t328 * t65;
-    const auto t356 = t232 * t64;
-    const auto t357 = -t350 * t356;
-    const auto t358 = t341 * t356;
-    const auto t359 = t38 * t64;
-    const auto t360 = t337 + t355 * t41;
-    const auto t361 = t232 * (-t359 - t360);
-    const auto t362 = t232 * (-t21 * t58 + t31 * t61 + t340 * t58);
-    const auto t363 = t232 * (t31 * t64 + t345 * t58 + t38 * t58);
-    const auto t364 = t232 * (-t21 * t64 + t340 * t64 + t38 * t61);
-    hess[0] = t68
-        * (-t41 * (t2 + t5) - t43 * t46 + std::pow(t45, 2) * t50 + t46 * t56);
-    hess[1] = t68 * (t56 * t72 + t80);
-    hess[2] = t68 * (t4 * t41 * t69 + t43 * t86 - t56 * t88 - t85 - t87);
-    hess[3] = t68 * (t41 * t94 + t43 * t91 - t56 * t97 - t90 - t96);
-    hess[4] = t67 * (-t104 + 2 * t41 * t55 * t98 * t99);
-    hess[5] = t67 * (2 * t105 * t41 * t55 * t99 - t110);
-    hess[6] = t68 * (2 * t111 * t41 * t49 * t55 - t116);
-    hess[7] = t67 * (2 * t117 * t41 * t55 * t99 - t121);
-    hess[8] = t67 * (-t123 * t124 * t55 - t127);
-    hess[9] = t137;
-    hess[10] = t140;
-    hess[11] = t141;
-    hess[12] = t143;
-    hess[13] = t145;
-    hess[14] = t147;
-    hess[15] = t68 * (t150 * t46 + t80);
-    hess[16] = t68
-        * (t150 * t72 - t41 * (t151 + t5) + t50 * std::pow(t71, 2) - t72 * t77);
-    hess[17] = t68 * (t1 * t4 * t41 - t150 * t88 - t152 - t154 + t77 * t86);
-    hess[18] = t67 * (-t155 * t97 - t159);
-    hess[19] = t68 * (t150 * t160 + t163);
-    hess[20] = t67 * (2 * t105 * t149 * t41 * t99 - t165);
-    hess[21] = t67 * (2 * t111 * t149 * t41 * t99 - t166);
-    hess[22] = t68 * (2 * t117 * t149 * t41 * t49 - t168);
-    hess[23] = t67 * (-t123 * t155 - t169);
-    hess[24] = t172;
-    hess[25] = t173;
-    hess[26] = t174;
-    hess[27] = t176;
-    hess[28] = t177;
-    hess[29] = t178;
-    hess[30] = t68 * (t180 * t46 + t4 * t78 + t43 * t86 - t85 - t87);
-    hess[31] = t68 * (t1 * t4 * t41 - t152 - t154 + t180 * t72 + t77 * t86);
-    hess[32] = t68
-        * (-t180 * t88 + t41 * t49 * std::pow(t86, 2) - t41 * (t151 + t2)
-           + 2 * t84 * t86);
-    hess[33] = t67 * (-t124 * t179 * t97 + t185);
-    hess[34] = t67 * (2 * t179 * t41 * t98 * t99 - t187);
-    hess[35] = t68 * (t108 * t86 + t180 * t190 - t188 + t189 * t41 - t192);
-    hess[36] = t67 * (2 * t111 * t179 * t41 * t99 - t193);
-    hess[37] = t67 * (2 * t117 * t179 * t41 * t99 - t194);
-    hess[38] = t68 * (-t123 * t180 + t196);
-    hess[39] = t199;
-    hess[40] = t200;
-    hess[41] = t201;
-    hess[42] = t205;
-    hess[43] = t206;
-    hess[44] = t208;
-    hess[45] = t68 * (t210 * t46 + t41 * t94 + t43 * t91 - t90 - t96);
-    hess[46] = t67 * (-t159 + 2 * t209 * t41 * t71 * t99);
-    hess[47] = t67 * (t185 - t211 * t88);
-    hess[48] = t68
-        * (-t210 * t97 - t41 * (t212 + t213) + t50 * std::pow(t91, 2)
-           + t89 * t97);
-    hess[49] = t68 * (t101 * t91 + t160 * t210 - t214 + t215 * t28 - t216);
-    hess[50] = t68 * (t108 * t91 + t11 * t215 + t190 * t210 - t217 - t218);
-    hess[51] = t68 * (t113 * t91 + t210 * t223 - t219 + t222 * t41 - t224);
-    hess[52] = t67 * (2 * t117 * t209 * t41 * t99 - t227);
-    hess[53] = t67 * (-t123 * t211 + t230);
-    hess[54] = t234;
-    hess[55] = t236;
-    hess[56] = t237;
-    hess[57] = t239;
-    hess[58] = t240;
-    hess[59] = t242;
-    hess[60] = t67 * (-t104 + 2 * t243 * t41 * t45 * t99);
-    hess[61] = t68 * (t163 + t244 * t72);
-    hess[62] = t67 * (-t187 - t245 * t88);
-    hess[63] = t68 * (t101 * t91 + t15 * t28 * t41 - t214 - t216 - t244 * t97);
-    hess[64] = t68
-        * (-t101 * t160 + 2 * t243 * t41 * t49 * t98
-           + t41 * t49 * std::pow(t98, 2) - t41 * (t213 + t246));
-    hess[65] = t68 * (t190 * t244 + t247);
-    hess[66] = t67 * (2 * t111 * t243 * t41 * t99 - t249);
-    hess[67] = t68 * (t244 * t250 + t252);
-    hess[68] = t67 * (-t123 * t245 - t253);
-    hess[69] = t256;
-    hess[70] = t257;
-    hess[71] = t258;
-    hess[72] = t259;
-    hess[73] = t260;
-    hess[74] = t262;
-    hess[75] = t67 * (-t110 + 2 * t263 * t41 * t45 * t99);
-    hess[76] = t67 * (-t165 + 2 * t263 * t41 * t71 * t99);
-    hess[77] = t68 * (t108 * t86 - t188 + t189 * t41 - t192 - t264 * t88);
-    hess[78] = t68 * (t108 * t91 + t11 * t15 * t41 - t217 - t218 - t264 * t97);
-    hess[79] = t68 * (t160 * t264 + t247);
-    hess[80] = t68
-        * (std::pow(t105, 2) * t41 * t49 + 2 * t105 * t263 * t41 * t49
-           - t108 * t190 - t41 * (t212 + t246));
-    hess[81] = t67 * (2 * t111 * t263 * t41 * t99 - t267);
-    hess[82] = t67 * (2 * t117 * t263 * t41 * t99 - t268);
-    hess[83] = t68 * (t108 * t122 - t123 * t264 - t269 + t270 * t41 - t271);
-    hess[84] = t274;
-    hess[85] = t275;
-    hess[86] = t276;
-    hess[87] = t277;
-    hess[88] = t278;
-    hess[89] = t279;
-    hess[90] = t68 * (-t116 + 2 * t280 * t41 * t45 * t49);
-    hess[91] = t67 * (-t166 + 2 * t280 * t41 * t71 * t99);
-    hess[92] = t67 * (-t124 * t280 * t88 - t193);
-    hess[93] = t68 * (t113 * t91 - t219 + t222 * t41 - t224 - t281 * t97);
-    hess[94] = t67 * (-t249 + 2 * t280 * t41 * t98 * t99);
-    hess[95] = t67 * (2 * t105 * t280 * t41 * t99 - t267);
-    hess[96] = t68
-        * (std::pow(t111, 2) * t41 * t49 + 2 * t111 * t280 * t41 * t49
-           - t113 * t223 - t41 * (t282 + t283));
-    hess[97] = t68 * (t250 * t281 + t286);
-    hess[98] = t68 * (t113 * t122 - t123 * t281 + t18 * t41 * t9 - t287 - t288);
-    hess[99] = t291;
-    hess[100] = t292;
-    hess[101] = t293;
-    hess[102] = t295;
-    hess[103] = t297;
-    hess[104] = t299;
-    hess[105] = t67 * (-t121 + 2 * t300 * t41 * t45 * t99);
-    hess[106] = t68 * (-t168 + 2 * t300 * t41 * t49 * t71);
-    hess[107] = t67 * (-t194 - t301 * t88);
-    hess[108] = t67 * (-t227 - t301 * t97);
-    hess[109] = t68 * (t160 * t302 + t252);
-    hess[110] = t67 * (2 * t105 * t300 * t41 * t99 - t268);
-    hess[111] = t68 * (t223 * t302 + t286);
-    hess[112] = t68
-        * (std::pow(t117, 2) * t41 * t49 + 2 * t117 * t300 * t41 * t49
-           - t120 * t250 - t41 * (t283 + t303));
-    hess[113] =
-        t68 * (t120 * t122 - t123 * t302 + t18 * t26 * t41 - t304 - t306);
-    hess[114] = t309;
-    hess[115] = t310;
-    hess[116] = t311;
-    hess[117] = t313;
-    hess[118] = t314;
-    hess[119] = t315;
-    hess[120] = t67 * (-t127 + 2 * t316 * t41 * t45 * t99);
-    hess[121] = t67 * (-t169 + 2 * t316 * t41 * t71 * t99);
-    hess[122] = t68 * (t196 - t317 * t88);
-    hess[123] = t67 * (-t124 * t316 * t97 + t230);
-    hess[124] = t67 * (-t253 + 2 * t316 * t41 * t98 * t99);
-    hess[125] = t68 * (t108 * t122 + t190 * t317 - t269 + t270 * t41 - t271);
-    hess[126] = t68 * (t113 * t122 + t18 * t284 + t223 * t317 - t287 - t288);
-    hess[127] =
-        t68 * (t120 * t122 + t18 * t26 * t41 + t250 * t317 - t304 - t306);
-    hess[128] = t68
-        * (std::pow(t122, 2) * t50 + t123 * t126 - t123 * t317
-           - t41 * (t282 + t303));
-    hess[129] = t320;
-    hess[130] = t321;
-    hess[131] = t322;
-    hess[132] = t325;
-    hess[133] = t326;
-    hess[134] = t327;
-    hess[135] = t137;
-    hess[136] = t172;
-    hess[137] = t199;
-    hess[138] = t234;
-    hess[139] = t256;
-    hess[140] = t274;
-    hess[141] = t291;
-    hess[142] = t309;
-    hess[143] = t320;
-    hess[144] = t330 * (t329 - 1);
-    hess[145] = t334;
-    hess[146] = t335;
-    hess[147] = t339;
-    hess[148] = t343;
-    hess[149] = t347;
-    hess[150] = t140;
-    hess[151] = t173;
-    hess[152] = t200;
-    hess[153] = t236;
-    hess[154] = t257;
-    hess[155] = t275;
-    hess[156] = t292;
-    hess[157] = t310;
-    hess[158] = t321;
-    hess[159] = t334;
-    hess[160] = t330 * (t348 - 1);
-    hess[161] = t349;
-    hess[162] = t351;
-    hess[163] = t232 * (t352 + t353);
-    hess[164] = t354;
-    hess[165] = t141;
-    hess[166] = t174;
-    hess[167] = t201;
-    hess[168] = t237;
-    hess[169] = t258;
-    hess[170] = t276;
-    hess[171] = t293;
-    hess[172] = t311;
-    hess[173] = t322;
-    hess[174] = t335;
-    hess[175] = t349;
-    hess[176] = t330 * (t355 - 1);
-    hess[177] = t357;
-    hess[178] = t358;
-    hess[179] = t361;
-    hess[180] = t143;
-    hess[181] = t176;
-    hess[182] = t205;
-    hess[183] = t239;
-    hess[184] = t259;
-    hess[185] = t277;
-    hess[186] = t295;
-    hess[187] = t313;
-    hess[188] = t325;
-    hess[189] = t339;
-    hess[190] = t351;
-    hess[191] = t357;
-    hess[192] = t232 * (2 * t336 + t338);
-    hess[193] = t362;
-    hess[194] = t363;
-    hess[195] = t145;
-    hess[196] = t177;
-    hess[197] = t206;
-    hess[198] = t240;
-    hess[199] = t260;
-    hess[200] = t278;
-    hess[201] = t297;
-    hess[202] = t314;
-    hess[203] = t326;
-    hess[204] = t343;
-    hess[205] = t135 * (-t130 * t348 + t130 + t61 * t76);
-    hess[206] = t358;
-    hess[207] = t362;
-    hess[208] = t232 * (-2 * t352 - t353);
-    hess[209] = t364;
-    hess[210] = t147;
-    hess[211] = t178;
-    hess[212] = t208;
-    hess[213] = t242;
-    hess[214] = t262;
-    hess[215] = t279;
-    hess[216] = t299;
-    hess[217] = t315;
-    hess[218] = t327;
-    hess[219] = t347;
-    hess[220] = t354;
-    hess[221] = t361;
-    hess[222] = t363;
-    hess[223] = t364;
-    hess[224] = t232 * (2 * t359 + t360);
-}
-
 // hess is (144×1) flattened in column-major order
 void triangle_closest_point_hessian_0(
     double p_x,

From 9d0903ae5be075c086440c14ecd3233628a8a42c Mon Sep 17 00:00:00 2001
From: Zizhou Huang <zizhouhuang@roblox.com>
Date: Tue, 21 Oct 2025 20:47:11 -0700
Subject: [PATCH 2/2] remove unreferenced code

---
 src/ipc/smooth_contact/distance/autogen.cpp | 590 --------------------
 1 file changed, 590 deletions(-)

diff --git a/src/ipc/smooth_contact/distance/autogen.cpp b/src/ipc/smooth_contact/distance/autogen.cpp
index 8472e4023..85f1b6237 100644
--- a/src/ipc/smooth_contact/distance/autogen.cpp
+++ b/src/ipc/smooth_contact/distance/autogen.cpp
@@ -3592,594 +3592,4 @@ void triangle_closest_point_hessian_1(
            + 4 * t238 * t28 * t580 * t76 - t239 * t579 + t248 * t41 * t579
            - t253 * t255 * t580);
 }
-
-void face_term_aux_fast_gradient(
-    double t0_x,
-    double t0_y,
-    double t0_z,
-    double t1_x,
-    double t1_y,
-    double t1_z,
-    double t2_x,
-    double t2_y,
-    double t2_z,
-    double p_x,
-    double p_y,
-    double p_z,
-    double d,
-    double grad[13])
-{
-    const auto t0 = t0_x - t1_x;
-    const auto t1 = -t2_z;
-    const auto t2 = t0_z + t1;
-    const auto t3 = -t2_x;
-    const auto t4 = t0_x + t3;
-    const auto t5 = t0_z - t1_z;
-    const auto t6 = t0 * t2 - t4 * t5;
-    const auto t7 = -t2_y;
-    const auto t8 = t0_y + t7;
-    const auto t9 = t0 * t8;
-    const auto t10 = t0_y - t1_y;
-    const auto t11 = t10 * t4;
-    const auto t12 = -t11 + t9;
-    const auto t13 = t10 * t2 - t5 * t8;
-    const auto t14 = std::pow(t12, 2) + std::pow(t13, 2);
-    const auto t15 = t14 + std::pow(t6, 2);
-    const auto t16 = std::pow(d * t15, -1.0 / 2.0);
-    const auto t17 = p_y - t0_y;
-    const auto t18 = t1 + t1_z;
-    const auto t19 = p_z - t0_z;
-    const auto t20 = t1_y + t7;
-    const auto t21 = p_x - t0_x;
-    const auto t22 = t12 * t19 + t13 * t21;
-    const auto t23 = (-t17 * t6 + t22) / t15;
-    const auto t24 = t1_x + t3;
-    const auto t25 = t4 * t5;
-    const auto t26 = t0 * t2 - t25;
-    const auto t27 = t14 + std::pow(t26, 2);
-    const auto t28 = std::pow(d * t27, -1.0 / 2.0);
-    const auto t29 = (-t17 * t26 + t22) / t27;
-    grad[0] =
-        t16 * (-t13 - t17 * t18 + t19 * t20 - t23 * (t12 * t20 + t18 * t6));
-    grad[1] = t16
-        * (t0 * t2 + t18 * t21 - t19 * t24 - t23 * (-t12 * t24 + t13 * t18)
-           - t25);
-    grad[2] =
-        t16 * (t11 + t17 * t24 - t20 * t21 + t23 * (t13 * t20 + t24 * t6) - t9);
-    grad[3] = t28 * (t17 * t2 - t19 * t8 + t29 * (t12 * t8 + t2 * t6));
-    grad[4] = -t28 * (-t19 * t4 + t2 * t21 + t29 * (t12 * t4 - t13 * t2));
-    grad[5] = t28 * (-t17 * t4 + t21 * t8 - t29 * (t13 * t8 + t26 * t4));
-    grad[6] = -t28 * (-t10 * t19 + t17 * t5 + t29 * (t10 * t12 + t26 * t5));
-    grad[7] = t28 * (-t0 * t19 + t21 * t5 + t29 * (t0 * t12 - t13 * t5));
-    grad[8] = t28 * (t0 * t17 - t10 * t21 + t29 * (t0 * t6 + t10 * t13));
-    grad[9] = t13 * t28;
-    grad[10] = -t26 * t28;
-    grad[11] = t12 * t28;
-    grad[12] = -1.0 / 2.0 * t15 * t28 * t29 / d;
-}
-
-// hess is (169×1) flattened in column-major order
-void face_term_aux_fast_hessian(
-    double t0_x,
-    double t0_y,
-    double t0_z,
-    double t1_x,
-    double t1_y,
-    double t1_z,
-    double t2_x,
-    double t2_y,
-    double t2_z,
-    double p_x,
-    double p_y,
-    double p_z,
-    double d,
-    double hess[169])
-{
-    const auto t0 = -t2_y;
-    const auto t1 = t0 + t1_y;
-    const auto t2 = std::pow(t1, 2);
-    const auto t3 = -t2_z;
-    const auto t4 = t1_z + t3;
-    const auto t5 = std::pow(t4, 2);
-    const auto t6 = -t0_y;
-    const auto t7 = p_y + t6;
-    const auto t8 = -t1_x;
-    const auto t9 = t0_x + t8;
-    const auto t10 = -t9;
-    const auto t11 = t0_z + t3;
-    const auto t12 = -t11;
-    const auto t13 = -t2_x;
-    const auto t14 = t0_x + t13;
-    const auto t15 = -t14;
-    const auto t16 = -t1_z;
-    const auto t17 = t0_z + t16;
-    const auto t18 = -t17;
-    const auto t19 = t10 * t12 - t15 * t18;
-    const auto t20 = -t0_x;
-    const auto t21 = p_x + t20;
-    const auto t22 = -t1_y;
-    const auto t23 = t0_y + t22;
-    const auto t24 = t11 * t23;
-    const auto t25 = t0 + t0_y;
-    const auto t26 = t17 * t25;
-    const auto t27 = -t26;
-    const auto t28 = t24 + t27;
-    const auto t29 = -t0_z;
-    const auto t30 = p_z + t29;
-    const auto t31 = t25 * t9;
-    const auto t32 = t14 * t23;
-    const auto t33 = -t32;
-    const auto t34 = t31 + t33;
-    const auto t35 = t21 * t28 + t30 * t34;
-    const auto t36 = -t19 * t7 + t35;
-    const auto t37 = -t1 * t30 + t28 + t4 * t7;
-    const auto t38 = -t37;
-    const auto t39 = t1 * t34;
-    const auto t40 = t19 * t4 + t39;
-    const auto t41 = 2 * t40;
-    const auto t42 = std::pow(t28, 2) + std::pow(t34, 2);
-    const auto t43 = std::pow(t19, 2) + t42;
-    const auto t44 = 1.0 / t43;
-    const auto t45 = t36 * t44;
-    const auto t46 = -t25;
-    const auto t47 = -t23;
-    const auto t48 = t10 * t46 - t15 * t47;
-    const auto t49 = -t19;
-    const auto t50 = t1 * t48 - t4 * t49;
-    const auto t51 = t36 * t50;
-    const auto t52 = t44 * t51;
-    const auto t53 = std::pow(d * t43, -1.0 / 2.0);
-    const auto t54 = t44 * t53;
-    const auto t55 = t13 + t1_x;
-    const auto t56 = -t28 * t4 + t34 * t55;
-    const auto t57 = -t56;
-    const auto t58 = 2 * t57;
-    const auto t59 = t21 * t4;
-    const auto t60 = t30 * t55;
-    const auto t61 = t11 * t9;
-    const auto t62 = t14 * t17;
-    const auto t63 = -t62;
-    const auto t64 = t61 + t63;
-    const auto t65 = t59 - t60 + t64;
-    const auto t66 = t36 * t55;
-    const auto t67 = t40 * t45;
-    const auto t68 = t1 * t66 - t38 * t57 - t40 * t65 + t57 * t67;
-    const auto t69 = t1 * t21;
-    const auto t70 = -t34 + t55 * t7 - t69;
-    const auto t71 = t40 * t70;
-    const auto t72 = t1 * t28 + t19 * t55;
-    const auto t73 = t67 * t72;
-    const auto t74 = 2 * t72;
-    const auto t75 = t11 * t7 - t25 * t30;
-    const auto t76 = t40 * t75;
-    const auto t77 = t11 * t19 + t25 * t34;
-    const auto t78 = t1 * t25;
-    const auto t79 = t11 * t4;
-    const auto t80 = t78 + t79;
-    const auto t81 = t67 * t77;
-    const auto t82 = 2 * t77;
-    const auto t83 = -t11 * t28 + t14 * t34;
-    const auto t84 = std::pow(t43, -2);
-    const auto t85 = -p_z;
-    const auto t86 = t11 * t21 - t14 * t30;
-    const auto t87 = -t86;
-    const auto t88 = t40 * t44;
-    const auto t89 = t38 * t44;
-    const auto t90 = t2_z - t36 * t40 * t83 * t84 + t45 * (t1 * t14 + t34)
-        + t83 * t89 + t85 + t87 * t88;
-    const auto t91 = t25 * t28;
-    const auto t92 = t14 * t19 + t91;
-    const auto t93 = t21 * t25;
-    const auto t94 = t14 * t7;
-    const auto t95 = t93 - t94;
-    const auto t96 = p_y + t0;
-    const auto t97 = -t36 * t40 * t84 * t92 + t45 * (t14 * t4 + t64) + t88 * t95
-        + t89 * t92 + t96;
-    const auto t98 = t23 * t34;
-    const auto t99 = t17 * t19 + t98;
-    const auto t100 = t17 * t7 - t23 * t30;
-    const auto t101 = -t100;
-    const auto t102 = t1 * t23;
-    const auto t103 = t17 * t4;
-    const auto t104 =
-        t101 * t40 - t36 * t40 * t44 * t99 + t36 * (t102 + t103) + t38 * t99;
-    const auto t105 = -t17 * t28 + t34 * t9;
-    const auto t106 = -t105;
-    const auto t107 = t17 * t21;
-    const auto t108 = t30 * t9;
-    const auto t109 = t107 - t108;
-    const auto t110 = p_z - t106 * t36 * t40 * t84 + t106 * t89 + t109 * t88
-        + t16 + t45 * (-t1 * t9 - t34);
-    const auto t111 = t19 * t9 + t23 * t28;
-    const auto t112 = 2 * t111;
-    const auto t113 = -p_y;
-    const auto t114 = t21 * t23;
-    const auto t115 = -t114 + t7 * t9;
-    const auto t116 = t36 * t84;
-    const auto t117 = t111 * t116 * t40 - t111 * t38 * t44 + t113 + t115 * t88
-        + t1_y + t45 * (-t4 * t9 - t64);
-    const auto t118 = t39 + t4 * t64;
-    const auto t119 = t42 + std::pow(t64, 2);
-    const auto t120 = std::pow(d * t119, -1.0 / 2.0);
-    const auto t121 = 1.0 / t119;
-    const auto t122 = t121 * t28;
-    const auto t123 = t120 * t122;
-    const auto t124 = -t118 * t123;
-    const auto t125 = t53 * (t16 + t19 * t88 + t2_z);
-    const auto t126 = t121 * t34;
-    const auto t127 = t22 + t2_y;
-    const auto t128 = t120 * (-t118 * t126 - t127);
-    const auto t129 = t12 * t47 - t18 * t46;
-    const auto t130 = std::pow(t129, 2) + std::pow(t48, 2) + std::pow(t49, 2);
-    const auto t131 = (1.0 / 2.0) * t130;
-    const auto t132 = t130 * t44;
-    const auto t133 = 1.0 / d;
-    const auto t134 = t133 * t54;
-    const auto t135 = t129 * t4 - t48 * t55;
-    const auto t136 = t135 * t36;
-    const auto t137 = t136 * t44;
-    const auto t138 = std::pow(t55, 2);
-    const auto t139 = t57 * t70;
-    const auto t140 = t45 * t57;
-    const auto t141 = t140 * t72;
-    const auto t142 = t136 * t84;
-    const auto t143 = t44 * t57;
-    const auto t144 = -t31 + t32;
-    const auto t145 = t116 * t57;
-    const auto t146 = p_z + t143 * t75 + t145 * t77 + t3 - t44 * t65 * t77
-        + t45 * (t144 + t25 * t55);
-    const auto t147 = 2 * t83;
-    const auto t148 = t14 * t55;
-    const auto t149 = t140 * t83 + t36 * (t148 + t79) - t57 * t87 - t65 * t83;
-    const auto t150 = -p_x;
-    const auto t151 = t44 * t65;
-    const auto t152 = t143 * t95 + t150 + t151 * t92 + t2_x
-        - t36 * t57 * t84 * t92 + t45 * (t25 * t4 + t28);
-    const auto t153 = t101 * t143 + t151 * t99 + t1_z - t36 * t57 * t84 * t99
-        + t45 * (-t144 - t23 * t55) + t85;
-    const auto t154 = t55 * t9;
-    const auto t155 =
-        -t106 * t36 * t44 * t57 + t106 * t65 + t109 * t57 + t36 * (t103 + t154);
-    const auto t156 = p_x + t111 * t145 - t111 * t44 * t65 + t115 * t143
-        + t45 * (-t23 * t4 - t28) + t8;
-    const auto t157 = t120 * (t122 * t56 + t4);
-    const auto t158 = t121 * t64;
-    const auto t159 = t120 * t158;
-    const auto t160 = -t159 * t56;
-    const auto t161 = -t53 * (t143 * t34 + t55);
-    const auto t162 = -t1 * t129 + t49 * t55;
-    const auto t163 = t162 * t36;
-    const auto t164 = t163 * t44;
-    const auto t165 = t44 * t72;
-    const auto t166 = t44 * t70;
-    const auto t167 = -t61 + t62;
-    const auto t168 = t116 * t72;
-    const auto t169 =
-        t165 * t75 + t166 * t77 + t168 * t77 - t45 * (t11 * t55 + t167) + t96;
-    const auto t170 = -t24 + t26;
-    const auto t171 = p_x + t13 + t166 * t83 + t168 * t83 - t44 * t72 * t87
-        + t45 * (t1 * t11 + t170);
-    const auto t172 = t70 * t92;
-    const auto t173 = t148 + t78;
-    const auto t174 = 2 * t92;
-    const auto t175 = t45 * t72;
-    const auto t176 = t175 * t92;
-    const auto t177 = p_y - t101 * t44 * t72 + t166 * t99 + t168 * t99 + t22
-        + t45 * (-t167 - t17 * t55);
-    const auto t178 = t106 * t166 + t106 * t168 - t109 * t44 * t72 + t150 + t1_x
-        + t45 * (-t1 * t17 - t170);
-    const auto t179 =
-        t111 * t175 + t111 * t70 + t115 * t72 - t36 * (t102 + t154);
-    const auto t180 = t53 * (t127 + t165 * t28);
-    const auto t181 = t120 * (-t158 * t72 - t2_x - t8);
-    const auto t182 = t120 * t126;
-    const auto t183 = t182 * t72;
-    const auto t184 = t11 * t49 + t46 * t48;
-    const auto t185 = t184 * t36;
-    const auto t186 = t185 * t44;
-    const auto t187 = t185 * t84;
-    const auto t188 = std::pow(t25, 2);
-    const auto t189 = std::pow(t11, 2);
-    const auto t190 = t75 * t83;
-    const auto t191 = t14 * t36;
-    const auto t192 = t45 * t77;
-    const auto t193 = t192 * t83;
-    const auto t194 = t75 * t92;
-    const auto t195 = t192 * t92;
-    const auto t196 = t75 * t99;
-    const auto t197 = t23 * t25;
-    const auto t198 = t11 * t17;
-    const auto t199 = t197 + t198;
-    const auto t200 = 2 * t99;
-    const auto t201 = t192 * t99;
-    const auto t202 = t44 * t75;
-    const auto t203 = t116 * t77;
-    const auto t204 = t0_z + t106 * t202 + t106 * t203 - t109 * t44 * t77
-        + t45 * (2 * t31 + t33) + t85;
-    const auto t205 = t44 * t77;
-    const auto t206 = t0_y + t113;
-    const auto t207 =
-        t111 * t202 + t111 * t203 + t115 * t205 + t206 - t45 * (2 * t61 + t63);
-    const auto t208 = t123 * t77;
-    const auto t209 = t120 * (-t158 * t77 - t29 - t2_z);
-    const auto t210 = t2_y + t6;
-    const auto t211 = t53 * (t205 * t34 + t210);
-    const auto t212 = t12 * t129 + t14 * t48;
-    const auto t213 = t212 * t36;
-    const auto t214 = t213 * t44;
-    const auto t215 = t213 * t84;
-    const auto t216 = std::pow(t14, 2);
-    const auto t217 = t45 * t83;
-    const auto t218 = t11 * t25 * t36 + t217 * t92 - t83 * t95 - t87 * t92;
-    const auto t219 = t44 * t83;
-    const auto t220 = t101 * t219 + t30 - t36 * t83 * t84 * t99
-        + t44 * t87 * t99 + t45 * (2 * t14 * t23 - t31);
-    const auto t221 = 2 * t106;
-    const auto t222 = t14 * t9;
-    const auto t223 =
-        t106 * t217 - t106 * t87 - t109 * t83 + t36 * (t198 + t222);
-    const auto t224 = t0_x + t111 * t116 * t83 - t111 * t44 * t87 + t115 * t219
-        + t150 + t45 * (2 * t24 + t27);
-    const auto t225 = -t53 * (t11 + t219 * t28);
-    const auto t226 = t159 * t83;
-    const auto t227 = t20 + t2_x;
-    const auto t228 = t120 * (-t126 * t83 - t227);
-    const auto t229 = t129 * t25 + t15 * t49;
-    const auto t230 = t229 * t36;
-    const auto t231 = t230 * t44;
-    const auto t232 = t44 * t95;
-    const auto t233 = t44 * t92;
-    const auto t234 = t101 * t233 + t206 + t232 * t99 - t36 * t84 * t92 * t99
-        + t45 * (2 * t14 * t17 - t61);
-    const auto t235 = t106 * t232 - t106 * t36 * t84 * t92 + t109 * t233 + t21
-        + t45 * (2 * t17 * t25 - t24);
-    const auto t236 = t115 * t92;
-    const auto t237 = t197 + t222;
-    const auto t238 = t45 * t92;
-    const auto t239 = t111 * t238;
-    const auto t240 = t14 * t64 + t91;
-    const auto t241 = t120 * (-t122 * t240 - t210);
-    const auto t242 = t53 * (t19 * t233 + t227);
-    const auto t243 = -t182 * t240;
-    const auto t244 = t18 * t49 + t23 * t48;
-    const auto t245 = t244 * t36;
-    const auto t246 = t245 * t44;
-    const auto t247 = std::pow(t23, 2);
-    const auto t248 = std::pow(t17, 2);
-    const auto t249 = t36 * t9;
-    const auto t250 = t45 * t99;
-    const auto t251 = -t101 * t106 + t106 * t250 - t109 * t99 + t23 * t249;
-    const auto t252 = t115 * t99;
-    const auto t253 = t111 * t250;
-    const auto t254 = t17 * t64 + t98;
-    const auto t255 = -t123 * t254;
-    const auto t256 = t53 * (t19 * t44 * t99 + t1_z + t29);
-    const auto t257 = t1_y + t6;
-    const auto t258 = t120 * (-t126 * t254 - t257);
-    const auto t259 = t10 * t48 + t129 * t17;
-    const auto t260 = t259 * t36;
-    const auto t261 = t260 * t84;
-    const auto t262 = t260 * t44;
-    const auto t263 = std::pow(t9, 2);
-    const auto t264 = t106 * t115;
-    const auto t265 = t106 * t45;
-    const auto t266 = t111 * t265;
-    const auto t267 = t120 * (t105 * t122 + t17);
-    const auto t268 = -t105 * t159;
-    const auto t269 = -t53 * (t106 * t34 * t44 + t9);
-    const auto t270 = t129 * t47 + t49 * t9;
-    const auto t271 = t270 * t36;
-    const auto t272 = t271 * t44;
-    const auto t273 = t53 * (t111 * t28 * t44 + t257);
-    const auto t274 = t120 * (-t111 * t158 - t1_x - t20);
-    const auto t275 = t111 * t182;
-    const auto t276 = t111 * t45;
-    const auto t277 = t121 * t43;
-    const auto t278 = t120 * t277;
-    const auto t279 = (1.0 / 2.0) * t133 * t278;
-    const auto t280 = -t279 * t28;
-    const auto t281 = t279 * t64;
-    const auto t282 = -t279 * t34;
-    const auto t283 = t131 * t134;
-    hess[0] = t54
-        * (-t36 * (t2 + t5) - t38 * t41 + std::pow(t40, 2) * t45 + t41 * t52);
-    hess[1] = t54 * (t52 * t58 + t68);
-    hess[2] = t54 * (t36 * t4 * t55 + t38 * t72 - t52 * t74 - t71 - t73);
-    hess[3] = t54 * (t36 * t80 + t38 * t77 - t52 * t82 - t76 - t81);
-    hess[4] = t53 * (2 * t36 * t50 * t83 * t84 - t90);
-    hess[5] = t53 * (2 * t36 * t50 * t84 * t92 - t97);
-    hess[6] = t54 * (-t104 + 2 * t36 * t44 * t50 * t99);
-    hess[7] = t53 * (2 * t106 * t36 * t50 * t84 - t110);
-    hess[8] = t53 * (-t112 * t51 * t84 - t117);
-    hess[9] = t124;
-    hess[10] = t125;
-    hess[11] = t128;
-    hess[12] = t134 * (-t131 * t38 + t131 * t67 + t132 * t51 - t51);
-    hess[13] = t54 * (t137 * t41 + t68);
-    hess[14] = t54
-        * (t137 * t58 - t36 * (t138 + t5) + t45 * std::pow(t57, 2) - t58 * t65);
-    hess[15] = t54 * (t1 * t36 * t4 - t137 * t74 - t139 - t141 + t65 * t72);
-    hess[16] = t53 * (-t142 * t82 - t146);
-    hess[17] = t54 * (t137 * t147 + t149);
-    hess[18] = t53 * (2 * t135 * t36 * t84 * t92 - t152);
-    hess[19] = t53 * (2 * t135 * t36 * t84 * t99 - t153);
-    hess[20] = t54 * (2 * t106 * t135 * t36 * t44 - t155);
-    hess[21] = t53 * (-t112 * t142 - t156);
-    hess[22] = t157;
-    hess[23] = t160;
-    hess[24] = t161;
-    hess[25] = t134 * (t131 * t140 - t131 * t65 + t132 * t136 - t136);
-    hess[26] = t54 * (t164 * t41 + t38 * t72 + t4 * t66 - t71 - t73);
-    hess[27] = t54 * (t1 * t36 * t4 - t139 - t141 + t164 * t58 + t65 * t72);
-    hess[28] = t54
-        * (-t164 * t74 + t36 * t44 * std::pow(t72, 2) - t36 * (t138 + t2)
-           + 2 * t70 * t72);
-    hess[29] = t53 * (-t163 * t82 * t84 + t169);
-    hess[30] = t53 * (2 * t162 * t36 * t83 * t84 - t171);
-    hess[31] = t54 * (t164 * t174 - t172 + t173 * t36 - t176 + t72 * t95);
-    hess[32] = t53 * (2 * t162 * t36 * t84 * t99 - t177);
-    hess[33] = t53 * (2 * t106 * t162 * t36 * t84 - t178);
-    hess[34] = t54 * (-t112 * t164 + t179);
-    hess[35] = t180;
-    hess[36] = t181;
-    hess[37] = t183;
-    hess[38] =
-        t134 * (t130 * t162 * t36 * t44 - t131 * t175 - t131 * t70 - t163);
-    hess[39] = t54 * (t186 * t41 + t36 * t80 + t38 * t77 - t76 - t81);
-    hess[40] = t53 * (-t146 + 2 * t184 * t36 * t57 * t84);
-    hess[41] = t53 * (t169 - t187 * t74);
-    hess[42] = t54
-        * (-t186 * t82 - t36 * (t188 + t189) + t45 * std::pow(t77, 2)
-           + t75 * t82);
-    hess[43] = t54 * (t147 * t186 - t190 + t191 * t25 - t193 + t77 * t87);
-    hess[44] = t54 * (t11 * t191 + t174 * t186 - t194 - t195 + t77 * t95);
-    hess[45] = t54 * (t101 * t77 + t186 * t200 - t196 + t199 * t36 - t201);
-    hess[46] = t53 * (2 * t106 * t184 * t36 * t84 - t204);
-    hess[47] = t53 * (-t112 * t187 + t207);
-    hess[48] = t208;
-    hess[49] = t209;
-    hess[50] = t211;
-    hess[51] =
-        t134 * (t130 * t184 * t36 * t44 - t131 * t192 - t131 * t75 - t185);
-    hess[52] = t53 * (2 * t212 * t36 * t40 * t84 - t90);
-    hess[53] = t54 * (t149 + t214 * t58);
-    hess[54] = t53 * (-t171 - t215 * t74);
-    hess[55] = t54 * (t14 * t25 * t36 - t190 - t193 - t214 * t82 + t77 * t87);
-    hess[56] = t54
-        * (-t147 * t87 + 2 * t212 * t36 * t44 * t83
-           + t36 * t44 * std::pow(t83, 2) - t36 * (t189 + t216));
-    hess[57] = t54 * (t174 * t214 + t218);
-    hess[58] = t53 * (2 * t212 * t36 * t84 * t99 - t220);
-    hess[59] = t54 * (t214 * t221 + t223);
-    hess[60] = t53 * (-t112 * t215 - t224);
-    hess[61] = t225;
-    hess[62] = t226;
-    hess[63] = t228;
-    hess[64] = t134
-        * (t130 * t212 * t36 * t44 + (1.0 / 2.0) * t130 * t36 * t44 * t83
-           - t131 * t87 - t213);
-    hess[65] = t53 * (2 * t229 * t36 * t40 * t84 - t97);
-    hess[66] = t53 * (-t152 + 2 * t229 * t36 * t57 * t84);
-    hess[67] = t54 * (-t172 + t173 * t36 - t176 - t231 * t74 + t72 * t95);
-    hess[68] = t54 * (t11 * t14 * t36 - t194 - t195 - t231 * t82 + t77 * t95);
-    hess[69] = t54 * (t147 * t231 + t218);
-    hess[70] = t54
-        * (-t174 * t95 + 2 * t229 * t36 * t44 * t92
-           + t36 * t44 * std::pow(t92, 2) - t36 * (t188 + t216));
-    hess[71] = t53 * (2 * t229 * t36 * t84 * t99 - t234);
-    hess[72] = t53 * (2 * t106 * t229 * t36 * t84 - t235);
-    hess[73] = t54 * (t111 * t95 - t112 * t231 - t236 + t237 * t36 - t239);
-    hess[74] = t241;
-    hess[75] = t242;
-    hess[76] = t243;
-    hess[77] = t134
-        * (t130 * t229 * t36 * t44 + (1.0 / 2.0) * t130 * t36 * t44 * t92
-           - t131 * t95 - t230);
-    hess[78] = t54 * (-t104 + 2 * t244 * t36 * t40 * t44);
-    hess[79] = t53 * (-t153 + 2 * t244 * t36 * t57 * t84);
-    hess[80] = t53 * (-t177 - t245 * t74 * t84);
-    hess[81] = t54 * (t101 * t77 - t196 + t199 * t36 - t201 - t246 * t82);
-    hess[82] = t53 * (-t220 + 2 * t244 * t36 * t83 * t84);
-    hess[83] = t53 * (-t234 + 2 * t244 * t36 * t84 * t92);
-    hess[84] = t54
-        * (-t101 * t200 + 2 * t244 * t36 * t44 * t99
-           + t36 * t44 * std::pow(t99, 2) - t36 * (t247 + t248));
-    hess[85] = t54 * (t221 * t246 + t251);
-    hess[86] = t54 * (t101 * t111 - t112 * t246 + t17 * t36 * t9 - t252 - t253);
-    hess[87] = t255;
-    hess[88] = t256;
-    hess[89] = t258;
-    hess[90] = t134
-        * (-t101 * t131 + t130 * t244 * t36 * t44
-           + (1.0 / 2.0) * t130 * t36 * t44 * t99 - t245);
-    hess[91] = t53 * (-t110 + 2 * t259 * t36 * t40 * t84);
-    hess[92] = t54 * (-t155 + 2 * t259 * t36 * t44 * t57);
-    hess[93] = t53 * (-t178 - t261 * t74);
-    hess[94] = t53 * (-t204 - t261 * t82);
-    hess[95] = t54 * (t147 * t262 + t223);
-    hess[96] = t53 * (-t235 + 2 * t259 * t36 * t84 * t92);
-    hess[97] = t54 * (t200 * t262 + t251);
-    hess[98] = t54
-        * (std::pow(t106, 2) * t36 * t44 + 2 * t106 * t259 * t36 * t44
-           - t109 * t221 - t36 * (t248 + t263));
-    hess[99] =
-        t54 * (t109 * t111 - t112 * t262 + t17 * t23 * t36 - t264 - t266);
-    hess[100] = t267;
-    hess[101] = t268;
-    hess[102] = t269;
-    hess[103] = t134
-        * ((1.0 / 2.0) * t106 * t130 * t36 * t44 - t109 * t131
-           + t130 * t259 * t36 * t44 - t260);
-    hess[104] = t53 * (-t117 + 2 * t270 * t36 * t40 * t84);
-    hess[105] = t53 * (-t156 + 2 * t270 * t36 * t57 * t84);
-    hess[106] = t54 * (t179 - t272 * t74);
-    hess[107] = t53 * (t207 - t271 * t82 * t84);
-    hess[108] = t53 * (-t224 + 2 * t270 * t36 * t83 * t84);
-    hess[109] = t54 * (t111 * t95 + t174 * t272 - t236 + t237 * t36 - t239);
-    hess[110] = t54 * (t101 * t111 + t17 * t249 + t200 * t272 - t252 - t253);
-    hess[111] =
-        t54 * (t109 * t111 + t17 * t23 * t36 + t221 * t272 - t264 - t266);
-    hess[112] = t54
-        * (std::pow(t111, 2) * t45 + t112 * t115 - t112 * t272
-           - t36 * (t247 + t263));
-    hess[113] = t273;
-    hess[114] = t274;
-    hess[115] = t275;
-    hess[116] =
-        t134 * (-t115 * t131 + t130 * t270 * t36 * t44 - t131 * t276 - t271);
-    hess[117] = t124;
-    hess[118] = t157;
-    hess[119] = t180;
-    hess[120] = t208;
-    hess[121] = t225;
-    hess[122] = t241;
-    hess[123] = t255;
-    hess[124] = t267;
-    hess[125] = t273;
-    hess[126] = 0;
-    hess[127] = 0;
-    hess[128] = 0;
-    hess[129] = t280;
-    hess[130] = t125;
-    hess[131] = t160;
-    hess[132] = t181;
-    hess[133] = t209;
-    hess[134] = t226;
-    hess[135] = t242;
-    hess[136] = t256;
-    hess[137] = t268;
-    hess[138] = t274;
-    hess[139] = 0;
-    hess[140] = 0;
-    hess[141] = 0;
-    hess[142] = t281;
-    hess[143] = t128;
-    hess[144] = t161;
-    hess[145] = t183;
-    hess[146] = t211;
-    hess[147] = t228;
-    hess[148] = t243;
-    hess[149] = t258;
-    hess[150] = t269;
-    hess[151] = t275;
-    hess[152] = 0;
-    hess[153] = 0;
-    hess[154] = 0;
-    hess[155] = t282;
-    hess[156] = t283 * (t37 + t67);
-    hess[157] = t283 * (t140 + t167 - t59 + t60);
-    hess[158] = -t283 * (t144 + t175 + t55 * t7 - t69);
-    hess[159] = -t283 * (t192 + t75);
-    hess[160] = t283 * (t217 + t86);
-    hess[161] = t283 * (t238 - t93 + t94);
-    hess[162] = t283 * (t100 + t250);
-    hess[163] = t283 * (-t107 + t108 + t265);
-    hess[164] = -t283 * (-t114 + t276 + t7 * t9);
-    hess[165] = t280;
-    hess[166] = t281;
-    hess[167] = t282;
-    hess[168] =
-        (1.0 / 4.0) * t278 * (t277 + 2) * (t35 - t64 * t7) / std::pow(d, 2);
-}
-
 } // namespace ipc::autogen