On Tue, Aug 28, 2012 at 4:58 PM, Chris Eldred <span dir="ltr"><<a href="mailto:chris.eldred@gmail.com" target="_blank">chris.eldred@gmail.com</a>></span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Building on the previous example code (for the Doublet mesh), I am<br>
getting an error when trying to set Cone orientations:<br>
<br>
[0]PETSC ERROR: --------------------- Error Message<br>
------------------------------------<br>
[0]PETSC ERROR: Argument out of range!<br>
[0]PETSC ERROR: Cone orientation 0 is not in the valid range [-1. 0)!<br>
[0]PETSC ERROR:<br>
------------------------------------------------------------------------<br>
[0]PETSC ERROR: Petsc Development HG revision:<br>
b6fe82991deee4a0d3f9a20654bc7750a6b1fe0f  HG Date: Mon Aug 27 13:17:10<br>
2012 -0500<br>
[0]PETSC ERROR: See docs/changes/index.html for recent updates.<br>
[0]PETSC ERROR: See docs/faq.html for hints about trouble shooting.<br>
[0]PETSC ERROR: See docs/index.html for manual pages.<br>
[0]PETSC ERROR:<br>
------------------------------------------------------------------------<br>
[0]PETSC ERROR: ./test on a arch-linu named Puget-101334 by user Tue<br>
Aug 28 15:54:36 2012<br>
[0]PETSC ERROR: Libraries linked from<br>
/home/user/Desktop/LIBRARIES/petsc-dev/arch-linux2-cxx-debug/lib<br>
[0]PETSC ERROR: Configure run at Mon Aug 27 12:55:14 2012<br>
[0]PETSC ERROR: Configure options --download-boost --download-chaco<br>
--download-ctetgen --download-f-blas-lapack --download-fiat<br>
--download-generator --download-metis --download-ml --download-mpich<br>
--download-parmetis --download-scientificpython --download-triangle<br>
--with-clanguage=cxx --with-dynamic-loading --with-shared-libraries<br>
--with-sieve PETSC_ARCH=arch-linux2-cxx-debug<br>
[0]PETSC ERROR:<br>
------------------------------------------------------------------------<br>
[0]PETSC ERROR: DMComplexSetConeOrientation() line 1207 in<br>
/home/user/Desktop/LIBRARIES/petsc-dev/src/dm/impls/complex/complex.c<br>
application called MPI_Abort(MPI_COMM_WORLD, 63) - process 0<br>
[unset]: aborting job:<br>
application called MPI_Abort(MPI_COMM_WORLD, 63) - process 0<br>
<br>
The code is:<br>
<br>
Orientations(1) = 0<br>
Orientations(2) = 1<br>
pOrientations => Orientations<br>
call DMComplexSetConeOrientation(dm,4,pOrientations,ierr)<br>
CHKERRQ(ierr)<br>
call DMComplexGetConeOrientation(dm,4,pOrientations,ierr)<br>
CHKERRQ(ierr)<br>
write(*,*) pOrientations<br>
<br>
Mesh point 4 is an edge (so it has a cone of size 2)- It is unclear to<br>
me why the cone orientation range is [-1,0). My understanding of cone<br></blockquote><div><br></div><div>This means that the number of dof on that point is 0.</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">

orientations is that they specified the order in which cone's were<br>
traversed for a point- ie for mesh point 4 (an edge) I could give it a<br>
cone Orientation of (1,0) which would mean look at point 1 of the cone<br>
first, then point 0.<br>
<span class="HOEnZb"><font color="#888888"><br>
--<br>
Chris Eldred<br>
DOE Computational Science Graduate Fellow<br>
Graduate Student, Atmospheric Science, Colorado State University<br>
B.S. Applied Computational Physics, Carnegie Mellon University, 2009<br>
<a href="mailto:chris.eldred@gmail.com">chris.eldred@gmail.com</a><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<br>