On Tue, Aug 28, 2012 at 4:58 PM, Chris Eldred <chris.eldred at gmail.com>wrote:
> Building on the previous example code (for the Doublet mesh), I am > getting an error when trying to set Cone orientations: > > [0]PETSC ERROR: --------------------- Error Message > ------------------------------------ > [0]PETSC ERROR: Argument out of range! > [0]PETSC ERROR: Cone orientation 0 is not in the valid range [-1. 0)! > [0]PETSC ERROR: > ------------------------------------------------------------------------ > [0]PETSC ERROR: Petsc Development HG revision: > b6fe82991deee4a0d3f9a20654bc7750a6b1fe0f HG Date: Mon Aug 27 13:17:10 > 2012 -0500 > [0]PETSC ERROR: See docs/changes/index.html for recent updates. > [0]PETSC ERROR: See docs/faq.html for hints about trouble shooting. > [0]PETSC ERROR: See docs/index.html for manual pages. > [0]PETSC ERROR: > ------------------------------------------------------------------------ > [0]PETSC ERROR: ./test on a arch-linu named Puget-101334 by user Tue > Aug 28 15:54:36 2012 > [0]PETSC ERROR: Libraries linked from > /home/user/Desktop/LIBRARIES/petsc-dev/arch-linux2-cxx-debug/lib > [0]PETSC ERROR: Configure run at Mon Aug 27 12:55:14 2012 > [0]PETSC ERROR: Configure options --download-boost --download-chaco > --download-ctetgen --download-f-blas-lapack --download-fiat > --download-generator --download-metis --download-ml --download-mpich > --download-parmetis --download-scientificpython --download-triangle > --with-clanguage=cxx --with-dynamic-loading --with-shared-libraries > --with-sieve PETSC_ARCH=arch-linux2-cxx-debug > [0]PETSC ERROR: > ------------------------------------------------------------------------ > [0]PETSC ERROR: DMComplexSetConeOrientation() line 1207 in > /home/user/Desktop/LIBRARIES/petsc-dev/src/dm/impls/complex/complex.c > application called MPI_Abort(MPI_COMM_WORLD, 63) - process 0 > [unset]: aborting job: > application called MPI_Abort(MPI_COMM_WORLD, 63) - process 0 > > The code is: > > Orientations(1) = 0 > Orientations(2) = 1 > pOrientations => Orientations > call DMComplexSetConeOrientation(dm,4,pOrientations,ierr) > CHKERRQ(ierr) > call DMComplexGetConeOrientation(dm,4,pOrientations,ierr) > CHKERRQ(ierr) > write(*,*) pOrientations > > Mesh point 4 is an edge (so it has a cone of size 2)- It is unclear to > me why the cone orientation range is [-1,0). My understanding of cone > This means that the number of dof on that point is 0. Matt > orientations is that they specified the order in which cone's were > traversed for a point- ie for mesh point 4 (an edge) I could give it a > cone Orientation of (1,0) which would mean look at point 1 of the cone > first, then point 0. > > -- > Chris Eldred > DOE Computational Science Graduate Fellow > Graduate Student, Atmospheric Science, Colorado State University > B.S. Applied Computational Physics, Carnegie Mellon University, 2009 > chris.eldred at gmail.com > -- 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-dev/attachments/20120828/735affb1/attachment-0001.html>