From d5adbde6b4aea2562cddbb35170ded83cb35a6a6 Mon Sep 17 00:00:00 2001 From: wolf Date: Mon, 16 Sep 2002 22:00:11 +0000 Subject: [PATCH] Implement some things with gradients in 3d for the Nedelec element. git-svn-id: https://svn.dealii.org/trunk@6436 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/deal.II/source/fe/fe_nedelec.cc | 77 +- tests/fe/non_primitive_2.cc | 8 +- tests/fe/non_primitive_2.checked | 1105 +++++++++++++++++++++++ 3 files changed, 1184 insertions(+), 6 deletions(-) create mode 100644 tests/fe/non_primitive_2.checked diff --git a/deal.II/deal.II/source/fe/fe_nedelec.cc b/deal.II/deal.II/source/fe/fe_nedelec.cc index 4ac1f20bee..2f283743cf 100644 --- a/deal.II/deal.II/source/fe/fe_nedelec.cc +++ b/deal.II/deal.II/source/fe/fe_nedelec.cc @@ -296,7 +296,7 @@ FE_Nedelec::shape_value_component (const unsigned int i, case 2: return (component == 0 ? p(1) : 0); // (0, 1-x) case 3: return (component == 0 ? 0 : 1-p(0)); - + // there are // only four // shape @@ -402,7 +402,7 @@ FE_Nedelec::shape_value_component (const unsigned int i, template Tensor<1,dim> FE_Nedelec::shape_grad_component (const unsigned int i, - const Point &/*p*/, + const Point &p, const unsigned int component) const { Assert (idofs_per_cell, ExcIndexRange(i,0,this->dofs_per_cell)); @@ -448,6 +448,79 @@ FE_Nedelec::shape_grad_component (const unsigned int i, }; }; + case 3: // 3d + { + switch (degree) + { + // first order Nedelec + // elements + case 1: + { + // on the unit cell, + // the gradients of + // these shape + // functions are + // linear. we pack them + // into an array, + // knowing that it may + // be expensive to + // recompute the whole + // array each + // time. maybe some + // clever compiler can + // optimize this out, + // seeing that except + // for one element all + // the other ones are + // dead stores... + // + // the format is: first + // index=shape function + // number; second + // index=vector + // component, thrid + // index=component + // within gradient + const double x = p(0), + y = p(1), + z = p(2); + static const double unit_gradients[12][3][3] + = { { {0,-(1-z), -(1-y)}, {0,0,0}, { 0, 0, 0} }, + { {0, 0, 0}, {0,0,0}, { (1-y), -x, 0} }, + { {0, -z, (1-y)}, {0,0,0}, { 0, 0, 0} }, + { {0, 0, 0}, {0,0,0}, {-(1-y), -(1-x), 0} }, + + { {0, (1-z), -y}, {0,0,0}, { 0, 0, 0} }, + { {0, 0, 0}, {0,0,0}, { y, x, 0} }, + { {0, z, y}, {0,0,0}, { 0, 0, 0} }, + { {0, 0, 0}, {0,0,0}, { -y, (1-x), 0} }, + + { {0, 0, 0}, {-(1-z), 0, -(1-x)}, {0, 0, 0} }, + { {0, 0, 0}, { (1-z), 0, -x}, {0, 0, 0} }, + { {0, 0, 0}, { z, 0, x}, {0, 0, 0} }, + { {0, 0, 0}, { -z, 0, (1-x)}, {0, 0, 0} } }; + // note: simple check + // whether this can at + // all be: build the + // sum over all these + // tensors. since the + // sum of the shape + // functions is a + // constant, the + // gradient must + // necessarily be + // zero. this is in + // fact the case here, + // so test successfull + return Tensor<1,dim>(unit_gradients[i][component]); + }; + + // no other degrees + // implemented + default: + Assert (false, ExcNotImplemented()); + }; + }; // presently no other space // dimension implemented default: diff --git a/tests/fe/non_primitive_2.cc b/tests/fe/non_primitive_2.cc index f8f0b88423..52533c176e 100644 --- a/tests/fe/non_primitive_2.cc +++ b/tests/fe/non_primitive_2.cc @@ -226,10 +226,10 @@ test () // for the stokes equation the // pressure does not couple to // itself - std::vector > mask (dim+1, std::vector (dim+1, true)); - mask[dim][dim] = false; +// std::vector > mask (dim+1, std::vector (dim+1, true)); +// mask[dim][dim] = false; - DoFTools::make_sparsity_pattern (dof_handler, mask, sparsity); + DoFTools::make_sparsity_pattern (dof_handler, /*mask,*/ sparsity); sparsity.compress (); SparseMatrix A2 (sparsity); @@ -257,7 +257,7 @@ test () int main () { - std::ofstream logfile("non_primitive_1.output"); + std::ofstream logfile("non_primitive_2.output"); deallog.attach(logfile); deallog.depth_console(0); diff --git a/tests/fe/non_primitive_2.checked b/tests/fe/non_primitive_2.checked new file mode 100644 index 0000000000..9f9714038c --- /dev/null +++ b/tests/fe/non_primitive_2.checked @@ -0,0 +1,1105 @@ + +DEAL::dim=2, n_cells=100 +DEAL::0 0.00000 +DEAL::8 0.00000 +DEAL::16 0.00000 +DEAL::24 0.00000 +DEAL::32 0.00000 +DEAL::40 0.00000 +DEAL::48 0.00000 +DEAL::56 0.00000 +DEAL::64 0.00000 +DEAL::72 0.00000 +DEAL::80 -2.01062 +DEAL::88 0.00000 +DEAL::96 0.00000 +DEAL::104 0.00000 +DEAL::112 0.00000 +DEAL::120 0.00000 +DEAL::128 0.00000 +DEAL::136 0.00000 +DEAL::144 -2.01062 +DEAL::152 0.00000 +DEAL::160 0.00000 +DEAL::168 0.00000 +DEAL::176 0.00000 +DEAL::184 0.00000 +DEAL::192 0.00000 +DEAL::200 2.01062 +DEAL::208 4.02124 +DEAL::216 0.00000 +DEAL::224 0.00000 +DEAL::232 0.00000 +DEAL::240 0.00000 +DEAL::248 0.00000 +DEAL::256 0.00000 +DEAL::264 0.00000 +DEAL::272 2.01062 +DEAL::280 0.00000 +DEAL::288 0.00000 +DEAL::296 0.00000 +DEAL::304 0.00000 +DEAL::312 0.00000 +DEAL::320 0.00000 +DEAL::328 0.00000 +DEAL::336 0.00000 +DEAL::344 0.00000 +DEAL::352 2.01062 +DEAL::360 0.00000 +DEAL::368 0.00000 +DEAL::376 0.00000 +DEAL::384 -2.01062 +DEAL::392 -2.01062 +DEAL::400 0.00000 +DEAL::408 0.00000 +DEAL::416 0.00000 +DEAL::424 0.00000 +DEAL::432 0.00000 +DEAL::440 0.00000 +DEAL::448 0.00000 +DEAL::456 0.00000 +DEAL::464 0.00000 +DEAL::472 0.00000 +DEAL::480 0.00000 +DEAL::488 0.00000 +DEAL::496 0.00000 +DEAL::504 0.00000 +DEAL::512 0.00000 +DEAL::520 0.00000 +DEAL::528 0.00000 +DEAL::536 0.00000 +DEAL::544 0.00000 +DEAL::552 0.00000 +DEAL::560 0.00000 +DEAL::568 0.00000 +DEAL::576 0.00000 +DEAL::584 0.00000 +DEAL::592 -2.01062 +DEAL::600 0.00 +DEAL::608 0.00000 +DEAL::616 0.00000 +DEAL::624 0.00000 +DEAL::632 0.00000 +DEAL::640 0.00000 +DEAL::648 0.00000 +DEAL::656 -2.01062 +DEAL::664 0.00 +DEAL::672 0.00000 +DEAL::680 0.00 +DEAL::688 0.00000 +DEAL::696 0.00000 +DEAL::704 0.00000 +DEAL::712 0.00000 +DEAL::720 0.00000 +DEAL::728 0.00 +DEAL::736 -2.01062 +DEAL::744 0.00 +DEAL::752 0.00000 +DEAL::760 0.00 +DEAL::768 0.00000 +DEAL::776 0.00000 +DEAL::784 0.00 +DEAL::792 0.00000 +DEAL::800 0.00000 +DEAL::808 0.00000 +DEAL::816 0.00000 +DEAL::824 0.00000 +DEAL::832 2.01062 +DEAL::840 4.02124 +DEAL::848 0.00000 +DEAL::856 0.00 +DEAL::864 0.00 +DEAL::872 0.00000 +DEAL::880 0.00000 +DEAL::888 0.00000 +DEAL::896 0.00000 +DEAL::904 0.00000 +DEAL::912 0.00000 +DEAL::920 0.00000 +DEAL::928 0.00000 +DEAL::936 0.00000 +DEAL::944 0.00000 +DEAL::952 0.00000 +DEAL::960 0.00000 +DEAL::968 0.00000 +DEAL::976 0.00 +DEAL::984 4.02124 +DEAL::992 0.00000 +DEAL::1000 0.00000 +DEAL::1008 0.00000 +DEAL::1016 0.00000 +DEAL::1024 0.00000 +DEAL::1032 0.00000 +DEAL::1040 0.00000 +DEAL::1048 -2.01062 +DEAL::1056 0.00000 +DEAL::1064 0.00 +DEAL::1072 0.00000 +DEAL::1080 0.00000 +DEAL::1088 0.00000 +DEAL::1096 -2.01062 +DEAL::1104 0.00 +DEAL::1112 0.00000 +DEAL::1120 0.00000 +DEAL::1128 0.00000 +DEAL::1136 0.00000 +DEAL::1144 0.00000 +DEAL::1152 0.00 +DEAL::1160 0.00 +DEAL::1168 0.00000 +DEAL::1176 0.00 +DEAL::1184 0.00000 +DEAL::1192 0.00000 +DEAL::1200 0.00000 +DEAL::1208 0.00 +DEAL::1216 0.00000 +DEAL::1224 0.00000 +DEAL::1232 0.00000 +DEAL::1240 0.00000 +DEAL::1248 0.00000 +DEAL::1256 0.00000 +DEAL::1264 0.00 +DEAL::1272 -2.01062 +DEAL::1280 0.00000 +DEAL::1288 0.00000 +DEAL::1296 0.00000 +DEAL::1304 0.00 +DEAL::1312 0.00000 +DEAL::1320 0.00000 +DEAL::1328 0.00000 +DEAL::1336 0.00 +DEAL::1344 0.00000 +DEAL::1352 0.00 +DEAL::1360 0.00000 +DEAL::1368 2.01062 +DEAL::1376 0.00000 +DEAL::1384 0.00000 +DEAL::1392 0.00000 +DEAL::1400 0.00000 +DEAL::1408 0.00000 +DEAL::1416 0.00000 +DEAL::1424 0.00000 +DEAL::1432 0.00000 +DEAL::1440 0.00000 +DEAL::1448 0.00000 +DEAL::1456 -2.01062 +DEAL::1464 0.00 +DEAL::1472 0.00000 +DEAL::1480 0.00 +DEAL::1488 0.00000 +DEAL::1496 0.00000 +DEAL::1504 0.00000 +DEAL::1512 0.00000 +DEAL::1520 -2.01062 +DEAL::1528 0.00000 +DEAL::1536 0.00000 +DEAL::1544 0.00 +DEAL::1552 0.00000 +DEAL::1560 0.00000 +DEAL::1568 0.00000 +DEAL::1576 0.00000 +DEAL::1584 0.00000 +DEAL::1592 0.00 +DEAL::1600 -2.01062 +DEAL::1608 0.00 +DEAL::1616 0.00000 +DEAL::1624 0.00 +DEAL::1632 0.00000 +DEAL::1640 0.00000 +DEAL::1648 0.00 +DEAL::1656 0.00000 +DEAL::1664 0.00000 +DEAL::1672 0.00000 +DEAL::1680 0.00000 +DEAL::1688 0.00 +DEAL::1696 4.02124 +DEAL::1704 0.00000 +DEAL::1712 0.00000 +DEAL::1720 0.00000 +DEAL::1728 0.00000 +DEAL::1736 0.00000 +DEAL::1744 0.00000 +DEAL::1752 0.00000 +DEAL::1760 0.00000 +DEAL::1768 0.00000 +DEAL::1776 0.00000 +DEAL::1784 0.00000 +DEAL::1792 0.00000 +DEAL::1800 0.00000 +DEAL::1808 0.00000 +DEAL::1816 0.00000 +DEAL::1824 4.02124 +DEAL::1832 0.00000 +DEAL::1840 0.00 +DEAL::1848 -2.01062 +DEAL::1856 0.00000 +DEAL::1864 0.00000 +DEAL::1872 0.00000 +DEAL::1880 0.00000 +DEAL::1888 0.00000 +DEAL::1896 0.00 +DEAL::1904 0.00 +DEAL::1912 0.00 +DEAL::1920 0.00000 +DEAL::1928 0.00000 +DEAL::1936 0.00000 +DEAL::1944 0.00000 +DEAL::1952 4.02124 +DEAL::1960 0.00000 +DEAL::1968 -2.01062 +DEAL::1976 -2.01062 +DEAL::1984 0.00000 +DEAL::1992 0.00000 +DEAL::2000 0.00000 +DEAL::2008 0.00000 +DEAL::2016 0.00 +DEAL::2024 0.00 +DEAL::2032 0.00000 +DEAL::2040 0.00000 +DEAL::2048 0.00000 +DEAL::2056 0.00 +DEAL::2064 -2.01062 +DEAL::2072 0.00000 +DEAL::2080 0.00000 +DEAL::2088 0.00000 +DEAL::2096 0.00000 +DEAL::2104 0.00000 +DEAL::2112 0.00000 +DEAL::2120 0.00000 +DEAL::2128 0.00 +DEAL::2136 0.00000 +DEAL::2144 0.00000 +DEAL::2152 0.00 +DEAL::2160 0.00000 +DEAL::2168 0.00000 +DEAL::2176 0.00000 +DEAL::2184 0.00000 +DEAL::2192 0.00000 +DEAL::2200 0.00 +DEAL::2208 -2.01062 +DEAL::2216 0.00000 +DEAL::2224 0.00000 +DEAL::2232 4.02124 +DEAL::2240 0.00000 +DEAL::2248 -2.01062 +DEAL::2256 0.00 +DEAL::2264 0.00000 +DEAL::2272 0.00000 +DEAL::2280 0.00000 +DEAL::2288 0.00 +DEAL::2296 0.00000 +DEAL::2304 0.00 +DEAL::2312 0.00000 +DEAL::2320 0.00000 +DEAL::2328 0.00000 +DEAL::2336 0.00000 +DEAL::2344 0.00000 +DEAL::2352 0.00000 +DEAL::2360 4.02124 +DEAL::2368 0.00000 +DEAL::2376 -2.01062 +DEAL::2384 0.00 +DEAL::2392 0.00000 +DEAL::2400 0.00000 +DEAL::2408 0.00000 +DEAL::2416 0.00000 +DEAL::2424 0.00000 +DEAL::2432 0.00 +DEAL::2440 0.00 +DEAL::2448 0.00000 +DEAL::2456 0.00000 +DEAL::2464 -2.01062 +DEAL::2472 0.00 +DEAL::2480 0.00000 +DEAL::2488 0.00000 +DEAL::2496 0.00000 +DEAL::2504 0.00000 +DEAL::2512 0.00000 +DEAL::2520 0.00 +DEAL::2528 0.00 +DEAL::2536 0.00000 +DEAL::2544 0.00000 +DEAL::2552 0.00000 +DEAL::2560 0.00 +DEAL::2568 0.00000 +DEAL::2576 0.00000 +DEAL::2584 0.00000 +DEAL::2592 0.00000 +DEAL::2600 0.00000 +DEAL::2608 0.00 +DEAL::2616 4.02124 +DEAL::2624 0.00000 +DEAL::2632 0.00 +DEAL::2640 -2.01062 +DEAL::2648 0.00000 +DEAL::2656 0.00000 +DEAL::2664 0.00000 +DEAL::2672 0.00000 +DEAL::2680 0.00000 +DEAL::2688 0.00000 +DEAL::2696 0.00000 +DEAL::2704 -2.01062 +DEAL::2712 0.00 +DEAL::2720 0.00000 +DEAL::2728 0.00000 +DEAL::2736 0.00000 +DEAL::2744 0.00 +DEAL::2752 -2.01062 +DEAL::2760 0.00000 +DEAL::2768 0.00000 +DEAL::2776 0.00000 +DEAL::2784 -2.01062 +DEAL::2792 0.00 +DEAL::2800 0.00000 +DEAL::2808 0.00000 +DEAL::2816 0.00 +DEAL::2824 0.00000 +DEAL::2832 0.00000 +DEAL::2840 0.00 +DEAL::2848 0.00000 +DEAL::2856 0.00000 +DEAL::2864 0.00000 +DEAL::2872 0.00000 +DEAL::2880 0.00000 +DEAL::2888 0.00000 +DEAL::2896 0.00000 +DEAL::2904 0.00000 +DEAL::2912 0.00000 +DEAL::2920 0.00000 +DEAL::2928 0.00000 +DEAL::2936 -8.04248 +DEAL::2944 0.00000 +DEAL::2952 0.00000 +DEAL::2960 0.00000 +DEAL::2968 0.00000 +DEAL::2976 0.00000 +DEAL::2984 0.00000 +DEAL::2992 0.00000 +DEAL::3000 -8.04248 +DEAL::3008 0.00000 +DEAL::3016 0.00000 +DEAL::3024 0.00000 +DEAL::3032 0.00000 +DEAL::3040 0.00000 +DEAL::3048 0.00000 +DEAL::3056 0.00000 +DEAL::3064 0.00000 +DEAL::3072 0.00000 +DEAL::3080 0.00000 +DEAL::3088 -8.04248 +DEAL::3096 0.00000 +DEAL::3104 0.00000 +DEAL::3112 0.00000 +DEAL::3120 0.00000 +DEAL::3128 0.00000 +DEAL::3136 0.00000 +DEAL::3144 0.00000 +DEAL::3152 0.00000 +DEAL::3160 0.00000 +DEAL::3168 0.00000 +DEAL::3176 8.04248 +DEAL::3184 16.0850 +DEAL::3192 0.00000 +DEAL::3200 0.00000 +DEAL::3208 0.00000 +DEAL::3216 0.00000 +DEAL::3224 0.00000 +DEAL::3232 0.00000 +DEAL::3240 0.00000 +DEAL::3248 0.00000 +DEAL::3256 0.00000 +DEAL::3264 0.00000 +DEAL::3272 0.00000 +DEAL::3280 -8.04248 +DEAL::3288 0.00000 +DEAL::3296 0.00000 +DEAL::3304 -8.04248 +DEAL::3312 -8.04248 +DEAL::3320 0.00000 +DEAL::3328 0.00000 +DEAL::3336 0.00000 +DEAL::3344 0.00000 +DEAL::3352 0.00000 +DEAL::3360 0.00000 +DEAL::3368 0.00000 +DEAL::3376 0.00000 +DEAL::3384 0.00000 +DEAL::3392 8.04248 +DEAL::3400 0.00000 +DEAL::3408 0.00000 +DEAL::3416 0.00000 +DEAL::3424 0.00000 +DEAL::3432 -8.04248 +DEAL::3440 0.00000 +DEAL::3448 0.00000 +DEAL::3456 0.00000 +DEAL::3464 0.00000 +DEAL::3472 0.00 +DEAL::3480 0.00000 +DEAL::3488 -8.04248 +DEAL::3496 0.00000 +DEAL::3504 0.00000 +DEAL::3512 0.00000 +DEAL::3520 0.00000 +DEAL::3528 0.00000 +DEAL::3536 0.00000 +DEAL::3544 0.00000 +DEAL::3552 0.00000 +DEAL::3560 0.00000 +DEAL::3568 0.00000 +DEAL::3576 0.00000 +DEAL::3584 0.00000 +DEAL::3592 0.00000 +DEAL::3600 0.00000 +DEAL::3608 0.00000 +DEAL::3616 0.00 +DEAL::3624 0.00 +DEAL::3632 0.00000 +DEAL::3640 0.00000 +DEAL::3648 0.00 +DEAL::3656 0.00000 +DEAL::3664 0.00000 +DEAL::3672 8.04248 +DEAL::3680 16.0850 +DEAL::3688 0.00 +DEAL::3696 0.00000 +DEAL::3704 -8.04248 +DEAL::3712 0.00000 +DEAL::3720 0.00000 +DEAL::3728 0.00000 +DEAL::3736 0.00000 +DEAL::3744 0.00000 +DEAL::3752 0.00000 +DEAL::3760 0.00 +DEAL::3768 0.00 +DEAL::3776 0.00000 +DEAL::3784 0.00000 +DEAL::3792 0.00000 +DEAL::3800 16.0850 +DEAL::3808 0.00000 +DEAL::3816 0.00 +DEAL::3824 0.00000 +DEAL::3832 0.00000 +DEAL::3840 0.00000 +DEAL::3848 0.00000 +DEAL::3856 0.00000 +DEAL::3864 0.00 +DEAL::3872 0.00000 +DEAL::3880 0.00000 +DEAL::3888 0.00 +DEAL::3896 0.00000 +DEAL::3904 0.00000 +DEAL::3912 0.00000 +DEAL::3920 0.00000 +DEAL::3928 0.00000 +DEAL::3936 0.00000 +DEAL::3944 16.0850 +DEAL::3952 0.00 +DEAL::3960 0.00 +DEAL::3968 -8.04248 +DEAL::3976 0.00000 +DEAL::3984 0.00000 +DEAL::3992 0.00000 +DEAL::4000 0.00000 +DEAL::4008 0.00000 +DEAL::4016 0.00 +DEAL::4024 0.00 +DEAL::4032 0.00 +DEAL::4040 0.00000 +DEAL::4048 0.00000 +DEAL::4056 0.00000 +DEAL::4064 8.04248 +DEAL::4072 16.0850 +DEAL::4080 0.00 +DEAL::4088 0.00000 +DEAL::4096 0.00000 +DEAL::4104 0.00 +DEAL::4112 0.00 +DEAL::4120 0.00000 +DEAL::4128 0.00000 +DEAL::4136 0.00000 +DEAL::4144 0.00 +DEAL::4152 0.00000 +DEAL::4160 0.00000 +DEAL::4168 0.00000 +DEAL::4176 0.00000 +DEAL::4184 0.00000 +DEAL::4192 0.00 +DEAL::4200 0.00000 +DEAL::4208 8.04248 +DEAL::4216 16.0850 +DEAL::4224 0.00000 +DEAL::4232 0.00000 +DEAL::4240 0.00000 +DEAL::4248 0.00000 +DEAL::4256 0.00000 +DEAL::4264 0.00000 +DEAL::4272 0.00000 +DEAL::4280 0.00000 +DEAL::4288 0.00000 +DEAL::4296 0.00000 +DEAL::4304 0.00000 +DEAL::4312 0.00000 +DEAL::4320 0.00000 +DEAL::4328 0.00000 +DEAL::4336 0.00000 +DEAL::4344 0.00 +DEAL::4352 -32.1699 +DEAL::4360 32.1699 +DEAL::4368 0.00000 +DEAL::4376 0.00000 +DEAL::4384 0.00000 +DEAL::4392 0.00000 +DEAL::4400 0.00000 +DEAL::4408 -32.1699 +DEAL::4416 -32.1699 +DEAL::4424 0.00000 +DEAL::4432 0.00000 +DEAL::4440 0.00000 +DEAL::4448 0.00000 +DEAL::4456 0.00000 +DEAL::4464 0.00000 +DEAL::4472 -32.1699 +DEAL::4480 0.00000 +DEAL::4488 0.00000 +DEAL::4496 0.00 +DEAL::4504 0.00000 +DEAL::4512 0.00000 +DEAL::4520 -32.1699 +DEAL::4528 -32.1699 +DEAL::4536 0.00000 +DEAL::4544 0.00000 +DEAL::4552 0.00000 +DEAL::4560 0.00000 +DEAL::4568 0.00000 +DEAL::4576 0.00000 +DEAL::4584 -32.1699 +DEAL::4592 0.00 +DEAL::4600 0.00000 +DEAL::4608 0.00000 +DEAL::4616 0.00000 +DEAL::4624 0.00000 +DEAL::4632 0.00000 +DEAL::4640 0.00 +DEAL::4648 0.00000 +DEAL::4656 0.00000 +DEAL::4664 0.00000 +DEAL::4672 0.00000 +DEAL::4680 0.00000 +DEAL::4688 0.00 +DEAL::4696 0.00 +DEAL::4704 0.00000 +DEAL::4712 0.00000 +DEAL::4720 0.00000 +DEAL::4728 0.00000 +DEAL::4736 0.00000 +DEAL::4744 0.00000 +DEAL::4752 0.00000 +DEAL::4760 0.00000 +DEAL::4768 0.00000 +DEAL::4776 0.00000 +DEAL::4784 64.3398 +DEAL::4792 -32.1699 +DEAL::4800 0.00000 +DEAL::4808 0.00000 +DEAL::4816 0.00000 +DEAL::4824 0.00 +DEAL::4832 0.00000 +DEAL::4840 0.00000 +DEAL::4848 0.00 +DEAL::4856 -32.1699 +DEAL::4864 0.00000 +DEAL::4872 -32.1699 +DEAL::4880 0.00 +DEAL::4888 0.00 +DEAL::4896 0.00000 +DEAL::4904 0.00000 +DEAL::4912 0.00000 +DEAL::4920 0.00000 +DEAL::4928 64.3398 +DEAL::4936 0.00000 +DEAL::4944 0.00000 +DEAL::4952 -32.1699 +DEAL::4960 -32.1699 +DEAL::4968 0.00 +DEAL::4976 0.00000 +DEAL::4984 0.00000 +DEAL::4992 0.00000 +DEAL::5000 -32.1699 +DEAL::5008 -32.1699 +DEAL::5016 0.00 +DEAL::5024 0.00 +DEAL::dim=3, n_cells=43 +DEAL::0 0.00000 +DEAL::27 0.00000 +DEAL::54 0.00000 +DEAL::81 0.00000 +DEAL::108 0.00000 +DEAL::135 0.00000 +DEAL::162 0.00000 +DEAL::189 0.00000 +DEAL::216 0.00000 +DEAL::243 0.00000 +DEAL::270 0.00000 +DEAL::297 0.00000 +DEAL::324 0.00000 +DEAL::351 0.00000 +DEAL::378 0.00000 +DEAL::405 0.00000 +DEAL::432 -0.0431842 +DEAL::459 0.00 +DEAL::486 -0.0125664 +DEAL::513 0.00000 +DEAL::540 -0.0376991 +DEAL::567 0.00 +DEAL::594 0.00 +DEAL::621 0.00 +DEAL::648 0.00000 +DEAL::675 0.00 +DEAL::702 0.00000 +DEAL::729 -0.0431842 +DEAL::756 0.00000 +DEAL::783 0.00000 +DEAL::810 0.00000 +DEAL::837 0.00000 +DEAL::864 0.00000 +DEAL::891 0.00000 +DEAL::918 0.0502655 +DEAL::945 0.00 +DEAL::972 0.00 +DEAL::999 0.00000 +DEAL::1026 0.00000 +DEAL::1053 0.00 +DEAL::1080 0.00 +DEAL::1107 0.00000 +DEAL::1134 0.00000 +DEAL::1161 0.00000 +DEAL::1188 0.00000 +DEAL::1215 0.00000 +DEAL::1242 0.00000 +DEAL::1269 0.00000 +DEAL::1296 0.00 +DEAL::1323 0.00000 +DEAL::1350 0.00000 +DEAL::1377 -0.0431842 +DEAL::1404 -0.0376991 +DEAL::1431 0.00 +DEAL::1458 0.00000 +DEAL::1485 0.00000 +DEAL::1512 0.00000 +DEAL::1539 0.00 +DEAL::1566 0.00 +DEAL::1593 0.00000 +DEAL::1620 0.00000 +DEAL::1647 0.00000 +DEAL::1674 0.00000 +DEAL::1701 0.00000 +DEAL::1728 0.00000 +DEAL::1755 0.00000 +DEAL::1782 0.00000 +DEAL::1809 0.00000 +DEAL::1836 0.00000 +DEAL::1863 0.00000 +DEAL::1890 0.00000 +DEAL::1917 0.00000 +DEAL::1944 0.00000 +DEAL::1971 0.00000 +DEAL::1998 0.00 +DEAL::2025 -0.0863685 +DEAL::2052 0.00000 +DEAL::2079 -0.0863685 +DEAL::2106 0.00000 +DEAL::2133 -0.0753982 +DEAL::2160 0.00 +DEAL::2187 0.00000 +DEAL::2214 0.00000 +DEAL::2241 0.00 +DEAL::2268 0.00000 +DEAL::2295 0.402124 +DEAL::2322 0.00 +DEAL::2349 0.00000 +DEAL::2376 0.00 +DEAL::2403 0.00000 +DEAL::2430 0.00000 +DEAL::2457 0.00000 +DEAL::2484 0.00000 +DEAL::2511 0.00 +DEAL::2538 0.00000 +DEAL::2565 0.298401 +DEAL::2592 -0.0863685 +DEAL::2619 0.00000 +DEAL::2646 0.00 +DEAL::2673 0.00 +DEAL::2700 0.103723 +DEAL::2727 0.00 +DEAL::2754 0.00000 +DEAL::2781 0.00000 +DEAL::2808 0.00000 +DEAL::2835 0.00000 +DEAL::2862 0.00000 +DEAL::2889 0.00000 +DEAL::2916 0.00000 +DEAL::2943 0.00 +DEAL::2970 0.00000 +DEAL::2997 0.00000 +DEAL::3024 0.00 +DEAL::3051 0.00 +DEAL::3078 0.00000 +DEAL::3105 0.00 +DEAL::3132 0.00000 +DEAL::3159 0.00 +DEAL::3186 0.00 +DEAL::3213 0.00 +DEAL::3240 0.00000 +DEAL::3267 0.00000 +DEAL::3294 0.00000 +DEAL::3321 0.00000 +DEAL::3348 0.00 +DEAL::3375 0.00000 +DEAL::3402 0.00 +DEAL::3429 0.00000 +DEAL::3456 0.00000 +DEAL::3483 0.00000 +DEAL::3510 0.00000 +DEAL::3537 0.00000 +DEAL::3564 0.00000 +DEAL::3591 0.00000 +DEAL::3618 0.00000 +DEAL::3645 0.00000 +DEAL::3672 0.00000 +DEAL::3699 0.00000 +DEAL::3726 0.00 +DEAL::3753 0.00 +DEAL::3780 -0.0251327 +DEAL::3807 0.00 +DEAL::3834 0.00 +DEAL::3861 -0.0251327 +DEAL::3888 0.298401 +DEAL::3915 0.00 +DEAL::3942 0.00000 +DEAL::3969 0.00 +DEAL::3996 0.00000 +DEAL::4023 0.00000 +DEAL::4050 0.00000 +DEAL::4077 0.00000 +DEAL::4104 0.00 +DEAL::4131 0.00000 +DEAL::4158 0.00 +DEAL::4185 0.00 +DEAL::4212 0.00000 +DEAL::4239 0.00000 +DEAL::4266 0.00000 +DEAL::4293 0.00000 +DEAL::4320 0.00000 +DEAL::4347 0.00 +DEAL::4374 0.00000 +DEAL::4401 0.00 +DEAL::4428 0.00000 +DEAL::4455 0.00 +DEAL::4482 0.00 +DEAL::4509 0.00000 +DEAL::4536 0.00000 +DEAL::4563 0.00000 +DEAL::4590 0.00000 +DEAL::4617 0.00000 +DEAL::4644 0.00000 +DEAL::4671 0.00000 +DEAL::4698 0.00000 +DEAL::4725 0.00000 +DEAL::4752 0.00000 +DEAL::4779 0.00 +DEAL::4806 0.00 +DEAL::4833 0.00 +DEAL::4860 0.00000 +DEAL::4887 0.00000 +DEAL::4914 0.00000 +DEAL::4941 0.00 +DEAL::4968 0.00000 +DEAL::4995 0.00000 +DEAL::5022 0.00000 +DEAL::5049 0.00000 +DEAL::5076 0.00000 +DEAL::5103 0.00000 +DEAL::5130 0.00000 +DEAL::5157 0.00000 +DEAL::5184 0.402124 +DEAL::5211 0.00000 +DEAL::5238 0.00000 +DEAL::5265 0.0109702 +DEAL::5292 0.00000 +DEAL::5319 0.00000 +DEAL::5346 0.00000 +DEAL::5373 0.00000 +DEAL::5400 0.00 +DEAL::5427 0.00000 +DEAL::5454 0.103723 +DEAL::5481 0.0109702 +DEAL::5508 0.00 +DEAL::5535 0.00000 +DEAL::5562 -0.0753982 +DEAL::5589 0.103723 +DEAL::5616 0.00 +DEAL::5643 0.00000 +DEAL::5670 0.00000 +DEAL::5697 0.00000 +DEAL::5724 0.00 +DEAL::5751 0.00 +DEAL::5778 0.00000 +DEAL::5805 0.00000 +DEAL::5832 0.00000 +DEAL::5859 0.00000 +DEAL::5886 0.00000 +DEAL::5913 0.00000 +DEAL::5940 0.00000 +DEAL::5967 0.00000 +DEAL::5994 0.00000 +DEAL::6021 0.00000 +DEAL::6048 0.00000 +DEAL::6075 0.00 +DEAL::6102 0.00 +DEAL::6129 0.00 +DEAL::6156 0.00 +DEAL::6183 0.00000 +DEAL::6210 0.00000 +DEAL::6237 0.00000 +DEAL::6264 -0.0863685 +DEAL::6291 0.00000 +DEAL::6318 0.00000 +DEAL::6345 0.00000 +DEAL::6372 0.00000 +DEAL::6399 0.00 +DEAL::6426 0.00 +DEAL::6453 0.00000 +DEAL::6480 0.00000 +DEAL::6507 -0.0753982 +DEAL::6534 0.00000 +DEAL::6561 0.00000 +DEAL::6588 0.00 +DEAL::6615 0.00000 +DEAL::6642 0.00 +DEAL::6669 0.00000 +DEAL::6696 0.0109702 +DEAL::6723 0.00000 +DEAL::6750 0.00 +DEAL::6777 0.00000 +DEAL::6804 0.00000 +DEAL::6831 0.00000 +DEAL::6858 0.00000 +DEAL::6885 0.00000 +DEAL::6912 0.00000 +DEAL::6939 0.00000 +DEAL::6966 0.00000 +DEAL::6993 -0.0863685 +DEAL::7020 0.402124 +DEAL::7047 0.00000 +DEAL::7074 0.00000 +DEAL::7101 0.00 +DEAL::7128 -0.0251327 +DEAL::7155 0.00000 +DEAL::7182 0.00 +DEAL::7209 0.00000 +DEAL::7236 -0.0753982 +DEAL::7263 0.00 +DEAL::7290 0.00000 +DEAL::7317 0.00000 +DEAL::7344 0.00000 +DEAL::7371 0.00000 +DEAL::7398 0.00000 +DEAL::7425 0.00000 +DEAL::7452 0.00000 +DEAL::7479 0.00 +DEAL::7506 0.0109702 +DEAL::7533 0.00000 +DEAL::7560 0.00000 +DEAL::7587 0.00000 +DEAL::7614 0.00000 +DEAL::7641 0.00000 +DEAL::7668 0.00000 +DEAL::7695 0.00 +DEAL::7722 0.00 +DEAL::7749 0.00 +DEAL::7776 0.00000 +DEAL::7803 -0.0251327 +DEAL::7830 0.00000 +DEAL::7857 0.00 +DEAL::7884 0.00000 +DEAL::7911 0.00000 +DEAL::7938 0.00000 +DEAL::7965 0.00000 +DEAL::7992 0.00000 +DEAL::8019 0.00000 +DEAL::8046 0.00000 +DEAL::8073 -0.0863685 +DEAL::8100 0.00000 +DEAL::8127 0.00000 +DEAL::8154 0.00000 +DEAL::8181 0.00000 +DEAL::8208 0.00000 +DEAL::8235 0.00 +DEAL::8262 0.00000 +DEAL::8289 0.00000 +DEAL::8316 0.00000 +DEAL::8343 0.00000 +DEAL::8370 0.00000 +DEAL::8397 -0.0251327 +DEAL::8424 0.00 +DEAL::8451 0.00000 +DEAL::8478 0.00000 +DEAL::8505 0.00000 +DEAL::8532 0.00000 +DEAL::8559 0.00 +DEAL::8586 0.00 +DEAL::8613 0.00000 +DEAL::8640 0.00000 +DEAL::8667 0.00000 +DEAL::8694 0.00 +DEAL::8721 -0.0251327 +DEAL::8748 0.00 +DEAL::8775 0.00 +DEAL::8802 -0.0753982 +DEAL::8829 0.00000 +DEAL::8856 0.00000 +DEAL::8883 0.00 +DEAL::8910 0.00000 +DEAL::8937 -0.0863685 +DEAL::8964 0.00000 +DEAL::8991 0.00000 +DEAL::9018 0.00000 +DEAL::9045 0.00000 +DEAL::9072 0.00000 +DEAL::9099 0.00 +DEAL::9126 -0.0251327 +DEAL::9153 0.00 +DEAL::9180 -0.0863685 +DEAL::9207 0.00000 +DEAL::9234 0.00000 +DEAL::9261 0.00 +DEAL::9288 0.00000 +DEAL::9315 0.00000 +DEAL::9342 0.00 +DEAL::9369 0.00 +DEAL::9396 0.00 +DEAL::9423 0.00000 +DEAL::9450 0.00 +DEAL::9477 0.00000 +DEAL::9504 0.00000 +DEAL::9531 0.00000 +DEAL::9558 0.00 +DEAL::9585 0.00 +DEAL::9612 0.00000 +DEAL::9639 0.00 +DEAL::9666 0.00000 +DEAL::9693 0.00000 +DEAL::9720 0.00000 +DEAL::9747 0.00 +DEAL::9774 0.00 +DEAL::9801 0.00000 +DEAL::9828 0.00000 +DEAL::9855 0.00000 +DEAL::9882 0.00000 +DEAL::9909 0.00000 +DEAL::9936 0.00000 +DEAL::9963 0.0109702 +DEAL::9990 0.00000 +DEAL::10017 0.00000 +DEAL::10044 0.00 +DEAL::10071 0.00000 +DEAL::10098 0.00000 +DEAL::10125 0.00000 +DEAL::10152 0.00 +DEAL::10179 0.00000 +DEAL::10206 0.00000 +DEAL::10233 0.00 +DEAL::10260 0.00000 +DEAL::10287 0.0109702 +DEAL::10314 0.00000 +DEAL::10341 0.00000 +DEAL::10368 0.0109702 +DEAL::10395 0.00000 +DEAL::10422 0.00000 +DEAL::10449 0.00000 +DEAL::10476 0.00000 +DEAL::10503 0.00000 +DEAL::10530 0.00000 +DEAL::10557 0.00000 +DEAL::10584 0.00000 +DEAL::10611 0.00000 +DEAL::10638 0.00000 +DEAL::10665 0.00000 +DEAL::10692 0.00000 +DEAL::10719 0.00000 +DEAL::10746 0.00000 +DEAL::10773 0.00000 +DEAL::10800 0.00 +DEAL::10827 0.00 +DEAL::10854 0.00 +DEAL::10881 -0.0502655 +DEAL::10908 0.00 +DEAL::10935 0.00 +DEAL::10962 0.00000 +DEAL::10989 0.596801 +DEAL::11016 -0.172737 +DEAL::11043 0.00000 +DEAL::11070 0.00 +DEAL::11097 0.00000 +DEAL::11124 0.596801 +DEAL::11151 0.00 +DEAL::11178 0.00000 +DEAL::11205 0.00000 +DEAL::11232 0.00000 +DEAL::11259 0.00000 +DEAL::11286 -0.172737 +DEAL::11313 0.00000 +DEAL::11340 0.00 +DEAL::11367 0.00 +DEAL::11394 0.00 +DEAL::11421 0.00000 +DEAL::11448 0.00 +DEAL::11475 0.00 +DEAL::11502 0.00 +DEAL::11529 0.00000 +DEAL::11556 0.00000 +DEAL::11583 0.00000 +DEAL::11610 0.00000 +DEAL::11637 0.00000 +DEAL::11664 0.00 +DEAL::11691 0.00000 +DEAL::11718 0.00000 +DEAL::11745 0.00 +DEAL::11772 0.00 +DEAL::11799 0.00 +DEAL::11826 0.00 +DEAL::11853 0.00 +DEAL::11880 0.00000 +DEAL::11907 0.00000 +DEAL::11934 -0.172737 +DEAL::11961 0.00 +DEAL::11988 0.00 +DEAL::12015 0.00000 +DEAL::12042 0.00000 +DEAL::12069 0.00000 +DEAL::12096 0.00000 +DEAL::12123 0.00000 +DEAL::12150 0.00 +DEAL::12177 0.00 +DEAL::12204 0.00000 +DEAL::12231 0.00000 +DEAL::12258 -0.172737 +DEAL::12285 0.00000 +DEAL::12312 0.00 +DEAL::12339 0.00 +DEAL::12366 -0.0502655 +DEAL::12393 0.00000 +DEAL::12420 0.00 +DEAL::12447 0.00000 +DEAL::12474 0.0219405 +DEAL::12501 0.00000 +DEAL::12528 0.0219405 +DEAL::12555 0.00000 +DEAL::12582 0.00638455 +DEAL::12609 0.00 +DEAL::12636 0.00000 +DEAL::12663 0.00000 +DEAL::12690 -0.0502655 +DEAL::12717 0.00000 +DEAL::12744 0.00 -- 2.39.5