From d9e6b29f46e0e608e97274a7c6c29d54405064ae Mon Sep 17 00:00:00 2001 From: Anton Khirnov Date: Fri, 26 Feb 2016 16:29:02 +0100 Subject: Add 5th order Largrange tensor product interpolation. Changes to other coeffs as well due to a different Maple version used. --- .../3d.cube.order2.smooth0/coeffs-dyz.compute.c | 276 +++++++++++---------- 1 file changed, 141 insertions(+), 135 deletions(-) (limited to 'src/Lagrange-tensor-product/3d.coeffs/3d.cube.order2.smooth0/coeffs-dyz.compute.c') diff --git a/src/Lagrange-tensor-product/3d.coeffs/3d.cube.order2.smooth0/coeffs-dyz.compute.c b/src/Lagrange-tensor-product/3d.coeffs/3d.cube.order2.smooth0/coeffs-dyz.compute.c index a17d581..8ed8838 100644 --- a/src/Lagrange-tensor-product/3d.coeffs/3d.cube.order2.smooth0/coeffs-dyz.compute.c +++ b/src/Lagrange-tensor-product/3d.coeffs/3d.cube.order2.smooth0/coeffs-dyz.compute.c @@ -1,164 +1,170 @@ -fp t716; -fp t720; -fp t714; -fp t737; -fp t694; -fp t682; -fp t756; -fp t711; -fp t740; -fp t693; -fp t710; -fp t742; -fp t696; -fp t755; fp t715; +fp t722; +fp t736; +fp t689; +fp t682; +fp t713; fp t754; -fp t705; +fp t716; fp t753; +fp t705; fp t691; +fp t683; +fp t714; +fp t704; +fp t685; fp t752; fp t717; -fp t738; -fp t692; -fp t681; -fp t708; +fp t709; +fp t693; fp t751; +fp t718; +fp t710; fp t699; -fp t750; -fp t704; -fp t739; fp t684; -fp t680; +fp t750; fp t749; -fp t697; -fp t701; -fp t688; +fp t686; +fp t703; fp t748; -fp t695; +fp t708; fp t747; fp t746; -fp t698; -fp t686; +fp t690; +fp t711; fp t745; -fp t735; -fp t687; +fp t694; fp t744; fp t743; -fp t719; +fp t735; +fp t697; +fp t742; +fp t698; +fp t706; fp t741; +fp t721; +fp t737; fp t733; -fp t679; -fp t690; -fp t736; +fp t740; +fp t739; +fp t738; +fp t720; +fp t712; +fp t681; +fp t701; fp t734; fp t732; fp t731; fp t730; fp t729; +fp t687; +fp t688; fp t728; +fp t696; fp t727; +fp t695; fp t726; fp t725; fp t724; fp t723; -fp t722; -fp t721; -fp t718; -fp t689; - t716 = RATIONAL(1.0,8.0); - t720 = x*x; - t714 = RATIONAL(1.0,2.0); - t737 = t714*t720; - t694 = y*t737; - t682 = z*t694; - t756 = t682+t716*t720; - t711 = RATIONAL(-1.0,4.0); - t740 = x*t711; - t693 = y*t740; - t710 = RATIONAL(1.0,4.0); - t742 = x*t710; - t696 = z*t742; - t755 = t693+t696; - t715 = RATIONAL(-1.0,8.0); - t754 = t715*t720+t682; - t705 = t711*t720; - t753 = t705+t710; - t691 = z*t740; - t752 = t693+t691; - t717 = RATIONAL(-1.0,2.0); - t738 = x*t717; - t692 = y*t738; - t681 = z*t692; - t708 = t715*x; - t751 = t681+t708; - t699 = y*t742; - t750 = t696+t699; - t704 = t716*x; - t739 = x*t714; - t684 = z*t739; - t680 = y*t684; - t749 = t704+t680; - t697 = y*t705; - t701 = t710*t720; - t688 = z*t701; - t748 = t697+t688; - t695 = z*t737; - t747 = t717*z+t695; - t746 = t711+t701; - t698 = z*t705; - t686 = y*t701; - t745 = t698+t686; - t735 = t717*t720; - t687 = z*t735; - t744 = t714*z+t687; - t743 = t699+t691; - t719 = y*z; - t741 = y*t720; - t733 = RATIONAL(2.0,1.0)*t719; - t679 = t720*t733; - t690 = RATIONAL(-2.0,1.0)*t719; - t736 = t679+t690; - t734 = z*t741; - t732 = (t720-1.0)*z; - t731 = -y+t741; - t730 = -t734+z*t738; - t729 = -t734+t684; - t728 = t694-t734; - t727 = y*t735-t734; - t726 = t704+t681+t754; - t725 = t686+t688+t756; - t724 = t697+t698+t756; - t723 = t680+t708+t754; - t722 = t719+t714*y+t727; - t721 = t719+t717*y+t728; - t718 = x*t719; - t689 = y*t739; - coeffs_dyz->coeff_m1_m1_m1 = t724+t750+t751; - coeffs_dyz->coeff_0_m1_m1 = t721+t747+t753; - coeffs_dyz->coeff_p1_m1_m1 = t724+t749+t752; - coeffs_dyz->coeff_m1_0_m1 = t692+t718+t728; - coeffs_dyz->coeff_0_0_m1 = -t731+t736; - coeffs_dyz->coeff_p1_0_m1 = t689-t718+t728; - coeffs_dyz->coeff_m1_p1_m1 = t726+t743+t748; - coeffs_dyz->coeff_0_p1_m1 = t721+t744+t746; - coeffs_dyz->coeff_p1_p1_m1 = t723+t748+t755; - coeffs_dyz->coeff_m1_m1_0 = t695+t718+t730; - coeffs_dyz->coeff_0_m1_0 = -t732+t736; - coeffs_dyz->coeff_p1_m1_0 = t695-t718+t729; - coeffs_dyz->coeff_m1_0_0 = x*t690+t679; - coeffs_dyz->coeff_0_0_0 = (RATIONAL(4.0,1.0)+RATIONAL(-4.0,1.0)*t720)* -t719; - coeffs_dyz->coeff_p1_0_0 = x*t733+t679; - coeffs_dyz->coeff_m1_p1_0 = t687+t718+t729; - coeffs_dyz->coeff_0_p1_0 = t732+t736; - coeffs_dyz->coeff_p1_p1_0 = t687-t718+t730; - coeffs_dyz->coeff_m1_m1_p1 = t726+t745+t755; - coeffs_dyz->coeff_0_m1_p1 = t722+t746+t747; - coeffs_dyz->coeff_p1_m1_p1 = t723+t743+t745; - coeffs_dyz->coeff_m1_0_p1 = t689+t718+t727; - coeffs_dyz->coeff_0_0_p1 = t731+t736; - coeffs_dyz->coeff_p1_0_p1 = t692-t718+t727; - coeffs_dyz->coeff_m1_p1_p1 = t725+t751+t752; - coeffs_dyz->coeff_0_p1_p1 = t722+t744+t753; - coeffs_dyz->coeff_p1_p1_p1 = t725+t749+t750; +fp t700; +fp t692; + t715 = RATIONAL(1.0,2.0); + t722 = x*x; + t736 = t715*t722; + t689 = z*t736; + t682 = y*t689; + t713 = RATIONAL(1.0,8.0); + t754 = t682+t713*t722; + t716 = RATIONAL(-1.0,8.0); + t753 = t682+t716*t722; + t705 = t715*z; + t691 = x*t705; + t683 = y*t691; + t714 = RATIONAL(1.0,4.0); + t704 = t714*t722; + t685 = z*t704; + t752 = t683+t685; + t717 = RATIONAL(-1.0,4.0); + t709 = t717*t722; + t693 = z*t709; + t751 = t683+t693; + t718 = RATIONAL(-1.0,2.0); + t710 = t718*z; + t699 = x*t710; + t684 = y*t699; + t750 = t684+t685; + t749 = t684+t693; + t686 = y*t704; + t703 = t713*x; + t748 = t686+t703; + t708 = t716*x; + t747 = t686+t708; + t746 = t689+t710; + t690 = y*t736; + t711 = t718*y; + t745 = t690+t711; + t694 = y*t709; + t744 = t694+t703; + t743 = t694+t708; + t735 = t718*t722; + t697 = z*t735; + t742 = t697+t705; + t698 = y*t735; + t706 = t715*y; + t741 = t698+t706; + t721 = y*z; + t737 = t722*z; + t733 = y*t737; + t740 = t721-t733; + t739 = t714*x; + t738 = t717*x; + t720 = x*t721; + t712 = RATIONAL(2.0,1.0); + t681 = t712*t733; + t701 = RATIONAL(-2.0,1.0)*t721; + t734 = t681+t701; + t732 = (-1.0+t722)*y; + t731 = t737-z; + t730 = -t720-t733; + t729 = t720-t733; + t687 = z*t739; + t688 = y*t739; + t728 = t687+t688+t754; + t696 = y*t738; + t727 = t687+t696+t753; + t695 = z*t738; + t726 = t688+t695+t753; + t725 = t695+t696+t754; + t724 = t704+t717+t740; + t723 = t709+t714+t740; + t700 = x*t711; + t692 = x*t706; + coeffs_dyz->coeff_m1_m1_m1 = t728+t743+t749; + coeffs_dyz->coeff_0_m1_m1 = t723+t745+t746; + coeffs_dyz->coeff_p1_m1_m1 = t725+t744+t751; + coeffs_dyz->coeff_m1_0_m1 = t700+t690+t729; + coeffs_dyz->coeff_0_0_m1 = -t732+t734; + coeffs_dyz->coeff_p1_0_m1 = t692+t690+t730; + coeffs_dyz->coeff_m1_p1_m1 = t726+t744+t750; + coeffs_dyz->coeff_0_p1_m1 = t724+t742+t745; + coeffs_dyz->coeff_p1_p1_m1 = t727+t743+t752; + coeffs_dyz->coeff_m1_m1_0 = t689+t699+t729; + coeffs_dyz->coeff_0_m1_0 = -t731+t734; + coeffs_dyz->coeff_p1_m1_0 = t691+t689+t730; + coeffs_dyz->coeff_m1_0_0 = t681+x*t701; + coeffs_dyz->coeff_0_0_0 = (RATIONAL(4.0,1.0)+RATIONAL(-4.0,1.0)*t722)* +t721; + coeffs_dyz->coeff_p1_0_0 = t712*t720+t681; + coeffs_dyz->coeff_m1_p1_0 = t691+t697+t729; + coeffs_dyz->coeff_0_p1_0 = t731+t734; + coeffs_dyz->coeff_p1_p1_0 = t699+t697+t730; + coeffs_dyz->coeff_m1_m1_p1 = t727+t748+t749; + coeffs_dyz->coeff_0_m1_p1 = t724+t741+t746; + coeffs_dyz->coeff_p1_m1_p1 = t726+t747+t751; + coeffs_dyz->coeff_m1_0_p1 = t692+t698+t729; + coeffs_dyz->coeff_0_0_p1 = t732+t734; + coeffs_dyz->coeff_p1_0_p1 = t700+t698+t730; + coeffs_dyz->coeff_m1_p1_p1 = t725+t747+t750; + coeffs_dyz->coeff_0_p1_p1 = t723+t741+t742; + coeffs_dyz->coeff_p1_p1_p1 = t728+t748+t752; -- cgit v1.2.3