1 files changed, 804 insertions, 867 deletions
diff --git a/src/Lagrange-tensor-product/3d.coeffs/3d.cube.order3.smooth0/coeffs-dxx.compute.c b/src/Lagrange-tensor-product/3d.coeffs/3d.cube.order3.smooth0/coeffs-dxx.compute.c
index 70ed1db..7f0c36d 100644
--- a/src/Lagrange-tensor-product/3d.coeffs/3d.cube.order3.smooth0/coeffs-dxx.compute.c
+++ b/src/Lagrange-tensor-product/3d.coeffs/3d.cube.order3.smooth0/coeffs-dxx.compute.c
@@ -1,891 +1,828 @@
-fp t4232;
-fp t4245;
-fp t4248;
-fp t4207;
-fp t4229;
-fp t4246;
-fp t4201;
-fp t4171;
-fp t4434;
-fp t4247;
-fp t4198;
-fp t4173;
-fp t4206;
-fp t4433;
-fp t4368;
-fp t4227;
-fp t4313;
-fp t4066;
-fp t4231;
-fp t4365;
-fp t4242;
-fp t4432;
-fp t4230;
-fp t4316;
-fp t4070;
-fp t4228;
-fp t4375;
-fp t4185;
-fp t4063;
-fp t4431;
-fp t4373;
-fp t4160;
-fp t4314;
-fp t4364;
-fp t4134;
-fp t4430;
+fp t3992;
+fp t4012;
+fp t4011;
+fp t4014;
+fp t4107;
+fp t3913;
+fp t3998;
+fp t3973;
+fp t4169;
+fp t3983;
+fp t4079;
+fp t3776;
+fp t4001;
+fp t4007;
+fp t4077;
+fp t3843;
 fp t4168;
-fp t4120;
-fp t4429;
-fp t4243;
-fp t4360;
-fp t4428;
-fp t4233;
-fp t4354;
-fp t4219;
-fp t4196;
-fp t4147;
-fp t4427;
-fp t4317;
-fp t4019;
-fp t4136;
-fp t4021;
-fp t4426;
-fp t4218;
-fp t4199;
-fp t4328;
-fp t4155;
-fp t4425;
-fp t4226;
-fp t4326;
-fp t4144;
-fp t4034;
-fp t4141;
-fp t4026;
-fp t4424;
-fp t4363;
+fp t4002;
+fp t4013;
+fp t4122;
+fp t3959;
+fp t3984;
+fp t4108;
+fp t4078;
+fp t4121;
+fp t4104;
+fp t4167;
+fp t4123;
+fp t3957;
+fp t4117;
+fp t4105;
+fp t4166;
+fp t3884;
+fp t4165;
+fp t4081;
+fp t4164;
+fp t3982;
+fp t3878;
+fp t3768;
+fp t3979;
+fp t4118;
+fp t4004;
 fp t4163;
-fp t4038;
-fp t4223;
-fp t4324;
+fp t4109;
+fp t4080;
+fp t3773;
 fp t4162;
-fp t4052;
-fp t4423;
-fp t4186;
-fp t4057;
+fp t3989;
+fp t3889;
+fp t3995;
+fp t3969;
+fp t3933;
+fp t4161;
+fp t3891;
+fp t3967;
+fp t4160;
+fp t3996;
+fp t3936;
+fp t3823;
+fp t3990;
+fp t4008;
+fp t3904;
+fp t4159;
+fp t3911;
+fp t3975;
+fp t3954;
+fp t4158;
+fp t4157;
+fp t3809;
+fp t4119;
+fp t3952;
+fp t4156;
+fp t3955;
+fp t4155;
+fp t3966;
+fp t3921;
+fp t3810;
 fp t4154;
-fp t4422;
-fp t4214;
-fp t4330;
-fp t4108;
-fp t3992;
-fp t4216;
-fp t4337;
-fp t4107;
-fp t4421;
-fp t4203;
-fp t4420;
-fp t4133;
-fp t4419;
-fp t4418;
-fp t4183;
+fp t3974;
+fp t4153;
+fp t3962;
+fp t3897;
+fp t3925;
+fp t4152;
+fp t3927;
+fp t3961;
 fp t4151;
-fp t4058;
-fp t4417;
-fp t4119;
-fp t4042;
+fp t3965;
+fp t3920;
+fp t3945;
+fp t4150;
+fp t3947;
+fp t3963;
 fp t4149;
-fp t4043;
-fp t4416;
-fp t4080;
-fp t4415;
-fp t4414;
-fp t4413;
-fp t4412;
-fp t4159;
-fp t4411;
-fp t4195;
-fp t4137;
-fp t4410;
-fp t4033;
-fp t4409;
-fp t4237;
-fp t4200;
-fp t4123;
-fp t4211;
-fp t4319;
-fp t4104;
-fp t4408;
-fp t4350;
-fp t4156;
-fp t4407;
-fp t4220;
-fp t4383;
-fp t4165;
-fp t4406;
+fp t3948;
+fp t3834;
+fp t3964;
+fp t4148;
+fp t4147;
+fp t3837;
+fp t4124;
+fp t4146;
+fp t3994;
+fp t3923;
+fp t3978;
+fp t4145;
+fp t3863;
+fp t3997;
+fp t3971;
+fp t3943;
+fp t4144;
+fp t4120;
+fp t3928;
+fp t4143;
+fp t3888;
+fp t3991;
+fp t3909;
 fp t4142;
-fp t4205;
-fp t4128;
-fp t4405;
-fp t4210;
-fp t4121;
-fp t4404;
+fp t4113;
+fp t4141;
+fp t4111;
 fp t4140;
-fp t4208;
-fp t4403;
-fp t4380;
-fp t4174;
-fp t4118;
-fp t4402;
-fp t4401;
-fp t4177;
-fp t4062;
-fp t4400;
-fp t4197;
+fp t3914;
+fp t3937;
+fp t4139;
+fp t3934;
+fp t3950;
+fp t4138;
+fp t3864;
+fp t3752;
+fp t4116;
+fp t3894;
+fp t3794;
+fp t3825;
+fp t4137;
+fp t3883;
+fp t3787;
+fp t3840;
+fp t4136;
+fp t3919;
+fp t3813;
+fp t4135;
+fp t3839;
+fp t3960;
+fp t3851;
+fp t3958;
+fp t3853;
+fp t4134;
+fp t4133;
+fp t3882;
+fp t4132;
+fp t3784;
+fp t3842;
+fp t4131;
+fp t4005;
+fp t4130;
+fp t3916;
+fp t3812;
+fp t3976;
+fp t4129;
+fp t4006;
+fp t4128;
+fp t4010;
+fp t4009;
 fp t4127;
-fp t4399;
-fp t4039;
-fp t4040;
-fp t4398;
-fp t4349;
-fp t4397;
-fp t4093;
-fp t3988;
-fp t4025;
-fp t4370;
-fp t4158;
-fp t4028;
-fp t4396;
-fp t4395;
-fp t4179;
-fp t4394;
-fp t4333;
-fp t4393;
-fp t4047;
-fp t4122;
-fp t4060;
-fp t4392;
-fp t4175;
-fp t4176;
-fp t4204;
-fp t4191;
-fp t4193;
-fp t4327;
-fp t4391;
-fp t4369;
+fp t3986;
+fp t4126;
 fp t4125;
-fp t4046;
-fp t4390;
-fp t4367;
-fp t4132;
-fp t4389;
-fp t4145;
-fp t4344;
-fp t4388;
-fp t4030;
-fp t4029;
-fp t4387;
-fp t4315;
-fp t4087;
-fp t4386;
-fp t4238;
-fp t4385;
-fp t4384;
-fp t4244;
-fp t4241;
-fp t4382;
-fp t4221;
-fp t4381;
-fp t4379;
-fp t4212;
-fp t4378;
-fp t4239;
-fp t4213;
-fp t4377;
-fp t4376;
-fp t4374;
-fp t4372;
-fp t4371;
-fp t4215;
-fp t4366;
-fp t4202;
-fp t4180;
-fp t4094;
-fp t4322;
-fp t4095;
-fp t4362;
-fp t4361;
-fp t4192;
-fp t4325;
+fp t3970;
+fp t3968;
+fp t4115;
+fp t4114;
+fp t4112;
 fp t4110;
