[petsc-users] DMPlexGetSupport

[Matthew Knepley]

On Mon, Apr 29, 2013 at 4:02 AM, Dharmendar Reddy
<dharmareddy84 at gmail.com>wrote:

> Hello,
>           If i extract the support and cone using
> DMPlexGetTransitiveclousre, I am getting correct answers. Please have a
> look at the attached code.
>

Yes, the interface definition for DMPlexRestoreCone/Support() was missing
from petscdmplex.h90, so Fortran
was just stomping through random memory without giving a warning. I pushed
the fix to 'next', and there is now
a test of it in DMPlex ex1f90.

  Thanks,

     Matt


> On Sun, Apr 28, 2013 at 9:30 PM, Dharmendar Reddy <dharmareddy84 at gmail.com
> > wrote:
>
>> Hello,
>>          Looks like the problem is with getCone and getSupport called
>> from F90 code.
>>
>> Thanks
>> Reddy
>>
>>
>> On Sun, Apr 28, 2013 at 7:36 PM, Dharmendar Reddy <
>> dharmareddy84 at gmail.com> wrote:
>>
>>> Hello,
>>>           I am seeing an issue using DMPlexGetSupport.
>>>
>>> I am trying to get the Ids of cells sharing a given facet in mesh.
>>>
>>> I have attached a test code for a small 3 x 3 square mesh. When i print
>>> the support for the 12 facets it has, i see 4 internal facets and 8
>>> boundary faces, as expected. But the cell Ids are wrong.
>>>
>>>  Surcessfully created DM From Cell List
>>> Mesh in 2 dimensions:
>>>   0-cells: 9
>>>   1-cells: 12
>>>   2-cells: 4
>>> Labels:
>>>   depth: 3 strata of sizes (9, 12, 4)
>>>  printing support for facet
>>>  edge id start=          13 edge id end=          24
>>>  f:          13 c:           0
>>>  f:          14 c:           0           1
>>>  f:          15 c:           0           2
>>>  f:          16 c:           0
>>>  f:          17 c:           0
>>>  f:          18 c:           0
>>>  f:          19 c:           0           0
>>>  f:          20 c:           0           0
>>>  f:          21 c:           0
>>>  f:          22 c:           0
>>>  f:          23 c:           0
>>>  f:          24 c:           0
>>>
>>> --
>>> -----------------------------------------------------
>>> Dharmendar Reddy Palle
>>> Graduate Student
>>> Microelectronics Research center,
>>> University of Texas at Austin,
>>> 10100 Burnet Road, Bldg. 160
>>> MER 2.608F, TX 78758-4445
>>> e-mail: dharmareddy84 at gmail.com
>>> Phone: +1-512-350-9082
>>> United States of America.
>>> Homepage: https://webspace.utexas.edu/~dpr342
>>>
>>
>>
>>
>> --
>> -----------------------------------------------------
>> Dharmendar Reddy Palle
>> Graduate Student
>> Microelectronics Research center,
>> University of Texas at Austin,
>> 10100 Burnet Road, Bldg. 160
>> MER 2.608F, TX 78758-4445
>> e-mail: dharmareddy84 at gmail.com
>> Phone: +1-512-350-9082
>> United States of America.
>> Homepage: https://webspace.utexas.edu/~dpr342
>>
>
>
>
> --
> -----------------------------------------------------
> Dharmendar Reddy Palle
> Graduate Student
> Microelectronics Research center,
> University of Texas at Austin,
> 10100 Burnet Road, Bldg. 160
> MER 2.608F, TX 78758-4445
> e-mail: dharmareddy84 at gmail.com
> Phone: +1-512-350-9082
> United States of America.
> Homepage: https://webspace.utexas.edu/~dpr342
>



-- 
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20130429/849544a1/attachment-0001.html>