<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Mon, Feb 17, 2014 at 1:23 PM, Dharmendar Reddy <span dir="ltr"><<a href="mailto:dharmareddy84@gmail.com" target="_blank">dharmareddy84@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hello,<br>
         My petsc code is at version e4acbc4 on next. I get undefined<br>
reference error:<br>
undefined reference to `petscdualspacecreatereferencecell_'<br>
<br>
Here is the code, i used.<br>
<br>
   print*,'Testing Reference Cell'<br>
   call PetscDualSpaceCreate(comm, Q, ierr)<br>
   Kdim=3<br>
   call PetscDualSpaceCreateReferenceCell(Q,Kdim, PETSC_TRUE, K, ierr)<br>
   call PetscViewerPushFormat(PETSC_VIEWER_STDOUT_WORLD,<br>
PETSC_VIEWER_ASCII_INFO_DETAIL, ierr)<br>
   call DMView(K, PETSC_VIEWER_STDOUT_WORLD, ierr)<br></blockquote><div><br></div><div>I am very sorry. This fell through the email cracks. I have put in the Fortran bindings, and docs</div><div>for all the PetscDualSpace functions into next.</div>
<div><br></div><div>  Thanks,</div><div><br></div><div>    Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Thanks<br>
Reddy<br>
<br>
On Mon, Feb 17, 2014 at 7:28 AM, Matthew Knepley <<a href="mailto:knepley@gmail.com">knepley@gmail.com</a>> wrote:<br>
> On Tue, Feb 11, 2014 at 7:50 AM, Dharmendar Reddy <<a href="mailto:dharmareddy84@gmail.com">dharmareddy84@gmail.com</a>><br>
> wrote:<br>
>><br>
>> Hello,<br>
>>        Where do i find the information on reference elements used by<br>
>> plex for local numbering.<br>
><br>
><br>
> You can use<br>
><br>
>   PetscDualSpaceCreateReferenceCell(sp, dim, isSimplex, &refCell)<br>
>   PetscViewerPushFormat(PETSC_VIEWER_STDOUT_WORLD,<br>
> PETSC_VIEWER_ASCII_INFO_DETAIL)<br>
>   DMView(refCell, PETSC_VIEWER_STDOUT_WORLD)<br>
><br>
>     Matt<br>
><br>
>> I currently use the following for the tetrahedron. Is this correct ?<br>
>> what is the corresponding information for a  triangle ?<br>
>><br>
>>         ! Vertex coordiantes<br>
>>        vert(1:3,1)=[0.0, 0.0, 0.0]<br>
>>         vert(1:3,2)=[1.0, 0.0, 0.0]<br>
>>         vert(1:3,3)=[0.0, 1.0, 0.0]<br>
>>         vert(1:3,4)=[0.0, 0.0, 1.0]<br>
>><br>
>>         ! vertex intdices<br>
>>         vs(1:4)= [1, 2, 3, 4]<br>
>><br>
>>         ! Edges id to node Id<br>
>>         edge(:,1) = [3, 4]<br>
>>         edge(:,2) = [2, 4]<br>
>>         edge(:,3) = [2, 3]<br>
>>         edge(:,4) = [1, 4]<br>
>>         edge(:,5) = [1, 3]<br>
>>         edge(:,6) = [1, 2]<br>
>><br>
>>         ! Faces id to node id<br>
>>         face(1:3,1)=[2,3,4]<br>
>>         face(1:3,2)=[1,3,4]<br>
>>         face(1:3,3)=[1,2,4]<br>
>>         face(1:3,4)=[1,2,3]<br>
>>         numFacet = 4<br>
>><br>
>>         ! Tets<br>
>>         tet(1:4,1) = [1,2,3,4]<br>
>> thanks<br>
>> Reddy<br>
><br>
><br>
><br>
<span class="HOEnZb"><font color="#888888">><br>
> --<br>
> What most experimenters take for granted before they begin their experiments<br>
> is infinitely more interesting than any results to which their experiments<br>
> lead.<br>
> -- Norbert Wiener<br>
</font></span></blockquote></div><br><br clear="all"><div><br></div>-- <br>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>