-fp t4359;
-fp t4146;
-fp t4055;
-fp t4358;
-fp t4357;
-fp t4356;
-fp t4355;
-fp t4153;
-fp t4077;
-fp t4017;
-fp t4353;
-fp t4352;
-fp t4167;
-fp t4054;
-fp t4157;
-fp t4011;
-fp t4351;
-fp t4031;
-fp t4348;
-fp t4089;
-fp t4092;
-fp t4347;
-fp t4338;
-fp t4014;
-fp t4045;
-fp t4346;
-fp t4345;
-fp t4090;
+fp t3977;
+fp t3858;
+fp t3862;
+fp t3868;
+fp t3981;
+fp t4106;
+fp t3988;
+fp t3887;
+fp t3895;
+fp t3931;
+fp t3855;
+fp t3745;
 fp t3980;
-fp t3987;
-fp t4343;
-fp t4161;
-fp t4098;
-fp t4342;
+fp t3870;
+fp t4103;
+fp t3869;
+fp t3759;
+fp t3857;
+fp t4102;
+fp t3751;
+fp t3877;
+fp t4101;
+fp t3881;
+fp t4100;
+fp t3867;
+fp t3756;
+fp t3875;
+fp t3765;
 fp t4099;
-fp t4105;
-fp t4341;
-fp t4318;
-fp t4106;
-fp t3981;
-fp t4112;
-fp t4340;
+fp t4098;
+fp t3893;
+fp t3793;
+fp t3820;
+fp t4097;
+fp t3898;
+fp t3796;
 fp t4096;
-fp t3989;
-fp t4339;
-fp t4103;
-fp t4336;
-fp t4182;
-fp t4036;
-fp t4181;
-fp t4015;
-fp t4335;
+fp t3924;
+fp t4095;
+fp t3824;
+fp t4094;
+fp t3949;
+fp t3836;
+fp t4093;
+fp t3860;
+fp t3861;
+fp t4092;
 fp t4091;
-fp t4334;
-fp t4240;
-fp t4332;
-fp t4234;
-fp t4331;
-fp t4329;
-fp t4114;
-fp t4217;
-fp t4323;
-fp t4321;
-fp t4005;
-fp t3994;
-fp t4320;
-fp t4130;
-fp t4131;
-fp t4135;
-fp t4097;
-fp t4101;
-fp t4109;
-fp t4113;
-fp t4115;
-fp t4184;
-fp t4164;
-fp t4126;
-fp t4150;
-fp t4312;
-fp t4311;
-fp t4129;
-fp t4102;
-fp t4310;
-fp t4309;
-fp t3996;
-fp t4022;
-fp t4308;
-fp t4044;
-fp t4307;
-fp t4194;
-fp t4306;
-fp t4305;
-fp t3985;
-fp t4304;
-fp t4041;
-fp t4303;
-fp t4302;
-fp t4301;
-fp t4300;
-fp t4299;
-fp t4012;
-fp t4298;
-fp t4297;
-fp t4296;
-fp t4224;
-fp t4295;
-fp t4294;
-fp t4293;
-fp t4292;
-fp t4291;
-fp t4290;
-fp t4079;
-fp t4289;
-fp t4288;
-fp t4287;
-fp t4286;
-fp t4008;
-fp t4064;
-fp t4285;
-fp t4284;
-fp t3997;
-fp t4283;
-fp t4282;
-fp t4006;
-fp t4281;
-fp t4280;
-fp t4279;
-fp t4277;
-fp t4276;
-fp t4048;
-fp t4275;
-fp t4274;
-fp t4273;
-fp t4272;
-fp t4271;
-fp t4074;
-fp t4078;
-fp t4270;
+fp t4090;
+fp t3865;
+fp t4089;
+fp t3866;
+fp t3879;
+fp t4088;
+fp t4087;
+fp t3892;
+fp t4086;
+fp t4085;
+fp t4084;
+fp t4083;
 fp t4082;
-fp t4081;
-fp t4269;
-fp t4268;
-fp t4267;
-fp t4266;
-fp t4236;
-fp t4265;
-fp t4010;
-fp t4264;
+fp t3856;
+fp t3871;
+fp t3873;
+fp t3874;
+fp t3880;
+fp t3885;
+fp t3900;
+fp t3903;
+fp t3915;
+fp t3929;
+fp t3930;
+fp t3938;
+fp t3939;
+fp t3940;
+fp t3876;
+fp t3899;
+fp t3901;
+fp t3902;
+fp t3935;
+fp t3985;
+fp t4000;
 fp t4076;
-fp t4263;
-fp t4088;
-fp t4262;
+fp t3987;
+fp t4003;
+fp t4075;
+fp t3766;
+fp t4074;
+fp t3815;
+fp t4073;
+fp t3790;
+fp t3822;
+fp t4072;
+fp t3754;
+fp t4071;
+fp t3797;
+fp t4070;
+fp t3799;
+fp t3818;
+fp t4069;
+fp t3800;
+fp t3817;
+fp t4068;
+fp t3805;
+fp t4067;
+fp t3832;
+fp t4066;
+fp t3807;
+fp t3821;
+fp t4065;
+fp t3811;
+fp t4064;
+fp t4063;
+fp t4062;
+fp t3830;
 fp t4061;
-fp t4261;
-fp t4260;
-fp t4259;
-fp t4258;
-fp t4257;
-fp t4256;
-fp t4255;
-fp t4254;
-fp t4253;
-fp t4252;
-fp t4251;
-fp t4235;
-fp t4222;
-fp t4189;
-fp t4172;
-fp t4170;
-fp t4166;
-fp t4100;
+fp t4060;
 fp t4059;
+fp t4058;
+fp t4057;
+fp t3770;
+fp t4056;
+fp t3753;
+fp t4055;
+fp t3827;
+fp t4054;
+fp t3798;
+fp t4053;
+fp t3802;
+fp t4052;
+fp t4051;
+fp t4050;
 fp t4049;
+fp t4048;
+fp t4047;
+fp t3785;
+fp t4045;
+fp t4044;
+fp t3803;
+fp t4043;
+fp t3795;
+fp t3826;
+fp t4042;
+fp t3828;
+fp t4041;
+fp t3835;
+fp t4040;
+fp t4039;
+fp t4038;
+fp t3838;
 fp t4037;
+fp t4036;
 fp t4035;
+fp t4034;
+fp t3890;
+fp t4033;
+fp t4032;
+fp t3774;
+fp t3850;
+fp t4031;
+fp t4030;
+fp t4029;
+fp t3829;
+fp t4028;
 fp t4027;
+fp t4026;
+fp t4025;
+fp t3926;
+fp t4024;
 fp t4023;
