<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
</head>
<body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">
<br class="">
<div><br class="">
<blockquote type="cite" class="">
<div class="">On 4 Jul 2019, at 15:36, Matthew Knepley <<a href="mailto:knepley@gmail.com" class="">knepley@gmail.com</a>> wrote:</div>
<br class="Apple-interchange-newline">
<div class="">
<div dir="ltr" class="">
<div dir="ltr" class="">On Wed, Jul 3, 2019 at 1:03 PM Fabian.Jakub <<a href="mailto:Fabian.Jakub@physik.uni-muenchen.de" class="">Fabian.Jakub@physik.uni-muenchen.de</a>> wrote:<br class="">
</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 class="">
it was in Commit 6d563f2afc7bcefff9dc7869193a6debf1d639a5<br class="">
Date: Wed Jun 26 12:32:01 2019 +0000<br class="">
<br class="">
remove remnants of non-existing<br class="">
DMPlexRestore{Cone,ConeOrientation,Support}<br class="">
</blockquote>
<div class=""><br class="">
</div>
<div class="">You are right!</div>
<div class=""><br class="">
</div>
<div class="">Vaclav, you cannot remove the Fortran RestoreCone() because it is actually doing something.</div>
<div class="">Fortran needs to create arrays, and this code deallocates them. We do not need that in C.</div>
<div class="">Can you put those back?</div>
</div>
</div>
</div>
</blockquote>
<div><br class="">
</div>
Sure, but then it should be documented somehow. Having that functions just for Fortran without explanation is very confusing.</div>
<div><br class="">
</div>
<div>But I only removed the following interface declarations and there are no other occurrences in the whole source tree:</div>
<div><br class="">
</div>
<div>
<div>git diff 6d563f2~ -- src/dm/f90-mod/petscdmplex.h90</div>
<div>diff --git a/src/dm/f90-mod/petscdmplex.h90 b/src/dm/f90-mod/petscdmplex.h90</div>
<div>index 7182cba..d20cdc4 100644</div>
<div>--- a/src/dm/f90-mod/petscdmplex.h90</div>
<div>+++ b/src/dm/f90-mod/petscdmplex.h90</div>
<div>@@ -10,16 +10,6 @@</div>
<div> End Interface</div>
<div> </div>
<div> Interface</div>
<div>- Subroutine DMPlexRestoreCone(m,p,cone,ierr)</div>
<div>- use petscdmdef</div>
<div>- PetscInt p</div>
<div>- PetscInt, pointer :: cone(:)</div>
<div>- PetscErrorCode ierr</div>
<div>- DM m</div>
<div>- End Subroutine</div>
<div>- End Interface</div>
<div>-</div>
<div>- Interface</div>
<div> Subroutine DMPlexGetConeOrientation(m,p,coneOrient,ierr)</div>
<div> use petscdmdef</div>
<div> PetscInt p</div>
<div>@@ -40,16 +30,6 @@</div>
<div> End Interface</div>
<div> </div>
<div> Interface</div>
<div>- Subroutine DMPlexRestoreSupport(m,p,support,ierr)</div>
<div>- use petscdmdef</div>
<div>- PetscInt p</div>
<div>- PetscInt, pointer :: support(:)</div>
<div>- PetscErrorCode ierr</div>
<div>- DM m</div>
<div>- End Subroutine</div>
<div>- End Interface</div>
<div>-</div>
<div>- Interface</div>
<div> Subroutine DMPlexGetTransitiveClosure(m,p,useCone,points,ierr)</div>
<div> use petscdmdef</div>
<div> PetscInt p</div>
<div class=""><br class="">
</div>
</div>
<div>So I wonder what these can actually do?</div>
<div><br class="">
</div>
<div>Where are the mentioned C functions? I don't see any.</div>
<div><br class="">
</div>
<div>
<div class="">When we are at it, why e.g. DMPlexGetCone() has the /*@C comment on top when no custom fortran stub is provided?</div>
<div class=""><br class="">
</div>
<div class="">Vaclav</div>
<div class=""><br class="">
</div>
<blockquote type="cite" class="">
<div class="">
<div dir="ltr" class="">
<div class="gmail_quote">
<div class=""><br class="">
</div>
<div class=""> Thanks,</div>
<div class=""><br class="">
</div>
<div class=""> Matt</div>
<div class=""> </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 class="">
In this case I guess that the fortran compiler reverts to the default<br class="">
calling convention and thus does not nullify the cone pointers but<br class="">
rather nulls the values in the cone.<br class="">
<br class="">
I suggest you remove the C functions to get a compiler error or revert<br class="">
the path.<br class="">
<br class="">
Many thanks,<br class="">
<br class="">
Fabian<br class="">
<br class="">
On 7/3/19 3:18 PM, Matthew Knepley wrote:<br class="">
> On Wed, Jul 3, 2019 at 7:38 AM Fabian.Jakub via petsc-users <<br class="">
> <a href="mailto:petsc-users@mcs.anl.gov" target="_blank" class="">petsc-users@mcs.anl.gov</a>> wrote:<br class="">
> <br class="">
>> Dear Matt, Dear Petsc<br class="">
>><br class="">
>> After my recent pull of petsc master I encounter test fails that use<br class="">
>> dmplex meshes.<br class="">
>><br class="">
>> I attached a simple script to show an issue where DMPlexRestoreCone<br class="">
>> mangles the cone relationships in the DM.<br class="">
>><br class="">
>> The behaviour happens with gcc 7.4.0 aswell as icc 19.0.3.199<br class="">
>> @ current petsc/master d3dc048301fdb3e4c282b528e7b28910551cfe13<br class="">
>><br class="">
>> We recently had a system update and I am not 100% sure if this is some<br class="">
>> weird stuff going on at my end or if this is indeed Plex related or if I<br class="">
>> should not call Restore in the inner loop.<br class="">
>><br class="">
>> Anyway, I would be most grateful if you could try to reproduce the error<br class="">
>> and hint me towards a solution.<br class="">
>><br class="">
> <br class="">
> Hi Fabian,<br class="">
> <br class="">
> I think it might be your system update. Seems to run fine for me:<br class="">
> <br class="">
> master *:~/Downloads/tmp$ ./plex_prism<br class="">
> <br class="">
> Cone 0 : 1 2 3 4<br class="">
> 5<br class="">
> <br class="">
> Cone 1 : 6 7 8<br class="">
> <br class="">
> Cone 2 : 12 13 14<br class="">
> <br class="">
> Cone 3 : 7 9 11 13<br class="">
> <br class="">
> Cone 4 : 6 9 10 12<br class="">
> <br class="">
> Cone 5 : 8 10 11 14<br class="">
> <br class="">
> Cone 6 : 15 16<br class="">
> <br class="">
> Cone 7 : 15 17<br class="">
> <br class="">
> Cone 8 : 16 17<br class="">
> <br class="">
> Cone 9 : 15 18<br class="">
> <br class="">
> Cone 10 : 16 19<br class="">
> <br class="">
> Cone 11 : 17 20<br class="">
> <br class="">
> Cone 12 : 18 19<br class="">
> <br class="">
> Cone 13 : 18 20<br class="">
> <br class="">
> Cone 14 : 19 20<br class="">
> <br class="">
> Cone 15 :<br class="">
> <br class="">
> Cone 16 :<br class="">
> <br class="">
> Cone 17 :<br class="">
> <br class="">
> Cone 18 :<br class="">
> <br class="">
> Cone 19 :<br class="">
> <br class="">
> Cone 20 :<br class="">
> <br class="">
> master *:~/Downloads/tmp$ ./plex_prism -show_plex ::ascii_info_detail<br class="">
> <br class="">
> DM Object: testplex 1 MPI processes<br class="">
> <br class="">
> type: plex<br class="">
> <br class="">
> testplex in 3 dimensions:<br class="">
> <br class="">
> Supports:<br class="">
> <br class="">
> [0] Max support size: 3<br class="">
> <br class="">
> [0]: 1 ----> 0<br class="">
> <br class="">
> [0]: 2 ----> 0<br class="">
> <br class="">
> [0]: 3 ----> 0<br class="">
> <br class="">
> [0]: 4 ----> 0<br class="">
> <br class="">
> [0]: 5 ----> 0<br class="">
> <br class="">
> [0]: 6 ----> 1<br class="">
> <br class="">
> [0]: 6 ----> 4<br class="">
> <br class="">
> [0]: 7 ----> 1<br class="">
> <br class="">
> [0]: 7 ----> 3<br class="">
> <br class="">
> [0]: 8 ----> 1<br class="">
> <br class="">
> [0]: 8 ----> 5<br class="">
> <br class="">
> [0]: 9 ----> 3<br class="">
> <br class="">
> [0]: 9 ----> 4<br class="">
> <br class="">
> [0]: 10 ----> 4<br class="">
> <br class="">
> [0]: 10 ----> 5<br class="">
> <br class="">
> [0]: 11 ----> 3<br class="">
> <br class="">
> [0]: 11 ----> 5<br class="">
> <br class="">
> [0]: 12 ----> 2<br class="">
> <br class="">
> [0]: 12 ----> 4<br class="">
> <br class="">
> [0]: 13 ----> 2<br class="">
> <br class="">
> [0]: 13 ----> 3<br class="">
> <br class="">
> [0]: 14 ----> 2<br class="">
> <br class="">
> [0]: 14 ----> 5<br class="">
> <br class="">
> [0]: 15 ----> 6<br class="">
> <br class="">
> [0]: 15 ----> 7<br class="">
> <br class="">
> [0]: 15 ----> 9<br class="">
> <br class="">
> [0]: 16 ----> 6<br class="">
> <br class="">
> [0]: 16 ----> 8<br class="">
> <br class="">
> [0]: 16 ----> 10<br class="">
> <br class="">
> [0]: 17 ----> 7<br class="">
> <br class="">
> [0]: 17 ----> 8<br class="">
> <br class="">
> [0]: 17 ----> 11<br class="">
> <br class="">
> [0]: 18 ----> 9<br class="">
> <br class="">
> [0]: 18 ----> 12<br class="">
> <br class="">
> [0]: 18 ----> 13<br class="">
> <br class="">
> [0]: 19 ----> 10<br class="">
> <br class="">
> [0]: 19 ----> 12<br class="">
> <br class="">
> [0]: 19 ----> 14<br class="">
> <br class="">
> [0]: 20 ----> 11<br class="">
> <br class="">
> [0]: 20 ----> 13<br class="">
> <br class="">
> [0]: 20 ----> 14<br class="">
> <br class="">
> Cones:<br class="">
> <br class="">
> [0] Max cone size: 5<br class="">
> <br class="">
> [0]: 0 <---- 1 (0)<br class="">
> <br class="">
> [0]: 0 <---- 2 (0)<br class="">
> <br class="">
> [0]: 0 <---- 3 (0)<br class="">
> <br class="">
> [0]: 0 <---- 4 (0)<br class="">
> <br class="">
> [0]: 0 <---- 5 (0)<br class="">
> <br class="">
> [0]: 1 <---- 6 (0)<br class="">
> <br class="">
> [0]: 1 <---- 7 (0)<br class="">
> <br class="">
> [0]: 1 <---- 8 (0)<br class="">
> <br class="">
> [0]: 2 <---- 12 (0)<br class="">
> <br class="">
> [0]: 2 <---- 13 (0)<br class="">
> <br class="">
> [0]: 2 <---- 14 (0)<br class="">
> <br class="">
> [0]: 3 <---- 7 (0)<br class="">
> <br class="">
> [0]: 3 <---- 9 (0)<br class="">
> <br class="">
> [0]: 3 <---- 11 (0)<br class="">
> <br class="">
> [0]: 3 <---- 13 (0)<br class="">
> <br class="">
> [0]: 4 <---- 6 (0)<br class="">
> <br class="">
> [0]: 4 <---- 9 (0)<br class="">
> <br class="">
> [0]: 4 <---- 10 (0)<br class="">
> <br class="">
> [0]: 4 <---- 12 (0)<br class="">
> <br class="">
> [0]: 5 <---- 8 (0)<br class="">
> <br class="">
> [0]: 5 <---- 10 (0)<br class="">
> <br class="">
> [0]: 5 <---- 11 (0)<br class="">
> <br class="">
> [0]: 5 <---- 14 (0)<br class="">
> <br class="">
> [0]: 6 <---- 15 (0)<br class="">
> <br class="">
> [0]: 6 <---- 16 (0)<br class="">
> <br class="">
> [0]: 7 <---- 15 (0)<br class="">
> <br class="">
> [0]: 7 <---- 17 (0)<br class="">
> <br class="">
> [0]: 8 <---- 16 (0)<br class="">
> <br class="">
> [0]: 8 <---- 17 (0)<br class="">
> <br class="">
> [0]: 9 <---- 15 (0)<br class="">
> <br class="">
> [0]: 9 <---- 18 (0)<br class="">
> <br class="">
> [0]: 10 <---- 16 (0)<br class="">
> <br class="">
> [0]: 10 <---- 19 (0)<br class="">
> <br class="">
> [0]: 11 <---- 17 (0)<br class="">
> <br class="">
> [0]: 11 <---- 20 (0)<br class="">
> <br class="">
> [0]: 12 <---- 18 (0)<br class="">
> <br class="">
> [0]: 12 <---- 19 (0)<br class="">
> <br class="">
> [0]: 13 <---- 18 (0)<br class="">
> <br class="">
> [0]: 13 <---- 20 (0)<br class="">
> <br class="">
> [0]: 14 <---- 19 (0)<br class="">
> <br class="">
> [0]: 14 <---- 20 (0)<br class="">
> <br class="">
> coordinates with 1 fields<br class="">
> <br class="">
> field 0 with 3 components<br class="">
> <br class="">
> Process 0:<br class="">
> <br class="">
> ( 0) dim 0 offset 0<br class="">
> <br class="">
> ( 1) dim 0 offset 0<br class="">
> <br class="">
> ( 2) dim 0 offset 0<br class="">
> <br class="">
> ( 3) dim 0 offset 0<br class="">
> <br class="">
> ( 4) dim 0 offset 0<br class="">
> <br class="">
> ( 5) dim 0 offset 0<br class="">
> <br class="">
> ( 6) dim 0 offset 0<br class="">
> <br class="">
> ( 7) dim 0 offset 0<br class="">
> <br class="">
> ( 8) dim 0 offset 0<br class="">
> <br class="">
> ( 9) dim 0 offset 0<br class="">
> <br class="">
> ( 10) dim 0 offset 0<br class="">
> <br class="">
> ( 11) dim 0 offset 0<br class="">
> <br class="">
> ( 12) dim 0 offset 0<br class="">
> <br class="">
> ( 13) dim 0 offset 0<br class="">
> <br class="">
> ( 14) dim 0 offset 0<br class="">
> <br class="">
> ( 15) dim 3 offset 0 0. 0. 0.<br class="">
> <br class="">
> ( 16) dim 3 offset 3 0. 2. 0.<br class="">
> <br class="">
> ( 17) dim 3 offset 6 2. 1. 0.<br class="">
> <br class="">
> ( 18) dim 3 offset 9 0. 0. 1.<br class="">
> <br class="">
> ( 19) dim 3 offset 12 0. 2. 1.<br class="">
> <br class="">
> ( 20) dim 3 offset 15 2. 1. 1.<br class="">
> <br class="">
> Cone 0 : 1 2 3 4<br class="">
> 5<br class="">
> <br class="">
> Cone 1 : 6 7 8<br class="">
> <br class="">
> Cone 2 : 12 13 14<br class="">
> <br class="">
> Cone 3 : 7 9 11 13<br class="">
> <br class="">
> Cone 4 : 6 9 10 12<br class="">
> <br class="">
> Cone 5 : 8 10 11 14<br class="">
> <br class="">
> Cone 6 : 15 16<br class="">
> <br class="">
> Cone 7 : 15 17<br class="">
> <br class="">
> Cone 8 : 16 17<br class="">
> <br class="">
> Cone 9 : 15 18<br class="">
> <br class="">
> Cone 10 : 16 19<br class="">
> <br class="">
> Cone 11 : 17 20<br class="">
> <br class="">
> Cone 12 : 18 19<br class="">
> <br class="">
> Cone 13 : 18 20<br class="">
> <br class="">
> Cone 14 : 19 20<br class="">
> <br class="">
> Cone 15 :<br class="">
> <br class="">
> Cone 16 :<br class="">
> <br class="">
> Cone 17 :<br class="">
> <br class="">
> Cone 18 :<br class="">
> <br class="">
> Cone 19 :<br class="">
> <br class="">
> Cone 20 :<br class="">
> <br class="">
> DM Object: testplex 1 MPI processes<br class="">
> <br class="">
> type: plex<br class="">
> <br class="">
> testplex in 3 dimensions:<br class="">
> <br class="">
> Supports:<br class="">
> <br class="">
> [0] Max support size: 3<br class="">
> <br class="">
> [0]: 1 ----> 0<br class="">
> <br class="">
> [0]: 2 ----> 0<br class="">
> <br class="">
> [0]: 3 ----> 0<br class="">
> <br class="">
> [0]: 4 ----> 0<br class="">
> <br class="">
> [0]: 5 ----> 0<br class="">
> <br class="">
> [0]: 6 ----> 1<br class="">
> <br class="">
> [0]: 6 ----> 4<br class="">
> <br class="">
> [0]: 7 ----> 1<br class="">
> <br class="">
> [0]: 7 ----> 3<br class="">
> <br class="">
> [0]: 8 ----> 1<br class="">
> <br class="">
> [0]: 8 ----> 5<br class="">
> <br class="">
> [0]: 9 ----> 3<br class="">
> <br class="">
> [0]: 9 ----> 4<br class="">
> <br class="">
> [0]: 10 ----> 4<br class="">
> <br class="">
> [0]: 10 ----> 5<br class="">
> <br class="">
> [0]: 11 ----> 3<br class="">
> <br class="">
> [0]: 11 ----> 5<br class="">
> <br class="">
> [0]: 12 ----> 2<br class="">
> <br class="">
> [0]: 12 ----> 4<br class="">
> <br class="">
> [0]: 13 ----> 2<br class="">
> <br class="">
> [0]: 13 ----> 3<br class="">
> <br class="">
> [0]: 14 ----> 2<br class="">
> <br class="">
> [0]: 14 ----> 5<br class="">
> <br class="">
> [0]: 15 ----> 6<br class="">
> <br class="">
> [0]: 15 ----> 7<br class="">
> <br class="">
> [0]: 15 ----> 9<br class="">
> <br class="">
> [0]: 16 ----> 6<br class="">
> <br class="">
> [0]: 16 ----> 8<br class="">
> <br class="">
> [0]: 16 ----> 10<br class="">
> <br class="">
> [0]: 17 ----> 7<br class="">
> <br class="">
> [0]: 17 ----> 8<br class="">
> <br class="">
> [0]: 17 ----> 11<br class="">
> <br class="">
> [0]: 18 ----> 9<br class="">
> <br class="">
> [0]: 18 ----> 12<br class="">
> <br class="">
> [0]: 18 ----> 13<br class="">
> <br class="">
> [0]: 19 ----> 10<br class="">
> <br class="">
> [0]: 19 ----> 12<br class="">
> <br class="">
> [0]: 19 ----> 14<br class="">
> <br class="">
> [0]: 20 ----> 11<br class="">
> <br class="">
> [0]: 20 ----> 13<br class="">
> <br class="">
> [0]: 20 ----> 14<br class="">
> <br class="">
> Cones:<br class="">
> <br class="">
> [0] Max cone size: 5<br class="">
> <br class="">
> [0]: 0 <---- 1 (0)<br class="">
> <br class="">
> [0]: 0 <---- 2 (0)<br class="">
> <br class="">
> [0]: 0 <---- 3 (0)<br class="">
> <br class="">
> [0]: 0 <---- 4 (0)<br class="">
> <br class="">
> [0]: 0 <---- 5 (0)<br class="">
> <br class="">
> [0]: 1 <---- 6 (0)<br class="">
> <br class="">
> [0]: 1 <---- 7 (0)<br class="">
> <br class="">
> [0]: 1 <---- 8 (0)<br class="">
> <br class="">
> [0]: 2 <---- 12 (0)<br class="">
> <br class="">
> [0]: 2 <---- 13 (0)<br class="">
> <br class="">
> [0]: 2 <---- 14 (0)<br class="">
> <br class="">
> [0]: 3 <---- 7 (0)<br class="">
> <br class="">
> [0]: 3 <---- 9 (0)<br class="">
> <br class="">
> [0]: 3 <---- 11 (0)<br class="">
> <br class="">
> [0]: 3 <---- 13 (0)<br class="">
> <br class="">
> [0]: 4 <---- 6 (0)<br class="">
> <br class="">
> [0]: 4 <---- 9 (0)<br class="">
> <br class="">
> [0]: 4 <---- 10 (0)<br class="">
> <br class="">
> [0]: 4 <---- 12 (0)<br class="">
> <br class="">
> [0]: 5 <---- 8 (0)<br class="">
> <br class="">
> [0]: 5 <---- 10 (0)<br class="">
> <br class="">
> [0]: 5 <---- 11 (0)<br class="">
> <br class="">
> [0]: 5 <---- 14 (0)<br class="">
> <br class="">
> [0]: 6 <---- 15 (0)<br class="">
> <br class="">
> [0]: 6 <---- 16 (0)<br class="">
> <br class="">
> [0]: 7 <---- 15 (0)<br class="">
> <br class="">
> [0]: 7 <---- 17 (0)<br class="">
> <br class="">
> [0]: 8 <---- 16 (0)<br class="">
> <br class="">
> [0]: 8 <---- 17 (0)<br class="">
> <br class="">
> [0]: 9 <---- 15 (0)<br class="">
> <br class="">
> [0]: 9 <---- 18 (0)<br class="">
> <br class="">
> [0]: 10 <---- 16 (0)<br class="">
> <br class="">
> [0]: 10 <---- 19 (0)<br class="">
> <br class="">
> [0]: 11 <---- 17 (0)<br class="">
> <br class="">
> [0]: 11 <---- 20 (0)<br class="">
> <br class="">
> [0]: 12 <---- 18 (0)<br class="">
> <br class="">
> [0]: 12 <---- 19 (0)<br class="">
> <br class="">
> [0]: 13 <---- 18 (0)<br class="">
> <br class="">
> [0]: 13 <---- 20 (0)<br class="">
> <br class="">
> [0]: 14 <---- 19 (0)<br class="">
> <br class="">
> [0]: 14 <---- 20 (0)<br class="">
> <br class="">
> coordinates with 1 fields<br class="">
> <br class="">
> field 0 with 3 components<br class="">
> <br class="">
> Process 0:<br class="">
> <br class="">
> ( 0) dim 0 offset 0<br class="">
> <br class="">
> ( 1) dim 0 offset 0<br class="">
> <br class="">
> ( 2) dim 0 offset 0<br class="">
> <br class="">
> ( 3) dim 0 offset 0<br class="">
> <br class="">
> ( 4) dim 0 offset 0<br class="">
> <br class="">
> ( 5) dim 0 offset 0<br class="">
> <br class="">
> ( 6) dim 0 offset 0<br class="">
> <br class="">
> ( 7) dim 0 offset 0<br class="">
> <br class="">
> ( 8) dim 0 offset 0<br class="">
> <br class="">
> ( 9) dim 0 offset 0<br class="">
> <br class="">
> ( 10) dim 0 offset 0<br class="">
> <br class="">
> ( 11) dim 0 offset 0<br class="">
> <br class="">
> ( 12) dim 0 offset 0<br class="">
> <br class="">
> ( 13) dim 0 offset 0<br class="">
> <br class="">
> ( 14) dim 0 offset 0<br class="">
> <br class="">
> ( 15) dim 3 offset 0 0. 0. 0.<br class="">
> <br class="">
> ( 16) dim 3 offset 3 0. 2. 0.<br class="">
> <br class="">
> ( 17) dim 3 offset 6 2. 1. 0.<br class="">
> <br class="">
> ( 18) dim 3 offset 9 0. 0. 1.<br class="">
> <br class="">
> ( 19) dim 3 offset 12 0. 2. 1.<br class="">
> <br class="">
> ( 20) dim 3 offset 15 2. 1. 1.<br class="">
> <br class="">
> <br class="">
> Thanks,<br class="">
> <br class="">
> <br class="">
> Matt<br class="">
> <br class="">
> <br class="">
>> As always, many many thanks!<br class="">
>><br class="">
>> Fabian<br class="">
>><br class="">
>> P.S. here the output of<br class="">
>> make run<br class="">
>><br class="">
>> and note the cones, the second time they are outputted...<br class="">
>><br class="">
>> ./plex_prism -show_plex ::ascii_info_detail<br class="">
>> DM Object: testplex 1 MPI processes<br class="">
>> type: plex<br class="">
>> testplex in 3 dimensions:<br class="">
>> Supports:<br class="">
>> [0] Max support size: 3<br class="">
>> [0]: 1 ----> 0<br class="">
>> [0]: 2 ----> 0<br class="">
>> [0]: 3 ----> 0<br class="">
>> [0]: 4 ----> 0<br class="">
>> [0]: 5 ----> 0<br class="">
>> [0]: 6 ----> 1<br class="">
>> [0]: 6 ----> 4<br class="">
>> [0]: 7 ----> 1<br class="">
>> [0]: 7 ----> 3<br class="">
>> [0]: 8 ----> 1<br class="">
>> [0]: 8 ----> 5<br class="">
>> [0]: 9 ----> 3<br class="">
>> [0]: 9 ----> 4<br class="">
>> [0]: 10 ----> 4<br class="">
>> [0]: 10 ----> 5<br class="">
>> [0]: 11 ----> 3<br class="">
>> [0]: 11 ----> 5<br class="">
>> [0]: 12 ----> 2<br class="">
>> [0]: 12 ----> 4<br class="">
>> [0]: 13 ----> 2<br class="">
>> [0]: 13 ----> 3<br class="">
>> [0]: 14 ----> 2<br class="">
>> [0]: 14 ----> 5<br class="">
>> [0]: 15 ----> 6<br class="">
>> [0]: 15 ----> 7<br class="">
>> [0]: 15 ----> 9<br class="">
>> [0]: 16 ----> 6<br class="">
>> [0]: 16 ----> 8<br class="">
>> [0]: 16 ----> 10<br class="">
>> [0]: 17 ----> 7<br class="">
>> [0]: 17 ----> 8<br class="">
>> [0]: 17 ----> 11<br class="">
>> [0]: 18 ----> 9<br class="">
>> [0]: 18 ----> 12<br class="">
>> [0]: 18 ----> 13<br class="">
>> [0]: 19 ----> 10<br class="">
>> [0]: 19 ----> 12<br class="">
>> [0]: 19 ----> 14<br class="">
>> [0]: 20 ----> 11<br class="">
>> [0]: 20 ----> 13<br class="">
>> [0]: 20 ----> 14<br class="">
>> Cones:<br class="">
>> [0] Max cone size: 5<br class="">
>> [0]: 0 <---- 1 (0)<br class="">
>> [0]: 0 <---- 2 (0)<br class="">
>> [0]: 0 <---- 3 (0)<br class="">
>> [0]: 0 <---- 4 (0)<br class="">
>> [0]: 0 <---- 5 (0)<br class="">
>> [0]: 1 <---- 6 (0)<br class="">
>> [0]: 1 <---- 7 (0)<br class="">
>> [0]: 1 <---- 8 (0)<br class="">
>> [0]: 2 <---- 12 (0)<br class="">
>> [0]: 2 <---- 13 (0)<br class="">
>> [0]: 2 <---- 14 (0)<br class="">
>> [0]: 3 <---- 7 (0)<br class="">
>> [0]: 3 <---- 9 (0)<br class="">
>> [0]: 3 <---- 11 (0)<br class="">
>> [0]: 3 <---- 13 (0)<br class="">
>> [0]: 4 <---- 6 (0)<br class="">
>> [0]: 4 <---- 9 (0)<br class="">
>> [0]: 4 <---- 10 (0)<br class="">
>> [0]: 4 <---- 12 (0)<br class="">
>> [0]: 5 <---- 8 (0)<br class="">
>> [0]: 5 <---- 10 (0)<br class="">
>> [0]: 5 <---- 11 (0)<br class="">
>> [0]: 5 <---- 14 (0)<br class="">
>> [0]: 6 <---- 15 (0)<br class="">
>> [0]: 6 <---- 16 (0)<br class="">
>> [0]: 7 <---- 15 (0)<br class="">
>> [0]: 7 <---- 17 (0)<br class="">
>> [0]: 8 <---- 16 (0)<br class="">
>> [0]: 8 <---- 17 (0)<br class="">
>> [0]: 9 <---- 15 (0)<br class="">
>> [0]: 9 <---- 18 (0)<br class="">
>> [0]: 10 <---- 16 (0)<br class="">
>> [0]: 10 <---- 19 (0)<br class="">
>> [0]: 11 <---- 17 (0)<br class="">
>> [0]: 11 <---- 20 (0)<br class="">
>> [0]: 12 <---- 18 (0)<br class="">
>> [0]: 12 <---- 19 (0)<br class="">
>> [0]: 13 <---- 18 (0)<br class="">
>> [0]: 13 <---- 20 (0)<br class="">
>> [0]: 14 <---- 19 (0)<br class="">
>> [0]: 14 <---- 20 (0)<br class="">
>> coordinates with 1 fields<br class="">
>> field 0 with 3 components<br class="">
>> Process 0:<br class="">
>> ( 0) dim 0 offset 0<br class="">
>> ( 1) dim 0 offset 0<br class="">
>> ( 2) dim 0 offset 0<br class="">
>> ( 3) dim 0 offset 0<br class="">
>> ( 4) dim 0 offset 0<br class="">
>> ( 5) dim 0 offset 0<br class="">
>> ( 6) dim 0 offset 0<br class="">
>> ( 7) dim 0 offset 0<br class="">
>> ( 8) dim 0 offset 0<br class="">
>> ( 9) dim 0 offset 0<br class="">
>> ( 10) dim 0 offset 0<br class="">
>> ( 11) dim 0 offset 0<br class="">
>> ( 12) dim 0 offset 0<br class="">
>> ( 13) dim 0 offset 0<br class="">
>> ( 14) dim 0 offset 0<br class="">
>> ( 15) dim 3 offset 0 0. 0. 0.<br class="">
>> ( 16) dim 3 offset 3 0. 2. 0.<br class="">
>> ( 17) dim 3 offset 6 2. 1. 0.<br class="">
>> ( 18) dim 3 offset 9 0. 0. 1.<br class="">
>> ( 19) dim 3 offset 12 0. 2. 1.<br class="">
>> ( 20) dim 3 offset 15 2. 1. 1.<br class="">
>><br class="">
>><br class="">
>><br class="">
>><br class="">
>><br class="">
>><br class="">
>><br class="">
>><br class="">
>><br class="">
>> DM Object: testplex 1 MPI processes<br class="">
>> type: plex<br class="">
>> testplex in 3 dimensions:<br class="">
>> Supports:<br class="">
>> [0] Max support size: 3<br class="">
>> [0]: 1 ----> 0<br class="">
>> [0]: 2 ----> 0<br class="">
>> [0]: 3 ----> 0<br class="">
>> [0]: 4 ----> 0<br class="">
>> [0]: 5 ----> 0<br class="">
>> [0]: 6 ----> 1<br class="">
>> [0]: 6 ----> 4<br class="">
>> [0]: 7 ----> 1<br class="">
>> [0]: 7 ----> 3<br class="">
>> [0]: 8 ----> 1<br class="">
>> [0]: 8 ----> 5<br class="">
>> [0]: 9 ----> 3<br class="">
>> [0]: 9 ----> 4<br class="">
>> [0]: 10 ----> 4<br class="">
>> [0]: 10 ----> 5<br class="">
>> [0]: 11 ----> 3<br class="">
>> [0]: 11 ----> 5<br class="">
>> [0]: 12 ----> 2<br class="">
>> [0]: 12 ----> 4<br class="">
>> [0]: 13 ----> 2<br class="">
>> [0]: 13 ----> 3<br class="">
>> [0]: 14 ----> 2<br class="">
>> [0]: 14 ----> 5<br class="">
>> [0]: 15 ----> 6<br class="">
>> [0]: 15 ----> 7<br class="">
>> [0]: 15 ----> 9<br class="">
>> [0]: 16 ----> 6<br class="">
>> [0]: 16 ----> 8<br class="">
>> [0]: 16 ----> 10<br class="">
>> [0]: 17 ----> 7<br class="">
>> [0]: 17 ----> 8<br class="">
>> [0]: 17 ----> 11<br class="">
>> [0]: 18 ----> 9<br class="">
>> [0]: 18 ----> 12<br class="">
>> [0]: 18 ----> 13<br class="">
>> [0]: 19 ----> 10<br class="">
>> [0]: 19 ----> 12<br class="">
>> [0]: 19 ----> 14<br class="">
>> [0]: 20 ----> 11<br class="">
>> [0]: 20 ----> 13<br class="">
>> [0]: 20 ----> 14<br class="">
>> Cones:<br class="">
>> [0] Max cone size: 5<br class="">
>> [0]: 0 <---- 0 (0)<br class="">
>> [0]: 0 <---- 0 (0)<br class="">
>> [0]: 0 <---- 3 (0)<br class="">
>> [0]: 0 <---- 4 (0)<br class="">
>> [0]: 0 <---- 5 (0)<br class="">
>> [0]: 1 <---- 0 (0)<br class="">
>> [0]: 1 <---- 0 (0)<br class="">
>> [0]: 1 <---- 8 (0)<br class="">
>> [0]: 2 <---- 0 (0)<br class="">
>> [0]: 2 <---- 0 (0)<br class="">
>> [0]: 2 <---- 14 (0)<br class="">
>> [0]: 3 <---- 0 (0)<br class="">
>> [0]: 3 <---- 0 (0)<br class="">
>> [0]: 3 <---- 11 (0)<br class="">
>> [0]: 3 <---- 13 (0)<br class="">
>> [0]: 4 <---- 0 (0)<br class="">
>> [0]: 4 <---- 0 (0)<br class="">
>> [0]: 4 <---- 10 (0)<br class="">
>> [0]: 4 <---- 12 (0)<br class="">
>> [0]: 5 <---- 0 (0)<br class="">
>> [0]: 5 <---- 0 (0)<br class="">
>> [0]: 5 <---- 11 (0)<br class="">
>> [0]: 5 <---- 14 (0)<br class="">
>> [0]: 6 <---- 0 (0)<br class="">
>> [0]: 6 <---- 0 (0)<br class="">
>> [0]: 7 <---- 0 (0)<br class="">
>> [0]: 7 <---- 0 (0)<br class="">
>> [0]: 8 <---- 0 (0)<br class="">
>> [0]: 8 <---- 0 (0)<br class="">
>> [0]: 9 <---- 0 (0)<br class="">
>> [0]: 9 <---- 0 (0)<br class="">
>> [0]: 10 <---- 0 (0)<br class="">
>> [0]: 10 <---- 0 (0)<br class="">
>> [0]: 11 <---- 0 (0)<br class="">
>> [0]: 11 <---- 0 (0)<br class="">
>> [0]: 12 <---- 0 (0)<br class="">
>> [0]: 12 <---- 0 (0)<br class="">
>> [0]: 13 <---- 0 (0)<br class="">
>> [0]: 13 <---- 0 (0)<br class="">
>> [0]: 14 <---- 0 (0)<br class="">
>> [0]: 14 <---- 0 (0)<br class="">
>> coordinates with 1 fields<br class="">
>> field 0 with 3 components<br class="">
>> Process 0:<br class="">
>> ( 0) dim 0 offset 0<br class="">
>> ( 1) dim 0 offset 0<br class="">
>> ( 2) dim 0 offset 0<br class="">
>> ( 3) dim 0 offset 0<br class="">
>> ( 4) dim 0 offset 0<br class="">
>> ( 5) dim 0 offset 0<br class="">
>> ( 6) dim 0 offset 0<br class="">
>> ( 7) dim 0 offset 0<br class="">
>> ( 8) dim 0 offset 0<br class="">
>> ( 9) dim 0 offset 0<br class="">
>> ( 10) dim 0 offset 0<br class="">
>> ( 11) dim 0 offset 0<br class="">
>> ( 12) dim 0 offset 0<br class="">
>> ( 13) dim 0 offset 0<br class="">
>> ( 14) dim 0 offset 0<br class="">
>> ( 15) dim 3 offset 0 0. 0. 0.<br class="">
>> ( 16) dim 3 offset 3 0. 2. 0.<br class="">
>> ( 17) dim 3 offset 6 2. 1. 0.<br class="">
>> ( 18) dim 3 offset 9 0. 0. 1.<br class="">
>> ( 19) dim 3 offset 12 0. 2. 1.<br class="">
>> ( 20) dim 3 offset 15 2. 1. 1.<br class="">
>><br class="">
>><br class="">
> <br class="">
<br class="">
</blockquote>
</div>
<br clear="all" class="">
<div class=""><br class="">
</div>
-- <br class="">
<div dir="ltr" class="gmail_signature">
<div dir="ltr" class="">
<div class="">
<div dir="ltr" class="">
<div class="">
<div dir="ltr" class="">
<div class="">What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br class="">
-- Norbert Wiener</div>
<div class=""><br class="">
</div>
<div class=""><a href="http://www.cse.buffalo.edu/~knepley/" target="_blank" class="">https://www.cse.buffalo.edu/~knepley/</a><br class="">
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</blockquote>
</div>
<br class="">
</body>
</html>