[petsc-users] Plex mesh reference element

[Matthew Knepley]

On Mon, Feb 17, 2014 at 1:23 PM, Dharmendar Reddy <dharmareddy84 at gmail.com>
wrote:

> Hello,
>          My petsc code is at version e4acbc4 on next. I get undefined
> reference error:
> undefined reference to `petscdualspacecreatereferencecell_'
>
> Here is the code, i used.
>
>    print*,'Testing Reference Cell'
>    call PetscDualSpaceCreate(comm, Q, ierr)
>    Kdim=3
>    call PetscDualSpaceCreateReferenceCell(Q,Kdim, PETSC_TRUE, K, ierr)
>    call PetscViewerPushFormat(PETSC_VIEWER_STDOUT_WORLD,
> PETSC_VIEWER_ASCII_INFO_DETAIL, ierr)
>    call DMView(K, PETSC_VIEWER_STDOUT_WORLD, ierr)
>

I am very sorry. This fell through the email cracks. I have put in the
Fortran bindings, and docs
for all the PetscDualSpace functions into next.

  Thanks,

    Matt


> Thanks
> Reddy
>
> On Mon, Feb 17, 2014 at 7:28 AM, Matthew Knepley <knepley at gmail.com>
> wrote:
> > On Tue, Feb 11, 2014 at 7:50 AM, Dharmendar Reddy <
> dharmareddy84 at gmail.com>
> > wrote:
> >>
> >> Hello,
> >>        Where do i find the information on reference elements used by
> >> plex for local numbering.
> >
> >
> > You can use
> >
> >   PetscDualSpaceCreateReferenceCell(sp, dim, isSimplex, &refCell)
> >   PetscViewerPushFormat(PETSC_VIEWER_STDOUT_WORLD,
> > PETSC_VIEWER_ASCII_INFO_DETAIL)
> >   DMView(refCell, PETSC_VIEWER_STDOUT_WORLD)
> >
> >     Matt
> >
> >> I currently use the following for the tetrahedron. Is this correct ?
> >> what is the corresponding information for a  triangle ?
> >>
> >>         ! Vertex coordiantes
> >>        vert(1:3,1)=[0.0, 0.0, 0.0]
> >>         vert(1:3,2)=[1.0, 0.0, 0.0]
> >>         vert(1:3,3)=[0.0, 1.0, 0.0]
> >>         vert(1:3,4)=[0.0, 0.0, 1.0]
> >>
> >>         ! vertex intdices
> >>         vs(1:4)= [1, 2, 3, 4]
> >>
> >>         ! Edges id to node Id
> >>         edge(:,1) = [3, 4]
> >>         edge(:,2) = [2, 4]
> >>         edge(:,3) = [2, 3]
> >>         edge(:,4) = [1, 4]
> >>         edge(:,5) = [1, 3]
> >>         edge(:,6) = [1, 2]
> >>
> >>         ! Faces id to node id
> >>         face(1:3,1)=[2,3,4]
> >>         face(1:3,2)=[1,3,4]
> >>         face(1:3,3)=[1,2,4]
> >>         face(1:3,4)=[1,2,3]
> >>         numFacet = 4
> >>
> >>         ! Tets
> >>         tet(1:4,1) = [1,2,3,4]
> >> thanks
> >> Reddy
> >
> >
> >
> >
> > --
> > 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
>



-- 
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/20140611/dd3d86bf/attachment.html>