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 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