-fp t4004;
-fp t4003;
-fp t4001;
-fp t3998;
-fp t3990;
-fp t3986;
-fp t3983;
-fp t3982;
-      t4232 = RATIONAL(-1.0,2.0);
-      t4245 = z*z;
-      t4248 = z*t4245;
-      t4207 = t4232*t4248;
-      t4229 = RATIONAL(1.0,2.0);
-      t4246 = y*y;
-      t4201 = t4229*t4246;
-      t4171 = t4248*t4201;
-      t4434 = t4207+t4171;
-      t4247 = t4246*y;
-      t4198 = t4229*t4245;
-      t4173 = t4247*t4198;
-      t4206 = t4232*t4247;
-      t4433 = t4173+t4206;
-      t4368 = x*t4247;
-      t4227 = RATIONAL(3.0,4.0);
-      t4313 = t4227*t4368;
-      t4066 = t4248*t4313;
-      t4231 = RATIONAL(3.0,2.0);
-      t4365 = x*t4231;
-      t4242 = y*z;
-      t4432 = t4365*t4242+t4066;
-      t4230 = RATIONAL(-3.0,4.0);
-      t4316 = t4230*t4368;
-      t4070 = t4248*t4316;
-      t4228 = RATIONAL(-3.0,2.0);
-      t4375 = x*t4228;
-      t4185 = t4245*t4375;
-      t4063 = y*t4185;
-      t4431 = t4070+t4063;
-      t4373 = x*t4248;
-      t4160 = t4228*t4373;
-      t4314 = t4227*t4373;
-      t4364 = t4246*x;
-      t4134 = t4231*t4364;
-      t4430 = y*t4314+t4248*t4134+t4160;
-      t4168 = t4228*t4368;
-      t4120 = t4245*t4365;
-      t4429 = t4168+z*t4313+t4247*t4120;
-      t4243 = t4246*t4245;
-      t4360 = t4246+t4245;
-      t4428 = (-t4243-1.0+t4360)*x;
-      t4233 = RATIONAL(1.0,3.0);
-      t4354 = t4233*t4245;
-      t4219 = RATIONAL(1.0,6.0);
-      t4196 = t4219*t4247;
-      t4147 = t4245*t4196;
-      t4427 = y*t4354+t4147;
-      t4317 = t4230*t4364;
-      t4019 = t4245*t4317;
-      t4136 = z*t4201;
-      t4021 = x*t4136;
-      t4426 = t4019+t4021;
-      t4218 = RATIONAL(-1.0,6.0);
-      t4199 = t4218*t4248;
-      t4328 = t4219*t4246;
-      t4155 = t4248*t4328;
-      t4425 = (t4155+t4199)*x;
-      t4226 = RATIONAL(-1.0,4.0);
-      t4326 = t4247*t4226;
-      t4144 = t4245*t4326;
-      t4034 = x*t4144;
-      t4141 = x*t4196;
-      t4026 = z*t4141;
-      t4424 = t4034+t4026;
-      t4363 = t4246*z;
-      t4163 = t4218*t4363;
-      t4038 = x*t4163;
-      t4223 = RATIONAL(1.0,4.0);
-      t4324 = t4223*t4246;
-      t4162 = t4245*t4324;
-      t4052 = x*t4162;
-      t4423 = t4038+t4052;
-      t4186 = x*t4201;
-      t4057 = t4245*t4186;
-      t4154 = t4232*t4364;
-      t4422 = t4057+t4154;
-      t4214 = RATIONAL(-1.0,12.0);
-      t4330 = t4248*t4214;
-      t4108 = t4246*t4330;
-      t3992 = x*t4108;
-      t4216 = RATIONAL(1.0,18.0);
-      t4337 = t4248*t4216;
-      t4107 = y*t4337;
-      t4421 = x*t4107+t3992;
-      t4203 = t4218*t4247;
-      t4420 = (t4203+t4147)*x;
-      t4133 = t4232*x*t4245;
-      t4419 = t4057+t4133;
-      t4418 = t4108+t4163;
-      t4183 = x*t4198;
-      t4151 = t4232*t4243;
-      t4058 = x*t4151;
-      t4417 = t4183+t4058;
-      t4119 = y*t4199;
-      t4042 = x*t4119;
-      t4149 = t4248*t4324;
-      t4043 = x*t4149;
-      t4416 = t4042+t4043;
-      t4080 = t4245*t4316;
-      t4415 = t4080+t4173;
-      t4414 = t4433*x;
-      t4413 = t4058+t4186;
-      t4412 = t4434*x;
-      t4159 = t4245*t4203;
-      t4411 = x*t4159+t4141;
-      t4195 = t4232*z;
-      t4137 = x*t4195;
-      t4410 = t4137+t4021;
-      t4033 = t4248*t4317;
-      t4409 = t4171+t4033;
-      t4237 = RATIONAL(-1.0,3.0);
-      t4200 = t4237*z;
-      t4123 = y*t4200;
-      t4211 = RATIONAL(1.0,12.0);
-      t4319 = t4247*t4211;
-      t4104 = t4248*t4319;
-      t4408 = t4123+t4104;
-      t4350 = t4246*t4248;
-      t4156 = t4226*t4350;
-      t4407 = t4136+t4156;
-      t4220 = RATIONAL(-1.0,9.0);
-      t4383 = z*t4220;
-      t4165 = y*t4383;
-      t4406 = t4247*t4337+t4165;
-      t4142 = t4246*t4199;
-      t4205 = t4219*t4248;
-      t4128 = x*t4205;
-      t4405 = t4128+x*t4142;
-      t4210 = t4229*t4247;
-      t4121 = t4245*t4206;
-      t4404 = (t4210+t4121)*x;
-      t4140 = t4246*t4207;
-      t4208 = t4229*t4248;
-      t4403 = (t4140+t4208)*x;
-      t4380 = t4233*z;
-      t4174 = y*t4380;
-      t4118 = t4248*t4196;
-      t4402 = t4174+t4118;
-      t4401 = t4123+t4118;
-      t4177 = y*t4208;
-      t4062 = x*t4177;
-      t4400 = t4043+t4062;
-      t4197 = t4232*y;
-      t4127 = z*t4197;
-      t4399 = t4127+t4177;
-      t4039 = x*t4123;
-      t4040 = y*t4183;
-      t4398 = t4039+t4040;
-      t4349 = t4233*t4246;
-      t4397 = t4155+z*t4349;
-      t4093 = t4245*t4319;
-      t3988 = x*t4093;
-      t4025 = t4226*x*t4243;
-      t4370 = y*t4245;
-      t4158 = t4219*t4370;
-      t4028 = x*t4158;
-      t4396 = t3988+t4025+t4028;
-      t4395 = t4136+t4140-t4246;
-      t4179 = t4247*t4365;
-      t4394 = z*t4316+t4179+t4247*t4185;
-      t4333 = t4158+t4159;
-      t4393 = t4333+t4218*y+t4196;
-      t4047 = t4246*t4314;
-      t4122 = t4228*t4364;
-      t4060 = z*t4122;
-      t4392 = t4047+t4060+t4019;
-      t4175 = t4248*t4203;
-      t4176 = z*t4196;
-      t4204 = t4218*z;
-      t4191 = y*t4204;
-      t4193 = y*t4205;
-      t4327 = t4191+t4193;
-      t4391 = t4175+t4176+t4327;
-      t4369 = t4223*z;
-      t4125 = t4247*t4369;
-      t4046 = x*t4125;
-      t4390 = t4046+t4414;
-      t4367 = y*t4248;
-      t4132 = t4248*t4365;
-      t4389 = t4230*x*t4367+t4248*t4122+t4132;
-      t4145 = z*t4328;
-      t4344 = t4142+t4145;
-      t4388 = t4205+t4344+t4204;
-      t4030 = x*t4156;
-      t4029 = z*t4154;
-      t4387 = t4030+t4029+t4397;
-      t4315 = t4227*t4364;
-      t4087 = t4245*t4134;
-      t4386 = t4122+t4087+z*t4315+t4033;
-      t4238 = x*t4370;
-      t4385 = t4080+t4227*t4238+t4087+t4185;
-      t4384 = t4247*t4208+t4375*t4242+t4070+t4242;
-      t4244 = y*x;
-      t4241 = x*z;
-      t4382 = t4219*y;
-      t4221 = RATIONAL(1.0,9.0);
-      t4381 = z*t4221;
-      t4379 = x*t4218;
-      t4212 = RATIONAL(-1.0,18.0);
-      t4378 = z*t4212;
-      t4239 = y*t4241;
-      t4213 = RATIONAL(-1.0,36.0);
-      t4377 = z*t4213;
-      t4376 = y-t4247;
-      t4374 = z*t4214;
-      t4372 = z*t4216;
-      t4371 = x*t4229;
-      t4215 = RATIONAL(1.0,36.0);
-      t4366 = y*t4215;
-      t4202 = t4237*y;
-      t4180 = y*t4369;
-      t4094 = t4247*t4330;
-      t4322 = t4248*t4211;
-      t4095 = y*t4322;
-      t4362 = t4094+t4095;
-      t4361 = -t4244+t4238;
-      t4192 = y*t4371;
-      t4325 = t4245*t4214;
-      t4110 = y*t4325;
-      t4359 = t3988+x*t4110;
-      t4146 = z*t4206;
-      t4055 = x*t4146;
-      t4358 = t4055+t4247*t4380;
-      t4357 = t4093+t4110;
-      t4356 = t4247*t4245;
-      t4355 = t4223*t4248;
-      t4153 = y*t4207;
-      t4077 = x*t4153;
-      t4017 = z*t4192;
-      t4353 = t4077+t4017;
-      t4352 = t4248*t4213;
-      t4167 = y*t4355;
-      t4054 = x*t4167;
-      t4157 = t4248*t4326;
-      t4011 = x*t4157;
-      t4351 = t4054+t4011;
-      t4031 = y*t4133;
-      t4348 = t4029+t4031;
-      t4089 = y*t4372;
-      t4092 = t4247*t4352;
-      t4347 = t4089+t4092;
-      t4338 = z*t4379;
-      t4014 = y*t4338;
-      t4045 = y*t4128;
-      t4346 = t4014+t4045;
-      t4345 = t4245*t4315+t4151;
-      t4090 = t4247*t4374;
-      t3980 = x*t4090;
-      t3987 = x*t4104;
-      t4343 = t3980+t3987;
-      t4161 = y*t4381;
-      t4098 = z*t4319;
-      t4342 = t4094+t4098;
-      t4099 = t4246*t4322;
-      t4105 = t4214*t4363;
-      t4341 = t4099+t4105;
-      t4318 = t4247*t4215;
-      t4106 = t4248*t4318;
-      t3981 = x*t4106;
-      t4112 = t4247*t4377;
-      t4340 = x*t4112+t3981;
-      t4096 = t4248*t4366;
-      t3989 = x*t4092;
-      t4339 = x*t4096+t3989;
-      t4103 = t4247*t4372;
-      t4336 = t4103+t4107;
-      t4182 = z*t4326;
-      t4036 = x*t4182;
-      t4181 = t4247*t4355;
-      t4015 = x*t4181;
-      t4335 = t4036+t4015;
-      t4091 = t4211*t4242;
-      t4334 = t4030+t4324*t4241;
-      t4240 = x*t4363;
-      t4332 = t4240-t4241;
-      t4234 = RATIONAL(-2.0,1.0);
-      t4331 = t4234-t4248;
-      t4329 = y*t4378;
-      t4114 = z*t4366;
-      t4217 = RATIONAL(-3.0,1.0);
-      t4323 = t4217*t4245;
-      t4321 = t4248*t4212;
-      t4005 = x*t4098;
-      t3994 = x*t4094;
-      t4320 = t4005+t3994;
-      t4130 = z*t4382;
-      t4131 = z*t4203;
-      t4135 = z*t4210;
-      t4097 = t4247*t4378;
-      t4101 = z*t4318;
-      t4109 = y*t4321;
-      t4113 = y*t4352;
-      t4115 = y*t4330;
-      t4184 = t4226*t4367;
-      t4164 = t4218*t4370;
-      t4126 = t4226*t4370;
-      t4150 = t4226*t4363;
-      t4312 = t3987+t4424;
-      t4311 = t4055+t4034+t4015;
-      t4129 = t4223*t4356;
-      t4102 = t4247*t4325;
-      t4310 = t4042+t3987+t3992;
-      t4309 = t3994+t4416;
-      t3996 = x*t4102;
-      t4022 = x*t4131;
-      t4308 = t3996+t4022+t3987;
-      t4044 = x*t4174;
-      t4307 = t4044+t4031+t4045;
-      t4194 = y*t4198;
-      t4306 = t4144+t4162+t4194;
-      t4305 = t3996+t3981+x*t4103;
-      t3985 = x*t4099;
-      t4304 = t3985+t3989+x*t4109;
-      t4041 = x*t4129;
-      t4303 = t4041+t4022+t3994;
-      t4302 = t4030+t4021+t4052;
-      t4301 = t4102+t4162+t4164;
-      t4300 = t4011+t4400;
-      t4299 = t4239+t4077+t4040;
-      t4012 = x*t4135;
-      t4298 = t4041+t4011+t4012;
-      t4297 = t4354+t4237;
-      t4296 = t4233+t4237*t4245;
-      t4224 = RATIONAL(3.0,1.0);
-      t4295 = t4224+t4217*t4246;
-      t4294 = t4237*t4246+t4233;
-      t4293 = t4237+t4349;
-      t4292 = t4217+t4224*t4246;
-      t4291 = t4243+t4194+t4121-t4245;
-      t4290 = t4201+t4149+t4150+t4151;
-      t4079 = x*t4184;
-      t4289 = t4015+t4079+t4403;
-      t4288 = t4129+t4126+t4198+t4151;
-      t4287 = t4115+t4104+t4090+t4091;
-      t4286 = t4351+t4412;
-      t4008 = y*t4120;
-      t4064 = z*t4134;
-      t4285 = t4066+t4248*t4206+t4008+t4064;
-      t4284 = t4335+t4404;
-      t3997 = x*t4115;
-      t4283 = x*t4091+t3997+t4343;
-      t4282 = t4113+t4114+t4112+t4106;
-      t4006 = x*t4095;
-      t4281 = t4006+t4374*t4244+t4320;
-      t4280 = t4294*z;
-      t4279 = t4296*y;
-      t4277 = t4096+t4097+t4347+t4357;
-      t4276 = t4191+t4176+t4357+t4362;
-      t4048 = t4223*t4238;
-      t4275 = t4048+t4017+t4079+t4311;
-      t4274 = t4048+t4014+t3997+t4312;
-      t4273 = t4327+t4341+t4342;
-      t4272 = t4101+t4109+t4341+t4347;
-      t4271 = t4334+t4343+t4346;
-      t4074 = x*t4150;
-      t4078 = x*t4130;
-      t4270 = t4074+t4078+t4005+t4309;
-      t4082 = x*t4126;
-      t4081 = y*t4137;
-      t4269 = t4054+t4082+t4081+t4298;
-      t4268 = t4006+t4082+t4078+t4303;
-      t4267 = t4046+t4081+t4074+t4300;
-      t4266 = t4334+t4335+t4353;
-      t4236 = RATIONAL(2.0,3.0);
-      t4265 = -t4239+t4175+t4015+t4236*t4242+t4345;
-      t4010 = t4245*t4122;
-      t4264 = t4243+t4047+z*t4317+t4134+t4010+t4395;
-      t4076 = t4226*t4239;
-      t4263 = t4046+t4076+t4287+t4351;
-      t4088 = t4245*t4313;
-      t4262 = t4230*t4238+t4088+t4120+t4010+t4291;
-      t4061 = x*t4180;
-      t4261 = x*t4197+t4061+t4040+t4079+t4284;
-      t4260 = t4281+t4287;
-      t4259 = t4153+t4146+t4181+t4242+t4306+t4407;
-      t4258 = t4193+t4131+t4301+t4407+t4408;
-      t4257 = t4161+t4106+t4301+t4336+t4418;
-      t4256 = t4176+t4119+t4306+t4408+t4418;
-      t4255 = t4167+t4127+t4135+t4157+t4288+t4434;
-      t4254 = t4125+t4157+t4290+t4399+t4433;
-      t4253 = t4199+t4200+t4130+t4131+t4288+t4362+t4397;
-      t4252 = t4202+t4130+t4203+t4119+t4290+t4342+t4427;
-      t4251 = t4184+t4208+RATIONAL(1.0,1.0)+t4210+t4195+t4180+t4197+t4181+t4182
-+t4291+t4395;
-      t4235 = RATIONAL(-2.0,3.0);
-      t4222 = RATIONAL(2.0,1.0);
-      t4189 = t4233*t4367;
-      t4172 = t4247*t4200;
-      t4170 = z*t4371;
-      t4166 = t4248*t4202;
-      t4100 = t4247*t4321;
-      t4059 = y*t4132;
-      t4049 = y*t4160;
-      t4037 = x*t4164;
-      t4035 = x*t4145;
-      t4027 = z*t4179;
-      t4023 = z*t4168;
-      t4004 = x*t4329;
-      t4003 = t4211*t4238;
-      t4001 = x*t4113;
-      t3998 = x*t4089;
-      t3990 = x*t4097;
-      t3986 = t4211*t4240;
-      t3983 = x*t4101;
-      t3982 = x*t4105;
-      coeffs_dxx->coeff_m1_m1_m1 = t4035+x*t4165+t4257+t3990+t4304+t4396;
-      coeffs_dxx->coeff_0_m1_m1 = t4100+t4247*t4383+(t4220*t4248+RATIONAL(-2.0,
-9.0)*z)*y+t4307+t4312+t4345+t4387+t4427;
-      coeffs_dxx->coeff_p1_m1_m1 = t4257+t4303+t4398+t4416+t4426;
-      coeffs_dxx->coeff_p2_m1_m1 = t4037+x*t4161+t4305+t4421+t4423;
-      coeffs_dxx->coeff_m1_0_m1 = x*t4280+t4253+t4274+t4405+t4419;
-      coeffs_dxx->coeff_0_0_m1 = t4262+(t4236+t4235*t4246)*z+t4294*t4248+t4119+
-t4017+t4286+t4332+t4358+t4401;
-      coeffs_dxx->coeff_p1_0_m1 = t4012+t4081+t4253+t4289-t4332+t4385;
-      coeffs_dxx->coeff_p2_0_m1 = t4293*t4241+t4268+t4417+t4425;
-      coeffs_dxx->coeff_m1_p1_m1 = t3985+t4256+t4035+t4025+t4303+t4307;
-      coeffs_dxx->coeff_0_p1_m1 = t4189+t4012+t4008+t4077-t4370+t4172+t4265+
-t4387+t4415;
-      coeffs_dxx->coeff_p1_p1_m1 = t4055+t4256+t4088+t4239+t4063+t4300+t4426;
-      coeffs_dxx->coeff_p2_p1_m1 = t4310+t4398+t4423+t4424;
-      coeffs_dxx->coeff_m1_p2_m1 = t4003+t4004+t4001+t4277+t4305;
-      coeffs_dxx->coeff_0_p2_m1 = t4109+t4247*t4381+t4268+t4333+t4406;
-      coeffs_dxx->coeff_p1_p2_m1 = t4274+t4277;
-      coeffs_dxx->coeff_p2_p2_m1 = t3998+t3990+t4339+t4359;
-      coeffs_dxx->coeff_m1_m1_0 = x*t4279+t4252+t4271+t4411+t4422;
-      coeffs_dxx->coeff_0_m1_0 = t4189+t4011+t4264+t4131+(t4236+t4235*t4245)*y+
-t4296*t4247+t4353+t4361+t4390+t4401;
-      coeffs_dxx->coeff_p1_m1_0 = t4252+t4081+t4062+t4284-t4361+t4386;
-      coeffs_dxx->coeff_p2_m1_0 = t4297*t4244+t4270+t4413+t4420;
-      coeffs_dxx->coeff_m1_0_0 = t4251+t4170+t4192+t4076+t4428+t4286+t4348+
-t4390;
-      coeffs_dxx->coeff_0_0_0 = z+t4285+(t4234*t4246+t4222-t4376)*t4245+(t4248+
-t4222-z)*t4246+(t4228*z+t4292*t4245+(t4228+t4227*z)*y+t4295)*x+t4135+t4331+
-t4376+t4389+t4394+t4399;
-      coeffs_dxx->coeff_p1_0_0 = (t4231*z+t4217+t4246*t4323+t4360*t4224+(t4230*
-z+t4231)*y)*x+t4251+t4060+t4429+t4430+t4431;
-      coeffs_dxx->coeff_p2_0_0 = -t4428+t4261+t4403+t4410;
-      coeffs_dxx->coeff_m1_p1_0 = y+t4254-t4370+t4266+t4361+t4404+t4422;
-      coeffs_dxx->coeff_0_p1_0 = t4247+t4059-t4356+t4264+t4146+(t4222*t4245+(
-t4323+t4224)*x+t4331)*y+t4384+t4429;
-      coeffs_dxx->coeff_p1_p1_0 = t4049+(-t4245+1.0+(t4217+t4224*t4245)*x)*y+
-t4254+t4386+t4394+t4432;
-      coeffs_dxx->coeff_p2_p1_0 = t4267-t4361+t4413+t4414;
-      coeffs_dxx->coeff_m1_p2_0 = x*t4382+t4037+t4260+t4393+t4420;
-      coeffs_dxx->coeff_0_p2_0 = t4297*t4247+t4279+t4261+t4391;
-      coeffs_dxx->coeff_p1_p2_0 = t4263+t4192+t4031+t4393+t4414;
-      coeffs_dxx->coeff_p2_p2_0 = y*t4379+t4028+t4283+t4411;
-      coeffs_dxx->coeff_m1_m1_p1 = t4044+t4029+t4026+t4258+t4309+t4396;
-      coeffs_dxx->coeff_0_m1_p1 = t4034+t4062+t4031+t4064+t4166+t4265-t4363+
-t4358+t4409+t4427;
-      coeffs_dxx->coeff_p1_m1_p1 = t4258+t4298+t4299+t4392;
-      coeffs_dxx->coeff_p2_m1_p1 = t4045+t4039+t4037+t4302+t4308;
-      coeffs_dxx->coeff_m1_0_p1 = z+t4255-t4363+t4275+t4332+t4403+t4419;
-      coeffs_dxx->coeff_0_0_p1 = t4262+(-t4247+t4234+t4222*t4246+t4295*x)*z+
-t4248+t4153-t4350+t4027+t4384+t4430;
-      coeffs_dxx->coeff_p1_0_p1 = t4023+(1.0-t4246+t4292*x)*z+t4255+t4385+t4389
-+t4432;
-      coeffs_dxx->coeff_p2_0_p1 = t4269-t4332+t4412+t4417;
-      coeffs_dxx->coeff_m1_p1_p1 = t4259+t4025-t4239+t4298+t4348+t4400;
-      coeffs_dxx->coeff_0_p1_p1 = t4049+t4023+(t4247-t4246)*z+(t4248-t4245+(
-t4224*x+t4234)*z)*y+t4285+t4345+t4409+t4415;
-      coeffs_dxx->coeff_p1_p1_p1 = t4217*t4239+t4059+t4259+t4088+t4027+t4392+
-t4431;
-      coeffs_dxx->coeff_p2_p1_p1 = t4299+t4302+t4311;
-      coeffs_dxx->coeff_m1_p2_p1 = t3997+t4003+t4078+t4276+t4308;
-      coeffs_dxx->coeff_0_p2_p1 = t4119+t4172+t4269+t4333+t4402;
-      coeffs_dxx->coeff_p1_p2_p1 = t4275+t4276;
-      coeffs_dxx->coeff_p2_p2_p1 = t3994+t4014+t4026+t4006+t4359;
-      coeffs_dxx->coeff_m1_m1_p2 = t4004+t3986+t4272+t4340+t4421;
-      coeffs_dxx->coeff_0_m1_p2 = t4221*t4367+t4097+t4270+t4344+t4406;
-      coeffs_dxx->coeff_p1_m1_p2 = t4271+t4272;
-      coeffs_dxx->coeff_p2_m1_p2 = t3982+t3998+t3983+t4304;
-      coeffs_dxx->coeff_m1_0_p2 = t4219*t4241+t4038+t4260+t4388+t4425;
-      coeffs_dxx->coeff_0_0_p2 = t4036+t4061+t4293*t4248+t4280+t4289+t4391+
-t4410;
-      coeffs_dxx->coeff_p1_0_p2 = t4170+t4029+t4263+t4388+t4412;
-      coeffs_dxx->coeff_p2_0_p2 = t4035+t4338+t4283+t4405;
-      coeffs_dxx->coeff_m1_p1_p2 = t3980+t3986+t4078+t4273+t4310;
-      coeffs_dxx->coeff_0_p1_p2 = t4166+t4131+t4267+t4344+t4402;
-      coeffs_dxx->coeff_p1_p1_p2 = t4266+t4273;
-      coeffs_dxx->coeff_p2_p1_p2 = t3985+t3982+t4320+t4346;
-      coeffs_dxx->coeff_m1_p2_p2 = t4377*t4244+t3983+t4282+t4339;
-      coeffs_dxx->coeff_0_p2_p2 = t4329+t4100+t4283+t4336;
-      coeffs_dxx->coeff_p1_p2_p2 = t4281+t4282;
-      coeffs_dxx->coeff_p2_p2_p2 = t4001+x*t4114+t4340;
+fp t3946;
+fp t4022;
+fp t3912;
+fp t4021;
+fp t4020;
+fp t4019;
+fp t4018;
+fp t4017;
+fp t3999;
+fp t3944;
+fp t3942;
+fp t3941;
+fp t3905;
+fp t3886;
+fp t3872;
+fp t3848;
+fp t3819;
+fp t3792;
+fp t3788;
+fp t3786;
+fp t3782;
+fp t3769;
+fp t3764;
+fp t3762;
+fp t3761;
+fp t3757;
+fp t3750;
+fp t3749;
+fp t3746;
+      t3992 = RATIONAL(1.0,2.0);
+      t4012 = z*z;
+      t4011 = z*t4012;
+      t4014 = y*y;
+      t4107 = t4014*t4011;
+      t3913 = t3992*t4107;
+      t3998 = RATIONAL(-1.0,2.0);
+      t3973 = t3998*t4011;
+      t4169 = t3913+t3973;
+      t3983 = RATIONAL(3.0,4.0);
+      t4079 = x*t4107;
+      t3776 = t3983*t4079;
+      t4001 = RATIONAL(-3.0,4.0);
+      t4007 = t4014*t4012;
+      t4077 = x*t4007;
+      t3843 = t4001*t4077;
+      t4168 = t3776+t3843;
+      t4002 = RATIONAL(-3.0,2.0);
+      t4013 = y*t4014;
+      t4122 = t4013*x;
+      t3959 = t4002*t4122;
+      t3984 = RATIONAL(3.0,2.0);
+      t4108 = t4013*t4012;
+      t4078 = x*t4108;
+      t4121 = t4013*z;
+      t4104 = x*t4121;
+      t4167 = t3983*t4104+t3984*t4078+t3959;
+      t4123 = t4011*x;
+      t3957 = t4002*t4123;
+      t4117 = y*t4011;
+      t4105 = x*t4117;
+      t4166 = t3983*t4105+t3984*t4079+t3957;
+      t3884 = t3984*t4122;
+      t4165 = t4001*t4104+t4002*t4078+t3884;
+      t4081 = t4012+t4014;
+      t4164 = (t4007+1.0-t4081)*x;
+      t3982 = RATIONAL(-1.0,12.0);
+      t3878 = t3982*t4107;
+      t3768 = x*t3878;
+      t3979 = RATIONAL(1.0,12.0);
+      t4118 = x*t4014;
+      t4004 = z*t4118;
+      t4163 = t3979*t4004+t3768;
+      t4109 = t4011*t4013;
+      t4080 = x*t4109;
+      t3773 = t3983*t4080;
+      t4162 = t3773+y*t3957;
+      t3989 = RATIONAL(1.0,6.0);
+      t3889 = t3989*t4108;
+      t3995 = RATIONAL(-1.0,6.0);
+      t3969 = t3995*t4013;
+      t3933 = x*t3969;
+      t4161 = x*t3889+t3933;
+      t3891 = t3989*t4107;
+      t3967 = t3995*t4011;
+      t4160 = (t3891+t3967)*x;
+      t3996 = RATIONAL(-1.0,4.0);
+      t3936 = t3996*t4121;
+      t3823 = x*t3936;
+      t3990 = RATIONAL(1.0,4.0);
+      t4008 = y*z;
+      t3904 = t3990*t4008;
+      t4159 = x*t3904+t3823;
+      t3911 = t3992*t4108;
+      t3975 = t3998*t4013;
+      t3954 = x*t3975;
+      t4158 = x*t3911+t3954;
+      t4157 = t4169*x;
+      t3809 = t3992*t4077;
+      t4119 = x*t4012;
+      t3952 = t3998*t4119;
+      t4156 = t3809+t3952;
+      t3955 = t3998*t4118;
+      t4155 = t3809+t3955;
+      t3966 = t3992*t4014;
+      t3921 = x*t3966;
+      t3810 = z*t3921;
+      t4154 = t3810+t3843;
+      t3974 = t3998*z;
+      t4153 = t3810+x*t3974;
+      t3962 = t3989*t4013;
+      t3897 = x*t3962;
+      t3925 = t3995*t4108;
+      t4152 = x*t3925+t3897;
+      t3927 = t3995*t4107;
+      t3961 = t3989*t4011;
+      t4151 = (t3927+t3961)*x;
+      t3965 = t3992*t4013;
+      t3920 = x*t3965;
+      t3945 = t3998*t4108;
+      t4150 = x*t3945+t3920;
+      t3947 = t3998*t4107;
+      t3963 = t3992*t4011;
+      t4149 = (t3947+t3963)*x;
+      t3948 = t3998*t4007;
+      t3834 = x*t3948;
+      t3964 = t3992*t4012;
+      t4148 = t3834+x*t3964;
+      t4147 = t3834+t3921;
+      t3837 = y*t3952;
+      t4124 = t3992*x;
+      t4146 = t3837+y*t4124;
+      t3994 = RATIONAL(-1.0,9.0);
+      t3923 = t3994*t4008;
+      t3978 = RATIONAL(1.0,18.0);
+      t4145 = t3978*t4109+t3923;
+      t3863 = t3979*t4109;
+      t3997 = RATIONAL(-1.0,3.0);
+      t3971 = t3997*z;
+      t3943 = y*t3971;
+      t4144 = t3863+t3943;
+      t4120 = t4014*z;
+      t3928 = t3995*t4120;
+      t4143 = t3878+t3928;
+      t3888 = t3989*t4109;
+      t3991 = RATIONAL(1.0,3.0);
+      t3909 = t3991*t4008;
+      t4142 = t3888+t3909;
+      t4113 = t3991*t4012;
+      t4141 = t3889+y*t4113;
+      t4111 = t3991*t4014;
+      t4140 = t3891+z*t4111;
+      t3914 = z*t3966;
+      t3937 = t3996*t4107;
+      t4139 = t3914+t3937;
+      t3934 = t3996*t4109;
+      t3950 = y*t3974;
+      t4138 = t3934+t3950;
+      t3864 = t3979*t4108;
+      t3752 = x*t3864;
+      t4116 = y*t4012;
+      t3894 = t3989*t4116;
+      t3794 = x*t3894;
+      t3825 = t3996*t4077;
+      t4137 = t3752+t3794+t3825;
+      t3883 = t3984*t4119;
+      t3787 = y*t3883;
+      t3840 = t4001*t4078;
+      t4136 = t3787+t3840+t3911;
+      t3919 = z*t4124;
+      t3813 = y*t3919;
+      t4135 = t3813+t3888+t3943;
+      t3839 = t4001*t4080;
+      t3960 = t4002*t4118;
+      t3851 = z*t3960;
+      t3958 = t4002*t4119;
+      t3853 = y*t3958;
+      t4134 = t3839+t3851+t3853;
+      t4133 = t3914+t3947-t4014;
+      t3882 = t3984*t4123;
+      t4132 = t3773+t4001*t4105+t4002*t4079+t3882;
+      t3784 = t3984*t4077;
+      t3842 = t4001*t4079;
+      t4131 = t3983*t4004+t3784+t3842+t3960;
+      t4005 = x*t4116;
+      t4130 = t3983*t4005+t3784+t3840+t3958;
+      t3916 = y*t3964;
+      t3812 = x*t3916;
+      t3976 = t3998*y;
+      t4129 = t3812+x*t3976+t4150;
+      t4006 = x*t4008;
+      t4128 = t3839+t4002*t4006+t3992*t4109+t4008;
+      t4010 = x*y;
+      t4009 = x*z;
+      t4127 = y+z;
+      t3986 = RATIONAL(2.0,3.0);
+      t4126 = t3986*y;
+      t4125 = t3989*z;
+      t3970 = t3995*y;
+      t3968 = t3995*z;
+      t4115 = t4011-z;
+      t4114 = t4013-y;
+      t4112 = t3991*t4013;
+      t4110 = t3997*t4013;
+      t3977 = RATIONAL(1.0,36.0);
+      t3858 = t3977*t4008;
+      t3862 = t3978*t4008;
+      t3868 = t3979*t4008;
+      t3981 = RATIONAL(-1.0,18.0);
+      t4106 = t3981*t4008;
+      t3988 = RATIONAL(1.0,9.0);
+      t3887 = t3988*t4008;
+      t3895 = y*t4125;
+      t3931 = y*t3968;
+      t3855 = t3977*t4109;
+      t3745 = x*t3855;
+      t3980 = RATIONAL(-1.0,36.0);
+      t3870 = t3980*t4121;
+      t4103 = t3745+x*t3870;
+      t3869 = t3980*t4109;
+      t3759 = x*t3869;
+      t3857 = t3977*t4117;
+      t4102 = x*t3857+t3759;
+      t3751 = x*t3863;
+      t3877 = t3982*t4121;
+      t4101 = t3751+x*t3877;
+      t3881 = t3982*t4116;
+      t4100 = t3752+x*t3881;
+      t3867 = t3979*t4117;
+      t3756 = x*t3867;
+      t3875 = t3982*t4109;
+      t3765 = x*t3875;
+      t4099 = t3756+t3765;
+      t4098 = t3983*t4077+t3948;
+      t3893 = y*t3961;
+      t3793 = x*t3893;
+      t3820 = x*t3931;
+      t4097 = t3793+t3820;
+      t3898 = t3990*t4109;
+      t3796 = x*t3898;
+      t4096 = t3796+t3823;
+      t3924 = t3995*t4109;
+      t4095 = t3796+t3924;
+      t3824 = x*t3937;
+      t4094 = t3990*t4004+t3824;
+      t3949 = y*t3973;
+      t3836 = x*t3949;
+      t4093 = t3836+t3991*t4117;
+      t3860 = t3978*t4121;
+      t3861 = t3978*t4117;
+      t4092 = t3860+t3861;
+      t4091 = t3862+t3869;
+      t4090 = t3864+t3881;
+      t3865 = t3979*t4121;
+      t4089 = t3865+t3875;
+      t3866 = t3979*t4107;
+      t3879 = t3982*t4120;
+      t4088 = t3866+t3879;
+      t4087 = t3867+t3875;
+      t3892 = t3989*t4120;
+      t4086 = t3892+t3927;
+      t4085 = t3893+t3931;
+      t4084 = t3894+t3925;
+      t4083 = t4009-t4004;
+      t4082 = t4010-t4005;
+      t3856 = t3977*t4121;
+      t3871 = t3980*t4117;
+      t3873 = t3981*t4121;
+      t3874 = t3981*t4117;
+      t3880 = t3982*t4117;
+      t3885 = t3984*t4118;
+      t3900 = t3990*t4121;
+      t3903 = t3990*t4117;
+      t3915 = y*t3963;
+      t3929 = y*t3967;
+      t3930 = t3995*t4116;
+      t3938 = t3996*t4120;
+      t3939 = t3996*t4117;
+      t3940 = t3996*t4116;
+      t3876 = t3982*t4108;
+      t3899 = t3990*t4108;
+      t3901 = t3990*t4107;
+      t3902 = t3990*t4007;
+      t3935 = t3996*t4108;
+      t3985 = RATIONAL(3.0,1.0);
+      t4000 = RATIONAL(-2.0,1.0);
+      t4076 = t3985*x+t4000;
+      t3987 = RATIONAL(2.0,1.0);
+      t4003 = RATIONAL(-3.0,1.0);
+      t4075 = t4003*x+t3987;
+      t3766 = x*t3876;
+      t4074 = t3745+x*t3860+t3766;
+      t3815 = z*t3933;
+      t4073 = t3751+t3766+t3815;
+      t3790 = z*t3897;
+      t3822 = x*t3935;
+      t4072 = t3751+t3790+t3822;
+      t3754 = x*t3866;
+      t4071 = t3754+t3759+x*t3874;
+      t3797 = x*t3899;
+      t4070 = t3765+t3797+t3815;
+      t3799 = x*t3901;
+      t3818 = x*t3929;
+      t4069 = t3765+t3799+t3818;
+      t3800 = x*t3902;
+      t3817 = x*t3928;
+      t4068 = t3768+t3800+t3817;
+      t3805 = x*t3909;
+      t4067 = t3793+t3805+t3837;
+      t3832 = z*t3954;
+      t4066 = t3796+t3822+t3832;
+      t3807 = z*t3920;
+      t3821 = x*t3934;
+      t4065 = t3797+t3807+t3821;
+      t3811 = x*t3915;
+      t4064 = t3799+t3811+t3821;
+      t4063 = t3800+t3810+t3824;
+      t4062 = t3811+t3837-t4006;
+      t3830 = x*t3943;
+      t4061 = t3812+t3818+t3830;
+      t4060 = t3812+t3836+t4006;
+      t4059 = t3876+t3902+t3930;
+      t4058 = t3902+t3916+t3935;
+      t4057 = t3991+t3997*t4014;
+      t3770 = x*t3880;
+      t4056 = x*t3868+t3770+t4101;
+      t3753 = x*t3865;
+      t4055 = t3753+t3982*t4006+t4099;
+      t3827 = x*t3939;
+      t4054 = t3796+t3827+t4149;
+      t3798 = x*t3900;
+      t4053 = t3798+t3821+t4158;
+      t3802 = x*t3903;
+      t4052 = t3802+t3821+t4157;
+      t4051 = t3855+t3858+t3870+t3871;
+      t4050 = t3863+t3868+t3877+t3880;
+      t4049 = t3899+t3940+t3948+t3964;
+      t4048 = t3901+t3938+t3948+t3966;
+      t4047 = t3916+t3945+t4007-t4012;
+      t3785 = z*t3885;
+      t4045 = t3785+t3842+t3913+t4098;
+      t4044 = t4094+t4097+t4101;
+      t3803 = t3990*t4005;
+      t4043 = t3770+t3803+t3820+t4072;
+      t3795 = x*t3895;
+      t3826 = x*t3938;
+      t4042 = t3753+t3795+t3826+t4069;
+      t3828 = x*t3940;
+      t4041 = t3756+t3795+t3828+t4070;
+      t3835 = z*t3955;
+      t4040 = t3824+t3835+t4098+t4140;
+      t4039 = t3813+t3836+t4094+t4096;
+      t4038 = t3803+t3813+t3827+t4066;
+      t3838 = x*t3950;
+      t4037 = t3802+t3828+t3838+t4065;
+      t4036 = t3798+t3826+t3838+t4064;
+      t4035 = t3856+t3874+t4088+t4091;
+      t4034 = t3857+t3873+t4090+t4091;
+      t3890 = z*t3962;
+      t4033 = t3890+t3931+t4087+t4090;
+      t4032 = t4085+t4088+t4089;
+      t3774 = t3983*t4078;
+      t3850 = t4002*t4077;
+      t4031 = t3774+t4001*t4005+t3850+t3883+t4047;
+      t4030 = t3776+t4001*t4004+t3850+t3885+t4007+t4133;
+      t4029 = t3827+t3890+t4085+t4095+t4159;
+      t3829 = t3996*t4006;
+      t4028 = t3798+t3829+t3835+t3919+t4052;
+      t4027 = t3961+t3968+t4050+t4086;
+      t4026 = t3962+t3970+t4050+t4084;
+      t4025 = t3855+t3887+t4059+t4092+t4143;
+      t3926 = z*t3969;
+      t4024 = t3893+t3926+t4059+t4139+t4144;
+      t4023 = t3890+t3929+t4058+t4143+t4144;
+      t3946 = z*t3975;
+      t4022 = t3898+t3946+t3949+t4008+t4058+t4139;
+      t3912 = z*t3965;
+      t4021 = t3903+t3912+t4049+t4138+t4169;
+      t4020 = t3900+t3911+t3915+t3975+t4048+t4138;
+      t4019 = t3895+t3929+t3969+t3997*y+t4048+t4089+t4141;
+      t4018 = t3895+t3926+t3967+t3971+t4049+t4087+t4140;
+      t4017 = t3898+t3904+t3936+t3939+t3963+t3965+t3974+t3976+RATIONAL(1.0,1.0)
++t4047+t4133;
+      t3999 = RATIONAL(-2.0,3.0);
+      t3944 = t3998*t4109;
+      t3942 = t3997*t4117;
+      t3941 = z*t4110;
+      t3905 = z*t4112;
+      t3886 = z*t4126;
+      t3872 = t3981*t4109;
+      t3848 = z*t3959;
+      t3819 = x*t3930;
+      t3792 = x*t3892;
+      t3788 = t3984*t4006;
+      t3786 = y*t3882;
+      t3782 = z*t3884;
+      t3769 = x*t3879;
+      t3764 = x*t4106;
+      t3762 = x*t3873;
+      t3761 = x*t3871;
+      t3757 = t3979*t4005;
+      t3750 = x*t3862;
+      t3749 = x*t3861;
+      t3746 = x*t3856;
+      coeffs_dxx->coeff_m1_m1_m1 = t3792+t4025+t3762+x*t3923+t4071+t4137;
+      coeffs_dxx->coeff_0_m1_m1 = RATIONAL(-2.0,9.0)*t4008+t3872+(t4121+t4117)*
+t3994+t4040+t4067+t4072+t4141;
+      coeffs_dxx->coeff_p1_m1_m1 = t3799+t4025+t4061+t4070+t4154;
+      coeffs_dxx->coeff_p2_m1_m1 = t3749+t3819+x*t3887+t4068+t4074;
+      coeffs_dxx->coeff_m1_0_m1 = t4057*t4009+t4018+t4043+t4151+t4156;
+      coeffs_dxx->coeff_0_0_m1 = (t3999*t4014+t3986)*z+t4031+t4057*t4011+t3832+
+t3929+t3905+t4052-t4083+t4135;
+      coeffs_dxx->coeff_p1_0_m1 = t4018+t3838+t3807+t4054+t4083+t4130;
+      coeffs_dxx->coeff_p2_0_m1 = (t3997+t4111)*t4009+t4041+t4148+t4160;
+      coeffs_dxx->coeff_m1_p1_m1 = t3792+t4023+t3825+t3754+t4067+t4070;
+      coeffs_dxx->coeff_0_p1_m1 = t3886+t4040-t4116+t3941-t4006+t3807+t4093+
+t4095+t4136;
+      coeffs_dxx->coeff_p1_p1_m1 = t3853+t3774+t4023+t4006+t3832+t4064+t4154;
+      coeffs_dxx->coeff_p2_p1_m1 = t4061+t4068+t4072;
+      coeffs_dxx->coeff_m1_p2_m1 = t3764+t3757+t3761+t4034+t4074;
+      coeffs_dxx->coeff_0_p2_m1 = t3874+t3988*t4121+t4041+t4084+t4145;
+      coeffs_dxx->coeff_p1_p2_m1 = t4034+t4043;
+      coeffs_dxx->coeff_p2_p2_m1 = t3762+t3750+t4100+t4102;
+      coeffs_dxx->coeff_m1_m1_0 = (t3991+t3997*t4012)*t4010+t4019+t4044+t4152+
+t4155;
+      coeffs_dxx->coeff_0_m1_0 = t4030+t4126+t4112+(t4110+t3999*y)*t4012+t3926+
+t4053-t4082+t4093+t4135;
+      coeffs_dxx->coeff_p1_m1_0 = t4019+t3811+t3838+t4082+t4096+t4131+t4150;
+      coeffs_dxx->coeff_p2_m1_0 = (t4113+t3997)*t4010+t4042+t4147+t4161;
+      coeffs_dxx->coeff_m1_0_0 = t4017-t4164+t4028+t4146+t4158;
+      coeffs_dxx->coeff_0_0_0 = t3912+t3785+t3915+t3944+t4000+t3787+t3950+(
+t3987+t4115)*t4014+(t3987+t4000*t4014+t4114)*t4012+(t3985*t4007+t3985+t3983*
+t4008+t4081*t4003+t4127*t4002)*x-t4114-t4115+t4132+t4165;
+      coeffs_dxx->coeff_p1_0_0 = t4017+(t4001*t4008+t4003*t4007+t4003+t4081*
+t3985+t4127*t3984)*x+t4134+t4166+t4167;
+      coeffs_dxx->coeff_p2_0_0 = t4164+t4054+t4129+t4153+t4159;
+      coeffs_dxx->coeff_m1_p1_0 = y+t4020-t4116+t4039-t4082+t4150+t4155;
+      coeffs_dxx->coeff_0_p1_0 = t4030-t4108+t4013+(-t4011+t4075*t4012+t4076)*y
++t3786+t3946+t4128+t4167;
+      coeffs_dxx->coeff_p1_p1_0 = t4020+t3788+(1.0-t4012+(t4003+t3985*t4012)*x)
+*y+t4131+t4162+t4165;
+      coeffs_dxx->coeff_p2_p1_0 = t4036+t4082+t4147+t4158;
+      coeffs_dxx->coeff_m1_p2_0 = t4026+t3819+t3989*t4010+t4055+t4161;
+      coeffs_dxx->coeff_0_p2_0 = (t4013+t4116)*t3997+(y+t4108)*t3991+t4029+
+t4129;
+      coeffs_dxx->coeff_p1_p2_0 = t4026+t3802+t3829+t4053+t4146;
+      coeffs_dxx->coeff_p2_p2_0 = x*t3970+t3794+t4056+t4152;
+      coeffs_dxx->coeff_m1_m1_p1 = t4024+t3835+t3790+t3805+t4069+t4137;
+      coeffs_dxx->coeff_0_m1_p1 = t3942+t3905+t3924+t3886-t4120+t4045+t4062+
+t4066+t4141;
+      coeffs_dxx->coeff_p1_m1_p1 = t4024+t3851+t4060+t4065+t4168;
+      coeffs_dxx->coeff_p2_m1_p1 = t3819+t3793+t3830+t4063+t4073;
+      coeffs_dxx->coeff_m1_0_p1 = z-t4120+t4021+t4038-t4083+t4149+t4156;
+      coeffs_dxx->coeff_0_0_p1 = t4011+t3782-t4107+t4031+(-t4013+t4075*t4014+
+t4076)*z+t3949+t4128+t4166;
+      coeffs_dxx->coeff_p1_0_p1 = t4021+t3788+t3848+(-t4014+1.0+(t3985*t4014+
+t4003)*x)*z+t4130+t4132;
+      coeffs_dxx->coeff_p2_0_p1 = t4037+t4083+t4148+t4157;
+      coeffs_dxx->coeff_m1_p1_p1 = t3799+t3835+t3825+t4022+t4062+t4065;
+      coeffs_dxx->coeff_0_p1_p1 = t3944+t3848+(-t4014+t4013)*z+(t4011-t4012+
+t4076*z)*y+t4045+t4136+t4162;
+      coeffs_dxx->coeff_p1_p1_p1 = t3782+t4003*t4006+t3774+t4022+t3786+t4134+
+t4168;
+      coeffs_dxx->coeff_p2_p1_p1 = t4060+t4063+t4066;
+      coeffs_dxx->coeff_m1_p2_p1 = t3757+t3770+t3795+t4033+t4073;
+      coeffs_dxx->coeff_0_p2_p1 = t3929+t3941+t4037+t4084+t4142;
+      coeffs_dxx->coeff_p1_p2_p1 = t4033+t4038;
+      coeffs_dxx->coeff_p2_p2_p1 = t3790+t3820+t4099+t4100;
+      coeffs_dxx->coeff_m1_m1_p2 = t3749+t3764+t4035+t4103+t4163;
+      coeffs_dxx->coeff_0_m1_p2 = t3988*t4117+t3873+t4042+t4086+t4145;
+      coeffs_dxx->coeff_p1_m1_p2 = t4035+t4044;
+      coeffs_dxx->coeff_p2_m1_p2 = t3769+t3750+t3746+t4071;
+      coeffs_dxx->coeff_m1_0_p2 = x*t4125+t3817+t4027+t4055+t4160;
+      coeffs_dxx->coeff_0_0_p2 = (t4011+t4120)*t3997+(z+t4107)*t3991+t4029+
+t4149+t4153;
+      coeffs_dxx->coeff_p1_0_p2 = t4027+t4028;
+      coeffs_dxx->coeff_p2_0_p2 = t3792+x*t3968+t4056+t4151;
+      coeffs_dxx->coeff_m1_p1_p2 = t3818+t3795+t4032+t4101+t4163;
+      coeffs_dxx->coeff_0_p1_p2 = t3926+t3942+t4036+t4086+t4142;
+      coeffs_dxx->coeff_p1_p1_p2 = t4032+t4039;
+      coeffs_dxx->coeff_p2_p1_p2 = t3753+t3765+t3754+t3769+t4097;
+      coeffs_dxx->coeff_m1_p2_p2 = t3980*t4006+t3746+t4051+t4102;
+      coeffs_dxx->coeff_0_p2_p2 = t4106+t3872+t4056+t4092;
+      coeffs_dxx->coeff_p1_p2_p2 = t4051+t4055;
+      coeffs_dxx->coeff_p2_p2_p2 = x*t3858+t3761+t4103;