<div dir="ltr"><div dir="ltr">On Wed, Jul 3, 2019 at 1:03 PM Fabian.Jakub <<a href="mailto:Fabian.Jakub@physik.uni-muenchen.de">Fabian.Jakub@physik.uni-muenchen.de</a>> wrote:<br></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Found the change...<br>
it was in Commit 6d563f2afc7bcefff9dc7869193a6debf1d639a5<br>
Date:   Wed Jun 26 12:32:01 2019 +0000<br>
<br>
    remove remnants of non-existing<br>
DMPlexRestore{Cone,ConeOrientation,Support}<br></blockquote><div><br></div><div>You are right!</div><div><br></div><div>Vaclav, you cannot remove the Fortran RestoreCone() because it is actually doing something.</div><div>Fortran needs to create arrays, and this code deallocates them. We do not need that in C.</div><div>Can you put those back?</div><div><br></div><div>   Thanks,</div><div><br></div><div>      Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
The patch removed the fortran stubs but not the C functions.<br>
In this case I guess that the fortran compiler reverts to the default<br>
calling convention and thus does not nullify the cone pointers but<br>
rather nulls the values in the cone.<br>
<br>
I suggest you remove the C functions to get a compiler error or revert<br>
the path.<br>
<br>
Many thanks,<br>
<br>
Fabian<br>
<br>
On 7/3/19 3:18 PM, Matthew Knepley wrote:<br>
> On Wed, Jul 3, 2019 at 7:38 AM Fabian.Jakub via petsc-users <<br>
> <a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a>> wrote:<br>
> <br>
>> Dear Matt, Dear Petsc<br>
>><br>
>> After my recent pull of petsc master I encounter test fails that use<br>
>> dmplex meshes.<br>
>><br>
>> I attached a simple script to show an issue where DMPlexRestoreCone<br>
>> mangles the cone relationships in the DM.<br>
>><br>
>> The behaviour happens with gcc 7.4.0 aswell as icc 19.0.3.199<br>
>> @ current petsc/master d3dc048301fdb3e4c282b528e7b28910551cfe13<br>
>><br>
>> We recently had a system update and I am not 100% sure if this is some<br>
>> weird stuff going on at my end or if this is indeed Plex related or if I<br>
>> should not call Restore in the inner loop.<br>
>><br>
>> Anyway, I would be most grateful if you could try to reproduce the error<br>
>> and hint me towards a solution.<br>
>><br>
> <br>
> Hi Fabian,<br>
> <br>
> I think it might be your system update. Seems to run fine for me:<br>
> <br>
> master *:~/Downloads/tmp$ ./plex_prism<br>
> <br>
>  Cone           0 :           1           2           3           4<br>
>   5<br>
> <br>
>  Cone           1 :           6           7           8<br>
> <br>
>  Cone           2 :          12          13          14<br>
> <br>
>  Cone           3 :           7           9          11          13<br>
> <br>
>  Cone           4 :           6           9          10          12<br>
> <br>
>  Cone           5 :           8          10          11          14<br>
> <br>
>  Cone           6 :          15          16<br>
> <br>
>  Cone           7 :          15          17<br>
> <br>
>  Cone           8 :          16          17<br>
> <br>
>  Cone           9 :          15          18<br>
> <br>
>  Cone          10 :          16          19<br>
> <br>
>  Cone          11 :          17          20<br>
> <br>
>  Cone          12 :          18          19<br>
> <br>
>  Cone          13 :          18          20<br>
> <br>
>  Cone          14 :          19          20<br>
> <br>
>  Cone          15 :<br>
> <br>
>  Cone          16 :<br>
> <br>
>  Cone          17 :<br>
> <br>
>  Cone          18 :<br>
> <br>
>  Cone          19 :<br>
> <br>
>  Cone          20 :<br>
> <br>
> master *:~/Downloads/tmp$ ./plex_prism -show_plex ::ascii_info_detail<br>
> <br>
> DM Object: testplex 1 MPI processes<br>
> <br>
>   type: plex<br>
> <br>
> testplex in 3 dimensions:<br>
> <br>
> Supports:<br>
> <br>
> [0] Max support size: 3<br>
> <br>
> [0]: 1 ----> 0<br>
> <br>
> [0]: 2 ----> 0<br>
> <br>
> [0]: 3 ----> 0<br>
> <br>
> [0]: 4 ----> 0<br>
> <br>
> [0]: 5 ----> 0<br>
> <br>
> [0]: 6 ----> 1<br>
> <br>
> [0]: 6 ----> 4<br>
> <br>
> [0]: 7 ----> 1<br>
> <br>
> [0]: 7 ----> 3<br>
> <br>
> [0]: 8 ----> 1<br>
> <br>
> [0]: 8 ----> 5<br>
> <br>
> [0]: 9 ----> 3<br>
> <br>
> [0]: 9 ----> 4<br>
> <br>
> [0]: 10 ----> 4<br>
> <br>
> [0]: 10 ----> 5<br>
> <br>
> [0]: 11 ----> 3<br>
> <br>
> [0]: 11 ----> 5<br>
> <br>
> [0]: 12 ----> 2<br>
> <br>
> [0]: 12 ----> 4<br>
> <br>
> [0]: 13 ----> 2<br>
> <br>
> [0]: 13 ----> 3<br>
> <br>
> [0]: 14 ----> 2<br>
> <br>
> [0]: 14 ----> 5<br>
> <br>
> [0]: 15 ----> 6<br>
> <br>
> [0]: 15 ----> 7<br>
> <br>
> [0]: 15 ----> 9<br>
> <br>
> [0]: 16 ----> 6<br>
> <br>
> [0]: 16 ----> 8<br>
> <br>
> [0]: 16 ----> 10<br>
> <br>
> [0]: 17 ----> 7<br>
> <br>
> [0]: 17 ----> 8<br>
> <br>
> [0]: 17 ----> 11<br>
> <br>
> [0]: 18 ----> 9<br>
> <br>
> [0]: 18 ----> 12<br>
> <br>
> [0]: 18 ----> 13<br>
> <br>
> [0]: 19 ----> 10<br>
> <br>
> [0]: 19 ----> 12<br>
> <br>
> [0]: 19 ----> 14<br>
> <br>
> [0]: 20 ----> 11<br>
> <br>
> [0]: 20 ----> 13<br>
> <br>
> [0]: 20 ----> 14<br>
> <br>
> Cones:<br>
> <br>
> [0] Max cone size: 5<br>
> <br>
> [0]: 0 <---- 1 (0)<br>
> <br>
> [0]: 0 <---- 2 (0)<br>
> <br>
> [0]: 0 <---- 3 (0)<br>
> <br>
> [0]: 0 <---- 4 (0)<br>
> <br>
> [0]: 0 <---- 5 (0)<br>
> <br>
> [0]: 1 <---- 6 (0)<br>
> <br>
> [0]: 1 <---- 7 (0)<br>
> <br>
> [0]: 1 <---- 8 (0)<br>
> <br>
> [0]: 2 <---- 12 (0)<br>
> <br>
> [0]: 2 <---- 13 (0)<br>
> <br>
> [0]: 2 <---- 14 (0)<br>
> <br>
> [0]: 3 <---- 7 (0)<br>
> <br>
> [0]: 3 <---- 9 (0)<br>
> <br>
> [0]: 3 <---- 11 (0)<br>
> <br>
> [0]: 3 <---- 13 (0)<br>
> <br>
> [0]: 4 <---- 6 (0)<br>
> <br>
> [0]: 4 <---- 9 (0)<br>
> <br>
> [0]: 4 <---- 10 (0)<br>
> <br>
> [0]: 4 <---- 12 (0)<br>
> <br>
> [0]: 5 <---- 8 (0)<br>
> <br>
> [0]: 5 <---- 10 (0)<br>
> <br>
> [0]: 5 <---- 11 (0)<br>
> <br>
> [0]: 5 <---- 14 (0)<br>
> <br>
> [0]: 6 <---- 15 (0)<br>
> <br>
> [0]: 6 <---- 16 (0)<br>
> <br>
> [0]: 7 <---- 15 (0)<br>
> <br>
> [0]: 7 <---- 17 (0)<br>
> <br>
> [0]: 8 <---- 16 (0)<br>
> <br>
> [0]: 8 <---- 17 (0)<br>
> <br>
> [0]: 9 <---- 15 (0)<br>
> <br>
> [0]: 9 <---- 18 (0)<br>
> <br>
> [0]: 10 <---- 16 (0)<br>
> <br>
> [0]: 10 <---- 19 (0)<br>
> <br>
> [0]: 11 <---- 17 (0)<br>
> <br>
> [0]: 11 <---- 20 (0)<br>
> <br>
> [0]: 12 <---- 18 (0)<br>
> <br>
> [0]: 12 <---- 19 (0)<br>
> <br>
> [0]: 13 <---- 18 (0)<br>
> <br>
> [0]: 13 <---- 20 (0)<br>
> <br>
> [0]: 14 <---- 19 (0)<br>
> <br>
> [0]: 14 <---- 20 (0)<br>
> <br>
> coordinates with 1 fields<br>
> <br>
>   field 0 with 3 components<br>
> <br>
> Process 0:<br>
> <br>
>   (   0) dim  0 offset   0<br>
> <br>
>   (   1) dim  0 offset   0<br>
> <br>
>   (   2) dim  0 offset   0<br>
> <br>
>   (   3) dim  0 offset   0<br>
> <br>
>   (   4) dim  0 offset   0<br>
> <br>
>   (   5) dim  0 offset   0<br>
> <br>
>   (   6) dim  0 offset   0<br>
> <br>
>   (   7) dim  0 offset   0<br>
> <br>
>   (   8) dim  0 offset   0<br>
> <br>
>   (   9) dim  0 offset   0<br>
> <br>
>   (  10) dim  0 offset   0<br>
> <br>
>   (  11) dim  0 offset   0<br>
> <br>
>   (  12) dim  0 offset   0<br>
> <br>
>   (  13) dim  0 offset   0<br>
> <br>
>   (  14) dim  0 offset   0<br>
> <br>
>   (  15) dim  3 offset   0 0. 0. 0.<br>
> <br>
>   (  16) dim  3 offset   3 0. 2. 0.<br>
> <br>
>   (  17) dim  3 offset   6 2. 1. 0.<br>
> <br>
>   (  18) dim  3 offset   9 0. 0. 1.<br>
> <br>
>   (  19) dim  3 offset  12 0. 2. 1.<br>
> <br>
>   (  20) dim  3 offset  15 2. 1. 1.<br>
> <br>
>  Cone           0 :           1           2           3           4<br>
>   5<br>
> <br>
>  Cone           1 :           6           7           8<br>
> <br>
>  Cone           2 :          12          13          14<br>
> <br>
>  Cone           3 :           7           9          11          13<br>
> <br>
>  Cone           4 :           6           9          10          12<br>
> <br>
>  Cone           5 :           8          10          11          14<br>
> <br>
>  Cone           6 :          15          16<br>
> <br>
>  Cone           7 :          15          17<br>
> <br>
>  Cone           8 :          16          17<br>
> <br>
>  Cone           9 :          15          18<br>
> <br>
>  Cone          10 :          16          19<br>
> <br>
>  Cone          11 :          17          20<br>
> <br>
>  Cone          12 :          18          19<br>
> <br>
>  Cone          13 :          18          20<br>
> <br>
>  Cone          14 :          19          20<br>
> <br>
>  Cone          15 :<br>
> <br>
>  Cone          16 :<br>
> <br>
>  Cone          17 :<br>
> <br>
>  Cone          18 :<br>
> <br>
>  Cone          19 :<br>
> <br>
>  Cone          20 :<br>
> <br>
> DM Object: testplex 1 MPI processes<br>
> <br>
>   type: plex<br>
> <br>
> testplex in 3 dimensions:<br>
> <br>
> Supports:<br>
> <br>
> [0] Max support size: 3<br>
> <br>
> [0]: 1 ----> 0<br>
> <br>
> [0]: 2 ----> 0<br>
> <br>
> [0]: 3 ----> 0<br>
> <br>
> [0]: 4 ----> 0<br>
> <br>
> [0]: 5 ----> 0<br>
> <br>
> [0]: 6 ----> 1<br>
> <br>
> [0]: 6 ----> 4<br>
> <br>
> [0]: 7 ----> 1<br>
> <br>
> [0]: 7 ----> 3<br>
> <br>
> [0]: 8 ----> 1<br>
> <br>
> [0]: 8 ----> 5<br>
> <br>
> [0]: 9 ----> 3<br>
> <br>
> [0]: 9 ----> 4<br>
> <br>
> [0]: 10 ----> 4<br>
> <br>
> [0]: 10 ----> 5<br>
> <br>
> [0]: 11 ----> 3<br>
> <br>
> [0]: 11 ----> 5<br>
> <br>
> [0]: 12 ----> 2<br>
> <br>
> [0]: 12 ----> 4<br>
> <br>
> [0]: 13 ----> 2<br>
> <br>
> [0]: 13 ----> 3<br>
> <br>
> [0]: 14 ----> 2<br>
> <br>
> [0]: 14 ----> 5<br>
> <br>
> [0]: 15 ----> 6<br>
> <br>
> [0]: 15 ----> 7<br>
> <br>
> [0]: 15 ----> 9<br>
> <br>
> [0]: 16 ----> 6<br>
> <br>
> [0]: 16 ----> 8<br>
> <br>
> [0]: 16 ----> 10<br>
> <br>
> [0]: 17 ----> 7<br>
> <br>
> [0]: 17 ----> 8<br>
> <br>
> [0]: 17 ----> 11<br>
> <br>
> [0]: 18 ----> 9<br>
> <br>
> [0]: 18 ----> 12<br>
> <br>
> [0]: 18 ----> 13<br>
> <br>
> [0]: 19 ----> 10<br>
> <br>
> [0]: 19 ----> 12<br>
> <br>
> [0]: 19 ----> 14<br>
> <br>
> [0]: 20 ----> 11<br>
> <br>
> [0]: 20 ----> 13<br>
> <br>
> [0]: 20 ----> 14<br>
> <br>
> Cones:<br>
> <br>
> [0] Max cone size: 5<br>
> <br>
> [0]: 0 <---- 1 (0)<br>
> <br>
> [0]: 0 <---- 2 (0)<br>
> <br>
> [0]: 0 <---- 3 (0)<br>
> <br>
> [0]: 0 <---- 4 (0)<br>
> <br>
> [0]: 0 <---- 5 (0)<br>
> <br>
> [0]: 1 <---- 6 (0)<br>
> <br>
> [0]: 1 <---- 7 (0)<br>
> <br>
> [0]: 1 <---- 8 (0)<br>
> <br>
> [0]: 2 <---- 12 (0)<br>
> <br>
> [0]: 2 <---- 13 (0)<br>
> <br>
> [0]: 2 <---- 14 (0)<br>
> <br>
> [0]: 3 <---- 7 (0)<br>
> <br>
> [0]: 3 <---- 9 (0)<br>
> <br>
> [0]: 3 <---- 11 (0)<br>
> <br>
> [0]: 3 <---- 13 (0)<br>
> <br>
> [0]: 4 <---- 6 (0)<br>
> <br>
> [0]: 4 <---- 9 (0)<br>
> <br>
> [0]: 4 <---- 10 (0)<br>
> <br>
> [0]: 4 <---- 12 (0)<br>
> <br>
> [0]: 5 <---- 8 (0)<br>
> <br>
> [0]: 5 <---- 10 (0)<br>
> <br>
> [0]: 5 <---- 11 (0)<br>
> <br>
> [0]: 5 <---- 14 (0)<br>
> <br>
> [0]: 6 <---- 15 (0)<br>
> <br>
> [0]: 6 <---- 16 (0)<br>
> <br>
> [0]: 7 <---- 15 (0)<br>
> <br>
> [0]: 7 <---- 17 (0)<br>
> <br>
> [0]: 8 <---- 16 (0)<br>
> <br>
> [0]: 8 <---- 17 (0)<br>
> <br>
> [0]: 9 <---- 15 (0)<br>
> <br>
> [0]: 9 <---- 18 (0)<br>
> <br>
> [0]: 10 <---- 16 (0)<br>
> <br>
> [0]: 10 <---- 19 (0)<br>
> <br>
> [0]: 11 <---- 17 (0)<br>
> <br>
> [0]: 11 <---- 20 (0)<br>
> <br>
> [0]: 12 <---- 18 (0)<br>
> <br>
> [0]: 12 <---- 19 (0)<br>
> <br>
> [0]: 13 <---- 18 (0)<br>
> <br>
> [0]: 13 <---- 20 (0)<br>
> <br>
> [0]: 14 <---- 19 (0)<br>
> <br>
> [0]: 14 <---- 20 (0)<br>
> <br>
> coordinates with 1 fields<br>
> <br>
>   field 0 with 3 components<br>
> <br>
> Process 0:<br>
> <br>
>   (   0) dim  0 offset   0<br>
> <br>
>   (   1) dim  0 offset   0<br>
> <br>
>   (   2) dim  0 offset   0<br>
> <br>
>   (   3) dim  0 offset   0<br>
> <br>
>   (   4) dim  0 offset   0<br>
> <br>
>   (   5) dim  0 offset   0<br>
> <br>
>   (   6) dim  0 offset   0<br>
> <br>
>   (   7) dim  0 offset   0<br>
> <br>
>   (   8) dim  0 offset   0<br>
> <br>
>   (   9) dim  0 offset   0<br>
> <br>
>   (  10) dim  0 offset   0<br>
> <br>
>   (  11) dim  0 offset   0<br>
> <br>
>   (  12) dim  0 offset   0<br>
> <br>
>   (  13) dim  0 offset   0<br>
> <br>
>   (  14) dim  0 offset   0<br>
> <br>
>   (  15) dim  3 offset   0 0. 0. 0.<br>
> <br>
>   (  16) dim  3 offset   3 0. 2. 0.<br>
> <br>
>   (  17) dim  3 offset   6 2. 1. 0.<br>
> <br>
>   (  18) dim  3 offset   9 0. 0. 1.<br>
> <br>
>   (  19) dim  3 offset  12 0. 2. 1.<br>
> <br>
>   (  20) dim  3 offset  15 2. 1. 1.<br>
> <br>
> <br>
>   Thanks,<br>
> <br>
> <br>
>      Matt<br>
> <br>
> <br>
>> As always, many many thanks!<br>
>><br>
>> Fabian<br>
>><br>
>> P.S. here the output of<br>
>> make run<br>
>><br>
>> and note the cones, the second time they are outputted...<br>
>><br>
>> ./plex_prism -show_plex ::ascii_info_detail<br>
>> DM Object: testplex 1 MPI processes<br>
>>   type: plex<br>
>> testplex in 3 dimensions:<br>
>> Supports:<br>
>> [0] Max support size: 3<br>
>> [0]: 1 ----> 0<br>
>> [0]: 2 ----> 0<br>
>> [0]: 3 ----> 0<br>
>> [0]: 4 ----> 0<br>
>> [0]: 5 ----> 0<br>
>> [0]: 6 ----> 1<br>
>> [0]: 6 ----> 4<br>
>> [0]: 7 ----> 1<br>
>> [0]: 7 ----> 3<br>
>> [0]: 8 ----> 1<br>
>> [0]: 8 ----> 5<br>
>> [0]: 9 ----> 3<br>
>> [0]: 9 ----> 4<br>
>> [0]: 10 ----> 4<br>
>> [0]: 10 ----> 5<br>
>> [0]: 11 ----> 3<br>
>> [0]: 11 ----> 5<br>
>> [0]: 12 ----> 2<br>
>> [0]: 12 ----> 4<br>
>> [0]: 13 ----> 2<br>
>> [0]: 13 ----> 3<br>
>> [0]: 14 ----> 2<br>
>> [0]: 14 ----> 5<br>
>> [0]: 15 ----> 6<br>
>> [0]: 15 ----> 7<br>
>> [0]: 15 ----> 9<br>
>> [0]: 16 ----> 6<br>
>> [0]: 16 ----> 8<br>
>> [0]: 16 ----> 10<br>
>> [0]: 17 ----> 7<br>
>> [0]: 17 ----> 8<br>
>> [0]: 17 ----> 11<br>
>> [0]: 18 ----> 9<br>
>> [0]: 18 ----> 12<br>
>> [0]: 18 ----> 13<br>
>> [0]: 19 ----> 10<br>
>> [0]: 19 ----> 12<br>
>> [0]: 19 ----> 14<br>
>> [0]: 20 ----> 11<br>
>> [0]: 20 ----> 13<br>
>> [0]: 20 ----> 14<br>
>> Cones:<br>
>> [0] Max cone size: 5<br>
>> [0]: 0 <---- 1 (0)<br>
>> [0]: 0 <---- 2 (0)<br>
>> [0]: 0 <---- 3 (0)<br>
>> [0]: 0 <---- 4 (0)<br>
>> [0]: 0 <---- 5 (0)<br>
>> [0]: 1 <---- 6 (0)<br>
>> [0]: 1 <---- 7 (0)<br>
>> [0]: 1 <---- 8 (0)<br>
>> [0]: 2 <---- 12 (0)<br>
>> [0]: 2 <---- 13 (0)<br>
>> [0]: 2 <---- 14 (0)<br>
>> [0]: 3 <---- 7 (0)<br>
>> [0]: 3 <---- 9 (0)<br>
>> [0]: 3 <---- 11 (0)<br>
>> [0]: 3 <---- 13 (0)<br>
>> [0]: 4 <---- 6 (0)<br>
>> [0]: 4 <---- 9 (0)<br>
>> [0]: 4 <---- 10 (0)<br>
>> [0]: 4 <---- 12 (0)<br>
>> [0]: 5 <---- 8 (0)<br>
>> [0]: 5 <---- 10 (0)<br>
>> [0]: 5 <---- 11 (0)<br>
>> [0]: 5 <---- 14 (0)<br>
>> [0]: 6 <---- 15 (0)<br>
>> [0]: 6 <---- 16 (0)<br>
>> [0]: 7 <---- 15 (0)<br>
>> [0]: 7 <---- 17 (0)<br>
>> [0]: 8 <---- 16 (0)<br>
>> [0]: 8 <---- 17 (0)<br>
>> [0]: 9 <---- 15 (0)<br>
>> [0]: 9 <---- 18 (0)<br>
>> [0]: 10 <---- 16 (0)<br>
>> [0]: 10 <---- 19 (0)<br>
>> [0]: 11 <---- 17 (0)<br>
>> [0]: 11 <---- 20 (0)<br>
>> [0]: 12 <---- 18 (0)<br>
>> [0]: 12 <---- 19 (0)<br>
>> [0]: 13 <---- 18 (0)<br>
>> [0]: 13 <---- 20 (0)<br>
>> [0]: 14 <---- 19 (0)<br>
>> [0]: 14 <---- 20 (0)<br>
>> coordinates with 1 fields<br>
>>   field 0 with 3 components<br>
>> Process 0:<br>
>>   (   0) dim  0 offset   0<br>
>>   (   1) dim  0 offset   0<br>
>>   (   2) dim  0 offset   0<br>
>>   (   3) dim  0 offset   0<br>
>>   (   4) dim  0 offset   0<br>
>>   (   5) dim  0 offset   0<br>
>>   (   6) dim  0 offset   0<br>
>>   (   7) dim  0 offset   0<br>
>>   (   8) dim  0 offset   0<br>
>>   (   9) dim  0 offset   0<br>
>>   (  10) dim  0 offset   0<br>
>>   (  11) dim  0 offset   0<br>
>>   (  12) dim  0 offset   0<br>
>>   (  13) dim  0 offset   0<br>
>>   (  14) dim  0 offset   0<br>
>>   (  15) dim  3 offset   0 0. 0. 0.<br>
>>   (  16) dim  3 offset   3 0. 2. 0.<br>
>>   (  17) dim  3 offset   6 2. 1. 0.<br>
>>   (  18) dim  3 offset   9 0. 0. 1.<br>
>>   (  19) dim  3 offset  12 0. 2. 1.<br>
>>   (  20) dim  3 offset  15 2. 1. 1.<br>
>><br>
>><br>
>><br>
>><br>
>><br>
>><br>
>><br>
>><br>
>><br>
>> DM Object: testplex 1 MPI processes<br>
>>   type: plex<br>
>> testplex in 3 dimensions:<br>
>> Supports:<br>
>> [0] Max support size: 3<br>
>> [0]: 1 ----> 0<br>
>> [0]: 2 ----> 0<br>
>> [0]: 3 ----> 0<br>
>> [0]: 4 ----> 0<br>
>> [0]: 5 ----> 0<br>
>> [0]: 6 ----> 1<br>
>> [0]: 6 ----> 4<br>
>> [0]: 7 ----> 1<br>
>> [0]: 7 ----> 3<br>
>> [0]: 8 ----> 1<br>
>> [0]: 8 ----> 5<br>
>> [0]: 9 ----> 3<br>
>> [0]: 9 ----> 4<br>
>> [0]: 10 ----> 4<br>
>> [0]: 10 ----> 5<br>
>> [0]: 11 ----> 3<br>
>> [0]: 11 ----> 5<br>
>> [0]: 12 ----> 2<br>
>> [0]: 12 ----> 4<br>
>> [0]: 13 ----> 2<br>
>> [0]: 13 ----> 3<br>
>> [0]: 14 ----> 2<br>
>> [0]: 14 ----> 5<br>
>> [0]: 15 ----> 6<br>
>> [0]: 15 ----> 7<br>
>> [0]: 15 ----> 9<br>
>> [0]: 16 ----> 6<br>
>> [0]: 16 ----> 8<br>
>> [0]: 16 ----> 10<br>
>> [0]: 17 ----> 7<br>
>> [0]: 17 ----> 8<br>
>> [0]: 17 ----> 11<br>
>> [0]: 18 ----> 9<br>
>> [0]: 18 ----> 12<br>
>> [0]: 18 ----> 13<br>
>> [0]: 19 ----> 10<br>
>> [0]: 19 ----> 12<br>
>> [0]: 19 ----> 14<br>
>> [0]: 20 ----> 11<br>
>> [0]: 20 ----> 13<br>
>> [0]: 20 ----> 14<br>
>> Cones:<br>
>> [0] Max cone size: 5<br>
>> [0]: 0 <---- 0 (0)<br>
>> [0]: 0 <---- 0 (0)<br>
>> [0]: 0 <---- 3 (0)<br>
>> [0]: 0 <---- 4 (0)<br>
>> [0]: 0 <---- 5 (0)<br>
>> [0]: 1 <---- 0 (0)<br>
>> [0]: 1 <---- 0 (0)<br>
>> [0]: 1 <---- 8 (0)<br>
>> [0]: 2 <---- 0 (0)<br>
>> [0]: 2 <---- 0 (0)<br>
>> [0]: 2 <---- 14 (0)<br>
>> [0]: 3 <---- 0 (0)<br>
>> [0]: 3 <---- 0 (0)<br>
>> [0]: 3 <---- 11 (0)<br>
>> [0]: 3 <---- 13 (0)<br>
>> [0]: 4 <---- 0 (0)<br>
>> [0]: 4 <---- 0 (0)<br>
>> [0]: 4 <---- 10 (0)<br>
>> [0]: 4 <---- 12 (0)<br>
>> [0]: 5 <---- 0 (0)<br>
>> [0]: 5 <---- 0 (0)<br>
>> [0]: 5 <---- 11 (0)<br>
>> [0]: 5 <---- 14 (0)<br>
>> [0]: 6 <---- 0 (0)<br>
>> [0]: 6 <---- 0 (0)<br>
>> [0]: 7 <---- 0 (0)<br>
>> [0]: 7 <---- 0 (0)<br>
>> [0]: 8 <---- 0 (0)<br>
>> [0]: 8 <---- 0 (0)<br>
>> [0]: 9 <---- 0 (0)<br>
>> [0]: 9 <---- 0 (0)<br>
>> [0]: 10 <---- 0 (0)<br>
>> [0]: 10 <---- 0 (0)<br>
>> [0]: 11 <---- 0 (0)<br>
>> [0]: 11 <---- 0 (0)<br>
>> [0]: 12 <---- 0 (0)<br>
>> [0]: 12 <---- 0 (0)<br>
>> [0]: 13 <---- 0 (0)<br>
>> [0]: 13 <---- 0 (0)<br>
>> [0]: 14 <---- 0 (0)<br>
>> [0]: 14 <---- 0 (0)<br>
>> coordinates with 1 fields<br>
>>   field 0 with 3 components<br>
>> Process 0:<br>
>>   (   0) dim  0 offset   0<br>
>>   (   1) dim  0 offset   0<br>
>>   (   2) dim  0 offset   0<br>
>>   (   3) dim  0 offset   0<br>
>>   (   4) dim  0 offset   0<br>
>>   (   5) dim  0 offset   0<br>
>>   (   6) dim  0 offset   0<br>
>>   (   7) dim  0 offset   0<br>
>>   (   8) dim  0 offset   0<br>
>>   (   9) dim  0 offset   0<br>
>>   (  10) dim  0 offset   0<br>
>>   (  11) dim  0 offset   0<br>
>>   (  12) dim  0 offset   0<br>
>>   (  13) dim  0 offset   0<br>
>>   (  14) dim  0 offset   0<br>
>>   (  15) dim  3 offset   0 0. 0. 0.<br>
>>   (  16) dim  3 offset   3 0. 2. 0.<br>
>>   (  17) dim  3 offset   6 2. 1. 0.<br>
>>   (  18) dim  3 offset   9 0. 0. 1.<br>
>>   (  19) dim  3 offset  12 0. 2. 1.<br>
>>   (  20) dim  3 offset  15 2. 1. 1.<br>
>><br>
>><br>
> <br>
<br>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div><div><br></div><div><a href="http://www.cse.buffalo.edu/~knepley/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div></div></div></div>